A Mathematical Model for Characterizing the Biomass and the Physiological/Biochemical Indicators of Salvia miltiorrhiza Based on Growth-Defense Tradeoff

Carbon(C) and nitrogen(N) metabolisms are important for plant growth and defense, and enzymes play a major role in these two metabolisms. Current studies show that the enzymes of N Metabolism, C Metabolism, and defense are correlated with biomass. Then, we conducted this research under the assumption that enzymes could characterize the relationship based on growth-defense tradeoff, and some of the enzymes could be used to represent the plant growth. From the mechanism model, we picked out 18 physiological/biochemical indicators and obtained the data from 24 tissue culture seedlings of Salvia miltiorrhiza (S.miltiorrhiza) which were grafted with 11 endophytic fungi. Then, the relationship between the biomass and the physiological/biochemical indicators was investigated by using statistical analysis, such as correlation analysis, variable screening, and regression analysis. The results showed that many physiological/biochemical indicators, especially enzyme activities, were related to biomass accumulation. Through a rigorous logical reasoning process, we established a mathematical model of the biomass and 6 key physiological/biochemical indicators, including glutamine synthetase (GS), glutamate synthase (GLS), glutamate dehydrogenase (GDH), peroxidase (POD), catalase (CAT), and soluble protein from Cobb-Douglas production function. This model had high prediction accuracy, and it could simplify the measurement of biomass. During the artificial cultivation of S.miltiorrhiza, we can monitor the biomass accumulation by scaling the key physiological/biochemical indicators in the leaves. Interestingly, the coefficients of Lasso regression during our analysis were consistent with the mechanism of growth-defense tradeoff. Perhaps, the key physiological/biochemical indicators obtained in the statistical analysis are related to the indicators affecting biomass accumulation in practice.

1. Introduction

As is known, Carbon (C) and nitrogen(N) are the main elements for the growth and defense of plants. Plants synthesize propane on the one hand through photosynthesis, which leads to starch, sucrose, and soluble sugars, and on the other hand they translate inorganic N into amino acids from ammonification, which leads to soluble proteins and are used for growth and defense in a certain C:N ratio (Lea and Morot-Gaudry, 2001). At the same time, plants consume the assimilated products of C/N metabolism through defense under biotic and abiotic stresses. Therefore, plants maintain a dynamic tradeoff between growth and defense by optimizing resource allocation through C/N metabolism and defense to enable plants to survive (Li et al., 2019).

Current studies show that some key enzyme activities of C/N metabolism have an available correlation with biomass accumulation (Noor et al., 2010). The activities of defensive enzymes show a significant correlation with adaptation to environmental stresses (Haddidi et al., 2020). Thus, the activities of sucrose phosphate synthase (SPS) and sucrose synthase (SS) in leaves, as well as the content of chlorophyll, soluble sugars, and starch are often used as important indicators to characterize the status of the C metabolism (Moriguchi et al., 1992; Verma et al., 2011). In N metabolism, nitrate reductase (NR) and nitrite reductase (NiR) translate N into NH${}_4^{+}$, formulate glutamate and glutamine through ammonia assimilation (Liang et al., 2018; Salehin et al., 2019). Glutaminase (GLS), glutamine synthetase (GS), glutamate synthase (GOGAT), and glutamate dehydrogenase (GDH) can promote amino acid accumulation. But with a large increase in the reduction of NH${}_4^{+}$ from NO${}_3^{-}$, a sustained raising in GDH amination activity will inhibit GLS, GS, and GOGAT activities. Thus, GLS, GS, and GOGAT have a positive influence on biomass accumulation, while GDH has a bidirectional effect on biomass accumulation (Liang et al., 2018; Salehin et al., 2019; Gonzíalez-Moro et al., 2021; Wei et al., 2021). Many soluble proteins are crucial components of enzymes. They are involved in the physiological/biochemical metabolic and are the key indicators of whether plants suffered from heavy metal stress.

In the competition of biotic and abiotic stresses, the activities of enzymes, such as GS, GDH, NR, GOGAT, NiR, and GLS, as well as the soluble protein content, are considerable indicators of the N metabolic (Salehin et al., 2019; Gonzíalez-Moro et al., 2021; Wei et al., 2021). The enzymatic activities of superoxide dismutase (SOD), peroxidase (POD), catalase(CAT), proline (Pro), phenylalanine ammonialyase (PAL), and malondialdehyde (MDA) are commonly used to characterize defense responses (Jaafar et al., 2012; Farooq et al., 2020; Sarker and Oba, 2020; Zaheer et al., 2020). The increase of the enzymes activities such as SOD, POD, CAT, and Pro can significantly improve plant growth, biomass, chlorophyll content, and gas exchange properties (Sarker and Oba, 2020; Zaheer et al., 2020), while PAL and MDA inhibit plant growth by reducing the activity of antioxidant enzymes through oxidative stress (Jaafar et al., 2012; Farooq et al., 2020). Although the relationship between C/N metabolism and plant growth defense is obvious, there is still a lack of systematic research.

