Spatiotemporal establishment of dense bacterial colonies growing on hard agar
Figures

Experimental observations of the growth and morphology of a bacterial colony.
(A–E) Confocal images of an E. coli colony harboring GFP expression growing on 1.5% agar (glucose minimal medium) taken at various time after seeding (). The center of the colony is indicated by the red dot. Single- and multi-layer regions are distinguished by red circles based on fluorescence intensity; see 'Experimental Methods'. (F) The radius of the first (red) and second layer (orange) of the colony, as well as their difference (green), versus time. (G) The cross-sectional profile of the growing bacterial colony at indicated time after single-cell inoculation. (H) After the buckling at around h, the colony radius (red symbols) increased at a constant speed (red line), while the colony height (blue symbols) increased linearly with speed (blue line). The latter slowed down some time after . (I) The dependence of the radial speed (red symbols) and the vertical speed (blue symbols) on cell growth rate (x-axis), for colonies grown in minimal medium with 8 different carbon sources (Supplementary file 1-Table S1): glucose (O); arabinose (□); mannitol (); maltose (◇); fructose (); melibiose (); sorbitol (); mannose (☆). The lines are best linear fit of the data.
-
Figure 1—source data 1
Experimental data for the temporal development of colony profiles and velocities.
- https://doi.org/10.7554/eLife.41093.005

Data for five repeats of E.coli EQ59 grown on 1.5% (w/v) agar in minimal medium with 0.2% glucose (11 mM), and incubated, covered, at 37°C for up to 3 days; cf. 'Experimental Methods'.
Their radii and heights were periodically monitored using a confocal microscope. To monitor the colony growth over long periods of time, we started with identical colonies at seed time separated by several hours. Growth curves extending over a period of multiple days were obtained by 'stitching together' the radii and height data at times where they overlapped; cf. section on Experimental Methods. The recorded radii data (A)-(E) and height data (G)-(K) clearly showed linear regimes. The data were fitted by straight lines over the linear regimes to obtain the speed of radial expansion (F) and that of the height increasing (L) for these five repeated experiments. has an average of and standard deviation of . has an average of and standard deviation of . Similar analyses were done for colonies grown in other carbon sources described in Figure 1I.
-
Figure 1—figure supplement 1—source data 1
Repetitions for the temporal development of colony height and radius.
- https://doi.org/10.7554/eLife.41093.004

Schematics of cell-cell, cell-agar, cell-fluid, and surface tension forces investigated in this study.
Green area indicates the colony (with cells in yellow). Blue area indicates the agar.

The simulated growth and morphology of a bacterial colony.
(A) A semi-log plot of number of cells vs. time showing the exponential growth of the population starting from a single cell at . Red line shows exponential growth rate of . (B) Cross sections of the growing colony at various times after seeding of a single cell. (C) Plots of the radius (red) and height (blue) vs. time up to , showing that, after an initial transient period of ~10 h, the growing colony increases linearly at the radial speed (red line) and vertical speed . The blue arrow at indicates the time when the colony height starts to increase. (D) The radius of the first (red) and second (orange) layer of the colony as well as their difference (green) vs. time. (E–M) Top view of the colony at various time. Cells in the bottom layer (blue) and upper layers (yellow) are fitted into red circles. The time evolution of buckling phenomenon is captured in detail in (J–M).

The spatially varying cell density (per unit volume of colony) is related to the spatially varying cell volume fraction by , where the volume fraction is defined as the volume of all cells in a unit volume of the colony and is the constant mass density of a typical mature cell; cf. Appendix 1.2 on nutrient update.
(A) and (B): The volume fraction in the cross section of the colony, calculated with the mean and standard deviation , with (A) , , and ; and (B) , , and . (C) and (D): The effect of coarse-graining size to the volume fraction . The average of (C) and the standard deviation of (D) in the coarse-graining with different sizes.

