A Bead Movement Based Computational Framework for 3-Dimensional Analysis of Biofilm Material Heterogeneity

Differences in the material properties of bacterial biofilms have been observed in biofilms of different bacterial species, within the same species under different growth conditions and after treatment with matrix modifying molecules. To better quantitate the material properties of 3D biofilms, an experimental and computational workflow was developed and applied to examine differences between Enterococcus faecalis, Salmonella enterica serotype Typhimurium and Escherichia coli biofilms as well as the role of the amyloid curli in confirming rigidity to Enterobacteriaceae biofilms. The spatio-temporal dynamics of 1 µm carboxylate beads in biofilms were tracked in 20 µm 3D biofilms over 20 minutes. The 4D image stacks were processed using the Mosaic plugin in ImageJ to produce 3D trajectory data of bead movement. This trajectory data was analyzed with a newly developed Bead Evaluator toolbox, where movement data, including trajectory lifespans, bead velocities, cell densities along trajectories, and bounding box information were computed and stored in csv-files. This paper presents the workflow from experimental setup and image recording to bead trajectory computation and analysis. The structure of curli-containing biofilms resulted in more stable bead interactions and less bead movement than in curli-mutant and Enterococcal biofilms. Bead movement did not appear strongly dependent on cell density when measuring the bead velocity and trajectory bounding box volume, supporting the hypothesis that other material properties of the biofilms control the bead dynamics. This technique is widely applicable to quantitating differences in biofilms of different matrix compositions as well as biofilms before and after matrix-modifying treatments.