Comprehensive comparison on the anti-inflammation and GC-MS-based metabolomics discrimination between Bupleuri chinense DC. and B. scorzonerifolium Willd

Bupleuri Radix (BR) is a traditional Chinese medicine and widely used for cold and fever, influenza, inflammation, hepatitis and menstrual diseases. Two authentic medicinal plants of Bupleuri chinense DC. (Beichaihu, BCH) and B. scorzonerifolium Willd. (Nanchiahu, NCH) are recommended by the current Chinese Pharmacopoeia for BR. In the present study, the comparative investigations on the anti-inflammatory effects and gas chromatography-mass spectrometry (GC-MS)-based metabolomics for the species discrimination of BCH and NCH were conducted and reported. The in vitro evaluations indicated that the supercritical fluid extracts (SFEs) (IC50 of 6.39 ± 0.52 and 1.32 ± 0.05 mg (herb)/mL for BCH and NCH) were determined to be more potent than those of the hydro-distillation extracts (HDEs) (IC50 of 203.90 ± 8.08 and 32.32 ± 2.27 mg (herb)/mL for BCH and NCH) against LPS-induced inflammation in RAW264.7 macrophages. The higher anti-inflammatory effects of NCH were associated to its different chemical compositions to the BCH as characterized by the GC-MS analysis. Furthermore, based on the metabolomics and deep chemometric approaches, a minimum combination containing 15 chemical markers was optimized from the identified components and successfully applied for the species discrimination of BCH and NCH. This study not only helps to comparative understand BCH and NCH both in phytochemistry and pharmacology, but also provides the potential chemical markers for improvement of methods for the quality control of BCH and NCH.


Introduction
Bupleuri Radix (BR, also named as Chaihu) is a traditional Chinese medicine and has been widely used in China for over 2000 years for its broad pharmacological activities and human health benefits (Yuan et al., 2017). Until now, the known plants of Bupleurum genus in China include 36 species, 17 varieties and seven variants, and 20 of which were/are practically used as BR in clinical applications Sun et al., 2019). Although most of these plants vary slightly in morphologic appearance, their appreciable difference in chemicals and bioactivities attracts increasing attention. For example, the plant of Bupleuri longiradiatum Turcz has been forbidden to be used as BR by the Chinese pharmacopeia in clinics due to its significant toxicities . Currently, only two species of B. chinense DC. (also called Beichaihu, BCH) ( Figure 1A) and B. scorzonerifolium Willd. (also called Nanchaihu, NCH). ( Figure 1B) are recommended by the Chinese pharmacopeia (2020 version) to serve as the plant sources for BR.
In clinical applications, BR has been widely used for treatment of fever, pain, and inflammation associated with influenza or the common cold . The modern pharmacological investigations reported that BR presented broad biological activities including anti-influenza (Yan et al., 2022), anti-tumor Zhou et al., 2021), anti-inflammation (Jiang et al., 2020), anti-depression (Guo et al., 2020;Chen et al., 2021), anti-stress (Wang P. et al., 2021;Wang Z. et al., 2021)and antipyretic , antiviral (Du et al., 2018), antimicrobial (Liu et al., 2016), immunomodulatory (Tang et al., 2021) and hepatoprotective effects (Ren et al., 2019). Although different types of chemicals have been reported to contain in BCH (saponins, volatile oils, flavonoids and polysaccharides)  and NCH (saponins, lignans, fatty acids, volatile oils and polysaccharides) (Jiang et al., 2020), the saponins and volatile oils have been investigated to be the most important parts for both BR species which are associated to their antipyretic Li et al., 2018) and anti-inflammatory effects (Kim et al., 2015;Ma et al., 2016). The pharmaceutical products (BR injection, major BR decoction and minor BR decoction) containing such ingredients have been well developed and widely used in clinics (Ni et al., 2009;Chen et al., 2010).
Phytochemically, BCH was reported to mainly contain saponins, volatile oils, flavonoids and polysaccharides ; while NCH was reported to be rich in saponins, lignans, fatty acids, volatile oils and polysaccharides (Jiang et al., 2020). The chemical difference between them could also be mentioned as their different description as aroma smelling for BCH and transmutative-oil-like smelling for NCH under the item of 'Properties' by the Chinese Pharmacopeia. All these evidences suggested the significant difference between BCH and NCH in chemistry despite both of them are used as BR in clinical applications. Nevertheless, up to date, the comparative investigation for BCH and NCH is still not reported.
In this study, we aim to provide a comparative comprehension on BCH and NCH in order to improve their quality control for product development and clinical applications. The anti-inflammatory effects were compared to find out the active extracts from these two herbs. Furthermore, a chemical combination with minimum numbers of components was optimized for the species discrimination of BCH and NCH using the GC-MS-based metabolomics and deep chemometric analyses. This study not only helps to comparatively understand BCH and NCH both in phytochemistry and pharmacology, but also provides the potential chemical markers for the further quality control of BR herbs.

