XJTU-China
The Growth Model
Table of Contents
Background and Assumptions
- We ignore other organisms and most environmental changes (such as rain, wind, CO2 and or so), aiming to simulate merely the interations between Bacillus Subtilis and Nostoc sp. Since the soil environment is really complicated, it is impossible to take everything into account. However, in regard of the nitrogen and phosphrus fixing, we may use the average rate considering all microorganisms in the soil.
- We assume that the two microorganisms and all related nutrients are distributed evenly in the soil. In other words, they are living in a solution-like system and spatial heterogeneity is not considered. Thus, we can use biomass concentration and nutrient concentrition to describe this system quantitively.
- “Nutrient” means matter that bacteria and cyanobacteria can directly and quickly, or “utilizable nutrients”. C, N, P exist in plenty of forms in the soil. Literature call them “dissolved nitrogen”, “soluable phosphorus” and so on. Only low molecular weight (LMW) or ions can be quickly utilized. For example, amino acids and soluable phosphate are “utilizable nutrients”, while humus and Ca3(PO4)2 are not, which need enzymatical decomposition.
- We agree that after Nostoc sp. die, the cell matter is quickly decomposed by bacteria and return to “utilizable nutrients”, while Bacillus Subtilis don’t do this[3].
- Another assumptions of established models we used.
- We are not considering spores, because they are not actually metabolizing. In other words, they neither produce or consume nutrients and EPS.
- We are not considering spore germination and metabolism of spores, because they are not actually metabolizing and the environment can not meet the suitable conditions for germination. In other words, they neither produce or consume nutrients and EPS.
Mathematical Model
Equations
- $N_i$ ($\mathrm{g/L}$): biomass concentration (cell dry weight per unit volume);
- $R_j$ ($\mathrm{g/L}$): nutrient concentrition (mass of the element per unit volume).
No. | Equation | Descriptions |
---|---|---|
1 | $\dfrac{\mathrm{d}N_i}{\mathrm{d}t}=Birth_i - Death_i$ | For both microorganisms, the net growth rate $\frac{\mathrm{d}N_i}{\mathrm{d}t}$ is defined by birth rate minus death rate. |
2 | $Birth_i=f_i r_i N_i$ | Birth rate $Birth_i$ is modified by a nutrient-constraint factor $f_i(R)$ which is not more than 1. |
3 | $Death_i=m_i N_i$ | Growth and death rates are proportional to current population $N_i$ quantity with a ratio. |
4 | $f_i(R)=\min\limits_j\{\dfrac{R_j}{R_j+K_{ij}N_i}\}$ | The constraint factor is built according to Monod growth kinetics with Liebig's law of the minimum. |
5 | $\dfrac{\mathrm{d}R_j}{\mathrm{d}t}=Income_j-\sum\limits_i Q_{ij} Grow_i$ or $\dfrac{\mathrm{d}N_i}{\mathrm{d}t}$ | The concentration of nutrients are dynamically regulated by production and consumption by microorganisms. |
Nutrient-constraint factor
-
$f_i(R)$ looks quite similar to Michaelis-Menten equation, which is used by many former iGem teams. When $R_j$ goes to infinity, $f_i(R)$ approaches to 1, thus have no effect on the birth rate; otherwise, the birth rate seriously decrease, even leading that net growth rate is less than 0.
-
$f_i(R)$ reflects how much each nutrient limits growth. Inspired by Liebig's law of the minimum, the overall effect is determined by the most lacking nutrient. Even if carbon and nitrogen is sufficient, when phosphrous is deficient, microorganisms cannot grow.
-
$K_{ij}$ is a half-saturation constant like Michaelis-Menten constant. However, most literatures ignored that the amount of nutrient needed obviously depends on biomass concentration. Therefore, we adopt half-saturation constant per unit biomass concentration, forming a variable and more realistic half-saturation constant. This is actually based on [4]'s research.
Dynamic change of nutrients
-
In order to monitor the nutrition limitation ($f_i$) encountered by microorganisms, we established the dynamic change equation of resources according to the physiological metabolism of microorganisms in the environment.
The detailed income and output is shown in equation 7:
- Based on assumption 4, the loss of nutrients is proportional to net growth rate ($\frac{\mathrm{d}N_2}{\mathrm{d}t}$) for Nostoc sp. and birth rate for *B. Subtilis* ($Grow_1$). So dead *B. Subtilis* (or spores) don't give its cell matter back to $R_c$, $R_n$ and $R_p$.
- Obviously, every resource will be consumed as a part of an organism when it grows. $Q_{ij}$ means consumption rate of nutrient $j$ per unit biomass concentration.
- Aerobic respiration is alway ongoing, so we have to use the ratio between C consumption and biomass increase as $Q$.
- The consumption of carbon resources includes the growth and reproduction of microorganisms and the dissipation of respiration; At the same time, photosynthesis from Nostoc sp. will supplement carbon source.
- $Re_1$ represents the consumption rate of carbon source per unit biomass B.S;
- $NPh$ means the net fixation rate of organic carbon obtained by photosynthesis rate per unit biomass vegetable Nostoc sp. ,which offsets its respiration.
- The formation of microorganisms is the main way to consume N resources, and the supplement comes from nitrogen mineralization mediated by Bacillus Subtilis (decomposition of organic matter by extracellular enzymes) and nitrogen fixation by Nostoc sp. using N~2~ in the atmosphere. We have not considered the nitrogen cycle process caused by other environmental factors, although they often play a potentially important role.
- $N_{sol}$ represents the nitrogen mineralization rate per unit biomass Bacillus Subtilis;
- $N_{fix}$ is the nitrogen fixation rate per unit biomass Nostoc sp.
- The consumption of P element is similar to that of N element, and Bacillus Subtilis takes the secretion of organic acids and extracellular enzymes to mediate the solubilization of phosphate minerals and organophosphorus in the environment as the supplement of resources.
- $P_{sol}$ represents the rate of producing available P per unit biomass Bacillus Subtilis.
Parameters' Values
All $K$ and $Q$ are dimensionless. Units of rates are all $\mathrm{g/(L\cdot min)}$ .
Name | Value | Descriptions | References |
---|---|---|---|
$K_{1c}, K_{1n}, K_{1p},K_{2n}$ | 1.1460,0.638,0.057,0.051 | The paper is using Contois model for Bacillus Subtilis and pseudomonas aeruginosa. | [5] |
$K_{2c},K_{2p}$ | 0,0.0187 | $K_{2c}$ given by the paper equals to 0; $K2p$ for Cyanothece sp. is given. | [6] |
$Q_{1c}$ | 0.6528 | Calculated from the slope of the “nutrient cosumption”-$\mu$ curve. It means the reciprocal of “growth yield”. | [7] |
$Q_{1n},Q_{1p}$ | 0.0734,0.00607 | Calculated from mass ratio between elements. $Q_{1n}$: P limited scenario; vice versa. | [8] |
$Q_{2c}$ | 0.2523 | We do’t have access to autotropical carbon consumption/biomass concentration. We just used elemental mass fraction. | [9] |
$Q_{2n}$ | 0.0846 | As for N and P, we use their proportion in cell dry weight because there is no extra loss. | [9] |
$Q_{2p}$ | 0.00248 | Ditto. Species: “blue-green algae” | [10] |
$Re_1$ | 1.34e-4 | We calculated it from “maintenance coefficient”, intercept of the “nutrient cosumption”-$\mu$ curve. | [7] |
$NPh$ | 9.27e-3 | Using the luminosity on a sunny day, we determined it based on a luminosity-net photosynthesis rate model $NPh=NPh_{max}\dfrac{I}{K_s+I+I^2/K_i}$ | [11] |
$N_{fix}$ | 7.38e-6 | Calculated from control group N-fixing rate and convert it to rate per unit biomass concentration. | [12] |
$N_{sol}$ | 2.60e-4 | Calculated from average N-fixing rate of mixed soil microorganisms under different biomass concentration. | [13] |
$P_{sol}$ | 2.15e-5 | Calculated from Bacillius megaterium P solubilizing rate of poultry bones. | [14] |
Result Analysis
Growth curve
We use the data obtained from the literature and the experimental analysis results to provide a set of parameters with high reliability for the growth model. After trial, we used $0.1 \mathrm{g/L}$ and $0.2 \mathrm{g/L}$ as the initial biomass of Bacillus subtilis and Nostoc sp., respectively. The following analysis uses this initial value. Using the parameters and initial values mentioned above, the following is the result of growth simulation by the model:
图。
Parameters sensitivity
Microorganism | 5 days | 10 days | 15 days | 20 days | 25 days |
---|---|---|---|---|---|
B. subtilis | 0.485 | 1.373 | 4.374 | 7.259 | 11.324 |
Nostoc. sp | 0.140 | 0.500 | 1.203 | 2.236 | 5.268 |
Expansion of parameters range
$m_1$
$P_{sol}$ and $N_{sol}$
$P_{sol}$
Explanation: A faster phosphorus dissolution rate could alleviate the P-limitation experienced by B. subtilis and Nostoc sp. at an early stage. After exceeding a certain threshold, due to the early arrival of the N-limiting stage, the fast-growing Nostoc sp. inhibits the growth of B. subtilis by competing for N, and this effect will gradually flatten out.
$N_{sol}$
Explanation: A higher nitrogen mineralization rate allowed B. subtilis and Nostoc sp. to enter the N-limiting stage later. B. subtilis without timely N-limiting would have a higher growth rate and inhibited the growth of Nostoc sp. by competing for P. When the nitrogen mineralization rate increases to a certain value, the N-limiting stage will be completely lifted and the growth condition will not be changed.
$K_{1N}$
$K_{1n}$
$K_{2n}$
$K_{1p}$
$K_{2p}$
$Q_{1n}$
$Q_{1p}$
The smaller the $Q_{1p}$, the more P will be accumulated in the early stage. Bacillus Subtilis will grow faster temporarily, and then the fast-growing Nostoc sp consumes a lot of N, which inhibits the growth of Bacillus Subtilis. The smaller $Q_{1p}$ is, the more obvious this phenomenon is, while the larger will not appear. Therefore, with the increase of $Q_{1p}$, the growth of both will first increase and then decrease. Therefore, we find that: $Q_{1p}$ is too small to cause Nostoc sp to grow too fast, which is not conducive to resource accumulation; Too much $Q_{1p}$ limits the growth of Bacillus Subtilis and is not conducive to resource accumulation.
polysaccharide production
Experiments proved that the polysaccharide yield of Bacillus subtilis can be up to 2.06 times as high as before by transferring different combinations of pgmA and galU. We predicted the real-time yield of polysaccharide by the growth model and the Luedeking-Piret equation, so as to predict whether the yield of exopolysaccharide can be significantly improved in the early stage. The yield of polysaccharide is very important for the stability and functionality of the system.
Under normal circumstances, there was no obvious difference in polysaccharide yield between them in the early stage. Comparing Fig. 45 and 46, we can find that the polysaccharide yield of B. subtilis is obviously improved in the early stage after the synthetic biological transformation, and the polysaccharide yield of B. subtilis will be more than that of Nostoc sp. in a long period of time. This indicates that the function of B. subtilis will be greatly improved in the early stage. This is in line with our project design. From the early stage when B. subtilisplayed a leading role, it will gradually transition to the period when Nostoc sp. played its role.
Conclusion and Prospects
- We found that Bacillus Subtilis and Nostoc sp. can successfully so-survive in the soil maily due to their high Complementarity in function. They can effectively increase the content of bioavailable nutrients in the soil while growing. This proves that our symbiotic system can improve the soil effectively. EPS can play an effective role when living things can survive stably.Through parameter sensitivity analysis, we found that the system was highly sensitive to the change of $m_1$, which indicated that the change of mortality rate of *Bacillus Subtilis* caused by external disturbance and environmental change had a large impact on the system.
- Though being a moderately good simulation of the growth of the two microorganisms, our growth model is not flawless due to neglecting various environmental factors and needs perfection. By extending the range of parameters, we think that the model can still be simulated for different environments and biological systems. Therefore, our model has great potential and prospects.
[2] He, Mingzhu & Hu, Rui & Jia, Rong-Liang. (2019). Biological soil crusts enhance the recovery of nutrient levels of surface dune soil in arid desert regions. Ecological Indicators. 10.1016/ j.ecolind.2019.105497.
[3] Xiao-jing, S. U. N., Boqiang, Q. I. N., Guangwei, Z. H. U., & Science, E. (2007). Release of colloidal phosphorus, nitrogen and organic carbon in the course of dying and decomposing of cyanobacteria. China Environmental Science, 27(3), 341–345.
[4] CONTOIS, D. E. (1959). Kinetics of Bacterial Growth: Relationship between Population Density and Specific Growth Rate of Continuous Cultures. Journal of General Microbiology, 21(1), 40– 50. doi: 10.1099/0022128721140
[5] Eswari, J. S. (2020). Therapeutic and Biodegradable Green Surfactant Kinetic Studies and Parameter Estimation with Multiple Substrates. AntiInfective Agents, 16(2), 114–120.
[6] Zhang, D., Dechatiwongse, P., DelRioChanona, E. A., Hellgardt, K., Maitland, G. C., & Vassil iadis, V. S. (2015). Analysis of the cyanobacterial hydrogen photoproduction process via model identification and process simulation. Chemical Engineering Science, 128, 130–146. Retrieved from http://dx.doi.org/10.1016/j.ces.2015.01.059 doi: 10.1016/j.ces.2015.01.059
[7] Zamboni, N., Mouncey, N., Hohmann, H. P., & Sauer, U. (2003). Reducing maintenance metabolism by metabolic engineering of respiration improves riboflavin production by Bacillus subtilis. Metabolic Engineering, 5(1), 49–55. doi: 10.1016/S10967176(03)000077
[8]M. Dauner, T. Storni, & U. Sauer. (2001). Bacillus subtilis metabolism and energetics in carbon-limited and excess-carbon chemostat culture. Journal of Bacteriology.
[9] Kimura, Y., Kimura, S., Sato, S., Katoh, H., Abe, T., Arai, M., & TomitaYokotani, K. (2015). Eval uation of a Cyanobacterium, Nostoc sp. HK01, as Food Material for Space Agriculture on Mars. Biological Sciences in Space, 29(0), 24–31. doi: 10.2187/bss.29.24
[10] Yan, R., & Li, Y. (2010). Determination of 14 elements in bluegreen algae of Chaohu Lake by ICPAES. Analytical Instrumentation(04), 39–41. doi: 10.1017/CBO9781107415324.004
[11] Dai, R., Chu, Z., & Yu, X. (2012). Effects of Gas Vesicle Collapse under Pressure on Growth and Photosynthesis of Three Planktonic Cyanobacteria Strains. Research of Environmental Sciences, 25(01), 30–35.
[12] Popović, S., Nikolić, N., Jovanović, J., Predojević, D., Trbojević, I., Manić, L., & Simić, G. S. (2019). Cyanobacterial and algal abundance and biomass in cave biofilms and relation to en vironmental and biofilm parameters. International Journal of Speleology, 48(1), 49–61. doi: 10.5038/1827806x.48.1.2224
[13] Song, X., Zhu, J., He, N., Huang, J., Tian, J., Zhao, X., … Wang, C. (2017). Asynchronous pulse responses of soil carbon and nitrogen mineralization to rewetting events at a shortterm: Regulation by microbes. Scientific Reports, 7(1), 1–10. Retrieved from http://dx.doi.org/ 10.1038/s41598-017-07744-1 doi: 10.1038/s41598017077441
[14] Wyciszkiewicz, M., Saeid, A., DobrowolskaIwanek, J., & Chojnacka, K. (2016). Utilization of mi croorganisms in the solubilization of lowquality phosphorus raw material. Ecological Engineer ing, 89, 109–113. Retrieved from http://dx.doi.org/10.1016/j.ecoleng.2016.01.065 doi: 10.1016/j.ecoleng.2016.01.065