Abstract:Exposure to environmental stressors, including certain antibiotics, induces stress responses in bacteria. Some of these responses increase mutagenesis and thus potentially accelerate resistance evolution. Many studies report increased mutation rates under stress, often using the standard experimental approach of fluctuation assays. However, single-cell studies have revealed that many stress responses are heterogeneously expressed in bacterial populations, which existing estimation methods have not yet addressed. We develop a population dynamic model that considers heterogeneous stress responses (subpopulations of cells with the response off or on ) that impact both mutation rate and cell division rate, inspired by the DNA-damage response in Escherichia coli (SOS response). We derive the mutant count distribution arising in fluctuation assays under this model and then implement maximum likelihood estimation of the mutation-rate increase specifically associated with the expression of the stress response. Using simulated mutant count data, we show that our inference method allows for accurate and precise estimation of the mutation-rate increase, provided that this increase is sufficiently large and the induction of the response also reduces the division rate. Moreover, we find that in many cases, either heterogeneity in stress responses or mutant fitness costs could explain similar patterns in fluctuation assay data, suggesting that separate experiments would be required to identify the true underlying process. In cases where stress responses and mutation rates are heterogeneous, current methods still correctly infer the effective increase in population mean mutation rate, but we provide a novel method to infer distinct stress-induced mutation rates, which could be important for parameterising evolutionary models. How does environmental stress, especially from antibiotics, affect mutation rates in bacteria? This question has often been examined by estimating mutation rates using fluctuation assays, an experiment dating back to Luria and Delbrück in the 1940s. In this study, we consider variation in stress responses within bacterial populations, as revealed by recent single-cell studies, which is neglected in currently available mutation-rate estimation methods. Our approach involves a population dynamic model inspired by the DNA-damage response in E. coli (SOS response). It accounts for a subpopulation with high expression of the stress response, which increases the mutation rate and decreases the division rate of a cell. We use computer simulations to generate synthetic fluctuation assay data. Notably, we find that over a wide range of scenarios, existing models and our heterogeneous-response model cannot be distinguished using fluctuation assay data alone. This emphasises the need for separate experiments to uncover the true underlying processes. Nevertheless, when stress responses are known to be heterogeneous, our study offers a novel method for accurately estimating mutation rates specifically associated with the high expression of the stress response. Uncovering the heterogeneity in stress-induced mutation rates could be important for predicting the evolution of antibiotic resistance.

Common regulatory mutation increases single-cell survival to antibiotic exposures in .

Evaluating the predictive power of combined gene expression dynamics from single cells on antibiotic survival

Mutations causing low level antibiotic resistance ensure bacterial survival in antibiotic-treated hosts

Phenotypic resistant single-cell characteristics under recurring ampicillin antibiotic exposure in Escherichia coli

Heteroresistance at the Single-Cell Level: Adapting to Antibiotic Stress Through a Population-Based Strategy and Growth-Controlled Interphenotypic Coordination

Measuring single-cell susceptibility to antibiotics within monoclonal bacterial populations

Dynamical model of antibiotic responses linking expression of resistance genes to metabolism explains emergence of heterogeneity during drug exposures

Nonoptimal Gene Expression Creates Latent Potential for Antibiotic Resistance

Antibiotic Treatment Enhances the Genome-Wide Mutation Rate of Target Cells.

Studies Of The Drug Resistance Response Of Sensitive And Drug-Resistant Strains In A Microfluidic System

Gene expression tradeoffs determine bacterial survival and adaptation to antibiotic stress

ppGpp and RNA-polymerase backtracking guide antibiotic-induced mutable gambler cells

Antibiotic resistant bacteria survive treatment by doubling while shrinking

Low-Level Cefepime Exposure Induces High-Level Resistance in Environmental Bacteria: Molecular Mechanism and Evolutionary Dynamics.

History-Dependent Physiological Adaptation to Lethal Genetic Modification under Antibiotic Exposure

Estimating mutation rates under heterogeneous stress responses

Phenotypic Plasticity of Antibiotic Resistance, Metabolism Byproduct Utilization and the Evolution of Mutually Beneficial Cooperation in Escherichia Coli

Single-cell Microfluidics Enabled Dynamic Evaluation of Drug Combinations on Antibiotic Resistance Bacteria

Integrating CRISPR-Enabled Trackable Genome Engineering and Transcriptomic Analysis of Global Regulators for Antibiotic Resistance Selection and Identification in Escherichia coli

Dynamics of drug delivery determines course of evolution of antibiotic responses in bacteria

Distinct Evolutionary Origins of Common Multi-Drug Resistance Phenotypes in Salmonella Typhimurium DT104: a Convergent Process for Adaptation under Stress.