Herbal materials
Total 18 batches of the BR samples including 10 batches of B. chinense DC (BCH) and eight batches of B. scorzonerifolium Willd (NCH) were collected from different regions in China (Table 1). All samples were authenticated by the corresponding authors and the voucher specimens were deposited at the Institute of Chinese Medical Sciences, University of Macau, Macao, China. Prior to experiments, the herbal materials were dried, powdered and passed through 65 mesh sieve. Herbal extraction Hydro-distillation Approximately 800 g of mixed powders of BCH (10 batches, 80 g for each) or NCH (8 batches, 100 g for each) were soaked with 8 L of distilled water for 3 h, and then extracted with the hydro-distillated method for another 12 h. The collected volatile oils (HDEs) of BCH (yield: 0.1%) and NCH (yield: 0.25%) were stored at -20°C in amber glass bottles for further experiments.

Supercritical fluid extraction
The SFE extraction of BCH and NCH was carried out on a supercritical fluid extractor coupled with fully automated pilot scale systems (Supercritical Fluid Technology, INC., United States, Model: SFT250). Briefly, 100 g powder of BCH or NCH was loaded into the cylindrical stainless-steel extraction tank and then extracted with CO 2 (purity 99.5%) fluid under the following parameters: extraction kettle pressure of 260 Pa; vessel temperature of 60°C; vessel oven temperature of 70°C; and the extraction time of 2 h. Finally, the SFE extracts (SFEs) were collected in a 50 ml centrifuge tube (average yields of 1.16 ± 0.06% and 3.57 ± 0.44% for BCH and NCH, respectively), and stored at −20°C for further experiments.

Comparison of anti-inflammatory activities between BCH and NCH
Cell culture and cytotoxicity RAW264.7 cells were obtained from the American Type Culture Collection (ATCC; Manassas, VA, United States). The cells were cultured in Dulbecco‫׳‬s modified eagle medium (DMEM) supplemented with 10% heat-inactivated fetal bovine serum (FBS) and 1% penicillin/streptomycin (P/S) (ThermoFisher Scientific Inc., United States) at 37°C in humidified 5% CO 2 atmosphere (Linghu et al., 2020a). When the cells reached 80% confluence, the cells were sub-cultured after scraping from a 25 cm 2 flask.
The cytotoxicity of all extracts on RAW264.7 cells was assessed using the MTT assay (Linghu et al., 2020a). Briefly, RAW264.7 cells (1×10 4 cells/well) were seeded in 96-well plates and allowed to adhere overnight. Prior to cell treatment, both HDEs and SFEs were completely dissolved in DMSO and diluted with culture medium for further cell treatment. The cells were treated with the indicated concentrations of HDEs or SFEs for 24 h, and then incubated for an additional 1 h with fresh culture medium containing 0.5 mg/ml MTT. Subsequently, the culture medium was removed and the absorbance of the dissolved precipitate (in 150 μL of DMSO) was measured at 490 nm using a microplate reader (FlexStation3; Molecular Devices, United States).

Anti-inflammatory effects
The anti-inflammatory effects of BCH and NCH extracts (HDEs and SFEs) were comparatively investigated by evaluating on the inhibition of nitric oxide (NO) production in LPS-induced RAW264.7 macrophages. Briefly, RAW264.7 cells (1×10 4 cells/ well) were seeded onto 96-well plates and allowed to adhere overnight. Subsequently, the cells were pre-treated with indicated concentrations of individual extract for 1 h, and then stimulated with LPS (1 μg/ml) for 24 h. NO production was determined by measuring the nitrite accumulated in the medium with Griess reagent.

