An efficient moments-based inference method for within-host bacterial infection dynamics

The Infectious Dose Shapes Vibrio cholerae Within-Host Dynamics

During infection, the rates of pathogen replication, death, and migration affect disease progression, dissemination, transmission, and resistance evolution. Here, we follow the population dynamics of Vibrio cholerae in a mouse model by labeling individual bacteria with one of >500 unique, fitness-neutral genomic tags. Using the changes in tag frequencies and CFU numbers, we inform a mathematical model that describes the within-host spatiotemporal bacterial dynamics. This allows us to disentangle growth, death, forward, and retrograde migration rates continuously during infection. Our model has robust predictive power across various experimental setups. The population dynamics of V. cholerae shows substantial spatiotemporal heterogeneity in replication, death, and migration. Importantly, we find that the niche available to V. cholerae in the host increases with inoculum size, suggesting cooperative effects during infection. Therefore, it is not enough to consider just the likelihood of exposure (50% infectious dose) but rather the magnitude of exposure to predict outbreaks. IMPORTANCE Determining the rates of bacterial migration, replication, and death during infection is important for understanding how infections progress. Separately measuring these rates is often difficult in systems where multiple processes happen simultaneously. Here, we use next-generation sequencing to measure V. cholerae migration, replication, death, and niche size along the mouse gastrointestinal tract. We show that the small intestine of the mouse is a heterogeneous environment, and the population dynamic characteristics change substantially between adjacent gut sections. Our approach also allows us to characterize the effect of inoculum size on these processes. We find that the niche size in mice increases with the infectious dose, hinting at cooperative effects in larger inocula. The dose-response relationship between inoculum size and final pathogen burden is important for the infected individual and is thought to influence the progression of V. cholerae epidemics.

RESTAMP - Rate estimates by sequence-tag analysis of microbial populations

Microbial division rates determine the speed of mutation accumulation and thus the emergence of antimicrobial resistance. Microbial death rates are affected by antibiotic action and the immune system. Therefore, measuring these rates has advanced our understanding of host-pathogen interactions and antibiotic action. Several methods based on marker-loss or few inheritable neutral markers exist that allow estimating microbial division and death rates, each of which has advantages and limitations. Technical bottlenecks, i.e., experimental sampling events, during the experiment can distort the rate estimates and are typically unaccounted for or require additional calibration experiments. In this work, we introduce RESTAMP (Rate Estimates by Sequence Tag Analysis of Microbial Populations) as a method for determining bacterial division and death rates. This method uses hundreds of fitness neutral sequence barcodes to measure the rates and account for experimental bottlenecks at the same time. We experimentally validate RESTAMP and compare it to established plasmid loss methods. We find that RESTAMP has a number of advantages over plasmid loss or previous marker based techniques. (i) It enables to correct the distortion of rate estimates by technical bottlenecks. (ii) Rate estimates are independent of the sequence tag distribution in the starting culture allowing the use of an arbitrary number of tags. (iii) It introduces a bottleneck sensitivity measure that can be used to maximize the accuracy of the experiment. RESTAMP allows studying microbial population dynamics with great resolution over a wide dynamic range and can thus advance our understanding of host-pathogen interactions or the mechanisms of antibiotic action.

An experimental design tool to optimize inference precision in data-driven mathematical models of bacterial infections in vivo

The management of bacterial diseases calls for a detailed knowledge about the dynamic changes in host-bacteria interactions. Biological insights are gained by integrating experimental data with mechanistic mathematical models to infer experimentally unobservable quantities. This inter-disciplinary field would benefit from experiments with maximal information content yielding high-precision inference. Here, we present a computationally efficient tool for optimizing experimental design in terms of parameter inference in studies using isogenic-tagged strains. We study the effect of three experimental design factors: number of biological replicates, sampling timepoint selection and number of copies per tagged strain. We conduct a simulation study to establish the relationship between our optimality criterion and the size of parameter estimate confidence intervals, and showcase its application in a range of biological scenarios reflecting different dynamics patterns observed in experimental infections. We show that in low-variance systems with low killing and replication rates, predicting high-precision experimental designs is consistently achieved; higher replicate sizes and strategic timepoint selection yield more precise estimates. Finally, we address the question of resource allocation under constraints; given a fixed number of host animals and a constraint on total inoculum size per host, infections with fewer strains at higher copies per strain lead to higher-precision inference.

