Next Article in Journal
Is the High Proportion of Males in a Population of the Self-Incompatible Fraxinus platypoda (Oleaceae) Indicative of True Androdioecy or Cryptic-Dioecy?
Next Article in Special Issue
Nutrient Concentrations Induced Abiotic Stresses to Sweet Pepper Seedlings in Hydroponic Culture
Previous Article in Journal
Features of Profiles of Biologically Active Compounds of Primary and Secondary Metabolism of Lines from VIR Flax Genetic Collection, Contrasting in Size and Color of Seeds
Previous Article in Special Issue
Bioactive Compounds and Antifungal Activity of Leaves and Fruits Methanolic Extracts of Ziziphus spina-christi L.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 11
Issue 6
10.3390/plants11060752
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metabolite Profiles Provide Insights into Underlying Mechanism in Bupleurum (Apiaceae) in Response to Three Levels of Phosphorus Fertilization

by
Jialin Sun
1,2,
Zejia Duan
1,
Ye Zhang
1,
Sisi Cao
3,
Zhonghua Tang
1,* and
Ann Abozeid
4,*
1
College of Chemistry, Chemical Engineering and Resource Utilization, Northeast Forestry University, Harbin 150040, China
2
Biological Science and Technology Department, Heilongjiang Vocational College for Nationalities, Harbin 150066, China
3
Medical Department, Harbin Vocational & Technical College, Harbin 150040, China
4
Botany and Microbiology Department, Faculty of Science, Menoufia University, Shebin Elkoom 32511, Egypt
*
Authors to whom correspondence should be addressed.
Plants 2022, 11(6), 752; https://doi.org/10.3390/plants11060752
Submission received: 27 January 2022 / Revised: 7 March 2022 / Accepted: 8 March 2022 / Published: 11 March 2022
(This article belongs to the Special Issue The Role of Fertilizers in Boosting Plant Adaptation to Biotic and Abiotic Stresses)
Browse Figures
Versions Notes

Abstract

:
Phosphorus (P) deficiency affects plant yield and quality, yet at the same time, excessive phosphorus application does not necessarily promote the growth of plants. How to maintain a balance between biomass accumulation and phosphorus application is a problem. Therefore, the purpose of this research was to explore the relationship between yield and quality of Bupleurum and phosphorus fertilization, based on three phosphorus fertilization levels (20 kg∙ha−1; 10 kg∙ha−1; and 0 kg∙ha−1). We adopted gas chromatography-mass spectrometry to assess the response of primary metabolites of different plant tissues (flowers, main shoots, lateral shoots and roots) to phosphorus fertilization. At the same time, high-performance liquid chromatography was used to quantify saikosaponin A and saikosaponin D, the main active ingredients of Bupleurum. Our research showed that low phosphorus level application has a positive impact on the yield and quality of Bupleurum, especially the above-ground parts increasing the fresh weight of flowers and lateral shoots and the length of main shoots, and moreover, increasing the saikosaponins content in all above-ground parts while decreasing the content in roots which show no significance increase in fresh weight and length. However, high phosphorus level showed a negative impact as it decreases the saikosaponins content significantly in flowers and roots. Furthermore, phosphorus application changed the proportion of saikosaponins, promoting the content of saikosaponin A and inhibiting the content of saikosaponin D in most organs of Bupleurum. Therefore, we can say that high phosphorus application is not preferable to the yield and quality of Bupleurum. To identify the metabolic pathways and special key metabolites, a total of 73 metabolites were discovered, and four differential metabolites—ether, glycerol, chlorogenic and L-rhamnose—were considered to be the key metabolites of Bupleurum’s response to phosphorus fertilization. Furthermore, Bupleurum’s response to phosphorus fertilization was mainly related to metabolic pathways, such as starch and sucrose metabolism and galactose metabolism. Under the phosphorus level, the content of sugars, organic acids and their derivatives, polyols and their derivatives and alkyl were upregulated in flowers. Furthermore, the contents of compounds in the main shoot and lateral shoots showed the same upward trend, except glycosides and polyols and their derivatives.

Graphical Abstract

1. Introduction

As it is known, the compositions and contents of the plant largely depend on genetics, growth stage, environment and cultivation, and so on [1]. Among numerous factors, the availability of plant nutrients has received more attention. Phosphorus (P) is an imperative nutrition element for plant growth and development, which affects energy metabolism, nucleic acid synthesis, glycolysis in plants and vital physiological processes [2,3,4]. As the available soil phosphorus is low, crop yield and growth is limited [5]. In order to cope with phosphorus deficiency, plants have developed flexible strategies through complex biological processes, such as changing metabolic processes, regulating root system structure, accumulating and secreting organic acids, making plant roots and mycorrhizas fungal symbiosis [6,7,8]. Changes in root system structure in reaction to insufficient phosphorus are usually involved in phosphorus uptake, while changes of metabolic processes are related to phosphorus use efficiency [9]. In addition, rational phosphorus fertilization and improved phosphorus use efficiency are effective methods to ensure sustainable development of agricultural production [10,11].
Chinese medicinal materials, as a typical complex chemical system, have multi-component, multi-target and synergistic effect characteristics. It is indispensable to comprehensively evaluate the quality of the chemical components of Chinese herbal medicines. Metabolomics has opened up a new way to qualitatively and quantitatively analyze all low-molecular-weight metabolites of a given organism or cell over a specific physiological period, and detect changes in metabolites in response to endogenous and exogenous stimuli [12,13,14]. It not only has the value of prediction, diagnosis and drug application in disease research, but has also been widely used in the characterization of natural products, evaluation of therapeutic effects and analysis of potential mechanisms [15]. During the last years, this technique has been successfully applied to focus on the mechanism of secondary metabolic pathways, such as identifying the potential biomarkers to control product quality and contributing to improve disease resistance [16,17]. In addition, it has been widely used in research on active components and quality evaluation of traditional Chinese medicine, which could help to obtain a complete picture of the metabolic reactions of medicinal plants [18,19]. Meantime, metabolomics analysis has been used to reveal the mechanism of low nutrient tolerance in plants, so as to supply a foundation for the effective development of plant quality control and germplasm improvement. The typical examples include the study of the effect of nitrogen and phosphorus fertilization on metabolome of goji fruit [20,21]. The purpose of metabolomics comprehensively analyze all metabolites in biological samples, and demonstrate huge potential in clarifying the metabolic processes of plants [22]. A few metabolomics studies have been conducted in medicinal plants, such as Panax ginseng [23] and Rehmannia glutinosa [24], as well as Astragalus mongholicus [25]. In contrast, the metabolomic profile of Bupleurum has not been investigated extensively.
Bupleurum is a perennial herb belonging to the family Apiaceae, with about 200 species all over the world [26]. It is grown diffusely worldwide, including in temperate zones of the Northern Hemisphere, Europe, North Africa, North America and Asia [27]. Bupleuri Radix, the root of Bupleurum scorzonerifolium Willd. or Bupleurum chinense DC., has a long history of use and important economic value in China [28,29]. Due to its remarkable and stable clinical efficacy, it is included in the 2020 Chinese Pharmacopoeia [30]. Active compounds isolated from the roots, such as volatile oils, flavonoids, saikosaponins and polysaccharides [31,32], possess various pharmacological functions, including common cough, fever and influenza, hepatitis, malaria menoxenia, thoracolumbar pain, and uterine as well as rectal prolapse [33,34,35,36]. Studies have shown that roots contain saikosaponin A, C and D, but saikosaponin A and D are the main active ingredients [37].
Recently, gas chromatography (GC) or liquid (LC) chromatography combined use with mass spectrometry were major technical methods for metabonomic analysis [38]. Due to sensitive selective mass detection, high separation ability and reproducible retention times, gas chromatography-mass spectrometry (GC-MS) is a very suitable technique for metabolomics analysis [39,40]. GC-MS could be used to determine a variety of primary metabolites, such as carbohydrates, organic acids, fatty acids, amino acids and so on [41]. Liquid chromatography-mass spectrometry (LC-MS) could be used to analyze secondary metabolites, such as phenylpropanoids, lignins, terpenes and phenolic compounds [42]. The combined analysis of GC-MS and LC-MS could significantly increase the range of metabolites, and then comprehensively summarized and analyzed the changes of metabolites in detail [43].
Despite many studies have confirmed biological activities and pharmacological properties of Bupleurum, the relationship between the metabolome and phosphorus fertilization of Bupleurum is still unclear. At this study, we applied gas chromatography-mass spectrometry to compare the metabolite profiles of different tissues of Bupleurum (roots, main shoots, lateral shoots, and flowers) under different levels of phosphorus fertilization. Our aim was to investigate the following: (1) effect of phosphorus fertilization on the mass accumulation and partition in Bupleurum; (2) changes in the metabolite profiles of different tissue under different level of phosphorus fertilization; and (3) key metabolites associated with phosphorus fertilization to explain key points of biosynthesis regulation. Therefore, studying overall metabolic changes of Bupleurum under phosphorus fertilization will provide a basis for nutritional regulation.