Discrimination of BCH and NCH with GC-MS-based metabolomics approach GC-MS analysis
The chemical profiles of the SFE extracts from BCH and NCH were characterized using a Thermo Trace 1,300 gaschromatography (GC) instrument equipped with a thermo ISQ LT single quadrupole mass spectrometer, a TRIPLUS RSH autosampler for liquid, static headspace and solid phase microextraction injections, a split/splitless injector and a Xcalibur chromatography processing system (Thermo Fisher, United States). Sample separation was preformed using an DB-5MS (30 m × 0.25 mm, 0.25 μm film thickness) capillary column (Agilent, United States). The electron-impact (EI) mode was selected and the ionization voltage was 70 eV. Both the interface and ionization source temperatures were 280°C. Mass spectral scan range was set at 50-550 m/z with a scan rate of 0.60 s per scan. The oven temperature programming was set as follows: 1 min at 140°C, 5°C/min to 180°C, holding for 2 min; 3°C/min to 280°C, holding for 2 min; 30°C/min to 300°C, holding for 1 min. The samples were injected in 30:1 split mode, and the injection temperature and volume were 280°C and 2 μL, respectively.

Data processing and markers identification
The raw data of SFE extracts for BCH and NCH from GC-MS analysis was extracted and introduced to Progenesis QI software (version 2.0, Waters Corporation, MA, United States) for data processing, including peak detection, alignment and normalization. The compounds were identified by comparing the mass spectra of the detected components in GC-MS with those records in National Institute of Standards and Technology (NIST) mass spectral library and linear retention indices (LRI) calculated relative to (C8-C20) n-alkanes with LRI database using the NIST Mass Spectral Library. To further evaluate comprehensively quality and explore differential markers between BCH and NCH, the processed data including peak intensity, m/z and retention time of identified compounds was Frontiers in Pharmacology frontiersin.org analyzed by a series of deep chemometric analysis including supervised regularized canonical correlation analysis (SRCCA), hierarchical clustering analysis (HCA) and partial least squares discriminant analysis (PLS-DA) using the Metaboanalyst website (https://www.metaboanalyst.ca).

Statistical analysis
Data were analyzed using GraphPad Prism 6.0 software. All data are presented as mean ± S.D., and each experiment was performed at least three times. Significant differences between groups were determined using one-way ANOVA with Dunnet's multiple comparisons test or unpaired t-test. p < 0.05 was considered difference significantly.

Anti-inflammatory activities between HDEs and SFEs of BCH and NCH
Excessive NO release is one of important indexes for activated macrophages to inflammatory response. The capability on inhibition of NO production could be used to assess the anti-inflammatory effects of drug candidates (Linghu et al., 2020a;Linghu et al., 2020b).
The cytotoxicity of HDEs and SFEs for BCH and NCH were illustrated in Figure 2 A and B. Under the nontoxic concentrations, both HDEs and SFEs of BCH and NCH presented the dose-dependent inhibition on NO production in LPS-induced RAW264.7 macrophages ( Figures 2C,D). The

FIGURE 2
The anti-inflammation activity of HDE and SFE from BCH and NCH on RAW264.7 cells. (A) Cytotoxicity of HDE from BCH and NCH on RAW264.7 cells; (B) cytotoxicity of SFE from BCH and NCH on RAW264.7 cells; (C) anti-inflammatory effects of HDE from BCH and NCH through the measurement of accumulated nitrite; (D) anti-inflammatory effects of SFE from BCH and NCH through the measurement of accumulated nitrite; (E) IC 50 of HDE from BCH and NCH on NO production in RAW264.7 cells; (F) IC 50 of SFE from BCH and NCH on NO production in RAW264.7 cells. *p < 0.05, **p < 0.01, ***p < 0.001 versus LPS group. ### p < 0.001 versus control group. Data are presented as mean ± S.D. (n = 3).
Frontiers in Pharmacology frontiersin.org Frontiers in Pharmacology frontiersin.org 06  Figure 2E), as well as those of 6.39 ± 0.52 and 1.32 ± 0.05 mg (herb)/mL for the SFEs of BCH and NCH, respectively ( Figure 2F). The results indicated the higher anti-inflammatory effects of SFEs than that of HDEs for BCH and NCH, as well as the better potency of NCH than that of BCH.