Based on the principles of metabolism, in this study we hypothesized that the correlation among the enzymes of C/N metabolic and defense could characterize the growth-defense tradeoff, and some of the enzymes could indicate the connection between the biomass and the physiological/biochemical indicators. Then, a mechanism model was established and 18 physiological/biochemical indicators were picked out (Figure 1). Since S.miltiorrhiza is considerable for the treatment of coronary heart disease and cerebrovascular disease (Su et al., 2015; Ma et al., 2016; Wang et al., 2017; Li, 2018; Shi et al., 2019), this study has practical significance to guide the production of medicinal plants. We cultured 24 tissue culture seedlings of Salvia miltiorrhiza (S.miltiorrhiza) which were grafted with 11 endophytic fungi. The 11 fungi are non-pathogenic and can intervene in physiological metabolism so that we can obtain data on different growth states. Through lasso screening variable (Efron et al., 2004) and regression analysis of plant physiological and biochemical indexes and biomass, the functional relationship between S.miltiorrhiza biomass and physiological /biochemical indexes was found. Thus, through the monitoring of relevant indicators of S.miltiorrhiza, we can understand the growth status of S.miltiorrhiza.

FIGURE 1

Figure 1. The path structure of biomass accumulation of S.miltiorrhiza. It shows the process of C/N metabolism and the role of the physiological/biochemical indicators in biomass accumulation.

2. Materials and Methods

2.1. Plant Materials

In this study, 24 tissue culture seedlings were inoculated from 11 strains of non-pathogenic endophytic fungi from 8 species of S.miltiorrhiza which our group obtained in previous studies (Ya-Li, 2018). The culture of tissue culture seedlings and endophytic fungal inoculation were obtained by the method in Lan et al. (2016) (daytime 25^oC/nighttime 20^oC, 14/10 h, light intensity 3,000 LX) where young leaves of S.miltiorrhiza were employed as explants. The seedlings with 3–5 roots were transplanted into artificial soil containing 50 ml of 1/2 MS medium, sealed with a sealing film, and incubated in vertical light for 3 days (Hesheng, 2000), then a small piece of fungal cake was taken with a hole punch and placed on the root of the seedlings in the culture flask, and the root of S.miltiorrhiza was gently punctured with a sterile needle. After 30 days of incubation, the plants were removed and the biomass and related physiological indicators were measured.

2.2. Determination of the Biomass

The S.miltiorrhiza seedlings were taken out and washed carefully to remove impurities attached to the roots. After absorbing the surface moisture with absorbent paper, weigh it and subtract the weight of seedlings when transplanting to get the net biomass.

2.3. Determination of the Physiological/Biochemical Indicators

The activities of SS, SPS, NR, GS, and GOGAT were determined through enzyme solution prepared from young leaves of S.miltiorrhiza (Deane-Drummonda et al., 1979; Islam et al., 1996; Li, 2000; Zhao et al., 2003). The contents of reducing sugar and soluble sugar were estimated by adopting 3,5-dinitrosalicylic acid method and anthrone method, respectively (Zhang, 1990; Li, 2000; Gao, 2006). The soluble protein content was assessed with the coomassie brilliant blue method (Elsharkawy et al., 2012). Chlorophyll content was determined by spectrophotometry (Strain and Svec, 1966). The activities of SOD, POD, CAT, PAL, and MDA and Pro in S.miltiorrhiza leaves were measured according to Zhu et al. (1983), Polle et al. (1994), and Gao (2006).

2.4. Statistical Analysis

2.4.1. Correlation Analysis

Because of the large difference in unit and quantity, the physiological/biochemical indicators were standardized by the following function:

$\begin{array}{ll}{x}_{ij}^{\prime }=\frac{x_{ij}-{\overline{x}}_j}{s_j},& (1)\end{array}$

where x_ij was the value of the jth physiological/biochemical indicator of the ith plant of S.miltiorrhiza, ${\overline{x}}_j$ denoted the average of jth indicator, and s_j denoted the SD of the jth indicator.

The correlation analysis between the physiological/biochemical indicators of S.miltiorrhiza and the biomass included linear correlation analysis and nonlinear correlation analysis, so the correlation coefficient matrix analysis was chosen to be used. It included three aspects: (1) analysis of the variation between the biomass and the physiological/biochemical indicators; (2) linear correlation between the biomass and the physiological/biochemical indicators; and (3) analysis of the nonlinear relationship between the biomass and the physiological/biochemical indicators.

2.4.2. Systematic Analysis Based on Lasso Algorithm

According to the mechanism analysis, we tried to find the physiological/biochemical indicators related to the biomass accumulation of S.miltiorrhiza, so there would be duplication or the introduction of irrelevant factors. However, there are multiple covariates among various physiological/biochemical indicators, in which case least squares and partial least squares work poorly. In order to select a concise set of physiological/biochemical indicators to effectively predict biomass and improve the prediction accuracy of the model (Efron et al., 2004), Lasso regression was introduced for variable screening. The regression function was as follow:

$\begin{array}{ll}{y}_i={\displaystyle \sum _{j=1}^p}{\beta }_j{x}_{ij}^{\prime }+{\epsilon }_i,(i=1,2,...,n),& (2)\end{array}$

