venturelli3

None

None

None

None

None

None

Deciphering microbial interactions in synthetic human gut microbiome communities.

  • Ophelia S Venturelli
  • Alex C Carr
  • Garth Fisher
  • Ryan H Hsu
  • Rebecca Lau
  • Benjamin P Bowen
  • Susan Hromada
  • Trent Northen
  • Adam P Arkin
Mol. Syst. Biol. 2018; 14 (6):
Abstract
The ecological forces that govern the assembly and stability of the human gut microbiota remain unresolved. We developed a generalizable model-guided framework to predict higher-dimensional consortia from time-resolved measurements of lower-order assemblages. This method was employed to decipher microbial interactions in a diverse human gut microbiome synthetic community. We show that pairwise interactions are major drivers of multi-species community dynamics, as opposed to higher-order interactions. The inferred ecological network exhibits a high proportion of negative and frequent positive interactions. Ecological drivers and responsive recipient species were discovered in the network. Our model demonstrated that a prevalent positive and negative interaction topology enables robust coexistence by implementing a negative feedback loop that balances disparities in monospecies fitness levels. We show that negative interactions could generate history-dependent responses of initial species proportions that frequently do not originate from bistability. Measurements of extracellular metabolites illuminated the metabolic capabilities of monospecies and potential molecular basis of microbial interactions. In sum, these methods defined the ecological roles of major human-associated intestinal species and illuminated design principles of microbial communities.

Unit definitions have no effect on the numerical analysis of the model. It remains the responsibility of the modeler to ensure the internal numerical consistency of the model. If units are provided, however, the consistency of the model units will be checked.

Name Definition
Id Name Spatial dimensions Size
default 0.0 1.0
Id Name Initial quantity Compartment Fixed
x1 1.0 default
x10 1.0 default
x11 1.0 default
x12 1.0 default
x2 1.0 default
x3 1.0 default
x4 1.0 default
x5 1.0 default
x6 1.0 default
x7 1.0 default
x8 1.0 default
x9 1.0 default
xBH 0.01 default
xBO 0.01 default
xBT 0.01 default
xBU 0.01 default
xBV 0.01 default
xCA 0.01 default
xCH 0.01 default
xDP 0.01 default
xEL 0.01 default
xER 0.01 default
xFP 0.01 default
xPC 0.01 default

Initial assignments are expressions that are evaluated at time=0. It is not recommended to create initial assignments for all model entities. Restrict the use of initial assignments to cases where a value is expressed in terms of values or sizes of other model entities. Note that it is not permitted to have both an initial assignment and an assignment rule for a single model entity.

Definition
Id Name Objective coefficient Reaction Equation and Kinetic Law Flux bounds
v1 x1 = xBH

xBH * (uBH + xBH * aBHBH + xCA * aBHCA + xBU * aBHBU + xPC * aBHPC + xBO * aBHBO + xBV * aBHBV + xBT * aBHBT + xEL * aBHEL + xFP * aBHFP + xCH * aBHCH + xDP * aBHDP + xER * aBHER)
v10 x10 = xCH

xCH * (uCH + xBH * aCHBH + xCA * aCHCA + xBU * aCHBU + xPC * aCHPC + xBO * aCHBO + xBV * aCHBV + xBT * aCHBT + xEL * aCHEL + xFP * aCHFP + xCH * aCHCH + xDP * aCHDP + xER * aCHER)
v11 x11 = xDP

xDP * (uDP + xBH * aDPBH + xCA * aDPCA + xBU * aDPBU + xPC * aDPPC + xBO * aDPBO + xBV * aDPBV + xBT * aDPBT + xEL * aDPEL + xFP * aDPFP + xCH * aDPCH + xDP * aDPDP + xER * aDPER)
v12 x12 = xER

xER * (uER + xBH * aERBH + xCA * aERCA + xBU * aERBU + xPC * aERPC + xBO * aERBO + xBV * aERBV + xBT * aERBT + xEL * aEREL + xFP * aERFP + xCH * aERCH + xDP * aERDP + xER * aERER)
v2 x2 = xCA

xCA * (uCA + xBH * aCABH + xCA * aCACA + xBU * aCABU + xPC * aCAPC + xBO * aCABO + xBV * aCABV + xBT * aCABT + xEL * aCAEL + xFP * aCAFP + xCH * aCACH + xDP * aCADP + xER * aCAER)
v3 x3 = xBU

xBU * (uBU + xBH * aBUBH + xCA * aBUCA + xBU * aBUBU + xPC * aBUPC + xBO * aBUBO + xBV * aBUBV + xBT * aBUBT + xEL * aBUEL + xFP * aBUFP + xCH * aBUCH + xDP * aBUDP + xER * aBUER)
v4 x4 = xPC