Within-host spatiotemporal dynamic of systemic salmonellosis: Ways to track infection, reaction to vaccination and antimicrobial treatment

During the last two decades our understanding of the complex in vivo host-pathogen interactions has increased due to technical improvements and new research tools. The rapid advancement of molecular biology, flow cytometry and microscopy techniques, combined with mathematical modelling, have empowered in-depth studies of systemic bacterial infections across scales from single molecules, to cells, to organs and systems to reach the whole organism level. By tracking subpopulations of bacteria in vivo using molecular or fluorescent tags, it has been possible to reconstruct the spread of infection within and between organs, allowing unprecedented quantification of the effects of antimicrobial treatment and vaccination. This review illustrates recent advances in the study of heterogeneous traits of the infection process and illustrate approaches to investigate the reciprocal interactions between antimicrobial treatments, bacterial growth/death as well as inter- and intra-organ spread. We also discuss how vaccines impact the in vivo behaviour of bacteria and how these findings can guide vaccine design and rational antimicrobial treatment.

A data-based mathematical modelling study to quantify the effects of ciprofloxacin and ampicillin on the within-host dynamics of Salmonella enterica during treatment and relapse

Antibiotic therapy has drastically reduced the mortality and sequelae of bacterial infections. From naturally occurring to chemically synthesized, different classes of antibiotics have been successfully used without detailed knowledge of how they affect bacterial dynamics in vivo. However, a proportion of patients receiving antimicrobial therapy develop recrudescent infections post-treatment. Relapsing infections are attributable to incomplete clearance of bacterial populations following antibiotic administration; the metabolic profile of this antibiotic-recalcitrant bacterial subpopulation, the spatio-temporal context of its emergence and the variance of antibiotic-bacterial interactions in vivo remain unclear. Here, we develop and apply a mechanistic mathematical model to data from a study comparing the effects of ciprofloxacin and ampicillin on the within-host dynamics of Salmonella enterica serovar Typhimurium in murine infections. Using the inferential capacity of our model, we show that the antibiotic-recalcitrant bacteria following ampicillin, but not ciprofloxacin, treatment belong to a non-replicating phenotype. Aligning with previous studies, we independently estimate that the lymphoid tissues and spleen are important reservoirs of non-replicating bacteria. Finally, we predict that post-treatment, the progenitors of the non-growing and growing bacterial populations replicate and die at different rates. Ultimately, the liver, spleen and mesenteric lymph nodes are all repopulated by progenitors of the previously non-growing phenotype in ampicillin-treated mice.

Integrating mathematical models with experimental data to investigate the within-host dynamics of bacterial infections

Bacterial infections still constitute a major cause of mortality and morbidity worldwide. The unavailability of therapeutics, antimicrobial resistance and the chronicity of infections due to incomplete clearance contribute to this phenomenon. Despite the progress in antimicrobial and vaccine development, knowledge about the effect that therapeutics have on the host–bacteria interactions remains incomplete. Insights into the characteristics of bacterial colonization and migration between tissues and the relationship between replication and host- or therapeutically induced killing can enable efficient design of treatment approaches. Recently, innovative experimental techniques have generated data enabling the qualitative characterization of aspects of bacterial dynamics. Here, we argue that mathematical modeling as an adjunct to experimental data can enrich the biological insight that these data provide. However, due to limited interdisciplinary training, efforts to combine the two remain limited. To promote this dialogue, we provide a categorization of modeling approaches highlighting their relationship to data generated by a range of experimental techniques in the area of in vivo bacterial dynamics. We outline common biological themes explored using mathematical models with case studies across all pathogen classes. Finally, this review advocates multidisciplinary integration to improve our mechanistic understanding of bacterial infections and guide the use of existing or new therapies.

Dose-response and transmission: the nexus between reservoir hosts, environment and recipient hosts