2. Results

2.1. Comparative Analysis of Quality and Traits of Bupleurum under Different Phosphorus Level

In this experiment, the whole Bupleurum individuals growing under three levels of phosphorus fertilization—low-phosphorus (LP), control-phosphorus (CP) and high-phosphorus (HP)—were divided into four parts (Figure 1), including flowers (F), main shoots (MS), lateral shoots (LS) and roots (R). We determined the fresh weight and length for the four parts to find the influence brought by different phosphorus fertilization levels. The result revealed that different levels of phosphorus fertilization showed a significant effect on the yield and biomass of Bupleurum (Figure 2). In details, the fresh weight of Bupleurum’s lateral shoots and flowers significantly increased (p ≤ 0.01) under LP level compared with the CP level, and a further weight increase was observed under HP level (Figure 2A). However, no significant increase in fresh weigh of main stem and roots was observed with phosphorus application. Considering the organs’ length (Figure 2B), the situation was reversed, and the increase in the length of lateral shoots and flowers was not significant (p ≥ 0.05), while a significant increase in the main stems length was observed under phosphorus fertilization.

2.2. Total Saikosaponins Content Accumulated in Bupleurum Different Tissues under Three Levels of Phosphorus Fertilization

Initially, this study exploited a HPLC method for simultaneous quantification of saikosaponin A and saikosaponin D. Simultaneously, we also measured the total content of saikosaponins in the different Bupleurum tissues under the three levels of phosphorus fertilization. We found that under control phosphorus level, the total saikosaponin content was the highest within roots, phosphorus application decreased the content in roots, while it increased in the above-ground organs. The decrease in roots’ total saikosaponin was very significant under a high phosphorus level. Moreover, comparing with low phosphorus level, high phosphorus shows an increasing in saikosaponin content in only the main stem and lateral stem, while the saikosaponins in flowers and roots decreased significantly, especially in roots (Figure 3A,B). When the total saikosaponins content was fixed and we compared among Saikosaponin A and Saikosaponin D, we found that when there was no phosphorus, and the Saikosaponin D percentage was higher than the percentage of saikosaponins A in all organs. With the increase of phosphorus fertilizer content, the proportion of Saikosaponin D gradually decreased, and Saikosaponins A increased instead (Figure 3C).

2.3. Overview of the Metabolites Profiles in Response to Three Levels of Phosphorus Fertilization

In order to understand more clearly the changes of metabolites in different phosphorus levels, GC-MS was used to the samples under phosphorus stress. A total of 73 primary metabolites were identified, which could be further divided into 9 groups, including 24 organic acids, 14 sugars, 12 alkyls, 5 polyols, 4 amino acids, 3 lipids, 3 glycosides, 2 esters and 6 others (Figure 4A). The clustering heat map was performed to assess the differences among the Bupleurum metabolites profiles in different tissues under three phosphorus levels. Two main clusters were obtained according to the relative differences of Bupleurum accumulation patterns (Figure 4B). The results indicates that different tissues have similar metabolic patterns. For example, with phosphorus application, metabolites such as esters, lipids, organic acids, sugar and glycosides accumulate preferentially within flowers and main shoots, while their overall levels within roots and lateral shoots are relatively low. However, roots and lateral shoots have substantially higher levels of sugar and glycosides under high phosphorus levels (Figure 4B). These specific metabolites likely reflect that the accumulation of metabolites during Bupleurum development is tissue-specific.
The PLS-DA could evaluate discrepancy between various phosphorus levels of multiple variables through some principal components (Figure 5). Obviously, the PLS-DA results showed that the samples of different phosphorus levels were clearly separated in all tissues, which indicated that the metabolic differences were significant, corresponding to the results of cluster analysis and heat map. The first partial least-square method (PLS-DA1) could account for approximately 27.8%, 20.4%, 28.1% and 20.7% of the variation in four tissues (flower, lateral shoot, main shoot, root) separately. In terms of composition, it can be found that for the flowers and main shoots, HP, LP and CP are completely separated. However, the lateral shoot and root showed that CP and LP groups were significantly separated. The second partial least-square method (PLS-DA2) could explain 24%, 19.2%, 15.6% and 13.2% of the features of the original data set.

2.4. Metabolite Profiling of Bupleurum under Phosphorus Fertilization in KEGG Enrichment Analysis and Volcanic Map

The metabolic profiles of Bupleurum were analyzed using a fold change value ≤ 0.05 and a p value ≥ 1. In most cases, the active ingredients of Bupleurum were closely related to the level of phosphorus fertilization (Table 1). Comparing HP with CP, significant changes occurred in 22 metabolites, including 18 up-regulated metabolites (8 F, 5 LS, 4 MS, 1 R) and 4 down-regulated metabolites (1 LS, 1 MS, 2 R) in the root. About 26 metabolites had changes when comparing LP to CP, with 20 (12 F, 4 LS, 4 MS) down-regulated and 6 up-regulated metabolites (1 F, 1 LS, 4 MS). Compared with HP, 17 metabolites were significantly changed in the LP group, including 8 (2 F, 2 LS, 3 MS, 1 R) up-regulated and 9 (5 F, 2 LS, 2 MS) down-regulated metabolites.
The Kyoto Encyclopedia of Genes and Genomes (KEGG) and enrichment results show that the effects of different levels of phosphorus fertilization (LP vs. CP, LP vs. HP, HP vs. CP) on metabolic pathways mainly involve the starch and sucrose metabolism and galactose metabolism (Figure 6). In HP vs. CP, various metabolites were mostly focused on glycerolipid metabolism, starch and sucrose metabolism, galactose metabolism, pentose phosphate pathway, fatty acid metabolism, glyoxylate and dicarboxylate metabolism, and even insulin signaling pathway (Figure 6 A,B). For HP vs. LP, the difference of KEGG enrichment classification was also related to the following pathways: starch and sucrose metabolism, glyoxylate and dicarboxylate metabolism, galactose metabolism, pentose phosphate pathway and insulin signaling pathway (Figure 6 C,D). For LP vs. CP, metabolic pathways of differential metabolites are involved with starch and sucrose metabolism, insulin signaling pathway, amino sugar and nucleotide sugar metabolism and galactose metabolism (Figure 6 E,F). The collective metabolic pathways are starch and sucrose metabolism and galactose metabolism, which are mainly related to sugars. The appraisal of potential metabolites involved into phosphorus stress might contribute to the growth and development of Bupleurum. We performed a Venn diagram (Figure 6G) to depict shared metabolites of different expressions among HP vs. CP, CP vs. LP and HP vs. CP. We found that the abundances of ether, glycerol, chlorogenic or L-rhamnose had an obvious change in the response to phosphorus fertilization; therefore, these four over-lapping metabolites could be regarded as key metabolites. Moreover, 9, 4 and 4 differential metabolites existed solely in (LP vs. CP) vs. (HP vs. LP), (LP vs. CP) vs. (HP vs. CP) and (HP vs. LP) vs. (HP vs. CP), respectively.