where $x_{ij}^{\prime }$ denoted the jth physiological/biochemical indicator of the ith sample, y_i denoted the biomass of S.miltiorrhiza of the ith sample, and p denoted the number of physiological/biochemical indicators screened.

Because $X_i={(x_{i1},\dots ,x_{ip})}^T,(i=1,2,...,p)$ were multicollinearity, the Lasso method was used to filter variables. The Lasso constructed a penalty function to obtain a more refined model, which made it compress some regression coefficients. Here, l₁-penalty function was used for regularization estimation parameter as in Efron et al. (2004), defined as

$\begin{array}{ll}\widehat{\beta }(\lambda )=\text{argmin}\frac{1}{N}\left|\right|Y-X\beta |{|}^2+\lambda {\displaystyle \sum _j}|{\beta }_j|,& (3)\end{array}$

where N denoted the number of samples, $Y={(y_1,y_2,...,y_N)}^T$ denoted the biomass, $X={(X_1,X_2,...,X_N)}^T$ denoted the physiological/biochemical indicators, β denoted the regression coefficient, and λ≥0 denoted the penalty parameter.

2.4.3. Nonlinear Regression Analysis Based on Cobb-Douglas Production Function

Based on the characteristics of important variables screened by Lasso, nonlinear regression analysis was utilized since correlations and interactions between variables still existed. Assuming that y denotes biomass, A_i(i = 1, 2, ⋯n) denotes factors with positive influence coefficients, and β_j(j = 1, 2, ⋯m) denotes factors with negative influence coefficients, similar to the Cobb-Douglas production model in Carter (2012), the models for the biomass and the physiological/biochemical indicators were established as follows:

$\begin{array}{ll}\begin{array}{c}Y=H·\frac{A_1^{a_1}{A}_2^{a_2}\cdots A_n^{a_n}}{B_1^{b_1}{B}_2^{b_2}\cdots B_n^{b_m}},\end{array}& (4)\end{array}$

Where H, a₁, a₂, ⋯ , a_n, b₁, b₂, ⋯ , b_m denoted parameters to be determined by the Equation (4). Then took common logarithms on both sides of the Equation (4) at the same time, that was

$\begin{array}{l}\ln Y=\ln H+a_1\ln A_1+a_2\ln A_2+\cdots +a_n\ln A_n-b_1\ln B_1\\ -b_2\ln B_2-\cdots -b_m\ln B_m.& (5)\end{array}$

The Equation (4) was transformed into a linear regression, from which the initial values of the nonlinear regression parameters could be obtained from the Equation (5). Then, the solution of Equation (4) could be optimized by performing a nonlinear regression.

2.4.4. Data Analysis Environment

All the data were analyzed in R4.0.3 + Rstudio (Chang, 2013; Lantz, 2013). glmnet() was used for filtering variables in Lasso, ls()/nls() was used for linear/nonlinear regression analysis, and ggplot() was used for graph plotting.

3. Results

3.1. Experimental Results and Data

To analyze the relationship between the biomass and the physiological/biochemical indicators, 24 samples of S.miltiorrhiza were collected according to the methods described in 2.1~2.3, and the data of the 18 physiological/biochemical indicators such as chlorophyll, SS, SPS, soluble sugar, and the values of corresponding biomass were displayed in Table 1.

TABLE 1

Table 1. Biomass and physiological/biochemical indicators of S. miltiorrhiza.

3.2. Descriptive Statistical Analysis of the Biomass and the Physiological/Biochemical Indicators

We defined biomass as the dependent variable(y) and chlorophyll(x₁), SS(x₂), SPS(x₃), soluble sugar(x₄), starch(x₅), GOGAT(x₆), GS(x₇), NR(x₈), NiR(x₉), GLS(x₁₀), GDH(x₁₁), soluble protein(x₁₂), POD(x₁₃), SOD(x₁₄), CAT(x₁₅), Pro(x₁₆), PAL(x₁₇), and MDA(x₁₈) as independent variables. To eliminate the difference in magnitude between variables, the criteria were standardized by using Equation (1) and then other corresponding analyses were performed.

3.2.1. Differences in Some Physiological/Biochemical Indicators and Lack of Consensus in the Biomass Accumulation

Some physiological/biochemical indicators were influenced by individual plants, especially SS(x₂) and GOGAT(x₆) were very different (Figure 2). It indicated that they lacked consensus in the biomass accumulation of S.miltiorrhiza and there was uncertainty among indicators.

FIGURE 2

Figure 2. The box line diagram of the physiological/biochemical indicators. The blue curve is the regression curve, the shaded part indicates the confidence interval, and the larger the shaded area the worse the regression effect of the curve. It shows that some physiological/biochemical indicators lack consensus in the biomass accumulation of S.miltiorrhiza.

3.2.2. Many Physiological/Biochemical Indicators Have a Linear (or Nonlinear) Correlation With the Biomass

From the result of correlation analysis, we can see there was a positive correlation of the biomass(y) with chlorophyll(x₁) and CAT(x₁₅) (Figure 3), and a negative correlation with GLS(x₁₀) and POD(x₁₃). Considering the relationship between the physiological/biochemical indicators, chlorophyll(x₁), NR(x₈), and CAT(x₁₅) were positively correlated, while chlorophyll(x₁) was negatively correlated with SPS(x₃), NiR(x₉), and POD(x₁₃). SPS(x₃) was positively correlated with CAT(x₁₅), while SPS(x₃) was negatively correlated with NiR(x₉), GLS(x₁₀) and POD(x₁₃). Similarly, these results indicated that linear correlations existed between the biomass and the physiological/biochemical indicators.