Chemical characterization of SFEs from BCH and NCH
Based on the above in vitro comparison on LPS-induced RAW264.7 macrophages, the SFEs of both BCH and NCH were determined to be much potent than the HDEs on antiinflammation. Therefore, the chemical compositions of the SFEs were further characterized by GC-MS analysis.
With the established GC-MS method, all SFE samples of BCH (10 batches) and NCH (8 batches) were chromatographically analyzed. The typical total ion chromatograms (TIC) of the representative samples for BCH (BCH001) and NCH (NCH001) were illustrated in Figure 3A. Furthermore, by comparing the mass spectra of the detected components in GC-MS with those records in National Institute of Standards and Technology (NIST) mass spectral library, as well as their retention indices with published literatures (Tykheev et al., 2018;Li et al., 2005), 28 components in BCHs (60.99-70.56% of the total peak area) and 44 components in NCHs (74.04-93.28% of the total peak area) were identified and summarized in Table 2. Among them, 24 components were found both in BCH and NCH, while 4 and 20 specie-specific components were found in BCH and NCH, respectively. The most abundant components in BCH were palmitic acid (25.4% on average), linoleic acid (15.8% on average), oleic Acid (6.1% on average) and octadecanoic acid (4.4% on average). And the most abundant components in NCH were linoleic acid (12.4% on average) and palmitic acid (8.8% on average), as well as other two specific components of mandenol (7.5% on average) and dodecanoic acid (4.7% on average). These results indicated the significant differences in chemical compositions between BCH and NCH both in quality and quantity.
According to classification based on the types of chemical structure ( Figure 3B), majorities of SFEs components in BCH and NCH belong to fatty acids (averaged relative contents of 51.07 and 50.36%), sterols (averaged relative contents of 4.53 and 1.06%) and, miscellaneous components (averaged relative contents of 2.86 and 9.98%). Nevertheless, the compositions and relative contents of the major components in SFEs of BCH and NCH varied greatly between these two species. As illustrated in Figure 3C, lactones, aldehydes, acetates and fatty alcohols were much more abundant in NCH, while BCH contained more sterols by contrast. These results further indicated the significant differences in chemical types of the major components between BCH and NCH.