2.5. Metabolic Network Diagram and Potential Metabolites in Bupleurm under Three Levels of Phosphorus Fertilization

On the basis of the KEGG annotation and enrichment data, significant metabolites were mapped to C-metabolism, lipids metabolism pathways to guarantee clear changes in the metabolic regulation of phosphorus (Figure 7). The overlap of differentially expressed metabolites was observed, suggesting the partial similarity of mechanisms of Bupleurum in response to phosphorus fertilization. The metabolic network diagram went on to verify the above hypothesis. Responses to phosphorus fertilization of Bupleurum were dynamic, and included intricate pathways. Therefore, phosphorus fertilization may activate some key physiological and metabolic activities that lead to the growth and development of Bupleurum.

3. Discussion

3.1. Influences of Phosphorus Fertilization on Bupleurum Yield, Quality and Saikosaponins Contents

Phosphorus deficiency severely restricts the growth and productivity of plants. Phosphorus can make stems and branches more rigid, cause early flowering and fruiting, and can improve the disease resistance and stress resistance of plants [44,45]. Plant Fresh weight and length are two important evaluation indexes of Bupleurum quality. Our results revealed that the fresh weight of flowers and lateral shoots showed significant increase under HP levels, compared with CP treatments. It is said that the long-term phosphorus deficiency could arrest proliferation of lateral shoots [46]. Moreover, similar results reported that high phosphorus levels increase flower number and weight [47,48,49]. The reason for this could be that the use of phosphorus fertilization encourages the accumulation of nutrients, such as sugars, organic acids and polyols. Plant roots play an important role in phosphorus absorption. In many studies, the rhizome ratio of phosphorus-stressed plants increased compared with phosphorus-sufficient plants due to the proliferation of lateral roots in phosphorus-stressed plants, in order to promote phosphorus absorption when phosphorus availability decreases [46,50,51,52]. This could explain our results, which show that low and high phosphorus application did not show a significant increase in root weight or length (Figure 2).
Recent studies indicated that saikosaponins A and D are the main active components to which the clinical efficacy evaluate the quality of Bupleurum in the Chinese Pharmacopoeia [37]. Therefore, we studied the relationship between phosphorus fertilizer application and the content of saikosaponins A and D in Bupleurum. Although phosphorus fertilization promoted the fullness of the shoots and leaves [53], high phosphorus levels (HP) have a limited or even negative impact on saikosaponins content within Bupleurum, decreasing the content in roots and flowers. Low phosphorus levels (LP) can promote plant growth and maintain the content of active ingredients.

3.2. Influences of Phosphorus Fertilization on Bupleurum Metabolites

3.2.1. Effect of Phosphorus Fertilization on Sugar Metabolism of Bupleurum

As an indispensable energy source, sugars have a significant effect not only on plant growth, but also on biomass production. Moreover, the stress response and growth process of plants are also regulated by sugars [54]. There are two main ways for plants to respond to low phosphorus, including changes in root system architecture (RSA) and the modification of metabolic processes [55]. We found that the variation of Bupleurum metabolites significantly affected the accumulation of chemical components in different parts under different levels of phosphorus fertilization. The specific performance is that the content of sugars, organic acids and derivatives, polyols and derivatives and alkyl were upregulated in flowers. Furthermore, the content of these compounds in main and lateral shoots showed the same upward trend, except glycosides and polyols and derivatives. Previous studies have confirmed the increased accumulation of sugars, such as d-glucose, d-fructose and d-mannose in both main and lateral shoots [56]. The accumulation of sugars was detected in barley [57] and maize [58], in the absence of phosphorus. Therefore, sugar accumulation under low phosphorus conditions might be a phosphorus efficiency mechanism. However, Muller treated up in with phosphorus deficiency [59], and observed that glucose, fructose and sucrose in shoots, after 14 days of treatment, were reduced, while there was no significant change in the amount of sugar levels in roots and shoots with 22 days of similar treatment. As for roots, we observed that all components presented the same variation: having first increased and then decreased. Comparing CP and HP, however, the final result was up-regulated. Overall, the least accumulation of compounds was in lateral shoots and roots, which were mainly concentrated in CP and LP treatment groups. This difference can be attributed to differences in farming methods and stages of harvest. Previous results confirmed that the accumulation of sugar began in the early stage, and the accumulation speed accelerated in the mid stage until the maturity stage, and then tended to be stable [60,61]. Furthermore, different levels of phosphorus deficiency lead to different results. In fruit, sugar accumulation was largely determined by sucrose metabolism [62,63,64], indicating that the application of phosphorus fertilization can enhance fruit’s sink strength, which may increase the amount of sugar that leaves deliver to the fruit [44].

3.2.2. Effect of Phosphorus Fertilization on Acid Metabolism of Bupleurum

Phosphoenolpyruvic acid (PEP) was closely related to the metabolic precursors of most organic acids. The effects of phosphorus deficiency on the secretion of organic acids in plant roots were as follows: short-term treatment increased the content of organic acids, while long-term phosphorus deficiency reduced the release of organic acids in the rhizosphere [51]. According to the current study, malic acid and citric acid content secreted by Bupleurum roots at the late stage of P deficiency increased. In plant tissues, the synthesis of organic acids was largely influenced by the citric acid cycle and glyoxylic acid cycle [56]. Comparing HP and LP, the level of malate and citrate were up-regulated in roots. However, both were down-regulated compared to the control group. Lin [65] supplied different concentrations of phosphate fertilizer on uniform tea, and the concentration of root malate and citrate did go down as the phosphorus supply increased from 0 to 40 μM, and then slightly increased with further increase the content of phosphate. The trend was same to the result we mentioned above. Additionally, the variation of content of malate and citrate showed in heat map had the same trends. Above all, we found that LP and HP treatment decreased the content of malate and citrate, which could be attributed to deceleration of TCA-cycle in roots. Inconsistent results suggest that the regulation of rhizosphere metabolism is not unique, and different phosphorus concentration may lead to contradictory results. The content of succinic acid and trans-aconitate in the rhizosphere increased when sufficient phosphorus fertilization was applied to Z. mays [66]. A similar reaction was found in G. max, which was manifested by an increase in the proportion of citric acid in root exudates [67], and the contents of fumaric acid and succinic acid in Hordeum vulgare roots increased significantly [57].
According to the content of phosphorus fertilization, a majority of research compared +P levels and -P levels to discuss the physiological and biochemical effects of phosphorus fertilizer on plants [10,44,46]. In addition, several phosphorus levels were set to discuss the optimal dosage. For example, Mohidin [68] set five phosphorus levels in combination with nitrogen and potassium to optimize oil palm seedlings in solution culture: 15 mg∙L−1, 30 mg∙L−1, 60 mg∙L−1, 90 mg∙L−1, 120 mg∙L−1. Additionally, according to the Dangi study [47] on the effect of phosphorus on the growth and flowering of Marigol, phosphorus levels were respectively 0 kg∙ha−1, 20 kg∙ha−1, 40 kg∙ha−1, 60 kg∙ha−1, 80 kg∙ha−1, 100 kg∙ha−1, 120 kg∙ha−1, and he also found that the highest plant length was obtained at 20 kg∙ha−1—the lowest concentration of phosphorous used. At the same time, in order to impact the phosphorus fertilizer level on the yield and metabolome of goji fruit, Feng [21] set three phosphorus levels—32.5 g∙per tree−1; 65 g∙per tree−1; and 97.5 g∙per tree−1—and he concluded that the yield of goji fruits had a strong negative correlation with phosphorus levels. Although these studies concluded that lower levels of phosphorous had better impact on plant yield, the lowest level used was 20 kg∙ha−1. In our study, we used lower levels of phosphorous 20 kg∙ha−1 and 10 kg∙ha−1, in trying to reach the optimal dosage.