FIGURE 3

Figure 3. Correlations between the biomass and the physiological/biochemical indicators of S.miltiorrhiza biomass. The blue rectangles represent a positive correlation, the red rectangles represent a negative correlation. The area of the circle represents the degree of correlation. According to correlograms shown in this figure, there were correlations among the biomass and the physiological/biochemical indicators.

There was a certain nonlinear relationship (Figure 4) between the biomass and the physiological/biochemical indicators. For instance, there was a nonlinear relationship between the biomass(y), and NiR(x₉), and soluble protein(x₁₂). Thereby, biomass could not be expressed by a single enzyme activity or physiological/biochemical indicator, but by a combination of some effective physiological/biochemical indicators (Figures 3, 4). Therefore, it was necessary to systematically analyze the relationship between the physiological/biochemical indicators and the biomass.

FIGURE 4

Figure 4. The scatter and regression curves of the biomass and the physiological/biochemical indicators of S.miltiorrhiza. The horizontal axis of x_i, is the physiological/biochemical indicator of i and the vertical axis is y, which is the biomass of S.miltiorrhiza. The blue curve is the regression curve, the shaded part indicates the confidence interval, and the larger the shaded area indicates the worse regression effect of the curve.

3.3. Some Physiological/Biochemical Indicators Play a Major Role in the Biomass Accumulation

We used the glmnet() of the R software to calculate the Lasso model (parameters set to default values) and optimized the model and function by cv.glmnet() to obtain the results shown in Figure 5 and Table 2.

FIGURE 5

Figure 5. Lasso regression analysis number of iterations vs. mean square error. It shows that the optimal λ value is obtained in 5th step.

TABLE 2

Table 2. Nonzero regression coefficient of Lasso.

The results showed that GS(x₇), GDH(x₁₁), and CAT(x₁₅) were positively correlated with biomass and acted as promoters when the whole system was considered. GLS(x₁₀), POD(x₁₅), and soluble protein(x₁₃) were negatively correlated with biomass accumulation and acted as inhibitors.

3.4. Some Key Physiological/Biochemical Indicators Are Able to Express the Biomass Well

A nonlinear regression analysis was performed by using the Douglas production model with the influencing factors of GS(x₁₅), GLS(x₁₅), GDH(x₁₅), soluble protein(x₁₅), POD(x₁₅), and CAT(x₁₅), and biomass. Took 5/6 of the data as the training set and 1/6 as the verification set, repeated the training 1,000 times, and selected the regression equation with the smallest verification set error. The model was as follows:

$\begin{array}{ll}Y=5.05853\times \frac{x_7^{0.24737}{x}_{11}^{0.09026}{x}_{15}^{0.22631}}{x_{10}^{0.14776}{x}_{12}^{0.64021}{x}_{13}^{0.44381}}.& (6)\end{array}$

The predicted values of biomass could be obtained from Equation (6) and the results were shown in Table 2. The test showed that the prediction effect of the model was good (Tables 3, 4). It showed that the regression equation was relatively successful and was able to predict the corresponding the biomass from the 6 important physiological/biochemical indicators.

TABLE 3

Table 3. The results of the training set of the biomass and the key physiological/biochemical indicators.

TABLE 4

Table 4. The results of the validation set of the biomass and the key physiological/biochemical indicators.

4. Discussion

Carbon/nitrogen metabolism are critical in growth and defense (Li et al., 2019). Adequate management of the plant in C/N metabolism and growth defense have a significant impact on crop productivity. In this study, we constructed a growth mechanism model by S.miltiorrhiza and picked out the relationship between the 18 physiological/biochemical indicators and the C/N metabolism, with the aim of finding important metabolic traits as an indication of the biomass (Table 2).

The variability of the physiological/biochemical indicators was analyzed by box-line plots (Figure 2), and it was observed that some indicators showed no consensus. That was, some physiological/biochemical indicators were greatly affected by individuals. In particular, SS(x₂) and GOGAT(x₆) were more discrete. SS(x₂) was an important enzyme catalyzing sucrose synthesis in plants, and sucrose was broken down into soluble sugars for drought resistance in plants on the one hand, and synthesized biomass with soluble proteins on the other hand, making SS showed uncertainty on the biomass (Kaur et al., 2006; Zhang et al., 2018). Similarly, GOGAT(x₆) consuming amino acids used by plants for the biomass accumulation during growth defense, formed aspartic acid, which was used for plant drought resistance under the action of Pro and producing glutamate, soluble protein, which leads to biomass accumulation. The uncertainty of the two conversion pathways, SS(x₂), and GOGAT(x₆), prevents them from being used as important indicators for characterizing the biomass (Wei et al., 2021).

There were linear(or nonlinear) correlations between the biomass and the physiological/biochemical indicators (Figures 2, 3), indicating interdependence in the growth-defense tradeoff of plants. For example, there was a positive correlation between biomass(y) and chlorophyll(x₁)/CAT(x₁₅), indicating that chlorophyll was the main substance for photosynthesis in plants, and it played an important role in C metabolism (Luo, 2018). Although CAT(x₁₅) was mainly used in the defense, it reduced the toxicity of H₂O₂ (Abogadallah, 2010; Sarker and Oba, 2020), and maintained normal growth and defense. The biomass was negatively correlated with some physiological/biochemical indicators, such as GLS(x₆) promoting amino acid accumulation, and POD(x₁₃) depleting C and N assimilation products, affecting the biomass accumulation. POD(x₁₃) reduces H₂O₂ and converts carbohydrates to lignin, improving the physical defense barrier and depleting assimilates required for growth. Considering the relationship between the physiological/biochemical indicators, chlorophyll(x₁), NR(x₁₅), and CAT(x₁₅) were positively correlated, while chlorophyll(x₁₅) was negatively correlated with SPS(x₃), NiR(x₉), and POD(x₁₃) (Figure 3). On the one hand, chlorophyll is used by plants for photosynthesis to convert CO₂ in the air into C required by plants, and on the other hand, about 90% of N in plants, coming from biological N fixation, is reduced from nitrate to nitrite in soil by NR. Plant cells can rapidly transfer nitrite produced by nitrate reduction from the stroma to the chloroplasts of leaf cells or the plastids of root cells and are reduced in the chloroplasts by NiR to ammonia, while the potential toxicity of nitrite is broken down by CAT(x₁₅). Chlorophyll(x₁), NR(x₈), NiR(x₉), and CAT(x₁₅) collaborate with each other in C/N metabolism as well as in the growth-defenses tradeoff (Zhou et al., 2021). SPS(x₃) and SS(x₂) are both important enzymes for catalyzing sucrose synthesis in plants, however, SPS exhibiting different characteristics from sucrose SPS SS(x₂) was positively correlated with CAT(x₁₅), indicating that it promoted the biomass accumulation (Figure 3). While SPS was negatively correlated with NiR(x₉), GLS(x₁₀), and POD(x₁₃)(Figure 3). It showed inhibition of biomass accumulation, but NiR(x₁₅) reduced nitrate in soil and then synthesized NH${}_4^{+}$ for the biomass formation, and GLS(x₁₀) formed NH${}_4^{+}$ into amino acids, which were against the previous analysis. All of these analyses indicated that biomass accumulation was a relatively complex process that required systematic analysis of the physiological/biochemical indicators affecting the biomass accumulation.

By Lasso analysis (Figure 5 and Table 2), GS(x₇) and GDH(x₁₁) were positively correlated with biomass (Table 3), indicating that GS(x₇) and GDH(x₁₁) were the main substrate in the seedling stage facilitate protein synthesis with amino acids and were conducive to coordinate C metabolism (Kaur et al., 2006), which was consistent with the rapid assimilation of N in early growth and laid the foundation for high rate C assimilation in later growth stages (Salehin et al., 2019; Gonzíalez-Moro et al., 2021; Wei et al., 2021). GS(x₇) was responsible for primary NH${}_4^{+}$ assimilation as well as germination, and elevating of GS(x₇) in leaves contributed to the conversion of inorganic NH${}_4^{+}$ to organic N (Lea and Morot-Gaudry, 2001). There was a positive correlation between GDH(x₁₁) and total amino acid content in the root system (Gonzíalez-Moro et al., 2021). When GDH(x₁₁) was sufficiently induced in leaves, it would indicated substrate specificity in plants with other enzymes. Although GLS facilitated amino acid accumulation, it was negatively correlated with the biomass (Table 3), possibly due to the lack of sufficient reducing agent(NADPH) or the unavailability of reducing agent(NADPH), making it possible that NH${}_3^{+}$ might not attack the ketone group, leading to its release into the environment (Yoneyama et al., 2015). Soluble proteins were negatively correlated with biomass (Table 3), indicating that soluble proteins would consume amino acids used for biomass accumulation. POD was negatively correlated with biomass (Table 3), indicating that glutamate synthesis, as well as the synthesis of metabolic enzymes represented by POD, consumed N and C assimilation products, thus affecting the biomass accumulation. POD reduced H₂O₂ and converted carbohydrates to lignin, improving the physical defense barrier and depleting assimilates required for growth. However, despite its synthesis of N and C assimilation products, CAT(x₁₅) specifically reduced the toxicity of H₂O₂ (Abogadallah, 2010; Sarker and Oba, 2020), maintained cell membrane stability, and promoted growth and defense. Meanwhile, it was a positive correlation between CAT(x₁₅) and biomass(y) (Table 3).

