In our design, we hope EcN will eventually colonize the ileum. To verify that the modified EcN can form an independent ecological niche at the ileum after erythritol consumption, we simulated the co-culture of the modified EcN with 24 species of ileal microbiota based on the dFBA algorithm, and designed a recipe for the average Chinese population using the Scientific Research Report on the Chinese Dietary Guidelines (2021). We simulated the intestinal environment of the average Chinese population based on the recipe, predicting the growth of multiple intestinal microbiota in co-culture, and verified that the modified EcN has better competitiveness in the microbiota. Meanwhile, the framework we developed can also be applied to other co-culture simulations, which can be promoted as a toolbox.
GEMs of 773 microbiota species in the gut have been constructed in the literature (Magnúsdóttir, S., et al., 2017), and we collected the GEMs of microbiota in the ileal terminus and obtained the abundance of each strain (Stolaki, M., et al., 2019). Since the samples used for the construction of these models are intestinal microorganisms, it is no longer necessary to specifically adjust them to the state in an anaerobic environment.
The total amount of microbiota in ileal is approximately to be 10^10 cfu/L. Hence, if we assume that one cfu is produced by one microbial, we can calculate the total biomass of the ileal micro-community is about 0.1g/L. The relevant abundance of the ileal micro-community is from Stolaki, M., et al., 2019.
Fig 1. relevant abundance of ileal microbial (EcN: 0.0673%(set), others microbials, Stolaki, M., et al., 2019)
Here we roughly set the relevant abundance of EcN to be 0.0678%. While the others are based on the paper mentioned above.
The environment of the human gut is closely related to a person's dietary habits and food intake and thus varies from person to person. To ensure that the model is both general and has Chinese characteristics, we entered the recipes recommended by Scientific Research Report on the Chinese Dietary Guidelines (2021) into the vmh website, which will return the fluxes that these foods can provide to the intestine, and use these data to approximate the concentrations of each substance in the human intestine.
Reaction | Value |
EX_starch1200(e) | 0.055320598 |
EX_strch1(e) | 5.975178884 |
EX_strch2[e] | 21.3399344 |
EX_ca2(e) | 9.553540057 |
EX_fe2(e) | 0.068240841 |
EX_fe3(e) | 0.068240841 |
EX_mg2(e) | 5.414733017 |
EX_pi(e) | 6.443149167 |
EX_k(e) | 35.78788526 |
EX_na1(e) | 20.43628026 |
EX_zn2(e) | 0.064912361 |
EX_cu2(e) | 0.006767984 |
EX_mn2(e) | 0.019159391 |
EX_h2o(e) | 28679.28387 |
EX_thm(e) | 0.00301431 |
EX_ribflv(e) | 0.002697518 |
EX_nac(e) | 0.042882923 |
EX_ncam(e) | 0.042874812 |
EX_pydam(e) | 0.002423892 |
EX_pydx(e) | 0.002453458 |
EX_pydxn(e) | 0.002424224 |
EX_5mthf(e) | 0.000103187 |
EX_fol(e) | 0.000107418 |
EX_thf(e) | 0.000106686 |
EX_ascb_L(e) | 0.349313552 |
EX_avite1(e) | 0.003242053 |
EX_but(e) | 1.133635357 |
EX_octa(e) | 0.689483605 |
EX_dca(e) | 0.623314774 |
EX_ddca(e) | 0.572278149 |
EX_ttdca(e) | 2.373634358 |
EX_hdca(e) | 13.95448135 |
EX_ocdca(e) | 5.292521997 |
EX_hdcea(e) | 2.173033524 |
EX_ocdcea(e) | 18.99451294 |
EX_CE2510(e) | 0.279557887 |
EX_doco13ac(e) | 0.398512927 |
EX_lnlc(e) | 7.784776551 |
EX_strdnc(e) | 0.070185018 |
EX_arachd(e) | 0.386771562 |
EX_tmndnc[e] | 0.270066543 |
EX_clpnd(e) | 0.091364251 |
EX_crvnc(e) | 0.163115003 |
EX_chsterol(e) | 0.342734031 |
EX_etoh(e) | 0 |
EX_trp_L(e) | 1.258446518 |
EX_thr_L(e) | 8.243926242 |
EX_ile_L(e) | 8.076639147 |
EX_leu_L(e) | 13.93394908 |
EX_lys_L(e) | 11.93422702 |
EX_met_L(e) | 3.805779321 |
EX_phe_L(e) | 5.943971455 |
EX_tyr_L(e) | 4.477115868 |
EX_val_L(e) | 10.69613606 |
EX_arg_L(e) | 7.03376007 |
EX_his_L(e) | 4.196621926 |
EX_ala_D[e] | 7.015297709 |
EX_ala_L(e) | 7.015297709 |
EX_asp_D[e] | 9.001163731 |
EX_asp_L(e) | 9.001163731 |
EX_glu_L(e) | 0 |
EX_gly(e) | 13.01341497 |
EX_pro_D[e] | 5.289697117 |
EX_pro_L(e) | 5.289697117 |
EX_ser_L(e) | 10.2893456 |
EX_pnto_R(e) | 0.010520324 |
EX_chol[e] | 1.483474813 |
EX_adpcbl(e) | 9.4134E-07 |
EX_retinol(e) | 0.000430404 |
EX_caro(e) | 0.009060913 |
EX_phyQ(e) | 5.86568E-05 |
EX_sucr(e) | 13.99658592 |
EX_glc_D(e) | 48 |
EX_fru(e) | 37.2261171 |
EX_lcts(e) | 19.26835772 |
EX_malt(e) | 0.222571789 |
EX_gal(e) | 0 |
EX_ptdca(e) | 0.006492597 |
EX_hpdca(e) | 0.015995702 |
EX_arach(e) | 0.001886596 |
EX_docosac(e) | 0.002307653 |
EX_lgnc(e) | 0 |
EX_ttdcea(e) | 0.006085456 |
EX_CE4843(e) | 0.011467997 |
EX_lnlnca(e) | 0.025421965 |
EX_vitd3(e) | 5.4328E-06 |
Table 1. Concentrations of some metabolites in ileal
1, microbiota is uniformly distributed at
the
ileum
2, the absorption coefficient of different substances by the same microbial is the same
3, the absorption rate of these substances by the body does not change with time
4, the model only simulates the existence of competing substances, that is, it is considered
that in addition to erythritol only one microbial can use the substance is sufficient, the
concentration will not be affected with time
5, there is no exchange of signal molecules between intestinal microbials, the competitive
relationship only through Competing substances in the form of expression
6, the intestinal microbiota dilution rate is the same
7, the goal of the growth of intestinal microbiota is to make the maximum Biomass
8, the intestinal material flow uniform rate
$I_j$: indexes of metabolites in organism j
$J$: indexes of the organism
$K_j$: the coupled metabolites of organism j
$X_j$: the biomass of organism j
$C_i$: the indexes of metabolites
$v_{lb j}^{i}$: uptake flux of an exchange reaction related to metabolites I in organism
j
$v_{max,j}^i$: max uptake flux of an exchange reaction related to metabolites I in organism
js
$k_{i,j}$: Monod constantss
$v_{s,j}^i$: the solution flux of an exchange reaction related to metabolites I in organism
js
$S_j$: Stoichiometric matrix of organism js
${\mu}_j$: Growth rate of organism js
$d_j$: dilution rate of organism js
$U_{human}^i$: uptake rate of ileals
\begin{equation} \frac{dX_j}{dt}={\mu}_jX_j=d_jX_j \label{(1)} \end{equation}
The biomass of each organism is determined by the growth rate and dilution rate. The growth rate is solved by the (3), the LP problem, while the dilution rate is constant with the specific microbial strain, won’t be affected by time or other factors.
\begin{equation} v_{lb j}^i=v_{max,j}^i\frac{C_i}{k_{i,j}+C_i} \label{(2)} \end{equation}
In constrain-based modelling(cobra), the lower bound of an exchange reaction means the ability for a microbial to uptake a specific metabolite i. And refers to the maximum ability to uptake a metabolite, the value of this term is set by the lower bound value in the GEMs accessed from AGORA. And the monod constant is set to 0.1mM, which is used in previous work(Popp, D., et al., 2020).
\begin{equation} max:{\mu}_j \label{(3)} \end{equation} \begin{equation} subject\ to: S_jv_j=0 \label{(4)} \end{equation} \begin{equation} v_{lb,j}^i \le v_{j}^i \le v_{ub,j}^i \label{(5)} \end{equation}
The lines above are the FBA. The object in each organism is set to be the biomass reaction, referring that the microbial will always try its best to grow. This object is widely used in cobra model. As the (4) shown, it is the expression of the constrain based on the conservation of mass theorem while the (5) is the constrain for the extent of reaction specific in the microbial.
\begin{equation} ({\exists}j' \in J, j' \neq j, i_j \in K_j \bigwedge i_j \in K_{j'})\longrightarrow i_j \in K_j \label{(6)} \end{equation} \begin{equation} \frac{dC_i}{dt}=\sum_{j=1,i\in K_j}^{n}v_{s,j}^iX_j-U_{human}^iC_i \label{(7)} \end{equation}
The (6) explains the concept of the coupled metabolites in math. By considering coupled metabolites, the algorithm can avoid some single metabolites that only utilized by single strain. We would rather assume these single metabolites to be adequate instead of the method used in (2) in order to keep the model stable. The single metabolites may arise from the inconsistency of the SBML files, and the coupled metabolites have an overwhelming amount compared to the single ones.
The (7) describes the time derivative of the coupled metabolites, the consumption of them is arisen from the microbial and human. Here, we assume the human uptake rate to be constant.(see code)
Fig 2. Biomass of 25 microbial in ileal in first 12 hours
We simulated the growth of the above-mentioned 25 species of microbial in the ileal terminal environment, and the change of substances in the starting 12 hours. The starting concentrations of microbiota concentrations were set according to the concentration of ileal microbiota and the relative abundance of various microbials, while the concentrations of various substances in the environment were set according to the data in Table1. The results showed that the growth trend of most microbials gradually leveled off within 12 hours, reflecting that the provided nutrients were going to be depleted within 12 hours, and the specific substance change trend is shown in Fig3. Within 12 hours, we can see the growth of EcN is negligibly low compared to other microbials, and the reason for this result is that the initial abundance of provided EcN is low, only 0.0678%. Therefore, it is obviously difficult to see the promoting effect of erythritol on EcN in terms of population density.
Fig 3. Concentration of 269 couple metabolites in ileal
But when considering the growth rates of ileal microbials in ileum, the advantage of EcN manifests. With erythritol, the extra-carbon source, existing, EcN maintains relatively low growth rate under the condition of nutrition-rich environment. However, given enough time, EcN is capable surpassed the grow-slow strain (such as ClostridiumsymbiosumATCC14940, etc.) under such condition.
The right half of the curve in FIG. 4 shows the growth rates of each microbial in the gut when glucose and other major carbon sources are no longer sufficient.
As you can see, the growth rate of Bacteroides vulgatus ATCC 8482 decreased greatly under the condition of insufficient nutrition. However, in contrast, under the same condition of insufficient main carbon source, EcN was supported by erythritol (as an independent carbon source, erythritol could only be absorbed by EcN, so it was relatively sufficient), and the decrease of growth rate was not so drastic as that of other microbials It was also noted that when the main carbon source was insufficient, the growth rate of EcN was above the medium level in the population, which further indicated the possibility of EcN colonization supported by erythritol as an independent carbon source.
Fig 4. (a) Growth rate of 25 microbial in ileal in first 12 hours
(b) emphasizing the EcN (iDK1463)
MICOM has been mentioned in the previous section. In this section, we still use MICOM to do a simple verification of the results of our model. The construction of growth medium in MICOM was realized by the method mentioned in this section (see Table1). For unmodified EcN, its growth rate in intestinal microbiota was the lowest (Fig5A), which also reflects the background on which our project is based -- EcN is often difficult in the process of colonization.
After the transformation, the growth rate of EcN that could utilize erythritol was increased to a certain extent (Figure 5B). Comparing the growth rate of EcN before and after being modified, it could be found that the growth rate of some microbials would change with the growth of EcN.
Fig 5. Growth rates of the ileal microbial community. (a), community growth rate with unmodified EcN. (b) community growth rate with modified EcN. (c) Difference between modified and unmodified EcN
The initial abundance of engineered EcN largely determines the outcome of the competition. If the initial abundance of EcN is too low, it will lead to a decline in the competitiveness of major carbon sources (such as glucose, etc.), resulting failure in competition. Only by ensuring that the initial abundance of engineering EcN reaches a specific proportion in the population, can EcN be competitive enough. (FIG. 7).
Fig 6. microbial community growth rate with modified EcN in different relevant abundance
This conclusion will also guide us to increase the abundance of EcN as much as possible in the actual design process, such as the use of wrapping pods or other wrapping materials.
1. Kumar, M., Ji, B., Babaei, P., Das, P., Lappa, D., Ramakrishnan, G., Fox, T.E., Haque, R., Petri, W.A., Bäckhed, F., Nielsen, J., 2018. Gut microbiota dysbiosis is associated with malnutrition and reduced plasma amino acid levels: Lessons from genome-scale metabolic modeling. Metabolic Engineering 49, 128–142. https://doi.org/10.1016/j.ymben.2018.07.018
2. Magnúsdóttir, S., Heinken, A., Kutt, L., Ravcheev, D.A., Bauer, E., Noronha, A., Greenhalgh, K., Jäger, C., Baginska, J., Wilmes, P., Fleming, R.M.T., Thiele, I., 2017. Generation of genome-scale metabolic reconstructions for 773 members of the human gut microbiota. Nat Biotechnol 35, 81–89. https://doi.org/10.1038/nbt.3703
3. Popp, D., Centler, F., 2020. μBialSim: Constraint-Based Dynamic Simulation of Complex Microbiomes. Front. Bioeng. Biotechnol. 8, 574. https://doi.org/10.3389/fbioe.2020.00574
4. Stolaki, M., Minekus, M., Venema, K., Lahti, L., Smid, E.J., Kleerebezem, M., Zoetendal, E.G., 2019. Microbial communities in a dynamic in vitro model for the human ileum resemble the human ileal microbiota. FEMS Microbiology Ecology 95, fiz096. https://doi.org/10.1093/femsec/fiz096
5. https://vmh.life