xPC * (uPC + xBH * aPCBH + xCA * aPCCA + xBU * aPCBU + xPC * aPCPC + xBO * aPCBO + xBV * aPCBV + xBT * aPCBT + xEL * aPCEL + xFP * aPCFP + xCH * aPCCH + xDP * aPCDP + xER * aPCER)
v5 x5 = xBO

xBO * (uBO + xBH * aBOBH + xCA * aBOCA + xBU * aBOBU + xPC * aBOPC + xBO * aBOBO + xBV * aBOBV + xBT * aBOBT + xEL * aBOEL + xFP * aBOFP + xCH * aBOCH + xDP * aBODP + xER * aBOER)
v6 x6 = xBV

xBV * (uBV + xBH * aBVBH + xCA * aBVCA + xBU * aBVBU + xPC * aBVPC + xBO * aBVBO + xBV * aBVBV + xBT * aBVBT + xEL * aBVEL + xFP * aBVFP + xCH * aBVCH + xDP * aBVDP + xER * aBVER)
v7 x7 = xBT

xBT * (uBT + xBH * aBTBH + xCA * aBTCA + xBU * aBTBU + xPC * aBTPC + xBO * aBTBO + xBV * aBTBV + xBT * aBTBT + xEL * aBTEL + xFP * aBTFP + xCH * aBTCH + xDP * aBTDP + xER * aBTER)
v8 x8 = xEL

xEL * (uEL + xBH * aELBH + xCA * aELCA + xBU * aELBU + xPC * aELPC + xBO * aELBO + xBV * aELBV + xBT * aELBT + xEL * aELEL + xFP * aELFP + xCH * aELCH + xDP * aELDP + xER * aELER)
v9 x9 = xFP

xFP * (uFP + xBH * aFPBH + xCA * aFPCA + xBU * aFPBU + xPC * aFPPC + xBO * aFPBO + xBV * aFPBV + xBT * aFPBT + xEL * aFPEL + xFP * aFPFP + xCH * aFPCH + xDP * aFPDP + xER * aFPER)

Global parameters

Id Value
aBHBH -0.9241
aBHBO -0.2199
aBHBT -0.2641
aBHBU -0.2284
aBHBV -0.131
aBHCA -0.312
aBHCH -0.3284
aBHDP -0.9182
aBHEL 0.1821
aBHER -0.5488
aBHFP -0.2367
aBHPC -0.5294
aBOBH 0.0
aBOBO -0.7134
aBOBT -0.5769
aBOBU -0.5064
aBOBV -0.5089
aBOCA 0.0
aBOCH -0.4668
aBODP 0.0
aBOEL 1.8096
aBOER -0.0391
aBOFP -0.1758
aBOPC -0.0001
aBTBH 0.0582
aBTBO -0.8237
aBTBT -0.9597
aBTBU -0.972
aBTBV -0.6819
aBTCA 0.0
aBTCH -0.0727
aBTDP 0.0
aBTEL 1.9668
aBTER 0.0
aBTFP -0.4498
aBTPC 0.0
aBUBH 0.0
aBUBO -0.7487
aBUBT -0.9377
aBUBU -0.9059
aBUBV -0.5581
aBUCA 0.0
aBUCH 0.0564
aBUDP 0.0
aBUEL 3.0816
aBUER 0.0
aBUFP -0.8309
aBUPC 0.0
aBVBH 0.1337
aBVBO -0.5768
aBVBT -0.6133
aBVBU -0.5773
aBVBV -0.6578
aBVCA -0.654
aBVCH -0.0464
aBVDP -0.108
aBVEL 1.2848
aBVER 0.0
aBVFP -0.6349
aBVPC 0.0
aCABH 0.4301
aCABO -0.283
aCABT -0.2737
aCABU -0.2612
aCABV -0.1701
aCACA -0.818
aCACH 0.3033
aCADP 0.0
aCAEL -0.4766
aCAER 0.0
aCAFP -1.1212
aCAPC -0.685
aCHBH -0.3846
aCHBO -0.4679
aCHBT -0.5906
aCHBU -0.2163
aCHBV 0.0
aCHCA -0.3186
aCHCH -1.242
aCHDP -2.3006
aCHEL 0.0
aCHER -0.5077
aCHFP -0.6943
aCHPC -0.7845
aDPBH 0.0
aDPBO -0.187
aDPBT -0.1915
aDPBU -0.1628
aDPBV -0.0172
aDPCA 0.0532
aDPCH -0.1472
aDPDP -1.3228
aDPEL 0.1989
aDPER 0.0
aDPFP 0.0
aDPPC -0.4407
aELBH 0.8285
aELBO -0.0932
aELBT -0.099
aELBU -0.0776
aELBV -0.0229
aELCA -0.9925
aELCH 0.0
aELDP -0.0003
aELEL -2.6202
aELER 0.0
aELFP 0.0
aELPC -1.091
aERBH 1.4712
aERBO 0.0
aERBT -0.0315
aERBU -0.0285
aERBV -0.0178
aERCA 0.0
aERCH 1.1905
aERDP 0.0
aEREL 0.0
aERER -1.3219
aERFP -0.0319
aERPC 0.0
aFPBH 0.9091
aFPBO -0.114
aFPBT -0.0742
aFPBU 0.1988
aFPBV 0.7253
aFPCA 0.0
aFPCH 0.4306
aFPDP 0.9762
aFPEL -0.8005
aFPER 0.0
aFPFP -0.9946
aFPPC -0.423
aPCBH 0.0
aPCBO -0.2718
aPCBT -0.3041
aPCBU -0.3241
aPCBV -0.2039
aPCCA -0.5827
aPCCH 0.2406
aPCDP -1.0009
aPCEL -0.9852
aPCER -0.9156
aPCFP -0.4513
aPCPC -0.6154
fBH 0.0
fBO 0.0
fBT 0.0
fBU 0.0
fBV 0.0
fCA 0.0
fCH 0.0
fDP 0.0
fEL 0.0
fER 0.0
fFP 0.0
fPC 0.0
uBH 0.2472
uBO 0.4624
uBT 0.6259
uBU 0.5991
uBV 0.4562
uCA 0.2464
uCH 0.4682
uDP 0.2361
uEL 0.4052
uER 0.1507
uFP 0.2164
uPC 0.2381

