Modelling the spatial dynamics of plasmid transfer and persistence

Abstract

Bacterial plasmids are extra-chromosomal genetic elements that code for a wide variety of phenotypes in their bacterial hosts and are maintained in bacterial communities through both vertical and horizontal transfer. Current mathematical models of plasmid-bacteria dynamics, based almost exclusively on mass-action differential equations that describe these interactions in completely mixed environments, fail to adequately explain phenomena such as the long-term persistence of plasmids in natural and clinical bacterial communities. This failure is, at least in part, due to the absence of any spatial structure in these models, whereas most bacterial populations are spatially structured in microcolonies and biofilms. To help bridge the gap between theoretical predictions and observed patterns of plasmid spread and persistence, an individual-based lattice model (interacting particle system) that provides a predictive framework for understanding the dynamics of plasmid-bacteria interactions in spatially structured populations is presented here. To assess the accuracy and flexibility of the model, a series of experiments that monitored plasmid loss and horizontal transfer of the IncP-1beta plasmid pB10 : : rfp in Escherichia coli K12 and other bacterial populations grown on agar surfaces were performed. The model-based visual patterns of plasmid loss and spread, as well as quantitative predictions of the effects of different initial parental strain densities and incubation time on densities of transconjugants formed on a 2D grid, were in agreement with this and previously published empirical data. These results include features of spatially structured populations that are not predicted by mass-action differential equation models.