Global dynamics of microbial communities emerge from local interaction rules

PLoS Comput Biol. 2022 Mar 4;18(3):e1009877. doi: 10.1371/journal.pcbi.1009877. eCollection 2022 Mar.

Abstract

Most microbes live in spatially structured communities (e.g., biofilms) in which they interact with their neighbors through the local exchange of diffusible molecules. To understand the functioning of these communities, it is essential to uncover how these local interactions shape community-level properties, such as the community composition, spatial arrangement, and growth rate. Here, we present a mathematical framework to derive community-level properties from the molecular mechanisms underlying the cell-cell interactions for systems consisting of two cell types. Our framework consists of two parts: a biophysical model to derive the local interaction rules (i.e. interaction range and strength) from the molecular parameters underlying the cell-cell interactions and a graph based model to derive the equilibrium properties of the community (i.e. composition, spatial arrangement, and growth rate) from these local interaction rules. Our framework shows that key molecular parameters underlying the cell-cell interactions (e.g., the uptake and leakage rates of molecules) determine community-level properties. We apply our model to mutualistic cross-feeding communities and show that spatial structure can be detrimental for these communities. Moreover, our model can qualitatively recapitulate the properties of an experimental microbial community. Our framework can be extended to a variety of systems of two interacting cell types, within and beyond the microbial world, and contributes to our understanding of how community-level properties emerge from microscopic interactions between cells.

Publication types

  • Research Support, Non-U.S. Gov't

MeSH terms

  • Biofilms
  • Microbial Interactions
  • Microbiota*
  • Symbiosis

Grants and funding

SVV was funded by the SNSF Postdoc Mobility fellowship nr. 175123, CH was funded by the Natural Sciences and Engineering Research Council of Canada (NSERC) grant nr. RGPIN-2015-05795 and RGPIN-2021-02608, MA and ADC were funded by SNSF grant nr. 31003A\_169978, Eawag, and ETH Zurich, and ADC and KF were funded by the Office of Naval Research grant nr. ONR N00014-17-1-3029. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.