The cross-sectional anatomy of a simulated colony.
(A) Snapshot of cross-sectional view of the colony at . Cyan represents horizontally oriented cells ( with z-axis); Yellow represents vertically oriented cells ( with z-axis). (B) Fraction of vertically oriented cells averaged over vs radius. (C) A side view of the azimuthally averaged director field, indicating the orientation of the rod-like cells. (D) A side view of the azimuthally averaged velocity field. (E) Vertical component of velocity, , at various values of along the center of the colony. Increase in vertical speed is seen only for the bottom 10 µm (F) A cross-sectional view of the colony, color representing the time since last division. Purple and blue represent cells that have not divided for the past 10 h, and red represents the actively dividing cells. (G) A cross-sectional view of the local growth rate in the colony, with the color bar showing the values of local growth rate. A disc-shaped 'growth zone' is revealed by the red color at the bottom of the colony.

A few top layers of cells in the colony visualized using simulation data.
(A) and experimental data (B). In (A), a cell is colored green (resp. red) if its angle with the vertical -axis is less (resp. greater) than .

Fraction of verticalized cells increase for large colonies.
(A) The bottom view of a simulated colony at , with the cells outlined by the black curves. Overlapping cells are of different layers. (B)-(D) show the fraction of vertically oriented cells (cells with an angle with axis) plotted versus the radius of the colony (B), the height of the colony up to 80 µm (C), and the height of the colony up to 10 µm (D).

Vertical penetration of nutrients.
(A) The profiles of nutrient concentration along the z-axis at different times. (B) The profile in the uniform scale vs. that in the rescaled z-axis, . (C) The numerical results for the height velocity vs. the batch culture growth rate with a fixed cell division length (open circles) and variable (asterisk), respectively. The square root fit for the open circles (solid line) is given by the expression ; the linear fit for circles with (dashed line) is given by the expression .

Semi-log plot of the steady-state nutrient profile : reconstructed from 3D simulations (green *); and the numerical solution to the 1D model (cf. Appendix 2.3 on nutrient penetration) with discretization (green circles) and (black line), respectively.
The linear part below the line of the Monod constant (dashed blue line) indicates the exponential decay of the concentration. Shown in the window is the quadratic curve (red dotted line), fitting the numerical solution to the 1D model with for the concentration value above the Monod constant.

Profiles of nutrient concentration versus for various values of the batch culture growth rate .
These nutrient profiles were reconstructed from 3D simulations. A larger value of the batch culture growth rate leads to a lower nutrient concentration, indicating a faster consumption of nutrients by cells.

Coarse-grained view of director, velocity and pressure fields in the bottom layer of colony.
(A) The bottom view of a simulated colony. Color scheme is the same as in Figure 4A. (B) Bars show the planar component of the coarse-grained director field at the bottom layer. (C) Arrows show the planar component of the coarse-grained velocity field at the bottom layer. Colors in (B) and (C) indicate the local pressure; see scale bar on the far right. The pressure is expressed in unit of where is the surface tension, the main force underlying pressure build-up in our model.

Spatiotemporal nutrient profiles.
(A) The xz cross sectional view of the nutrient concentration inside the colony and in the agar, at time . (B) The nutrient profile at the agar surface for different times. The nutrient concentration vs. time at the center (C) and at the periphery (D) of the colony at the agar surface.

Physical characteristics near the outer periphery of the colony.
(A) The azimuthally averaged local speed of radial expansion vs. the signed distance from the edge of the colony (cf. Equations (A1.5.1) in Appendix 1) at various times. (B) The azimuthally averaged local height vs. at various times. (C) The azimuthally averaged local speed of vertical expansion vs., averaging over hr. (D) The azimuthally averaged local pressure vs. at time As in Figure 6, pressure is expressed in unit of . (E) The azimuthally averaged height vs. at growth rate and . (F) The simulated colony horizontal expansion speed vs. the batch culture growth rate with a fixed (red open circles) and a growth-rate dependent (red closed circles) using dynamic friction. The dashed line that fits the open circles is given by ; the solid line that fits the closed circles for is ; the dash dotted line fits the closed circles with . In these expressions, the speeds are in unit of and growth rate in . For comparison, we also include simulated vs. with a growth-rate dependent , for a model with static friction alone (see Equations (A1.4.6) of Appendix 1) between cell and agar (blue triangles).