4. Materials and Methods

4.1. Experimental Design and Plant Source

The medicinal plant used in our study was cultivar Bupleurum scorzonerifolium Willd., which was widely cultivated in the northeast of China. The dried roots of Bupleurum chinense DC. and Bupleurum scorzonerifolium Willd. were the only two authentic sources of Chaihu [69]. Two-year-old Bupleurum seedlings were cultivated in natural environmental conditions at the research site of Lin Dian, Da Qing, Heilongjiang province, China. The chemical prosperities of the soil [70] were pH 7.78, electrical conductivity was 134 us∙cm−1; and organic matter was 4.86 g∙kg−1. In the 20182019 growing season, the average annual temperature was 4 °C, and the average annual precipitation was 417.2 mm [71].
The three phosphorus levels were distributed in a randomized-block design, and each plot measured 3 m × 1.2 m. Based on previous experience, (Ca(H2PO4)2·H2O) was used as phosphorus fertilization. Bupleurum was processed for control phosphorus (CP), low phosphorus (LP) and high phosphorus (HP) treatment, by adjusting the phosphorus concentration in the current study. The three phosphorus fertilization levels were: CP (0 kg∙ha−1); LP (10 kg∙ha−1); and HP (20 kg∙ha−1) (Table 1). Phosphorus fertilization was applied in July 2018. After the 30th day (August 2018), the Bupleurum were harvested and divided into flowers, main shoots, lateral shoots and roots, and quickly transferred to the laboratory within three hours. Each group obtained three technical replicates, and the fresh weight and length were then measured. Simultaneously, each sample was divided into two batches: one batch was dried to a constant weight in a blast oven at 42 °C and used for saikosaponins extraction, and the other batch was stored at −80 °C for sample processing of GC-MS.

4.2. Determination of the Saikosaponins Content

500 mg dried powder of Bupleurum (roots, main shoots, lateral shoots and flowers) was dissolved into 25 mL of methanol solution containing 8% ammonia water. Samples were ultrasonicated at 30 °C for 30 min, and then filtered. The filter residue was rinsed twice with 10 mL methanol, and the filtrates were combined and recovered to dryness. The residue was dissolved in methanol to a 10 mL volumetric flask. The solution was passed through a 0.45 μm filter membrane to use. An aliquot of 20 μL was injected for HPLC analysis. The content of saikosaponin A and D (Figure 8) were measured separately. A HPLC (Hitachi, Tokyo, Japan) system was equipped with L-2000 High Performance Liquid Chromatograph and L-2200 Autosampler, and a reversed phase column was adopted. According to the test requirements of saikosaponins in Chinese Pharmacopoeia, acetonitrile and pure water were used as solvents [69], and the ratio is as follows in Table 2. A Diamonsil C18 (4.6 × 250 mm, 5 μm) chromatographic column was selected, the column temperature was 25 °C, the flow rate was 0.8 mL∙min−1 and the detector was set at 210 nm.

4.3. Preparation and Extraction

60 mg various organs of Bupleurum (roots, main shoots, lateral shoots and flowers) were moved to a 1.5 mL centrifuge tube. Samples were extracted by 540 μL cold methanol (pre-cooled at −20 °C for use), and 60 μL of internal standard (L-2-chloro-phenylalanine, 0.3 mg∙mL−1, prepared in methanol) was added sequentially. After vortexing (60 Hz, 2 min) in a grinder, ultrasonic extraction was performed for 30 min. 300 μL of chloroform were added and vortexed in the grinder (20 Hz, 2 min), and then 600 μL of deionized water was again added, vortexed and extracted by ultrasonic for 30 min. At 4 °C, samples were centrifuged at 10,000 rpm for 10 min. 700 μL of supernatant was moved to a centrifuge tube, and evaporated to dryness with a rapid centrifugal concentrator. The dried residue was dissolved in 200 μL of methoxyamine pyridine solution (15 mg∙mL−1), and incubated in a shaking incubator at 37 °C for 90 min. Subsequently, 200 μL N,O−Bis (trimethylsilyl) trifluoroacetamide (BSTFA) and 40 μL Hexane were added, vortexed for 2 min and derivatized 60 min at 70 °C [23,25,72]. The solution was centrifuged at 12,000 rpm for 5 min, in order to obtain the supernatant for GC-MS analysis.

4.4. GC−MS Analysis

In this experiment, the analytical instrument was Agilent’s 7890A-5975C gas chromatography-mass spectrometer (Agilent, Santa Clara, FL, USA). 1 μL of the derivatized extract was injected into the instrument, and the sample was passed through a non-polar DB-5MS capillary. The column (30 m × 250 μm ID, J&W Scientific, Folsom, CA, USA) was separated and entered into mass spectrometry detection. High-purity helium was used as the carrier gas at a flow rate of 1.0 mL∙min−1. Program temperature [23,25,72]: 8 °C∙min−1, from 0 °C to125 °C; 4 °C∙min−1, from 125 to 210 °C; 5 °C∙min−1, from 210 to 270 °C; 10 °C∙min−1, from 270 to 305 °C; and a final maintenance at 305 °C for 3 min. The temperature of the injection port and EI source were both 260 °C, and the voltage was -70 V. The mass scanning range was m/z 50–600, and the acquisition rate was 20 spectrum∙s−1.

4.5. Multivariate Statistical Analysis

The partial least squares discriminant analysis (PLS-DA) method was used to comprehensively compare the identified metabolites. According to the definition of fold change being less than 0.05 and p value being greater than 1 in Student’s t-test analysis, the differential metabolites were selected. The Kyoto Genome and Genome Encyclopedia (KEGG) were used to analyze the pathways of metabolites. The difference in the biomass and metabolite among CP, LP and HP were calculated using a one-way analysis of variance. The Graph Pad Prism 9.0 was used for drawing related graphs in article. Using hierarchical cluster analysis, the content value of each compound was standardized to obtain a heat map of relative differences.

5. Conclusions

The findings of the current research highlighted the effect of phosphorus fertilizer on the traditional Chinese medicinal plant Bupleurum. By comparing the indexes and metabolomics analysis results of Bupleurum under three phosphorus fertilizer levels, it can be concluded that phosphorus plays a pivotal role in the growth of Bupleurum. The application of high level of phosphorus could not significantly prompt the increase in fresh weight, length and saikosaponins active metabolites, but enriched the plants with primary nutrients. Thus, we can say that although phosphorus is essential for plant growth, high phosphorus level is not preferable for maintaining the yield and quality of Bupleurum. Meanwhile, a proper amount of phosphorus could promote the yield and quality characteristics, especially in the above-ground parts of Bupleurum. Moreover, metabolomics analysis identified 73 kinds of metabolites of Bupleurum under phosphorus stress (CP, LP and HP). Under high phosphorus levels, the content of sugars, organic acids and derivatives, polyols and derivatives and alkyl were upregulated in flowers. Additionally, the contents of compound in main shoot and lateral shoots showed the same upward trend, except for glycosides, polyols and derivatives. In addition, 4 differentially accumulated metabolites in Bupleurum—ether, glycerol, chlorogenic and l-rhamnose—were identified as key metabolites under phosphorus fertilization. The identified metabolic pathways were mostly starch and sucrose metabolism and galactose metabolism, which set up a primary metabolic network diagram associated with Bupleurum. Above all, the impact of phosphorus fertilizer on Bupleurum is obvious, and it is meaningful for us to further explore the influence of phosphorus fertilizer on related Traditional Chinese Medicinal plants.

Author Contributions