Species discrimination of BCH and NCH with metabolomic-based chemometrics
Chemometric analysis including SRCCA, PLS-DA and HCA were employed to visualize the discrepancy and screen the potential makers for species discrimination of BCH and NCH. The GC-MS chromatograms of all BCH (BCH001-010) and NCH (NCH001-008) samples were illustrated in Figure 4A. Prior to statistical analysis, features with low repeatability in all samples were removed by setting the value of relative standard deviation (RSD) at 25%. Moreover, data of all samples were normalized by sum then were transformed to a binary logarithm (log 2 X, where X represents the peak area) to avoid the instrument error.
As illustrated in Figure 4B, the results of SRCCA indicated a significant low inter-species correlation while strong intraspecies correlation between these two species. The cross validation [R 2 Y (cum) = 0.93, Q 2 (cum) = 0.90] also suggested the good predictive capability and the significant explanatory power of PLS-DA model for effective species discrimination of BCH and NCH ( Figure 4C). In addition, the grouping trend was confirmed by HCA, which could be seen from the well clustered groups of BCH and NCH ( Figure 4D). All these results further suggested the significant difference in chemical compositions between BCH and NCH.
To further verify the accuracy and feasibility of the optimized combination of 15 chemical components for specie discrimination of BCH and NCH, the relative content of above 15 markers were normalized to build new Chemometric models. As illustrated in Figure 5, individual BCH and NCH species could be well discriminated under multiple chemometric models (including PLS-DA, HCA) with above 15 markers, and the new models presented a comparably satisfied performance for discrimination of two BR species (R 2 Y (cum) = 0.940, Q 2 (cum) = 0.891).
Frontiers in Pharmacology frontiersin.org 09 scorzonerifolium Willd. (NCH) are officially authenticated to be used for decoction pieces of BR, the difference in chemical compositions between them could be seen from the previous publications of chemical analysis for individual species (Yang et al., 2010;Tykheev et al., 2018). However, the comparative investigation on the similarity and difference between BCH and NCH is still not reported. Therefore, in this study, the chemical compositions of the anti-inflammatory extracts from BCH and NCH were firstly comprehensively compared for further selection of the potential chemical markers for species discrimination.
As reported by previous publications, the volatile components are one of the most important parts of chemicals in BR (Meng et al., 2014;Li et al., 2015) and associated to its pharmacological activities for antipyretic and anti-inflammation (Ni et al., 2009;Chen et al., 2010). Considering the low contents of essential oils (HDEs) in the BR (yield of 0.1 and 0.25% for BCH and NCH by hydro-distillation extraction in this study), the extracts of BCH (average yield: 1.16 ± 0.06%) and NCH (average yield: 3.57 ± 0.44%) prepared by supercritical fluid extraction (SFEs) were also involved for comparison. Beside, as summarized in Table 2, compounds including palmitic acid, linoleic acid, octadecanoic acid, 1-dodecanol, lauryl acetate, mandenol and dodecenyl succinic anhydride were common components of NCH independence on habitats; while falcarinol, palmitic acid, linoleic acid, oleic acid, octadecanoic acid, eicosanoic acid and docosanoic acid were consistent compounds of BCH independence on habitats. However, their content varied from traces to significant amounts in individual specie, and some components were not detected in samples from other regions, suggesting certain difference in intra-species for BCH and NCH from different regions. As shown in Figures 4B,D, these intra-species differences allow samples divided into groups almost according to the geographical regions, indicating habitats of individual specie also contributed to the chemical differences.
Apart from the higher potent of the NCH than the BCH, the SFEs (IC 50 of 6.39 ± 0.52 and 1.32 ± 0.05 mg (herb)/mL for BCH and NCH) were observed to present much better activities against LPS-stimulated inflammation (NO release) in RAW264.7 macrophages than those of the HDEs (IC 50 of 203.90 ± 8.08 and 32.32 ± 2.27 mg (herb)/mL for BCH and NCH). The obtained results suggested the contribution of some other hydrophobic components in BR to anti-inflammation.
The compositional analysis of the SFEs of BCH and NCH were preformed using the GC-MS system. BCH exhibited a significant difference to NCH in both the chemical constituents and the contents of the main components containing in the SFEs. A total of 48 components of different chemical types have been identified in the SFEs of BCH and NCH, including 24 common components as well as 4 and 20 specific components for BCH and NCH, respectively. Among them, lactones, aldehydes, acetates and fatty alcohols were significantly abundant in NCH, while sterols were more enriched in BCH. The characteristic distribution of the volatile (aromatic) components reveals the different smell for each decoction pieces as specifically described in the Chinese pharmacopeia for BCH (aroma) and NCH (transmutative-oil-like).
Chromatography-based untargeted metabolomics analyses are rapid and reliable approaches which have been quickly developed and applied for the quality research of Chinese herbal medicines in recent years (Tao et al., 2018;Lu, et al., 2022;Xiao, et al., 2022). However, appropriate screening and Accuracy evaluation of potential marker compounds. (A) BCH and NCH were well discriminated at specie level by PLS-DA using 15 optimized potential makers; (B) BCH and NCH were well separated by HCA using 15 optimized potential makers.
Frontiers in Pharmacology frontiersin.org selection of the chemical markers with minimum combination for effective species discrimination is very important for such studies. In this study, by combining various chemometric models (HCA, PLS-DA, and SRCCA) as well as the VIP values and intensity changes evaluated by the t-test, 15 makers from the identified components in the SFEs of BCH and NCH were discovered. As shown in Figure 4F, the selected components could be divided into three groups including the non-specific components (C2, C3, C6, C7, C8 and C10), BCH-specific components (C1 and C11) and NCH-specific components (C4, C5, C9, C12, C13 and C14).
Among the non-specific components, the contents of C2, C6, C7 and C8 in NCH were significantly higher in NCH than those in BCH. These higher-content and specific components in NCH might be mainly responsible for its significantly potency on anti-inflammation. Furthermore, the good species discrimination with the optimized components combination was verified under PLS-DA and hierarchical cluster analyses. The results suggested the potential selection of the chemical markers from this combination for the quality control of BCH and NCH.

Conclusion
In summary, the similarity and difference between two BR plants including B. chinenses DC (BCH) and B. scorzonerifolium Willd (NCH) were firstly comprehensively compared and reported both in chemistry and pharmacology. The higher anti-inflammatory activities of the SFEs than the HDEs suggested the involvement of some hydrophobic components apart from the essential oils on the functions of BR against fever and inflammation. Moreover, the developed minimum chemical combination including 15 components from both BCH and NCH was successfully applied for the discrimination of individual species. The contribution of this study is not only help us to better understand the activity related chemical difference between BCH and NCH, but also provide the potential selection of chemical markers for further improvement of the quality control standard for BCH and NCH.

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 authors.

Author contributions
MZ and HY conceived and designed the study. MZ, LX, K-GL, GZ, QC, LS, PD, and YH conducted the experiments and data analysis. MC and JZ provided the technical support and advices for the study. MZ wrote the manuscript and HY revised the manuscript. All authors contributed to the review and the approval of the final manuscript.