Local parameters

Id Value Reaction

Assignment rules

Definition
fEL = xEL / (xBH + xCA + xBU + xPC + xBO + xBV + xBT + xEL + xFP + xCH + xDP + xER)
fFP = xFP / (xBH + xCA + xBU + xPC + xBO + xBV + xBT + xEL + xFP + xCH + xDP + xER)
fDP = xDP / (xBH + xCA + xBU + xPC + xBO + xBV + xBT + xEL + xFP + xCH + xDP + xER)
fER = xER / (xBH + xCA + xBU + xPC + xBO + xBV + xBT + xEL + xFP + xCH + xDP + xER)
fBO = xBO / (xBH + xCA + xBU + xPC + xBO + xBV + xBT + xEL + xFP + xCH + xDP + xER)
fCH = xCH / (xBH + xCA + xBU + xPC + xBO + xBV + xBT + xEL + xFP + xCH + xDP + xER)
fBT = xBT / (xBH + xCA + xBU + xPC + xBO + xBV + xBT + xEL + xFP + xCH + xDP + xER)
fBH = xBH / (xBH + xCA + xBU + xPC + xBO + xBV + xBT + xEL + xFP + xCH + xDP + xER)
fBU = xBU / (xBH + xCA + xBU + xPC + xBO + xBV + xBT + xEL + xFP + xCH + xDP + xER)
fBV = xBV / (xBH + xCA + xBU + xPC + xBO + xBV + xBT + xEL + xFP + xCH + xDP + xER)
fPC = xPC / (xBH + xCA + xBU + xPC + xBO + xBV + xBT + xEL + xFP + xCH + xDP + xER)
fCA = xCA / (xBH + xCA + xBU + xPC + xBO + xBV + xBT + xEL + xFP + xCH + xDP + xER)

Rate rules

Definition

Algebraic rules

Definition
Trigger Assignments
gt(time, 24) xER = 0.05 * xER; xDP = 0.05 * xDP; xCH = 0.05 * xCH; xEL = 0.05 * xEL; xBT = 0.05 * xBT; xBV = 0.05 * xBV; xPC = 0.05 * xPC; xBU = 0.05 * xBU; xBH = 0.05 * xBH; xCA = 0.05 * xCA; xBO = 0.05 * xBO; xFP = 0.05 * xFP
gt(time, 48) xER = 0.05 * xER; xDP = 0.05 * xDP; xFP = 0.05 * xFP; xEL = 0.05 * xEL; xBT = 0.05 * xBT; xBO = 0.05 * xBO; xPC = 0.05 * xPC; xCA = 0.05 * xCA; xBH = 0.05 * xBH; xBU = 0.05 * xBU; xBV = 0.05 * xBV; xCH = 0.05 * xCH