Data curation, J.S.; formal analysis, J.S. and Z.D.; funding acquisition, Z.T.; investigation, Z.D. and S.C.; methodology, J.S. and A.A.; project administration, Z.T.; resources, Y.Z.; supervision, A.A. and Z.T.; writing—original draft, J.S. and A.A.; and writing—review and editing, J.S., Z.D., Y.Z., S.C., Z.D., Z.T. and A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Fundamental Research Funds for the Central Universities (2572020DY04) and the Forestry Science and Technology Promotion Project ([2019]12, 2020133122).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Thanks to Zejia Duan and Yanan Li for processing the massive data of this research, as well as the drawing and preliminary idea of the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Larbat, R.; Olsen, K.M.; Slimestad, R.; Løvdal, T.; Bénard, C.; Verheul, M.; Bourgaud, F.; Robin, C.; Lillo, C. Influence of repeated short-term nitrogen limitations on leaf phenolics metabolism in tomato. Phytochemistry 2012, 77, 119–128. [Google Scholar] [CrossRef] [PubMed]
  2. Yousaf, M.; Li, X.K.; Ren, T.; Cong, R.; Ata-Ul-Karim, S.; Shah, A.; Jamil Khan, M.; Zhang, Z.; Fahad, S.; Lu, J.W. Response of nitrogen, phosphorus and potassium fertilization on productivity and quality of winter rapeseed in Central China. Int. J. Agric. Biol. 2016, 18, 1137–1142. [Google Scholar] [CrossRef]
  3. Vance, C.P.; Uhde-Stone, C.; Allan, D.L. Phosphorus acquisition and use: Critical adaptations by plants for securing a nonrenewable resource. New Phytol. 2003, 157, 423–447. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Sulieman, S.; Tran, L.S. Phosphorus homeostasis in legume nodules as an adaptive strategy to phosphorus deficiency. Plant Sci. 2015, 239, 36–43. [Google Scholar] [CrossRef]
  5. Ferrol, N.; Azcón-Aguilar, C.; Pérez-Tienda, J. Review: Arbuscular mycorrhizas as key players in sustainable plant phosphorus acquisition: An overview on the mechanisms involved. Plant Sci. 2019, 280, 441–447. [Google Scholar] [CrossRef] [Green Version]
  6. Chiou, T.J.; Lin, S.I. Signaling network in sensing phosphate availability in plants. Annu. Rev. Plant Biol. 2011, 62, 185–206. [Google Scholar] [CrossRef] [Green Version]
  7. Crombez, H.; Motte, H.; Beeckman, T. Tackling plant phosphorus starvation by the roots. Dev. Cell 2019, 48, 599–615. [Google Scholar] [CrossRef] [Green Version]
  8. Zhang, Y.; Chen, H.; Liang, Y.; Lu, T.; Liu, Z.; Jin, X.; Hou, L.; Xu, J.; Zhao, H.; Shi, Y.; et al. Comparative transcriptomic and metabolomic analyses reveal the protective effects of silicon against low phosphorus stress in tomato plants. Plant Physiol. Biochem. 2021, 166, 78–87. [Google Scholar] [CrossRef]
  9. Niu, Y.F.; Chai, R.S.; Jin, G.L.; Wang, H.; Tang, C.X.; Zhang, Y.S. Responses of root architecture development to low phosphorus availability: A review. Ann. Bot. 2013, 112, 391–408. [Google Scholar] [CrossRef]
  10. Nguyen, V.L.; Palmer, L.; Roessner, U.; Stangoulis, J. Genotypic variation in the root and shoot metabolite profiles of wheat (Triticum aestivum L.) indicate sustained, preferential carbon allocation as a potential mechanism in phosphorus efficiency. Front. Plant Sci. 2019, 10, 995. [Google Scholar] [CrossRef] [Green Version]
  11. Mazlouzi, M.E.; Morel, C.; Robert, T.; Yan, B.F.; Mollier, A. Phosphorus uptake and partitioning in two durum wheat cultivars with contrasting biomass allocation as affected by different p supply during grain filling. Plant Soil 2020, 449, 179–192. [Google Scholar] [CrossRef]
  12. Wang, Y.; Zhou, S.; Wang, M.; Liu, S.; Hu, Y.; He, C.; Li, P.; Wan, J.B. UHPLC/Q-TOFMS-based metabolomics for the characterization of cold and hot properties of Chinese materiamedica. J. Ethnopharmacol. 2016, 179, 234–242. [Google Scholar] [CrossRef] [PubMed]
  13. Chu, H.; Zhang, A.; Han, Y.; Lu, S.; Kong, L.; Han, J.; Liu, Z.; Sun, H.; Wang, X. Metabolomics approach to explore the effects of Kai-Xin-San on Alzheimer’s disease using UPLC/ESI-Q-TOF mass spectrometry. J. Chromatogr. B Analyt. Technol. Biomed. Life Sci. 2016, 1015–1016, 50–61. [Google Scholar] [CrossRef]
  14. Goodacre, R.; Vaidyanathan, S.; Dunn, W.B.; Harrigan, G.G.; Kell, D.B. Metabolomics by numbers: Acquiring and understanding global metabolite data. Trends Biotechnol. 2004, 22, 245–252. [Google Scholar] [CrossRef] [PubMed]
  15. Zhang, A.H.; Yu, J.B.; Sun, H.; Kong, L.; Wang, X.Q.; Zhang, Q.Y.; Wang, X.J. Identifying quality-markers from Shengmai San protects against transgenic mouse model of Alzheimer’s disease using chinmedomics approach. Phytomedicine 2018, 45, 84–92. [Google Scholar] [CrossRef] [PubMed]
  16. Yang, L.; Fountain, J.C.; Ji, P.; Ni, X.; Chen, S.; Lee, R.D.; Kemerait, R.C.; Guo, B. Deciphering drought-induced metabolic responses and regulation in developing maize kernels. Plant Biotechnol. J. 2018, 16, 1616–1628. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Lu, Y.; Lam, H.; Pi, E.; Zhan, Q.; Tsai, S.; Wang, C.; Kwan, Y.; Ngai, S. Comparative metabolomics in Glycine max and Glycine soja under salt stress to reveal the phenotypes of their offspring. J. Agric. Food Chem. 2013, 61, 8711–8721. [Google Scholar] [CrossRef] [PubMed]
  18. Cao, Y.W.; Qu, R.J.; Tang, X.Q.; Sun, L.Q.; Chen, Q.Q.; Miao, Y.J. UPLC-Triple TOF-MS/MS based metabolomics approach to reveal the influence of nitrogen levels on Isatis indigotica seedling leaf. Sci. Hortic. 2020, 266, 109280. [Google Scholar] [CrossRef]
  19. Chang, X.; Zhang, J.; Li, D.; Zhou, D.; Zhang, Y.; Wang, J.; Hu, B.; Ju, A.; Ye, Z. Nontargeted metabolomics approach for the differentiation of cultivation ages of mountain cultivated ginseng leaves using UHPLC/QTOF-MS. J. Pharm. Biomed. Anal. 2017, 141, 108–122. [Google Scholar] [CrossRef]
  20. Shi, Z.; Wei, F.; Wan, R.; Li, Y.; Wang, Y.; An, W.; Qin, K.; Dai, G.; Cao, Y.; Feng, J. Impact of nitrogenfertilizer levels on metabolite profiling of the Lyciumbarbarum L. fruit. Molecules 2019, 24, 3879. [Google Scholar] [CrossRef] [Green Version]
  21. Wei, F.; Shi, Z.; Wan, R.; Li, Y.; Wang, Y.; An, W.; Qin, K.; Cao, Y.; Chen, X.; Wang, X.; et al. Impact of phosphorus fertilizer level on the yield and metabolome of goji fruit. Sci. Rep. 2020, 10, 14656. [Google Scholar] [CrossRef] [PubMed]
  22. Wang, S.; Huang, S.; Yang, J.; Li, Z.; Zhang, M.; Fang, Y.; Yang, Q.; Jin, W. Metabolite profiling of violet, white and pink flowers revealing flavonoids composition patterns in Rhododendron pulchrum Sweet. J. Biosci. 2021, 46, 3. [Google Scholar] [CrossRef] [PubMed]
  23. Liu, J.; Liu, Y.; Wang, Y.; Abozeid, A.; Zu, Y.G.; Zhang, X.N.; Tang, Z.H. GC-MS metabolomic analysis to reveal the metabolites and biological pathways involved in the developmental stages and tissue response of Panax ginseng. Molecules 2017, 22, 496. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Zhou, Y.; Shao, L.; Zhu, J.; Li, H.; Duan, H. Comparative analysis of tuberous root metabolites between cultivated and wild varieties of Rehmanniaglutinosa by widely targeted metabolomics. Sci. Rep. 2021, 11, 11460. [Google Scholar] [CrossRef]
  25. Liu, Y.; Liu, J.; Abozeid, A.; Wu, K.X.; Guo, X.R.; Mu, L.Q.; Tang, Z.H. UV-B radiation largely promoted the transformation of primary metabolites to phenols in Astragalusmongholicus seedlings. Biomolecules 2020, 10, 504. [Google Scholar] [CrossRef] [Green Version]
  26. Pimenov, M.G.; Leonov, M.V. The Genera of the Umbelliferae; Royal Botanic Gardens: Kew, UK, 1993. [Google Scholar]
  27. Neves, S.S.; Weakley, A.S.; Cox, P.B. Bupleurum gerardii All. (Apiaceae), an addition to the North American Flora, with comments onthe treatment of aliens in Floras. Castanea 2009, 73, 424–433. [Google Scholar] [CrossRef]
  28. Ren, S.; Liu, J.; Xue, Y.; Zhang, M.; Liu, Q.; Xu, J.; Zhang, Z.; Song, R. Comparative permeability of three saikosaponinsand corresponding saikogenins in Caco-2 model by a validated UHPLC-MS/MS method. J. Pharm. Anal. 2021, 11, 435–443. [Google Scholar] [CrossRef]
  29. Yang, F.; Dong, X.; Yin, X.; Wang, W.; You, L.; Ni, J. Radix Bupleuri: A review of traditional uses, botany, phytochemistry, pharmacology, and toxicology. Biomed. Res. Int. 2017, 2017, 7597596. [Google Scholar] [CrossRef] [Green Version]
  30. Xia, Z.; Liu, X.; Tong, L.; Wang, H.; Feng, M.; Xi, X.; He, P.; Qin, X. Comparison of chemical constituents of Bupleurummarginatum var. stenophyllum and Bupleurumchinense DC. using UHPLC-Q-TOF-MS based on a metabonomics approach. Biomed. Chromatogr. 2021, 35, e5133. [Google Scholar] [CrossRef]
  31. Li, D.Q.; Yue, D.D.; Liu, D.; Liu, X.G.; Song, S.J. Chemical constituents from Bupleurumchinese and their chemotaxonomic significance. Biochem. Syst. Ecol. 2019, 86, 103929. [Google Scholar] [CrossRef]
  32. Liu, X.J.; Hu, J.; Li, Z.Y.; Qin, X.M.; Zhang, L.Z.; Guo, X.Q. Species classification and quality assessment of Chaihu (Radix Bupleuri) based on high-performance liquid chromatographic fingerprint and combined chemometrics methods. Arch. Pharm. Res. 2011, 34, 961–969. [Google Scholar] [CrossRef] [PubMed]
  33. Song, X.; Ren, T.; Zheng, Z.; Lu, T.; Wang, Z.; Du, F.; Tong, H. Anti-tumor and immunomodulatoryactivitiesinduced by an alkali-extracted polysaccharide BCAP-1 from Bupleurumchinense via NF-κB signaling pathway. Int. J. Biol. Macromol. 2017, 95, 357–362. [Google Scholar] [CrossRef] [PubMed]
  34. Li, D.; Zhou, L.; Wang, D.; Wu, J.; Li, L.; Huang, X.; Liu, Q.; Wu, Y.; Lin, S.; Yang, J. Neuroprotectiveoleananetriterpenes from the roots of Bupleurumchinense. Bioorg. Med. Chem. Lett. 2016, 26, 1594–1598. [Google Scholar] [CrossRef] [PubMed]
  35. Yuan, B.; Yang, R.; Ma, Y.; Zhou, S.; Zhang, X.; Liu, Y. A systematic review of the active saikosaponinsandextracts isolated from Radix Bupleuri and their applications. Pharm. Biol. 2017, 55, 620–635. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Shin, M.R.; Park, H.J.; Seo, B.I.; Roh, S.S. New approach of medicinal herbs and sulfasalazi.ne mixture on ulcerative colitis induced by dextran sodium sulfate. World J. Gastroenterol. 2020, 26, 5272–5286. [Google Scholar] [CrossRef] [PubMed]
  37. Liang, Z.H.; Zhang, J.; Yang, G.Y.; Chen, H.B.; Zhao, Z.Z. Chemical profiling and histochemical analysis of Bupleurummarginatum roots from different growing areas of Hubei province. Acta Pharm. Sin. B 2013, 3, 193–204. [Google Scholar] [CrossRef] [Green Version]
  38. Galili, G. The aspartate-family pathway of plants: Linking production of essential amino acids with energy and stress regulation. Plant Signal. Behav. 2011, 6, 192–195. [Google Scholar] [CrossRef] [Green Version]
  39. Cui, S.; Wang, J.; Yang, L.; Wu, J.; Wang, X. Qualitative and quantitative analysis on aroma characteristics of ginseng at different ages using E-nose and GC-MS combined with chemometrics. J. Pharm. Biomed. Anal. 2015, 102, 64–77. [Google Scholar] [CrossRef]
  40. Wang, Y.; Xu, L.; Shen, H.; Wang, J.; Liu, W.; Zhu, X.; Wang, R.; Sun, X.; Liu, L. Metabolomic analysis withGC-MS to reveal potential metabolites and biological pathways involved in Pb& Cd stress response ofradish roots. Sci. Rep. 2015, 5, 18296. [Google Scholar] [CrossRef] [Green Version]
  41. Lisec, J.; Schauer, N.; Kopka, J.; Willmitzer, L.; Fernie, A.R. Corrigendum: Gas chromatography massspectrometry-based metabolite profiling in plants. Nat. Protoc. 2015, 10, 1457. [Google Scholar] [CrossRef]
  42. Farag, M.A.; Huhman, D.V.; Lei, Z.; Sumner, L.W. Metabolic profiling and systematic identification of flavonoids and isoflavonoids in roots and cell suspension cultures of Medicagotruncatula using HPLC-UV-ESI-MS and GC-MS. Phytochemistry 2007, 68, 342–354. [Google Scholar] [CrossRef] [PubMed]
  43. Doerfler, H.; Lyon, D.; Nägele, T.; Sun, X.; Fragner, L.; Hadacek, F.; Egelhofer, V.; Weckwerth, W. Grangercausality in integrated GC–MS and LC–MS metabolomics data reveals the interface of primary and secondarymetabolism. Metabolomics 2013, 9, 564–574. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Wu, S.; Li, M.; Zhang, C.; Tan, Q.; Hu, C. Effects of phosphorus on fruit soluble sugar and citric acid accumulations in citrus. Plant Physiol. Biochem. 2021, 160, 73–81. [Google Scholar] [CrossRef] [PubMed]
  45. Xu, Y.; Shi, Q.H.; Gong, B. Characterization of COMT1-mediated low phosphorus resistance mechanism by metabolomics in tomato plants. Environ. Exp. Bot. 2020, 179, 104187. [Google Scholar] [CrossRef]
  46. Hernandez, G.; Ramirez, M.; Valdes-Lopez, O.; Tesfaye, M.; Graham, M.A.; Czechowski, T. Phosphorus stress in common bean: Root transcript and metabolic responses. Plant Physiol. 2007, 144, 752–767. [Google Scholar] [CrossRef] [Green Version]
  47. Dangi, S.P.; Aryal, K.; Magar, P.S.; Bhattarai, S.; Shrestha, D.; Gyawali, S. Study on effect of phosphorus on growth and flowering of marigold (tageteserecta). JOJ Wildl. Biodivers. 2019, 1, 108–112. [Google Scholar] [CrossRef]
  48. Sharma, A.K.; Chaudhary, S.V.S.; Gupta, Y.C. Effect of nitrogen andphosphorus on flowering and yield of African marigold (Tagetes erecta Linn.). Progress. Agric. 2010, 10, 158–160. [Google Scholar]
  49. Ahirwar, M.K.; Ahirwar, K.; Shukla, M. Effect of plant density, nitrogen and phosphorus levels on growth, yield and quality of AfricanMarigold. Ann. Plant Soil Res. 2012, 14, 153–155. [Google Scholar]
  50. Lynch, J.P. Root phenes for enhanced soil exploration and phosphorus acquisition: Tools for future crops. Plant Physiol. 2011, 156, 1041–1049. [Google Scholar] [CrossRef] [Green Version]
  51. Yang, P.; Chen, H.C.; Fan, H.Y.; Li, Q.S.; Gao, Q.; Wang, D.S.; Wang, L.L.; Zhou, C.; Zeng, E.Y. Phosphorus supply alters the root metabolism of Chinese flowering cabbage (Brassica campestris L. ssp. chinensis var. utilis Tsenet Lee) and the mobilization of Cd bound to lepidocrocite in soil. Environ. Exp. Bot. 2019, 167, 103827. [Google Scholar] [CrossRef]
  52. Wang, Y.; Krogstad, T.; Clarke, J.L.; Hallama, M.; Øgaard, A.F.; Eich-Greatorex, S.; Kandeler, E.; Clarke, N. Rhizosphereorganicanions play a minor role in improving crop species’ ability to take up residual phosphorus (P) in agricultural soils low in Pavailability. Front. Plant Sci. 2016, 7, 1664. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Tawaraya, K.; Horie, R.; Wagatsuma, T.; Saito, K.; Oikawa, A. Metabolite profiling of shoot extract, root extract, and root exudate of rice under nitrogen and phosphorus deficiency. Soil Sci. Plant Nutr. 2018, 64, 1–11. [Google Scholar] [CrossRef]
  54. Ruan, Y.L. Sucrose metabolism: Gateway to diverse carbon use and sugar signaling. Annu. Rev. Plant Biol. 2014, 65, 33–67. [Google Scholar] [CrossRef] [PubMed]
  55. Wang, Y.L.; Almvik, M.; Clarke, N.; Eich-Greatorex, S.; Øgaard, A.F.; Krogstad, T.; Lambers, H.; Clarke, J.L. Contrasting responses of root morphology and root-exuded organic acids to low phosphorus availability in three important food crops with divergent root traits. AoB Plants 2015, 7, plv097. [Google Scholar] [CrossRef] [PubMed]
  56. Wang, Y.; Lysøe, E.; Armarego-Marriott, T.; Erban, A.; Paruch, L.; van Eerde, A.; Bock, R.; Liu-Clarke, J. Transcriptome and metabolome analyses provide insights into root and root-released organic anion responses to phosphorus deficiency in oat. J. Exp. Bot. 2018, 69, 3759–3771. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Huang, C.Y.; Roessner, U.; Eickmeier, I.; Genc, Y.; Callahan, D.L.; Shirley, N.; Langridge, P.; Bacic, A. Metabolite profiling reveals distinct changes in carbon and nitrogenmetabolism in phosphate-deficient barley plants (Hordeum vulgare L.). Plant Cell Physiol. 2008, 49, 691–703. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Ganie, A.H.; Ahmad, A.; Pandey, R.; Aref, I.M.; Yousuf, P.Y.; Ahmad, S.; Iqbal, M. Metabolite profiling of low-P tolerant and low-P sensitive maize genotypes under phosphorus starvation and restoration conditions. PLoS ONE 2015, 10, e0129520. [Google Scholar] [CrossRef]
  59. Müller, J.; Gödde, V.; Niehaus, K.; Zörb, C. Metabolic adaptations of white lupinroots and shoots under phosphorus deficiency. Front. Plant Sci. 2015, 6, 1014. [Google Scholar] [CrossRef] [Green Version]
  60. Antoine, S.; Pailly, O.; Gibon, Y.; Luro, F.; Santini, J.; Giannettini, J.; Berti, L. Short-and long-term effects of carbohydrate limitation on sugar and organic acid accumulation during mandarin fruit growth. J. Sci. Food Agric. 2016, 96, 3906–3914. [Google Scholar] [CrossRef]
  61. Zhou, Y.; He, W.; Zheng, W.; Tan, Q.; Xie, Z.; Zheng, C.; Hu, C. Fruit sugar andorganic acid were significantly related to fruit Mg of six citrus cultivars. Food Chem. 2018, 259, 278–285. [Google Scholar] [CrossRef]
  62. Zheng, C.; Lan, X.; Tan, Q.; Zhang, Y.; Gui, H.; Hu, C. Soil application of calciumand magnesium fertilizer influences the fruit pulp mastication characteristics ofNanfeng tangerine (Citrus reticulata Blanco cv. Kinokuni). Sci. Hortic. 2015, 191, 121–126. [Google Scholar] [CrossRef]
  63. Zhang, Y.; Hu, C.; Tan, Q.; Nie, Z.; Zheng, C.; Gui, H.; Sun, X.; Zhao, X. Soil application of Boron and Zinc influence fruit yield and quality of Satsuma Mandarin in acidic soils. Agron. J. 2015, 107, 1–8. [Google Scholar] [CrossRef]
  64. Li, X.; Li, C.; Sun, J.; Jackson, A. Dynamic changes of enzymes involved in sugar and organic acid level modification during blueberry fruit maturation. Food Chem. 2020, 309, 125617. [Google Scholar] [CrossRef] [PubMed]
  65. Lin, Z.H.; Chen, L.S.; Chen, R.B.; Zhang, F.Z.; Jiang, H.X.; Tang, N.; Smith, B.R. Root release and metabolism of organic acids in tea plants in response to phosphorus supply. J. Plant Physiol. 2011, 168, 644–652. [Google Scholar] [CrossRef] [PubMed]
  66. Wen, Z.; Li, H.; Shen, J.; Rengel, Z. Maize responds to low shoot P concentration by altering rootmorphology rather than increasing root exudation. Plant Soil. 2017, 416, 377–389. [Google Scholar] [CrossRef]
  67. Lyu, Y.; Tang, H.; Li, H.; Zhang, F.; Rengel, Z.; Whalley, W.R.; Shen, J. Major crop species show differential balance between root morphological and physiological responses to variable phosphorus supply. Front. Plant Sci. 2016, 7, 1939. [Google Scholar] [CrossRef] [PubMed]
  68. Mohidin, H.; Hanafi, M.M.; Rafii, Y.M.; Abdullah, S.; Idris, A.S.; Man, S. Determination of optimum levels of nitrogen, phosphorus and potassium of oil palm seedlings in solution culture. Bragantia 2015, 74, 247–254. [Google Scholar] [CrossRef] [Green Version]
  69. Chinese Pharmacopoeia Commission. Pharmacopoeia of People’s Republic of China; China Medical Sciences and Technology Press: Beijing, China, 2020; p. 293. (In Chinese) [Google Scholar]
  70. Xiao, J.; Gu, C.; He, S.; Zhu, D.; Huang, Y.; Zhou, Q. Widely targeted metabolomics analysis reveals new biomarkers and mechanistic insights on chestnut (Castaneamollissima Bl.) calcification process. Food Res. Int. 2021, 141, 110128. [Google Scholar] [CrossRef]
  71. Tilsner, J.; Kassner, N.; Struck, C.; Lohaus, G. Amino acid contents and transport in oilseed rape (Brassica napus L.) under different nitrogen conditions. Planta 2005, 221, 328–338. [Google Scholar] [CrossRef]
  72. Sun, J.; Li, W.; Zhang, Y.; Guo, Y.; Duan, Z.; Tang, Z.; Abozeid, A. Metabolomics analysis reveals potential mechanisms in Bupleurum L. (Apiaceae) induced by three levels of Nitrogen fertilization. Agronomy 2021, 11, 2291. [Google Scholar] [CrossRef]
