Stochastic microbiome assembly depends on context.

Observational studies reveal substantial variability in microbiome composition across individuals. Targeted studies in gnotobiotic animals underscore this variability by showing that some bacterial strains colonize deterministically, while others colonize stochastically. While some of this variability can be explained by external factors like environmental, dietary, and genetic differences between individuals, in this paper we show that for the model organism Drosophila melanogaster, interactions between bacteria can affect the microbiome assembly process, contributing to a baseline level of microbiome variability even among isogenic organisms that are identically reared, housed, and fed. In germ-free flies fed known combinations of bacterial species, we find that some species colonize more frequently than others even when fed at the same high concentration. We develop an ecological technique that infers the presence of interactions between bacterial species based on their colonization odds in different contexts, requiring only presence/absence data from two-species experiments. We use a progressive sequence of probabilistic models, in which the colonization of each bacterial species is treated as an independent stochastic process, to reproduce the empirical distributions of colonization outcomes across experiments. We find that incorporating context-dependent interactions substantially improves the performance of the models. Stochastic, context-dependent microbiome assembly underlies clinical therapies like fecal microbiota transplantation and probiotic administration and should inform the design of synthetic fecal transplants and dosing regimes.

Microbiome interactions shape host fitness.

Gut bacteria can affect key aspects of host fitness, such as development, fecundity, and lifespan, while the host, in turn, shapes the gut microbiome. However, it is unclear to what extent individual species versus community interactions within the microbiome are linked to host fitness. Here, we combinatorially dissect the natural microbiome of Drosophila melanogaster and reveal that interactions between bacteria shape host fitness through life history tradeoffs. Empirically, we made germ-free flies colonized with each possible combination of the five core species of fly gut bacteria. We measured the resulting bacterial community abundances and fly fitness traits, including development, reproduction, and lifespan. The fly gut promoted bacterial diversity, which, in turn, accelerated development, reproduction, and aging: Flies that reproduced more died sooner. From these measurements, we calculated the impact of bacterial interactions on fly fitness by adapting the mathematics of genetic epistasis to the microbiome. Development and fecundity converged with higher diversity, suggesting minimal dependence on interactions. However, host lifespan and microbiome abundances were highly dependent on interactions between bacterial species. Higher-order interactions (involving three, four, and five species) occurred in 13-44% of possible cases depending on the trait, with the same interactions affecting multiple traits, a reflection of the life history tradeoff. Overall, we found these interactions were frequently context-dependent and often had the same magnitude as individual species themselves, indicating that the interactions can be as important as the individual species in gut microbiomes.

In silico analysis of antibiotic-induced Clostridium difficile infection: Remediation techniques and biological adaptations.

In this paper we study antibiotic-induced C. difficile infection (CDI), caused by the toxin-producing C. difficile (CD), and implement clinically-inspired simulated treatments in a computational framework that synthesizes a generalized Lotka-Volterra (gLV) model with SIR modeling techniques. The gLV model uses parameters derived from an experimental mouse model, in which the mice are administered antibiotics and subsequently dosed with CD. We numerically identify which of the experimentally measured initial conditions are vulnerable to CD colonization, then formalize the notion of CD susceptibility analytically. We simulate fecal transplantation, a clinically successful treatment for CDI, and discover that both the transplant timing and transplant donor are relevant to the the efficacy of the treatment, a result which has clinical implications. We incorporate two nongeneric yet dangerous attributes of CD into the gLV model, sporulation and antibiotic-resistant mutation, and for each identify relevant SIR techniques that describe the desired attribute. Finally, we rely on the results of our framework to analyze an experimental study of fecal transplants in mice, and are able to explain observed experimental results, validate our simulated results, and suggest model-motivated experiments.

Pre-Guideline Trends Associated with Chest Radiograph Usage for Pediatric Community-Acquired Pneumonia in Emergency Department Settings.

BACKGROUND: Health care spending in the U.S. totaled $3 trillion in 2014 and continues to increase rapidly. Minimizing waste through clinical guidelines is a promising strategy to reduce spending without compromising patient care. In 2011, clinical guidelines recommended against the use of chest X-ray (CXR) for diagnosis of community-acquired pneumonia (CAP) in pediatric ambulatory settings. However, use of CXR has not changed post-guideline. Thus, understanding the drivers of CXR utilization prior to guideline implementation could improve guideline adherence. METHODS: Retrospective study using 2009 Nationwide Emergency Department Sample data set consisting of a representative sample of all emergency room admissions. Inclusion criteria consisted of: 18 years of age or younger and the diagnosis of outpatient CAP. Population was segmented by the presence of a CXR obtained during the visit. Socioeconomic status was determined by quartile classification of the estimated median household income based on patient ZIP code. RESULTS: In 2009, children living in wealthier ZIP codes presenting to the emergency department (ED) who were diagnosed with CAP were more likely to receive diagnostic CXR. The use of chest radiograph was not statistically correlated to gender, weekday versus weekend admission, number of diagnoses at discharge, or total ED charges. CONCLUSION: The research demonstrates a strong correlation between socioeconomic status of the pediatric patient and use of chest radiograph for CAP in the ED setting prior to 2011 guideline publication. Further research to determine the reason for this correlation could give rise to focused efforts to successfully encourage adherence to clinical practice guidelines.

First-passage-time statistics of growing microbial populations carry an imprint of initial conditions.

In exponential population growth, variability in the timing of individual division events and environmental factors (including stochastic inoculation) compound to produce variable growth trajectories. In several stochastic models of exponential growth we show power-law relationships that relate variability in the time required to reach a threshold population size to growth rate and inoculum size. Population-growth experiments in E. coli and S. aureus with inoculum sizes ranging between 1 and 100 are consistent with these relationships. We quantify how noise accumulates over time, finding that it encodes-and can be used to deduce-information about the early growth rate of a population.

A symbiotic physical niche in Drosophila melanogaster regulates stable association of a multi-species gut microbiota.

The gut is continuously invaded by diverse bacteria from the diet and the environment, yet microbiome composition is relatively stable over time for host species ranging from mammals to insects, suggesting host-specific factors may selectively maintain key species of bacteria. To investigate host specificity, we used gnotobiotic Drosophila, microbial pulse-chase protocols, and microscopy to investigate the stability of different strains of bacteria in the fly gut. We show that a host-constructed physical niche in the foregut selectively binds bacteria with strain-level specificity, stabilizing their colonization. Primary colonizers saturate the niche and exclude secondary colonizers of the same strain, but initial colonization by Lactobacillus species physically remodels the niche through production of a glycan-rich secretion to favor secondary colonization by unrelated commensals in the Acetobacter genus. Our results provide a mechanistic framework for understanding the establishment and stability of a multi-species intestinal microbiome.

Control of ecological outcomes through deliberate parameter changes in a model of the gut microbiome.

The generalized Lotka-Volterra (gLV) equations are a mathematical proxy for ecological dynamics. We focus on a gLV model of the gut microbiome, in which the evolution of the gut microbial state is determined in part by pairwise interspecies interaction parameters that encode environmentally mediated resource competition between microbes. We develop an in silico method that controls the steady-state outcome of the system by adjusting these interaction parameters. This approach is confined to a bistable region of the gLV model. In this method, a dimensionality reduction technique called steady-state reduction (SSR) is first used to generate a two-dimensional (2D) gLV model that approximates the high-dimensional dynamics on the 2D subspace spanned by the two steady states. Then a bifurcation analysis of the 2D model analytically determines parameter modifications that drive an initial condition to a target steady state. This parameter modification of the reduced 2D model guides parameter modifications of the original high-dimensional model, resulting in a change of steady-state outcome in the high-dimensional model. This control method, called SSR-guided parameter change (SPARC), bypasses the computational challenge of directly determining parameter modifications in the original high-dimensional system. SPARC could guide the development of indirect bacteriotherapies, which seek to change microbial compositions by deliberately modifying gut environmental variables such as gut acidity or macronutrient availability.

Steady-state reduction of generalized Lotka-Volterra systems in the microbiome.

The generalized Lotka-Volterra (gLV) equations, a classic model from theoretical ecology, describe the population dynamics of a set of interacting species. As the number of species in these systems grow in number, their dynamics become increasingly complex and intractable. We introduce steady-state reduction (SSR), a method that reduces a gLV system of many ecological species into two-dimensional subsystems that each obey gLV dynamics and whose basis vectors are steady states of the high-dimensional model. We apply this method to an experimentally-derived model of the gut microbiome in order to observe the transition between "healthy" and "diseased" microbial states. Specifically, we use SSR to investigate how fecal microbiota transplantation, a promising clinical treatment for dysbiosis, can revert a diseased microbial state to health.