The azimuthally averaged and rescaled local growth rate as function of the signed distance to the colony rim with various values of the batch culture growth rate .
Precise definition of this signed distance is given in Equations (A1.5.1) in Appendix 1.5 on coarse-grained variables. At the rim , the nutrient concentration is close to , the constant concentration value in the boundary condition, and is close to , which is with our choice of parameters and . (A) The average is taken over the entire colony, where the variation of the (rescaled) local growth rate is within . (B) The average is taken over of the colony near the rim, where the variation of the (rescaled) local growth rate is within .

A zoomed in view of the periphery of the colony shown in Figure 4A, overlaid with coarse-grained velocity field (zoomed in view of the same periphery region in Figure 4D).
is the signed distance from the edge of the colony. is the buckling width.

Parameter dependence of colony growth characteristics.
Simulation results using the full model with 2x increase in glucose concentration (panel A) and 4x decrease in all frictional parameters (panel B) for (red bars) and (blue bars). Specifically, in panel (A) we fix the friction at a high level (with , and ), and use as the lower glucose concentration, as the higher glucose concentration. In panel (B), we fix the glucose concentration at the lower level (), and vary the friction, from the higher value of , used in (A) to the lower value of , . The corresponding experimental results are shown in panels C and D: In (C), glucose concentration was varied with agar density fixed at 1.5%. In (D), agar density was varied with glucose fixed at 0.2% (w/v). The data for in panel C is consistent with a square root dependence on nutrient concentration (blue line) expected from the basic analysis in Figure 5.
-
Figure 9—source data 1
Experimental data on the horizontal and vertical colony expansion speeds at various glucose and agar concentrations.
- https://doi.org/10.7554/eLife.41093.021

A schematic summary of key mechanisms in the growth of an E.coli colony.
After an initial, exponential monolayer growth, buckling occurs at the center of the colony. Cells then grow actively only in the bottom layers (red vertical arrows) whose thickness () is determined by the nutrient penetration level (dashed blue line). Cells lying above them are passively pushed up. Throughout this yellow triangular region, cells are oriented vertically. Near the colony edge (cyan region), the cells are oriented planarly and grow outward (horizontal red arrow) continuously in a spread mode to expand the colony in the radial direction. The width of this annulus () is determined by mechanical effects arising from the surface tension which pulls the thin layer of cells into the agar, and cell-agar friction which builds up the pressure from the outer edge of the layer, eventually causing buckling at an inner radius where cells transition to the vertical orientation (the green region). These two characteristic parameters, and , set the speeds of radial and vertical expansions, and , respectively, as shown in red. The growth rate dependence of these parameters is shown in blue.

Schematic of the computational box and different regions in the model of simulation.
The computational box is , where all , , and are positive numbers in the units of length. This box is divided into the air region , colony region , and agar region , respectively. See Supplementary file 1-Table S5 for typical values of , , and used in our simulations. The colony surface or colony-air interface separates the colony from air. The plane in the computational box is divided into two parts. One is the interface that separates the colony from agar, and is denoted by . The remaining part, denoted , separates the air from agar. Note that, since the bacterial colony grows with time , all the air region , the colony region , the colony-air interface , and the colony-agar interface depend on time .

Schematic view of cell growth and division.
(A) A sphero-cylinder model of an E.coli cell. Here, is the diameter of each of the hemispheres, and are the centers of these hemispheres, is the cylindrical length of the cell, is the unit vector along the cylindrical axis of the cell, and is the center of the cell. (B) Cell division. Once the cylindrical length of a cell reaches a critical value , the cell divides into two daughter cells. The two centers of hemispheres of the mother cell become the centers of hemispheres of the daughter cells. Each of these two daughter cells has the cylindrical length with fluctuations, where is a constant cylindrical length for any new born cell and any initial cell in the simulation. Fluctuations of angular velocities are also introduced for the daughter cells; cf. Appendix 1.3.