Dose is the nexus between exposure and all upstream processes that determine pathogen pressure, and is thereby an important element underlying disease dynamics. Understanding the relationship between dose and disease is particularly important in the context of spillover, where nonlinearities in the dose-response could determine the likelihood of transmission. There is a need to explore dose-response models for directly transmitted and zoonotic pathogens, and how these interactions integrate within-host factors to consider, for example, heterogeneity in host susceptibility and dose-dependent antagonism. Here, we review the dose-response literature and discuss the unique role dose-response models have to play in understanding and predicting spillover events. We present a re-analysis of dose-response experiments for two important zoonotic pathogens (Middle East respiratory syndrome coronavirus and Nipah virus), to exemplify potential difficulties in differentiating between appropriate models with small exposure experiment datasets. We also discuss the data requirements needed for robust selection between dose-response models. We then suggest how these processes could be modelled to gain more realistic predictions of zoonotic transmission outcomes and highlight the exciting opportunities that could arise with increased collaboration between the virology and epidemiology disciplines. This article is part of the theme issue 'Dynamic and integrative approaches to understanding pathogen spillover'.

Computational Health Engineering Applied to Model Infectious Diseases and Antimicrobial Resistance Spread

Infectious diseases are the primary cause of mortality worldwide. The dangers of infectious disease are compounded with antimicrobial resistance, which remains the greatest concern for human health. Although novel approaches are under investigation, the World Health Organization predicts that by 2050, septicaemia caused by antimicrobial resistant bacteria could result in 10 million deaths per year. One of the main challenges in medical microbiology is to develop novel experimental approaches, which enable a better understanding of bacterial infections and antimicrobial resistance. After the introduction of whole genome sequencing, there was a great improvement in bacterial detection and identification, which also enabled the characterization of virulence factors and antimicrobial resistance genes. Today, the use of in silico experiments jointly with computational and machine learning offer an in depth understanding of systems biology, allowing us to use this knowledge for the prevention, prediction, and control of infectious disease. Herein, the aim of this review is to discuss the latest advances in human health engineering and their applicability in the control of infectious diseases. An in-depth knowledge of host–pathogen–protein interactions, combined with a better understanding of a host’s immune response and bacterial fitness, are key determinants for halting infectious diseases and antimicrobial resistance dissemination.

Death and population dynamics affect mutation rate estimates and evolvability under stress in bacteria

The stress-induced mutagenesis hypothesis postulates that in response to stress, bacteria increase their genome-wide mutation rate, in turn increasing the chances that a descendant is able to better withstand the stress. This has implications for antibiotic treatment: exposure to subinhibitory doses of antibiotics has been reported to increase bacterial mutation rates and thus probably the rate at which resistance mutations appear and lead to treatment failure. More generally, the hypothesis posits that stress increases evolvability (the ability of a population to generate adaptive genetic diversity) and thus accelerates evolution. Measuring mutation rates under stress, however, is problematic, because existing methods assume there is no death. Yet subinhibitory stress levels may induce a substantial death rate. Death events need to be compensated by extra replication to reach a given population size, thus providing more opportunities to acquire mutations. We show that ignoring death leads to a systematic overestimation of mutation rates under stress. We developed a system based on plasmid segregation that allows us to measure death and division rates simultaneously in bacterial populations. Using this system, we found that a substantial death rate occurs at the tested subinhibitory concentrations previously reported to increase mutation rate. Taking this death rate into account lowers and sometimes removes the signal for stress-induced mutagenesis. Moreover, even when antibiotics increase mutation rate, we show that subinhibitory treatments do not increase genetic diversity and evolvability, again because of effects of the antibiotics on population dynamics. We conclude that antibiotic-induced mutagenesis is overestimated because of death and that understanding evolvability under stress requires accounting for the effects of stress on population dynamics as much as on mutation rate. Our goal here is dual: we show that population dynamics and, in particular, the numbers of cell divisions are crucial but neglected parameters in the evolvability of a population, and we provide experimental and computational tools and methods to study evolvability under stress, leading to a reassessment of the magnitude and significance of the stress-induced mutagenesis paradigm.

The effect of environmental stress on bacterial mutagenesis has been a paradigm-shift discovery. Recent developments include evidence that various antibiotics increase mutation rates in bacteria when used at subinhibitory concentrations. It is therefore suggested that such treatments promote resistance evolution because they increase the generation of genetic variation on which natural selection can act. However, existing methods to compute mutation rate neglect the effect of stress on death and population dynamics. Developing new experimental and computational tools, we find that taking death into account significantly lowers the signal for stress-induced mutagenesis. Moreover, we show that treatments that increase mutation rate do not always lead to increased genetic diversity, which questions the standard paradigm of increased evolvability under stress.