The relationship between the physiological/biochemical indicators in the biomass accumulation mechanism model was consistent with the results of the Lasso model analysis. Our analysis not only illustrated that the biomass accumulation was related to some key physiological/biochemical indicators, but also the characteristics of these enzymes were the same as those of Lasso analysis. For example, GS(x₇) and GDH(x₁₀), which had positive coefficients in the Lasso model, were also present in the routes of C/N metabolism to promote biomass accumulation, and GLS(x₁₁) and POD(x₁₃), which had negative coefficients in the Lasso model, were also present in the routes of defense to consume the C and N and affect the biomass accumulation (Lea and Morot-Gaudry, 2001). Furthermore, although CAT(x₁₅) consumed C and N on the defense route, its negative synergistic effect with POD(x₁₃), or competition for substrates, could explain this phenomenon. Thereby, perhaps increasing GS(x₇), GDH(x₁₀), and CAT(x₁₅) activities, or decreasing GLS(x₁₁) and POD(x₁₃) activities, would be beneficial to promote biomass accumulation of S.miltiorrhiza.

Finally, we established a functional equation between the 6 key physiological/biochemical indicators and the biomass based on the Cobb-Douglas economic model. Although Cobb-Douglas economic model is an economic mathematical model used to predict the production of national and regional industrial systems or large enterprises and analyze the ways to develop production, this study skillfully uses it to establish the relationship between physiological/biochemical indicators and biomass. Interestingly, the correlation between the physiological/biochemical indicators and the biomass in this equation was the same as that obtained by Lasso regression. That is GS(x₇), GDH(x₁₁), and CAT(x₁₅) were positively correlated with the biomass, and GLS(x₁₀), POD(x₁₂), and soluble protein(x₁₃) were negatively correlated with biomass accumulation. The biomass predicted from the physiological/biochemical indicators of this equation had less error. In the artificial cultivation of S.miltiorrhiza, only the values of the key physiological/biochemical indicators of the above-ground parts are needed to obtain the corresponding biomass.

5. Conclusion

This study provides a rigorous logical reasoning process in terms of the selection of factors affecting biomass accumulation, the screening of key factors, and the establishment and validation of regression models. The method is applicable not only to metabolic engineering but also to phenomena with similar mechanistic features, such as the relationship between the activities and soil environment, and the self-organization of microbial communities. In addition, the generalized application of Lasso regression and the Cobb-Douglas production model used in this study provides a powerful tool for a comprehensive and systematic study of growth and active ingredient synthesis.

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Author Contributions

ZY and KW obtained funding. BL, HJ, NL, and ZW designed the experiments and analyzed the data, with the physiology support from YM and LQ generated the structure diagram of biomass accumulation mechanism, WW modified Lasso and nonlinear regression algorithm. KW completed all the program design and graphics drawing. KW and ZY wrote this study with the contributions of all the authors, who also approved the final version prior to submission.

Funding

This study was supported by Sichuan Science and Technology Program (2021 YFS0045), Special Scientific and Technological Research Project of Sichuan Province Administration of Traditional Chinese Medicine (Project 2018QN013), State Key Laboratory of Characteristic Chinese Medicine Resources in Southwest China, School of Pharmacy, Chengdu University of Traditional Chinese Medicine, Chengdu, China.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Abbreviations

C, carbon; CAT, catalase; GDH, glutamate dehydrogenase; GLS, glutaminase; GOGAT, glutamate synthase; GS, glutamine synthetase; MDA, malondialdehyde; N, nitrogen; NiR, nitrite reductase; NR, nitrate reductase; PAL, phenylalanine ammonialyase; PGA, 3- phosphoglyceric acid; POD, peroxidases; Pro, proline; NADPH, reducing agent; S.miltiorrhiza, Salvia miltiorrhiza; SOD, superoxide dismutase; SPS, synthase; SS, sucrose synthase.

References

Abogadallah, G.. (2010). Antioxidative defense under salt stress. Plant Signal. Behav. 5, 369–374. doi: 10.4161/psb.5.4.10873

PubMed Abstract | CrossRef Full Text | Google Scholar

Carter, S.. (2012). “On the Cobb-Douglas and all that... ”: the solow-simon correspondence over the aggregate neoclassical production function. J. Post Keynesian Econ. 34, 255–273. doi: 10.2753/PKE0160-3477340204

CrossRef Full Text | Google Scholar

Chang, W.. (2013). R Graphics Cookjournal. Sebastopol: O'Reilly Media, Inc.

Deane-Drummonda, C. E., Clarksona, D. T., and Johnsonb, C. B. (1979). Effect of shoot removal and malate on the activity of nitrate reductase assayed in vivo in barley roots (Hordeum Vulgare cv. midas). Plant Physiol. 64, 660–662. doi: 10.1104/pp.64.4.660

PubMed Abstract | CrossRef Full Text | Google Scholar

Efron, B., Hastie, T., Johnstone, I., Tibshirani, R., Ishwaran, H., Knight, K., et al. (2004). Least angle regression. Ann. Stat. 32, 407–499. doi: 10.1214/009053604000000067

PubMed Abstract | CrossRef Full Text | Google Scholar

Elsharkawy, M., Shimizu, M., Takahashi, H., and Hyakumachi, M. (2012). The plant growth-promoting fungus Fusarium equiseti and the arbuscular mycorrhizal fungus Glomus mosseae induce systemic resistance against Cucumber mosaic virus in cucumber plants. Plant Soil 361, 397–409. doi: 10.1007/s11104-012-1255-y

CrossRef Full Text | Google Scholar

