Linking complex microbial interactions and dysbiosis through a disordered Lotka–Volterra model

  1. Jacopo Pasqualini
  2. Amos Maritan
  3. Andrea Rinaldo
  4. Sonia Facchin
  5. Edoardo Vincenzo Savarino
  6. Ada Altieri  Is a corresponding author
  7. Samir Suweis  Is a corresponding author
  1. Dipartimento di Fisica “G. Galilei” e INFN sezione di Padova, Università di Padova, Italy
  2. Dipartimento di Ingegneria Civile, Edile e Ambientale (ICEA), Università di Padova, Italy
  3. École Polytechnique Fédérale de Lausanne, Switzerland
  4. Dipartimento di Scienze Chirurgiche, Oncologiche e Gastroenterologiche (DiSCOG), Università di Padova, Italy
  5. Laboratoire Matière et Systèmes Complexes (MSC), Université Paris Cité, CNRS, France
  6. Padova Neuroscience Center, University of Padova, Italy
5 figures and 1 additional file

Figures

The inference protocol of the dgLV generative model is performed by a moment matching optimization procedure.

We aim to infer the free parameters θ=(μ,σ,T,λ) – shown on the right – so that to match the mean abundance, higher-order correlations between species abundances, and the average carrying capacity for the two cohorts, that is, (h,q0,qd,K) – on the left. This procedure enables us to extract relevant information about the ecological dynamics from cross-sectional data of healthy (blue) and diseased (red) microbiomes, which are treated as independent disordered realizations.

Order parameters inferred from the data using Equation 4.

(a) shows h, (b) qd, and (c) q0 in healthy (blue) and unhealthy (red) cohorts. Circles denote the mean value of each order parameter, and error bars indicate the standard deviation across bootstrap realizations. Bootstrap estimates were obtained from 5000 iterations, each retaining 90% of samples within each cohort. Gray symbols show the corresponding null-model values obtained by randomizing cohort labels prior to estimation.

Distinct ecological organization in healthy vs diseased microbiomes.

(a) Inferred T (demographic noise strength) and σ (interactions heterogeneity) for healthy (blue) and diseased (red) microbiomes are clustered. Darker dots correspond to better solutions (i.e., solutions with a lower value of the cost function C), while the two points with hexagonal markers correspond to the best two (healthy and diseased, respectively) solutions. In the first panel inset, we also show (in log–log scale) the species abundance distributions (SADs) corresponding to each solution. To have a more concise representation, we present each SAD fixing the disorder to its average ζ¯=Kμh. (b) The probability density function of the inferred interactions αi,j for healthy (blue) and diseased (red) microbiomes. Dysbiosis reduces the heterogeneity of the interaction strengths. The quantities reported in the legend are the average and standard deviation of αi,j. They are calculated as μXα=μX/SX and σXα=σX/SX, where SX is the species pool size, estimated as the set of all observed species in a dataset, X can denote healthy (H) or diseased (D) individuals.

Stability of healthy vs diseased microbiomes.

(a) The replicon eigenvalue corresponding to each solution of our optimization procedure (shaded dots). The solid hexagon represents the replicon corresponding to the best solutions that minimize the error in predicting the order parameters of the theory (minimum C). The two investigated microbiome phenotypes (healthy in blue, diseased in red) are significantly different. In particular, diseased microbiomes are closer to the marginal stability of replica-symmetric ansatz (gray horizontal line). (b) Solutions of the moment-matching objective function are shown as a function of ψ and m, which in turn depend on the species abundance distribution (SAD) parameters (see main text). Healthy (blue) and diseased (red) microbiomes appear to be clustered. Therefore, distinct ecological organization scenarios (strong neutrality/emergent neutrality) take place. Darker dots correspond to solutions with lower values of the cost function, while hexagonal markers correspond to the two best solutions.

Author response image 1

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

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)

  1. Jacopo Pasqualini
  2. Amos Maritan
  3. Andrea Rinaldo
  4. Sonia Facchin
  5. Edoardo Vincenzo Savarino
  6. Ada Altieri
  7. Samir Suweis
(2026)
Linking complex microbial interactions and dysbiosis through a disordered Lotka–Volterra model
eLife 14:RP105948.
https://doi.org/10.7554/eLife.105948.3