Schematic view of cell-cell and cell-agar interactions.
(A) Cell-cell interaction. Two cells, centered at and , respectively, are in contact with each other. The shortest distance between the central cylindrical lines of the two cells is with and two points on the cylindrical central lines of these two cells, respectively. The amount of the overlap of these two cells is . We shall denote by the center of the line segment connecting and The total interaction force, exerted at center , is the sum of the normal force in the normal direction and the tangential force in a direction orthogonal to that is determined by the relative velocities of these two cells; see the details in Appendix 1.4 on force calculations. (B) Cell-agar interaction. A cell, centered at , touches the agar surface that has the mean position at and the roughness (the maximum fluctuation around the mean), with the amount of the overlap of the cell and agar. If the center of the hemispherical cap of the cell corresponding to the end that dips into the agar is , then . Denote , which is the midpoint of the line segment along the vertical line passing through the point between and . The total cell-agar interaction force , exerted at the center , is the superposition of a normal force in the vertical direction and the tangential force in a direction along the plane that is determined by the velocity of the cell; cf. Appendix 1.4.

Schematic view of surface tension acting on cells at the colony boundary.
(A) a snapshot of part of a growing colony from a typical computer simulation. The red curve defines the macroscopic colony-air interface which is used in updating the nutrient profile; cf. in Figure 11 and Appendix 1.1. Cells on the top of the colony are held down by the surface tension force. (B) A sequence of configurations of colony during the early stage of growth. Top: Initial cells grow exponentially and form a monolayer. All cells in this layer are held down by the surface tension force. Middle: As more cells are born in the monolayer, the frictions between these cells and the rough agar surface increase. The competition between such frictions and the surface tension force that pulls down cells leads to an accumulation of the lateral pressure in the cells. The monolayer buckles up once that the surface tension can no longer hold down all the cells in the monolayer. Such buckling occurs at certain radial distance from the center of colony; cf. Figure 8D. Bottom: Once the monolayer buckles, the colony starts to grow vertically. In the meantime, the buckling region moves outward as the colony expands radially. (C) The parameters used in the definition of surface tension force. The parameter is the height of a cell that sticks out of the macroscopic colony surface; it is measured from the mean height of the agar surface to the 'tallest point' in the cell. The blue and red lines describe the macroscopic water level (indicated by ) and colony height (indicated by ), respectively. The parameter is used to control how tightly the surface tension holds back those cells on the top of the colony. See more details in Appendix 1.4c.

Schematic view of a nested finite difference grid for the agar region .
https://doi.org/10.7554/eLife.41093.029
The standard (A) and modified (B) model for friction between two objects, as a function of tangential relative velocity between the two objects, where and are the static and dynamic friction constants, respectively.
https://doi.org/10.7554/eLife.41093.030
Schematic of the derivation of the surface tension force.
https://doi.org/10.7554/eLife.41093.031
Experimental measurement on and with various batch culture growth rates.
Data are fitted with the stright line and the squre root curve for the redial and vertical speeds of expansion and respectively. (B) Simulation results on and with various batch culture growth rates and with a variable diviting length . data are fitted for using the straight line the squre-root curve respectively.

Statistics of cells undergoing growth transition.
(A) The histogram of the number of generations of the cells that experienced high-to-low growth rate transition.(B) Curves of growth rates vs. shifted time for 100 cells randomly selected from those 80% cells.

Lack of verticalization for thin colonies.
(A) Reproduction of Figure 2A from (Su et al., 2012). (B) Bottom view of a simulated colony at t=7h. (C) The fraction of vertically oriented cells (defined by the angle with the z-axis less than 45 degrees) vs. height from the simulation.

Experimental results on the confocal images of a bacterial colony at different positions, from the center to the edge.
The bacterial strain and the growth medium are the same as described in “Experimental Methods” in the main text, except that the thickness of agar dish is ~0.3mm. The picture was taken at roughly 24 hours after the initial inoculation.

The dynamic and static friction models distinguished by the growth-rate dependence of radial expansion speed.
(A) From main text Figure 1F. Experimental results on the radial speed VR and the vertical speed VH with respect to various batch culture growth rates. (B) From Figure 8F in the main text (new version). Simulation results of radial speed VR of the colony, where both static and dynamic frictions are included in simulations. (C) New simulation results of radial speed VR of the colony, where only the static friction is included.
Additional files
-
Supplementary file 1
Supplemental tables S1-5.
- https://doi.org/10.7554/eLife.41093.027