Farooq, A., Bukhari, S. A., Akram, N. A., Ashraf, M., Wijaya, L., Alyemeni, M. N., et al. (2020). Exogenously applied ascorbic acid-mediated changes in osmoprotection and oxidative defense system enhanced water stress tolerance in different cultivars of saower (carthamus tinctorious l.). Plants 9:104. doi: 10.3390/plants9010104

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, J. (2006). The Experimental Guidance of Plant Physiology. Beijing: Higher Education Press.

Gonzíalez-Moro, M., Gonzíalez-Moro, I., Estavillo, J. M., Aparicio-Tejo, P. M., Marino, D., et al. (2021). A multi-species analysis defines anaplerotic enzymes and amides as metabolic markers for ammonium nutrition. Fronti. Plant Sci. 11:2250. doi: 10.3389/fpls.2020.632285

PubMed Abstract | CrossRef Full Text | Google Scholar

Haddidi, I., Duc, N. H., Tonk, S., Rápó, E., and Posta, K. (2020). Defense enzymes in mycorrhizal tomato plants exposed to combined drought and heat stresses. Agronomy 10:1657. doi: 10.3390/agronomy10111657

CrossRef Full Text | Google Scholar

Hesheng, L. (2000). Principles and Techniques of Plant Physiological and Biochemical Experiments. Beijing: Higher Education Press.

Islam, M., Matsui, T., and Yoshida, Y. (1996). Carbohydrate content and the activities of sucrose synthase, sucrose phosphate synthase and acid invertase in different tomato cultivars during fruit development. Sci. Horticult. 65, 125–136. doi: 10.1016/0304-4238(95)00871-3

CrossRef Full Text | Google Scholar

Jaafar, H. Z. E., Ibrahim, M. H., and Mohamad Fakri, N. F. (2012). Impact of soil field water capacity on secondary metabolites, phenylalanine ammonia-lyase (PAL), Maliondialdehyde (MDA) and photosynthetic responses of malaysian kacip fatimah (labisia pumila benth). Molecules 17, 7305–7322. doi: 10.3390/molecules17067305

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaur, S., Gupta, A. K., and Kaur, N. (2006). Effect of hydro- and osmopriming of chickpea (Cicer arietinum L.) seeds on enzymes of sucrose and nitrogen metabolism in nodules. Plant Growth Regulat. 49, 177–182. doi: 10.1007/s10725-006-9103-9

CrossRef Full Text | Google Scholar

Lan, Y., Liu, M., Yan, Z., Xie, H., Shen, X., Zhang, L., and Lin, C. (2016). Comparison of fast breeding and regeneration seedling traits of salvia miltiorrhiza divinorum from different geographical sources. Jiangsu Agric. Sci. 44, 103–107. doi: 10.21273/HORTSCI13072-18

CrossRef Full Text

Lantz, B. (2013). Machine Learning With R. Birmingham: Packt Publishing Ltd.

Google Scholar

Lea, P. J., and Morot-Gaudry, J.-F. (2001). Plant Nitrogen. Berlin; Heidelberg: Springer.

Google Scholar

Li, H. S. (2000). Principles and Techniques of Plant Physiological Biochemical Experiment. Beijing: Higher Education Press.

Li, Y., Yang, Y., Hu, Y., Liu, H., He, M., Yang, Z., et al. (2019). DELLA and EDS1 form a feedback regulatory module to fine-tune plant growth-defense tradeoff in arabidopsis. Mol. Plant 12, 1485–1498. doi: 10.1016/j.molp.2019.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Z. (2018). Salvia miltiorrhiza Burge (Danshen): a golden herbal medicine in cardiovascular therapeutics. Acta Pharmacol Sin. 39, 802–824. doi: 10.1038/aps.2017.193

PubMed Abstract | CrossRef Full Text | Google Scholar

Liang, B., Ma, C., Zhang, Z., Wei, Z., Gao, T., Zhao, Q., et al. (2018). Long-term exogenous application of melatonin improves nutrient uptake fluxes in apple plants under moderate drought stress. Environ. Exp. Bot. 155, 650–661. doi: 10.1016/j.envexpbot.2018.08.016

CrossRef Full Text | Google Scholar

Luo, H. (2018). Biomass accumulation, photosynthetic traits and root development of cotton as affected by Irrigation and nitrogen-fertilization. Front. Plant Sci. 9:14. doi: 10.3389/fpls.2018.00173

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, S., Zhang, D., Lou, H., Sun, L., and Ji, J. (2016). Evaluation of the anti-inflammatory activities of tanshinones isolated from < i>salvia miltiorrhiza < /i> var. alba roots in thp-1 macrophages. J. Ethnopharmacol. 199, 193–199. doi: 10.1016/j.jep.2016.05.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Moriguchi, T., Abe, K., Sanada, T., and Shohei, Y. (1992). Levels and role of sucrose synthase, sucrose-phosphate synthase, and acid invertase in sucrose accumulation in fruit of asian pear. J. Am. Soc. Hortic. Sci. 117, 274–278. doi: 10.21273/JASHS.117.2.274

CrossRef Full Text | Google Scholar

