AbstractIn this study, we present a mathematical model for plasmid spread in a growing biofilm, formulated as a nonlocal system of partial differential equations in a 1-D free boundary domain. Plasmids are mobile genetic elements able to transfer to different phylotypes, posing a global health problem when they carry antibiotic resistance factors. We model gene transfer regulation influenced by nearby potential receptors to account for recipient-sensing. We also introduce a promotion function to account for trace metal effects on conjugation, based on literature data. The model qualitatively matches experimental results, showing that contaminants like toxic metals and antibiotics promote plasmid persistence by favoring plasmid carriers and stimulating conjugation. Even at higher contaminant concentrations inhibiting conjugation, plasmid spread persists by strongly inhibiting plasmid-free cells. The model also replicates higher plasmid density in biofilm’s most active regions.