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.
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
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
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.
© 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.
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Glutamine synthetases (GS) are central enzymes essential for the nitrogen metabolism across all domains of life. Consequently, they have been extensively studied for more than half a century. Based on the ATP-dependent ammonium assimilation generating glutamine, GS expression and activity are strictly regulated in all organisms. In the methanogenic archaeon Methanosarcina mazei, it has been shown that the metabolite 2-oxoglutarate (2-OG) directly induces the GS activity. Besides, modulation of the activity by interaction with small proteins (GlnK1 and sP26) has been reported. Here, we show that the strong activation of M. mazei GS (GlnA1) by 2-OG is based on the 2-OG dependent dodecamer assembly of GlnA1 by using mass photometry (MP) and single particle cryo-electron microscopy (cryo-EM) analysis of purified strep-tagged GlnA1. The dodecamer assembly from dimers occurred without any detectable intermediate oligomeric state and was not affected in the presence of GlnK1. The 2.39 Å cryo-EM structure of the dodecameric complex in the presence of 12.5 mM 2-OG demonstrated that 2-OG is binding between two monomers. Thereby, 2-OG appears to induce the dodecameric assembly in a cooperative way. Furthermore, the active site is primed by an allosteric interaction cascade caused by 2-OG-binding towards an adaption of an open active state conformation. In the presence of additional glutamine, strong feedback inhibition of GS activity was observed. Since glutamine dependent disassembly of the dodecamer was excluded by MP, feedback inhibition most likely relies on the binding of glutamine to the catalytic site. Based on our findings, we propose that under nitrogen limitation the induction of M. mazei GS into a catalytically active dodecamer is not affected by GlnK1 and crucially depends on the presence of 2-OG.
The members of the Mycobacterium tuberculosis complex (MTBC) causing human tuberculosis comprise 10 phylogenetic lineages that differ in their geographical distribution. The human consequences of this phylogenetic diversity remain poorly understood. Here, we assessed the phenotypic properties at the host-pathogen interface of 14 clinical strains representing five major MTBC lineages. Using a human in vitro granuloma model combined with bacterial load assessment, microscopy, flow cytometry, and multiplexed-bead arrays, we observed considerable intra-lineage diversity. Yet, modern lineages were overall associated with increased growth rate and more pronounced granulomatous responses. MTBC lineages exhibited distinct propensities to accumulate triglyceride lipid droplets—a phenotype associated with dormancy—that was particularly pronounced in lineage 2 and reduced in lineage 3 strains. The most favorable granuloma responses were associated with strong CD4 and CD8 T cell activation as well as inflammatory responses mediated by CXCL9, granzyme B, and TNF. Both of which showed consistent negative correlation with bacterial proliferation across genetically distant MTBC strains of different lineages. Taken together, our data indicate that different virulence strategies and protective immune traits associate with MTBC genetic diversity at lineage and strain level.