Noor, E., Eden, E., Milo, R., and Alon, U. (2010). Central carbon metabolism as a minimal biochemical walk between precursors for biomass and energy. Mol. Cell 39, 809–820. doi: 10.1016/j.molcel.2010.08.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Polle, A., Otter, T., and Seifert, F. (1994). Apoplastic peroxidases and lignification in needles of norway spruce (Picea abies L.). Plant Physiol. 106, 53–60. doi: 10.1104/pp.106.1.53

PubMed Abstract | CrossRef Full Text | Google Scholar

Salehin, M., Li, B., Tang, M., Katz, E., Song, L., Ecker, J. R., et al. (2019). Auxin-sensitive aux/iaa proteins mediate drought tolerance in arabidopsis by regulating glucosinolate levels. Nat. Commun. 10, 1–9. doi: 10.1038/s41467-019-12002-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Sarker, U., and Oba, S. (2020). The response of salinity stress-induced A. tricolor to growth, anatomy, physiology, non-enzymatic and enzymatic antioxidants. Front. Plant Sci. 11:559876. doi: 10.3389/fpls.2020.559876

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, M., Huang, F., Deng, C., Wang, Y., and Kai, G. (2019). Bioactivities, biosynthesis and biotechnological production of phenolic acids in salvia miltiorrhiza. Crit. Rev. Food Sci. Nutr. 59, 953–964. doi: 10.1080/10408398.2018.1474170

PubMed Abstract | CrossRef Full Text | Google Scholar

Strain, H. H., and Svec, W. A. (1966). “Extraction separation, estimation and isolation of the chlorophylls,” in and GR (eds) The Chlorophylls, eds V. P. Vernon, and G. R. Svec (New York, NY: Academic Press). p. 21–66.

Google Scholar

Su, C., Ming, Q.-L., Rahman, K., Han, T., and Qin, L.-P. (2015). Salvia miltiorrhiza: traditional medicinal uses, chemistry, and pharmacology. Chin. J. Natural Med. 13, 163–182. doi: 10.1016/S1875-5364(15)30002-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Verma, A. K., Upadhyay, S. K., Verma, P. C., Solomon, S., and Singh, S. B. (2011). Functional analysis of sucrose phosphate synthase (sps) and sucrose synthase (ss) in sugarcane (saccharum) cultivars. Plant Biol. 13, 325–332. doi: 10.1111/j.1438-8677.2010.00379.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Ma, R., Liu, C., Liu, H., Zhu, R., Guo, S., et al. (2017). Salvia miltiorrhiza: a potential red light to the development of cardiovascular diseases. Curr. Pharm. Des. 23, 1077–1097. doi: 10.2174/1381612822666161010105242

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, Y., Xiong, S., Zhang, Z., Meng, X., Wang, L., Zhang, X., et al. (2021). Localization, gene expression, and functions of glutamine synthetase isozymes in wheat grain (triticum aestivum l.). Front. Plant Sci. 12:37. doi: 10.3389/fpls.2021.580405

PubMed Abstract | CrossRef Full Text | Google Scholar

Ya-Li, L. (2018). Effects of endophytic fungi on plant growth and metabolism of salvia miltiorrhiza bge (MA thesis). Chengdu University of Traditional Chinese Medicine.

Yoneyama, Tadakatsu, Fujimori, Tamaki, Yanagisawa, Shuichi, Hase, Toshiharu, Suzuki, and Akira (2015). N-15 tracing studies on in vitro reactions of ferredoxin-dependent nitrite reductase and glutamate synthase using reconstituted electron donation systems. Plant Cell Physiol. 56, 1154–1161. doi: 10.1093/pcp/pcv039

PubMed Abstract | CrossRef Full Text | Google Scholar

Zaheer, I. E., Ali, S., Saleem, M. H., Imran, M., Alnusairi, G. S. H., Alharbi, B. M., et al. (2020). Role of iron-lysine on morpho-physiological traits and combating chromium toxicity in rapeseed (brassica napus l.) plants irrigated with different levels of tannery wastewater. Plant Physiol. Biochem. 155, 70–84. doi: 10.1016/j.plaphy.2020.07.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Khan, S., Sun, M., Gao, Z. Q., Liang, Y. F., Yang, Q. S., et al. (2018). Coordinated improvement of grain yield and protein content in dryland wheat by subsoiling and optimum planting density. Appl. Ecol. Environ. Res. 16, 7847–7866. doi: 10.15666/aeer/1606_78477866

CrossRef Full Text | Google Scholar

Zhang, Z. (1990). The Experimental Guidance of Plant Physiology. Beijing: Higher Education Press.

Zhao, Y., Wei, Z., and Ma, F. (2003). Influence of ammoniacal nitrogen on sucrose synthase and sucrose phosphate synthase in sugar beet. Sugar Crops China 3, 1–5.

Google Scholar

Zhou, B., Zhao, G., Zhu, Y., Chen, X., Zhang, N., Yang, J., et al. (2021). Protective effects of nicotinamide riboside on H₂O₂-induced oxidative damage in lens epithelial cells. Curr. Eye Res. 46, 961–970. doi: 10.1080/02713683.2020.1855662

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, G., Deng, X., and Zuo, W. (1983). Determination of free proline in plants. Plant Physiol. J. 1, 37–39.

Google Scholar