Plants 11 00752 g001 550
Figure 1. The whole Bupleurum was divided into four organs, which are the root, main shoot, lateral shoot and flower.
Figure 1. The whole Bupleurum was divided into four organs, which are the root, main shoot, lateral shoot and flower.
Plants 11 00752 g001
Plants 11 00752 g002 550
Figure 2. Growth performances of three concentrations of phosphorus fertilization. (A) Various organs fresh weight (n = 3); (B) Various organs length (n = 3). Three phosphorus fertilization levels: low phosphorus (LP), control group (CP) and high phosphorus (HP); and Four organs: root (R), main shoot (MS), lateral shoot (LS) and flower (F). FW: fresh weight; PL: plant organ length. ** represent p < 0.01.
Figure 2. Growth performances of three concentrations of phosphorus fertilization. (A) Various organs fresh weight (n = 3); (B) Various organs length (n = 3). Three phosphorus fertilization levels: low phosphorus (LP), control group (CP) and high phosphorus (HP); and Four organs: root (R), main shoot (MS), lateral shoot (LS) and flower (F). FW: fresh weight; PL: plant organ length. ** represent p < 0.01.
Plants 11 00752 g002
Plants 11 00752 g003 550
Figure 3. (A) Saikosaponin A analysis at three phosphorus fertilization levels and four organs. (B) Saikosaponin D analysis at three phosphorus fertilization levels and four organs. (C) The percentage of saikosaponins content, including saikosaponin A and D. Three phosphorus fertilization levels: low phosphorus (LP), control phosphorus (CP) and high phosphorus (HP); and Four organs: flower (F), main shoot (MS), lateral shoot (LS) and root (R). ** represent p < 0.01; *** represent p < 0.001.
Figure 3. (A) Saikosaponin A analysis at three phosphorus fertilization levels and four organs. (B) Saikosaponin D analysis at three phosphorus fertilization levels and four organs. (C) The percentage of saikosaponins content, including saikosaponin A and D. Three phosphorus fertilization levels: low phosphorus (LP), control phosphorus (CP) and high phosphorus (HP); and Four organs: flower (F), main shoot (MS), lateral shoot (LS) and root (R). ** represent p < 0.01; *** represent p < 0.001.
Plants 11 00752 g003
Plants 11 00752 g004 550
Figure 4. (A) Classification of 73 metabolites of Bupleurum. (B) The cluster analysis and heat map for Bupleurum plant parts under different levels of phosphorus fertilization. CP: control phosphorus; LP: low phosphorus; HP: high phosphorus; R: root; MS: main shoots; LS: lateral shoots; F: flower. The down-regulated and up-regulated and metabolites are represented as blue and yellow, respectively. The blue color means that abundance is 0, as shown in the lower right corner.
Figure 4. (A) Classification of 73 metabolites of Bupleurum. (B) The cluster analysis and heat map for Bupleurum plant parts under different levels of phosphorus fertilization. CP: control phosphorus; LP: low phosphorus; HP: high phosphorus; R: root; MS: main shoots; LS: lateral shoots; F: flower. The down-regulated and up-regulated and metabolites are represented as blue and yellow, respectively. The blue color means that abundance is 0, as shown in the lower right corner.
Plants 11 00752 g004
Plants 11 00752 g005 550
Figure 5. Score plots of PLS−DA for metabolites variation in flowers (A), lateral shoots (B), main shoots (C) and roots (D). Different markers and colors represent different meanings; ●, ◆, ▲, ★ represent roots (R), flowers (F), lateral shoots (LS), main shoots (MS), respectively; and red, green and blue represent high phosphorus (HP), low phosphorus (LP) and control phosphorus (CP), respectively.
Figure 5. Score plots of PLS−DA for metabolites variation in flowers (A), lateral shoots (B), main shoots (C) and roots (D). Different markers and colors represent different meanings; ●, ◆, ▲, ★ represent roots (R), flowers (F), lateral shoots (LS), main shoots (MS), respectively; and red, green and blue represent high phosphorus (HP), low phosphorus (LP) and control phosphorus (CP), respectively.
Plants 11 00752 g005
Plants 11 00752 g006 550
Figure 6. (A–F) The volcanic map of differential metabolites and enrichment analysis of KEGG. The red and green colors represent the content of metabolic that was up-regulated and down-regulated, respectively; ●, ▇, ▲, ▼ represent flowers, lateral shoots, main shoots and roots, respectively. The p-value indicates the degree of enrichment, and the closer the p-value is to 0, the more remarkable the enrichment. The point stand for the quantity of differential metabolites. (G) The Venn diagram of differential metabolites.
Figure 6. (A–F) The volcanic map of differential metabolites and enrichment analysis of KEGG. The red and green colors represent the content of metabolic that was up-regulated and down-regulated, respectively; ●, ▇, ▲, ▼ represent flowers, lateral shoots, main shoots and roots, respectively. The p-value indicates the degree of enrichment, and the closer the p-value is to 0, the more remarkable the enrichment. The point stand for the quantity of differential metabolites. (G) The Venn diagram of differential metabolites.
Plants 11 00752 g006
Plants 11 00752 g007 550
Figure 7. Visualization of primary metabolites in biochemical pathways, and the cluster analysis and heat map for metabolite classification in phosphorus stress.
Figure 7. Visualization of primary metabolites in biochemical pathways, and the cluster analysis and heat map for metabolite classification in phosphorus stress.
Plants 11 00752 g007
Plants 11 00752 g008 550
Figure 8. Molecular formula and molecular weight of saikosaponin A and D.
Figure 8. Molecular formula and molecular weight of saikosaponin A and D.
Plants 11 00752 g008
Table
Table 1. Amount of phosphorus fertilization applied to Bupleurum.
Table 1. Amount of phosphorus fertilization applied to Bupleurum.
TreatmentCPLPHP
Phosphorus fertilization0 kg∙ha−110 kg∙ha−120 kg∙ha−1
Table
Table 2. Elution gradient.
Table 2. Elution gradient.
Times (min)Acetonitrile (%)Pure Water (%)
0–5025–90%75–10%
50–5590%10%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sun, J.; Duan, Z.; Zhang, Y.; Cao, S.; Tang, Z.; Abozeid, A. Metabolite Profiles Provide Insights into Underlying Mechanism in Bupleurum (Apiaceae) in Response to Three Levels of Phosphorus Fertilization. Plants 2022, 11, 752. https://doi.org/10.3390/plants11060752

AMA Style

Sun J, Duan Z, Zhang Y, Cao S, Tang Z, Abozeid A. Metabolite Profiles Provide Insights into Underlying Mechanism in Bupleurum (Apiaceae) in Response to Three Levels of Phosphorus Fertilization. Plants. 2022; 11(6):752. https://doi.org/10.3390/plants11060752

Chicago/Turabian Style

Sun, Jialin, Zejia Duan, Ye Zhang, Sisi Cao, Zhonghua Tang, and Ann Abozeid. 2022. "Metabolite Profiles Provide Insights into Underlying Mechanism in Bupleurum (Apiaceae) in Response to Three Levels of Phosphorus Fertilization" Plants 11, no. 6: 752. https://doi.org/10.3390/plants11060752

APA Style

Sun, J., Duan, Z., Zhang, Y., Cao, S., Tang, Z., & Abozeid, A. (2022). Metabolite Profiles Provide Insights into Underlying Mechanism in Bupleurum (Apiaceae) in Response to Three Levels of Phosphorus Fertilization. Plants, 11(6), 752. https://doi.org/10.3390/plants11060752

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop