Gut bacterial aggregates as living gels
Abstract
The spatial organization of gut microbiota influences both microbial abundances and host-microbe interactions, but the underlying rules relating bacterial dynamics to large-scale structure remain unclear. To this end we studied experimentally and theoretically the formation of three-dimensional bacterial clusters, a key parameter controlling susceptibility to intestinal transport and access to the epithelium. Inspired by models of structure formation in soft materials, we sought to understand how the distribution of gut bacterial cluster sizes emerges from bacterial-scale kinetics. Analyzing imaging-derived data on cluster sizes for eight different bacterial strains in the larval zebrafish gut, we find a common family of size distributions that decay approximately as power laws with exponents close to -2, becoming shallower for large clusters in a strain-dependent manner. We show that this type of distribution arises naturally from a Yule-Simons-type process in which bacteria grow within clusters and can escape from them, coupled to an aggregation process that tends to condense the system toward a single massive cluster, reminiscent of gel formation. Together, these results point to the existence of general, biophysical principles governing the spatial organization of the gut microbiome that may be useful for inferring fast-timescale dynamics that are experimentally inaccessible.
Data availability
A table of all bacterial cluster sizes analysed in this study is included in the Supplementary Data File. MATLAB code for simulating the models described in the study is available at https://github.com/rplab/cluster_kinetics
Article and author information
Author details
Funding
National Institutes of Health (P50GM09891)
- Brandon H Schlomann
- Raghuveer Parthasarathy
National Institutes of Health (P01GM125576)
- Brandon H Schlomann
- Raghuveer Parthasarathy
National Institutes of Health (F32AI112094)
- Brandon H Schlomann
- Raghuveer Parthasarathy
National Institutes of Health (T32GM007759)
- Raghuveer Parthasarathy
National Science Foundation (1427957)
- Brandon H Schlomann
- Raghuveer Parthasarathy
James S. McDonnell Foundation
- Brandon H Schlomann
Kavli Foundation (Kavli Microbiome Ideas Challenge)
- Brandon H Schlomann
- Raghuveer Parthasarathy
National Institutes of Health (P01HD22486)
- Brandon H Schlomann
- Raghuveer Parthasarathy
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Ethics
Animal experimentation: The studies that generated the data analyzed in this paper (see cited references) were done in strict accordance with protocols approved by the University of Oregon Institutional Animal Care and Use Committee and following standard protocols.
Copyright
© 2021, Schlomann & Parthasarathy
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,761
- views
-
- 261
- downloads
-
- 15
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Microbiology and Infectious Disease
The Staphylococcus aureus clonal complex 8 (CC8) is made up of several subtypes with varying levels of clinical burden; from community-associated methicillin-resistant S. aureus USA300 strains to hospital-associated (HA-MRSA) USA500 strains and ancestral methicillin-susceptible (MSSA) strains. This phenotypic distribution within a single clonal complex makes CC8 an ideal clade to study the emergence of mutations important for antibiotic resistance and community spread. Gene-level analysis comparing USA300 against MSSA and HA-MRSA strains have revealed key horizontally acquired genes important for its rapid spread in the community. However, efforts to define the contributions of point mutations and indels have been confounded by strong linkage disequilibrium resulting from clonal propagation. To break down this confounding effect, we combined genetic association testing with a model of the transcriptional regulatory network (TRN) to find candidate mutations that may have led to changes in gene regulation. First, we used a De Bruijn graph genome-wide association study to enrich mutations unique to the USA300 lineages within CC8. Next, we reconstructed the TRN by using independent component analysis on 670 RNA-sequencing samples from USA300 and non-USA300 CC8 strains which predicted several genes with strain-specific altered expression patterns. Examination of the regulatory region of one of the genes enriched by both approaches, isdH, revealed a 38-bp deletion containing a Fur-binding site and a conserved single-nucleotide polymorphism which likely led to the altered expression levels in USA300 strains. Taken together, our results demonstrate the utility of reconstructed TRNs to address the limits of genetic approaches when studying emerging pathogenic strains.
-
- Immunology and Inflammation
- Microbiology and Infectious Disease
Pseudomonas aeruginosa (PA) is an opportunistic, frequently multidrug-resistant pathogen that can cause severe infections in hospitalized patients. Antibodies against the PA virulence factor, PcrV, protect from death and disease in a variety of animal models. However, clinical trials of PcrV-binding antibody-based products have thus far failed to demonstrate benefit. Prior candidates were derivations of antibodies identified using protein-immunized animal systems and required extensive engineering to optimize binding and/or reduce immunogenicity. Of note, PA infections are common in people with cystic fibrosis (pwCF), who are generally believed to mount normal adaptive immune responses. Here, we utilized a tetramer reagent to detect and isolate PcrV-specific B cells in pwCF and, via single-cell sorting and paired-chain sequencing, identified the B cell receptor (BCR) variable region sequences that confer PcrV-specificity. We derived multiple high affinity anti-PcrV monoclonal antibodies (mAbs) from PcrV-specific B cells across three donors, including mAbs that exhibit potent anti-PA activity in a murine pneumonia model. This robust strategy for mAb discovery expands what is known about PA-specific B cells in pwCF and yields novel mAbs with potential for future clinical use.