Next Article in Journal
Identification and Bioaccessibility of Maillard Reaction Products and Phenolic Compounds in Buckwheat Biscuits Formulated from Flour Fermented by Rhizopus oligosporus 2710
Previous Article in Journal
MSF Enhances Human Antimicrobial Peptide β-Defensin (HBD2 and HBD3) Expression and Attenuates Inflammation via the NF-κB and p38 Signaling Pathways
Previous Article in Special Issue
Attomole-per Cell Atomic Mass Spectrometry Measurement of Platinum and Gold Drugs in Cultured Lung Cancer Cells
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 28
Issue 6
10.3390/molecules28062745
molecules-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

Untargeted Metabolomic Analysis and Chemometrics to Identify Potential Marker Compounds for the Chemical Differentiation of Panax ginseng, P. quinquefolius, P. notoginseng, P. japonicus, and P. japonicus var. major

by
Ruifeng Ji
1,2,†,
Thomas Avery Garran
1,†,
Yilu Luo
2,
Meng Cheng
1,
Mengyue Ren
2,* and
Xiuteng Zhou
1,*
1
State Key Laboratory of Dao-di Herbs, National Resource Center for Chinese Materia Medica, China Academy of Chinese Medical Sciences, Beijing 100700, China
2
School of Chinese Materia Medica, Guangdong Pharmaceutical University, Guangzhou 510006, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Molecules 2023, 28(6), 2745; https://doi.org/10.3390/molecules28062745
Submission received: 14 January 2023 / Revised: 14 March 2023 / Accepted: 15 March 2023 / Published: 18 March 2023
(This article belongs to the Special Issue State-of-the-Art Analytical Technologies for Metabolomics Analysis)
Browse Figures
Versions Notes

Abstract

:
The Panax L. genus is well-known for many positive physiological effects on humans, with major species including P. ginseng, P. quinquefolius, P. notoginseng, P. japonicus, and P. japonicus var. major, the first three of which are globally popular. The combination of UPLC-QTOF-MS and chemometrics were developed to profile “identification markers” enabling their differentiation. The establishment of reliable biomarkers that embody the intrinsic metabolites differentiating species within the same genus is a key in the modernization of traditional Chinese medicine. In this work, the metabolomic differences among these five species were shown, which is critical to ensure their appropriate use. Consequently, 49 compounds were characterized, including 38 identified robust biomarkers, which were mainly composed of saponins and contained small amounts of amino acids and fatty acids. VIP (projection variable importance) was used to identify these five kinds of ginseng. In conclusion, by illustrating the similarities and differences between the five species of ginseng with the use of an integrated strategy of combining UPLC-QTOF-MS and multivariate analysis, we provided a more efficient and more intelligent manner for explaining how the species differ and how their secondary metabolites affect this difference. The most important biomarkers that distinguished the five species included Notoginsenoside-R1, Majonoside R1, Vinaginsenoside R14, Ginsenoside-Rf, and Ginsenoside-Rd.

1. Introduction

Many plants from the ginseng genus (Panax) are regarded as valuable traditional medicine resources, as well as valuable herbs in the herb and supplement industry globally. They are well known for many positive physiological effects on humans and widely used in traditional medicine systems, as well as in the health food and supplement industry [1]. In China, the genus consists of six native and one introduced species (P. quinquefolius). Of these seven species, five are recognized as medicinal and officially recorded in the Pharmacopoeia of the People’s Republic of China (2020 Edition): P. ginseng C.A. Meyer (ginseng, Chinese ginseng, renshen, RS), P. quinquefolius L. (American ginseng, xiyangshen, XYS), P. notoginseng F. H. Chen (notoginseng, sanqi, SQ), P. japonicus C.A. Meyer (Japanese ginseng, zhujieshen, ZJS), P. japonicus var. major (Burkill) C.Y. Wu and K.M. Feng (zhuzishen, ZZS). In the global ginseng market, RS has the largest market share, followed by XYS, SQ, ZJS, and ZZS in the descending order of market share [2].
As medicinal materials from the homologous species, they basically contain ginsenosides, polysaccharides, volatile oils, proteins, amino acids, organic acids, flavonoids, vitamins, and trace elements and other active ingredients. At the same time, their traditional efficacy and pharmacological activity are relatively similar. However, there are also some differences. For example, both RS and XYS have the effect of replenishing “qi” (a concept of traditional Chinese medicine, which is the energy of movement that reaches the skin externally and the organs internally, maintaining vitality) and producing fluid. SQ, ZJS, and ZZS all have the effect of removing blood stasis, stopping bleeding, and relieving pain. RS has antioxidant effects, can improve immunity, is anti-fatigue and anti-tumor, and can regulate the nervous system and immune system. XYS can improve blood circulation, improve the nervous system, and regulate immune function. SQ has the function of protecting the cardio-cerebrovascular and nervous systems and being anti-tumor, anti-bacterial, and anti-inflammatory. ZJS can protect the liver and heart, reduce blood fat, resist fatigue, enhance immunity, and is anti-tumor. ZZS can be anti-inflammatory, is an analgesic, is anti-tumor, and can regulate the immune function of the cardiovascular and cerebrovascular systems. According to traditional Chinese medicine (TCM) theory, RS and XYS are herbs for replenishing qi, and RS is considered slightly warm, while XYS is considered cool [3]. SQ is sweet and slightly bitter, warm in nature, and has a special tropism to the liver and stomach. ZJS and ZZS have a sweet and bitter taste similar to that of SQ, and they show tropism activity in the liver and stomach similar to SQ; however, they have additional tropism activity in the lung and thus can be used to treat phlegm and stop coughing. In addition, ZZS is slightly cold in nature and thus has additional efficacy in nourishing the lung-yin [4]. By integrating the information about the properties (nature, taste, and organ-specific channel tropism), as well as indications of RS, XYS, SQ, ZJS, and ZZS into a Venn diagram (Figure 1), we created a visual means of identifying the commonalities and differences in these herbs according to their traditional properties [5].
The phytochemistry and biological activities of Panax has been widely investigated. Although different species of Panax have diverse properties and indications in the TCM system, most literature presume that saponins (also known as ginsenosides) are the major active ingredients [6]. The current literature available regarding the metabolic difference studies of the same parts among different Panax species [7,8,9], the different parts [10] or ages [11,12] of the same species, and their combination demonstrates that ginsenosides have great potential when differentiating various species within the ginseng genus. The content and composition of ginsenosides that belong to the protopanaxadiol (PPD)-, protopanaxatriol (PPT)-, ocotillol (OCT)-, and oleanane (OA)-types; the C-17 side-chain-varied; the malonylated; and the others saponins vary widely in different ginseng species and in different parts of the plants [13]. Discrimination of the differences in the chemical composition among the five most commonly used species of Panax could inspire new research into potentially novel clinical applications of these species [14,15]. In contrast to the three most important Panax species (P. ginseng, P. quinquefolius, and P. notoginseng) covering very extensive studies, P. japonicus and P. japonicus var. major have far less literature support, particularly for the differences in their metabolome [16,17]. Characterizing the chemical components can ensure the correct use [18].
A key segment of the modernization research for TCM is the clarification of the chemical compositions for each species used [19]. The development of powerful and feasible analytical methods capable of the comprehensive deconvolution of plant metabolites that are from similar species is currently an important topic in the field of analytical chemistry [20,21]. Ultra-high performance liquid chromatography (UPLC) combined with quadrupole time-of-flight tandem mass spectrometry (QTOF-MS) is the most popular approach for the systematic multi-component characterization of an herb or a bio-sample, and UPLC-QTOF-MS profiling combined with chemometrics can render the untargeted metabolomics analysis suitable for discovering potential chemical markers for the authentication and quality evaluation of easily confused herb [11,22,23]. Yoon et al. (2022) reported the characterization of a total of 62 saponins from the geographical origin discrimination of P. ginseng by UPLC-QTOF-MS with multivariate analysis [24]. The establishment of a more powerful untargeted metabolomics platform can offer useful insights into the evaluation of the properties of medicinal plants, including easily confused species with similar pharmacological activity.
Untargeted metabolomics can detect the dynamic changes of all small-molecule metabolites before and after stimulation or disturbance in cells, tissues, organs, or organisms without bias, screen differential metabolites through bioinformatics analysis, and conduct pathway analysis of differential metabolites to reveal the physiological mechanism of their changes. Targeted metabolomics is the study of a specific class of metabolites, used for the discovery and quantization of differential metabolites, and the in-depth research and analysis of subsequent metabolic molecular markers, which play an important role in disease research, animal model validation, biomarker discovery, disease diagnosis, drug development, plant metabolism, and other studies. Targeted metabolomics usually have several advantages over untargeted metabolomics, both in terms of sensitivity and specificity, as well as quantitative data processing.
In our case study, we presented an integral strategy of combining untargeted metabolomics with multivariate analysis and applied it to simultaneously differentiate five Panax species. PCA, PLS-DA, and OPLS-DA are all dimensionality reduction methods. They are used to look for a small number of principal components and to explain major changes in the data. We found that magnitudes were measured by the VIP of the chosen principal components. Our data showed that the newly established method was precise and rapid for distinguishing the five chosen species of ginseng.

2. Results and Discussion

2.1. Optimization of a UPLC-QTOF-MS Approach for the Enhanced Profiling and Characterization of Metabolites Simultaneously from RS, XYS, SQ, ZJS, and ZZS

The effects of the mobile phase system on the chromatographic peak and gradient elution on sample peak separation were investigated. The aqueous solution of formic acid containing 0 mol/L, 0.1 mol/L, 0.1 mol/L (containing 10 mmol/L ammonium formate), and 0.2 mol/L (containing 10 mmol/L ammonium formate) were investigated. The results showed that when gradient elution was applied with 0.2 mol/L formic acid solution (containing 10 mmol/L ammonium formate)–acetonitrile solution, both the peak shape and the separation effect were acceptable. Figure 2 shows the base peak chromatograms (BPC) at the optimized analysis conditions of the representative samples of five Panax species, which illustrate the metabolite differences among the five species. The similarity in the spectrograms of RS and XYS are obvious, while ZJS and ZZS have a similar composition. The essential characteristics are consistent with other literature [25,26]. The positive and negative ion data were interpreted so as to comprehensively and simultaneously characterize the five Panax species. After processing the chromatographic data, those listed as “Identified Compounds” were putatively identified using databases such as Pubchem and Massbank.

2.2. Multivariate Analysis for UPLC-QTOF-MS Results and Selection of Target Ion

To obtain more information on the components of the five species, UPLC-QTOF-MS data was used for untargeted component analysis. Multivariate statistical analysis methods, such as unsupervised principal component analysis (PCA), supervised principal component analysis (PLS-DA), and orthogonal partial least squares discrimination analysis (OPLS-DA), were performed to identify the differences in metabolic profiles among the species. The quality of the model is dependent on the values of R2Y and Q2; when the values are higher (>0.9 and close to 1), the model is more reliable [27]. R2Y and Q2 represent the explanatory rate model and the forecast rate, respectively. The MS data of samples were statistically analyzed by PCA, PLS-DA, and OPLS-DA. The score plots of each revealed the factors that accounted for the largest variations and grouping tendencies. We found that the clustering degree of the samples in the PCA model indicated that the instrument was stable during this experiment in the positive and negative ion modes (Figure 3A,B). We also found signs of partial separation among the five groups. Using the PLS-DA method to analyze the metabolite profile of the Panax samples, the samples from RS, XYS, SQ, ZJS, and ZZS were clearly clustered and segregated into different groups scattering in different quadrants up to the 95% Hotelling T2 ellipse. There were significant differences among the five groups (ZZS is a variant of ZJS that considers partial overlap of results acceptable): [R2X (cum) =0.592, R2Y (cum) = 0.992, Q2 (cum) = 0.985] in positive ion mode, and [R2X (cum) =0.671, R2Y (cum) = 0.992, Q2 (cum) = 0.986] in negative ion mode (Figure 3C,D). It can be seen that R2Y and Q2 are both above 0.5 and close to 1, indicating the good reliability, good predictability, and no over-fitting for the PLS-DA model. These findings indicate that the PLS-DA model can be used to distinguish the five species. To obtain the greatest separation of metabolites, OPLS-DA was performed (Figure 3E,F). In the positive and negative ion modes, there was a clear separation among the groups. (ZZS is a variant of ZJS that considers partial overlap of results acceptable.) The samples from different species tended to cluster in a concentrated manner, with a high degree of aggregation and without any obvious intragroup differences. We found that the result of the OPLS-DA model was consistent with the PLS-DA model, and the Q2 of PLS-DA model was higher than the OPLS-DA model ([R2X (cum) = 0.592, R2Y (cum) = 0.992, Q2 (cum) = 0.977] in positive ion mode, and [R2X (cum) = 0.671, R2Y (cum) = 0.992, Q2 (cum) = 0.981] in negative ion mode). These findings indicated that the PLS-DA model had the highest discrimination and prediction rates (p < 0.05).
Variable importance in projection (VIP) is a weighted sum of the squares of the PLS loadings. To distinguish the most important metabolites among the species, p-values and VIP scores were used to screen for differential components (Figure 4). Given VIP > 4.0 and p < 0.05, a total of 49 potential biomarkers were preliminarily screened. The heatmaps of the component of significant difference among the species, as detected in the positive and negative modes, are illustrated in Figure 5.
According to the screening results of potential markers, there were some significant differences in chemical composition between the five Panax species, including saponins, amino acids, and fatty acids. The primary differences were saponins (Table 1), such as tR 6.80 and m/z 817.4781, identified as ginsenoside Re5 or ginsenjilinol, and tR 8.45 and m/z 945.5336, identified as ginsenoside Rd or ginsenoside Re [28,29]. The levels of 12 ginsenosides showed significant differences in the five species. The spectral intensities of these 12 ginsenosides are presented as bar plots (Figure 6). These bar charts show the peak intensities for the target ion compounds, which vary significantly between species.

3. Materials and Methods

3.1. Reagents and Material

Acetonitrile, methanol (Merck, Darmstadt, Germany), leucine-enkephalin (LE, Sigma-Aldrich, St. Louis, MO, USA, No. L9133-50MG), formic acid (FA; CNW, Shanghai, China), and ammonium formate (AF; Sigma-Aldrich, MO, USA, No. 41470) were of HPLC grade. Deionized water was purchased from Guangzhou Watson’s Food & Beverage Co., Ltd. (Guangzhou, China). A total of five species of the officinal PSA herb genus Panax, including Panax ginseng, P. quinquefolius, P. japonicus, P. japonicus var. major, and P. notoginseng, were collected from wild or cultivated sources in China and the United States in 2018 (Figure 7).

3.2. Sample Preparation

The Panax ginseng, P. quinquefolius, P. japonicus, P. japonicus var. major, and P. notoginseng roots were powdered to a homogeneous size and sieved through a No. 60 mesh. The study involved 10 batches of RS, 8 XYS, 10 SQ, 8 ZJS, 10 ZZS. Ultrasonic extraction with 60% methanol at a solid-to-liquid ratio of 1:5 g·mL−1 for 1 h was used to extract the medicinal materials. Dried crude powder (0.3 g) was accurately weighed in a centrifuge tube, and then 1.5 mL of 70% methanol/water (v/v) solution was added. Then, the sample was processed using ultrasonic extraction for 45 min and cooled to room temperature. The sample solution was centrifuged at 12,000 rpm for 5 min at 20 °C, and the supernatant was stored at 4 °C and filtered through a 0.22 μm filter membrane before injection for UPLC analysis.

3.3. UPLC/QTOF-MS Analysis

The ultrahigh-performance liquid chromatography by an Acquity UPLC system (Waters, Milford, MA, USA) was coupled with high-resolution MS analysis by a Micromass QTOF mass spectrometer (Waters, Manchester, U.K.). An Acquity UPLC BEH C18 column (2.1 mm × 50 mm, 1.7 μm) was used to perform the metabolite profiling, and 10 μL of each sample was injected into a gradient system at a flow rate of 0.3 mL/min. The mobile phase consisted of 0.2 mol/L formic acid and 10 mmol/L ammonium formate in water (A) and acetonitrile (B). The starting eluent was 2% B, and its proportion was held constant for 1 min, increased linearly to 10% from 1.0 to 2.0 min, to 25% from 2.0 to 5.1 min, to 32% from 5.1 to 6.1 min, to 38% from 6.1 to 8.1 min, held constant at 38% for 0.5 min, increased to 50% from 8.6 to 9.1 min, to 55% from 9.1 to 9.6 min, held constant at 55% until 12.5 min, increased to 70% from 12.5 to 14.0 min, to 100% from 14.0 to 16.1 min, held constant at 100% until 18.5 min, returned to 2% B at 18.5 min, and held constant until 20 min to equilibrate the column. This UPLC elution condition was optimized to detect the maximal number of metabolites in P. ginseng, especially to separate ginsenosides for identifying markers. The column was maintained at 40 °C. The mass spectrometer was equipped with an ESI source and operated in positive and negative ion modes. The MS conditions were as follows: capillary and cone voltages were adjusted to 2800 V and 20 V in positive ion mode but 3000 V and 20 V in negative ion mode, separately. The source and desolvation temperatures were maintained at 100 °C and 350 °C, respectively. The desolvation gas used was N2. The flow rate of desolvation gas and cone gas were at 500L/h and 50L/h, respectively. The scanning m/z range was 50 to 1500. Mass resolution > 10,000 FWHM (standard mode). To ensure that mass was measured accurately, leucine-enkephalin (10 μg/mL) was infused as the reference lockmass compound for real-time correction, and the [M+H]+ ion at 556.2771 Da and the [M-H]− ion at 554.2615 Da were detected in the analysis.

3.4. Multivariate Analysis

To evaluate the potential characteristic components of the five species of Panax in our study, the raw data of all samples were analyzed with the MassLynx application manager version 4.1 (Waters MS Technologies) [30]. The method parameters were as follows: retention time range, 0.3–17 min; mass range, 100–1500 Da; mass tolerance, 50 mDa; and noise elimination level, 6.00. For further analysis, a combination of retention time (tR) and mass data (m/z) from the detected peaks was assigned as temporary ID (tR-m/z), and the identifier of each peak for data adjustment was based on their chromatographic elution order of UPLC [31]. The identities of the variables chosen as biomarkers were based on m/z. The list was then applied for PCA, PLS-DA, and OPLS-DA using the Ezinfo software (Waters). Markers that differentiated the five groups were selected according to the variable importance in the projection values (VIP). The data matrix involving tR, m/z, and normalized peak area was exported into the SIMCA-P 14.0 software (Umetrics, Umea, Sweden) for chemometric analysis [32]. The compound with a VIP>4 was evaluated as the potential characteristic component from the positive and negative chromatogram [33]. The ANOVA value was lower than 0.05. Internal standard standardization was used.

4. Conclusions

Our study used a metabolomic fingerprinting approach based on discrimination by combining UPLC-QTOF-MS and multivariate analysis with the aim of discovering robust biomarkers for the authentication of the five most commonly used Panax species. The target compounds were selected from UPLC-MS screening results using multivariate analysis. Samples of the five Panax species were analyzed to create PCA, PLS-DA, and OPLS-DA models, which were compared with the UPLC-MS multivariate models. The resulting score plots showed good separation between the five species. In addition, the intensities of the target markers were compared with one another using bar plots to illustrate differences between the species. This is the first report to systematically compare the metabolome differences among these five important Panax species. The method elucidates highly efficient discrimination, especially considering the similarly of the samples, as required for authenticity testing. Moreover, with further optimization, this method could be used to analyze other plant materials to facilitate “chemical fingerprinting” and could be used as a new and important tool for the quality control of natural products derived from the Panax genus.

Author Contributions

R.J. collected the samples, designed the experiments, and wrote the original draft; T.A.G. collected the samples and polished the language; Y.L. performed the experiments and analyzed the data; M.C. performed the formal analysis; M.R. and X.Z. edited and revised the manuscript. 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 Public Welfare Research Institute (grant number ZZ13-YQ-088).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the research are available from the authors.

Abbreviations List

VIPvariable importance in projection
RSP. ginseng
XYSP. quinquefolius
SQP. notoginseng
ZJSP. japonicus
ZZSP. japonicus var. major
TCMtraditional Chinese medicine
PPDprotopanaxadiol
PPTprotopanaxatriol
OCTocotillol
OAoleanane
PCAprincipal component analysis
PLS-DAsupervised principal component analysis
OPLS-DAorthogonal partial least squares discrimination analysis
HCAhierarchical clustering analysis

References

  1. Jee, H.S.; Chang, K.H.; Park, S.H.; Kim, K.T.; Paik, H.D. Morphological characterization, chemical components, and biofunctional activities of Panax ginseng, Panax quinquefolium, and Panax notoginseng roots: A comparative study. Food Rev. Int. 2014, 30, 91–111. [Google Scholar] [CrossRef]
  2. Baeg, I.H.; So, S.H. The world ginseng market and the ginseng (Korea). J. Ginseng Res. 2013, 37, 1–7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Li, M.; Hua, S.; Huang, X.; Yue, H.; Chen, C.; Liu, S. Non-targeted metabonomics to investigate the differences in the properties of ginseng and American ginseng based on rapid resolution liquid chromatography coupled with quadrupole-time-of-flight mass spectrometry. J. Sep. Sci. 2021, 44, 3497–3505. [Google Scholar] [CrossRef] [PubMed]
  4. Yang, B.R.; Yuen, S.C.; Fan, G.Y.; Cong, W.H.; Leung, S.W.; Lee, S.M. Identification of certain Panax species to be potential substitutes for Panax notoginseng in hemostatic treatments. Pharmacol. Res. 2018, 134, 1–15. [Google Scholar] [CrossRef]
  5. Fu, J.; Pang, J.; Zhao, X.; Han, J. The Quantitative Ideas and Methods in Assessment of Four Properties of Chinese Medicinal Herbs. Cell Biochem. Biophys. 2015, 71, 1307–1312. [Google Scholar] [CrossRef]
  6. Liu, L.; Xu, F.R.; Wang, Y.Z. Traditional uses, chemical diversity and biological activities of Panax L. (Araliaceae): A review. J. Ethnopharmacol. 2020, 263, 112792. [Google Scholar] [CrossRef]
  7. Kim, S.; Shin, B.K.; Lim, D.K.; Yang, T.J.; Lim, J.; Park, J.H.; Kwon, S.W. Expeditious discrimination of four species of the Panax genus using direct infusion-MS/MS combined with multivariate statistical analysis. J. Chromatogr. B Analyt. Technol. Biomed. Life Sci. 2015, 1002, 329–336. [Google Scholar] [CrossRef]
  8. Liu, H.; Lu, X.; Hu, Y.; Fan, X. Chemical constituents of Panax ginseng and Panax notoginseng explain why they differ in therapeutic efficacy. Pharmacol. Res. 2020, 161, 105263. [Google Scholar] [CrossRef]
  9. Yang, F.; Chen, B.; Jiang, M.; Wang, H.; Hu, Y.; Wang, H.; Xu, X.; Gao, X.; Yang, W. Integrating Enhanced Profiling and Chemometrics to Unveil the Potential Markers for Differentiating among the Leaves of Panax ginseng, P. quinquefolius, and P. notoginseng by Ultra-High Performance Liquid Chromatography/Ion Mobility-Quadrupole Time-of-Flight Mass Spectrometry. Molecules 2022, 27, 5549. [Google Scholar]
  10. In, G.; Seo, H.K.; Park, H.W.; Jang, K.H. A Metabolomic Approach for the Discrimination of Red Ginseng Root Parts and Targeted Validation. Molecules 2017, 22, 471. [Google Scholar] [CrossRef] [Green Version]
  11. Zhu, H.; Lin, H.; Tan, J.; Wang, C.; Wang, H.; Wu, F.; Dong, Q.; Liu, Y.; Li, P.; Liu, J. UPLC-QTOF/MS-Based Nontargeted Metabolomic Analysis of Mountain- and Garden-Cultivated Ginseng of Different Ages in Northeast China. Molecules 2018, 24, 33. [Google Scholar] [CrossRef] [Green Version]
  12. Shin, J.S.; Park, H.W.; In, G.; Seo, H.K.; Won, T.H.; Jang, K.H.; Cho, B.G.; Han, C.K.; Shin, J. Metabolomic Approach for Discrimination of Four- and Six-Year-Old Red Ginseng (Panax ginseng) Using UPLC-QToF-MS. Chem. Pharm. Bull. 2016, 64, 1298–1303. [Google Scholar] [CrossRef] [Green Version]
  13. Hou, M.; Wang, R.; Zhao, S.; Wang, Z. Ginsenosides in Panax genus and their biosynthesis. Acta Pharm. Sin. B 2021, 11, 1813–1834. [Google Scholar] [PubMed]
  14. Shi, Z.Y.; Zeng, J.Z.; Wong, A.S.T. Chemical Structures and Pharmacological Profiles of Ginseng Saponins. Molecules 2019, 24, 2443. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Yue, J.; Zuo, Z.; Huang, H.; Wang, Y. Application of Identification and Evaluation Techniques for Ethnobotanical Medicinal Plant of Genus Panax: A Review. Crit. Rev. Anal. Chem. 2021, 51, 373–398. [Google Scholar] [CrossRef] [PubMed]
  16. Wang, X.J.; Xie, Q.; Liu, Y.; Jiang, S.; Li, W.; Li, B.; Wang, W.; Liu, C.X. Panax japonicus and chikusetsusaponins: A review of diverse biological activities and pharmacology mechanism. Chin. Herb. Med. 2021, 13, 64–77. [Google Scholar]
  17. Du, Z.; Li, J.; Zhang, X.; Pei, J.; Huang, L. An Integrated LC-MS-Based Strategy for the Quality Assessment and Discrimination of Three Panax Species. Molecules 2018, 23, 2988. [Google Scholar] [CrossRef] [Green Version]
  18. Liu, J.; Wang, H.D.; Yang, F.F.; Chen, B.X.; Li, X.; Huang, Q.X.; Li, J.; Li, X.Y.; Li, Z.; Yu, H.S.; et al. Multi-level fingerprinting and cardiomyocyte protection evaluation for comparing polysaccharides from six Panax herbal medicines. Carbohydr. Polym. 2022, 277, 118867. [Google Scholar] [CrossRef]
  19. Guo, D.; Wu, W.; Ye, M.; Liu, X.; Cordell, G. A holistic approach to the quality control of traditional Chinese medicines. Science 2015, 347, S29–S31. [Google Scholar]
  20. Alseekh, S.; Aharoni, A.; Brotman, Y.; Contrepois, K.; D’Auria, J.; Ewald, J.C.; Fraser, P.D.; Giavalisco, P.; Hall, R.D.; Heinemann, M.; et al. Mass spectrometry-based metabolomics: A guide for annotation, quantification and best reporting practices. Nat. Methods 2021, 18, 747–756. [Google Scholar]
  21. Ren, J.L.; Zhang, A.H.; Kong, L.; Han, Y.; Yan, G.L.; Sun, H.; Wang, X.J. Analytical strategies for the discovery and validation of quality-markers of traditional Chinese medicine. Phytomedicine 2020, 67, 153165. [Google Scholar] [CrossRef] [PubMed]
  22. Li, M.M.; Wang, Y.K.; Sheng, Y.H.; Huang, X.; Yue, H.; Liu, S.Y. Application of UPLC-QTOF-MS non-targeted metabonomics in mechanism study of property differences of ginseng and American ginseng. Zhongguo Zhong Yao Za Zhi 2021, 46, 5930–5935. [Google Scholar] [PubMed]
  23. Fan, M.; Fan, Y.; Rao, Z.; Li, Y.; Qian, H.; Zhang, H.; Wu, G.; Qi, X.; Wang, L. Comparative investigation on metabolite changes in ‘wu mi’ production by Vaccinium bracteatum Thunb. leaves based on multivariate data analysis using UPLC-QToF-MS. Food Chem. 2019, 286, 146–153. [Google Scholar] [CrossRef] [PubMed]
  24. Yoon, D.; Shin, W.C.; Oh, S.M.; Choi, B.R.; Young Lee, D. Integration of multiplatform metabolomics and multivariate analysis for geographical origin discrimination of Panax ginseng. Food Res. Int. 2022, 159, 111610. [Google Scholar] [CrossRef]
  25. Wang, L.; Huang, Y.; Yin, G.; Wang, J.; Wang, P.; Chen, Z.Y.; Wang, T.; Ren, G. Antimicrobial activities of Asian ginseng, American ginseng, and notoginseng. Phytother. Res. 2020, 34, 1226–1236. [Google Scholar] [CrossRef]
  26. Shim, H.; Waminal, N.E.; Kim, H.H.; Yang, T.J. Dynamic evolution of Panax species. Genes Genom. 2021, 43, 209–215. [Google Scholar] [CrossRef] [PubMed]
  27. Du, Y.; Fan, P.; Zou, L.; Jiang, Y.; Gu, X.; Yu, J.; Zhang, C. Serum Metabolomics Study of Papillary Thyroid Carcinoma Based on HPLC-Q-TOF-MS/MS. Front. Cell Dev. Biol. 2021, 9, 593510. [Google Scholar] [CrossRef]
  28. Meng, Q.; Zhong, Y.M.; Guo, X.L.; FENG, Y.F. Rapid Identification and Preliminary Study of Fragmentation Regularity of Saponins Components of Ginseng. J. Chin. Med. Mater. 2013, 36, 240–245. [Google Scholar]
  29. Wang, Q.; Rui, Y.F.; Zhong, Y.M.; Guo, J. Rapid Analysis of Saponins in Panax notoginseng by UPLC/Q-TOF MS. Nat. Prod. Res. Dev. 2014, 26, 1233–1239. [Google Scholar]
  30. Fang, L.X.; Xiong, A.Z.; Wang, R.; Ji, S.; Yang, L.; Wang, Z.T. A strategy for screening and identifying mycotoxins in herbal medicine using ultra-performance liquid chromatography with tandem quadrupole time-of-flight mass spectrometry. J. Sep. Sci. 2013, 36, 3115–3122. [Google Scholar] [CrossRef]
  31. Wang, F.X.; Yuan, J.C.; Kang, L.P.; Pang, X.; Yan, R.Y.; Zhao, Y.; Zhang, J.; Sun, X.G.; Ma, B.P. An efficient approach to identify different chemical markers between fibrous root and rhizome of Anemarrhena asphodeloides by ultra high-performance liquid chromatography quadrupole time-of-flight tandem mass spectrometry with multivariate statistical analysis. J. Pharm. Biomed. Anal. 2016, 129, 105–116. [Google Scholar] [PubMed]
  32. Robotti, E.; Marengo, E. Chemometric Multivariate Tools for Candidate Biomarker Identification: LDA, PLS-DA, SIMCA, Ranking-PCA. Methods Mol. Biol. 2016, 1384, 237–267. [Google Scholar] [PubMed]
  33. Qiu, S.; Yang, W.Z.; Yao, C.L.; Qiu, Z.D.; Shi, X.J.; Zhang, J.X.; Hou, J.J.; Wang, Q.R.; Wu, W.Y.; Guo, D.A. Nontargeted metabolomic analysis and “commercial-homophyletic” comparison-induced biomarkers verification for the systematic chemical differentiation of five different parts of Panax ginseng. J. Chromatogr. A 2016, 1453, 78–87. [Google Scholar] [CrossRef] [PubMed]
Molecules 28 02745 g001 550
Figure 1. (A) Properties and (B) efficacy based on TCM theory for ginseng (RS), American ginseng (XYS), sanqi (SQ), zhujieshen (ZJS), and zhuzishen (ZZS).
Figure 1. (A) Properties and (B) efficacy based on TCM theory for ginseng (RS), American ginseng (XYS), sanqi (SQ), zhujieshen (ZJS), and zhuzishen (ZZS).
Molecules 28 02745 g001
Molecules 28 02745 g002 550
Figure 2. Base peak chromatograms for the representative samples of five Panax species. (A) In positive ion mode; (B) in negative ion mode.
Figure 2. Base peak chromatograms for the representative samples of five Panax species. (A) In positive ion mode; (B) in negative ion mode.
Molecules 28 02745 g002
Molecules 28 02745 g003 550
Figure 3. The PCA (A,B), PLS-DA(C,D), and OPLS-DA(E,F) score plots of five Panax species (RS, XYS, SQ, ZJS, and ZZS). (A,C,E) represents the score plots in positive ion mode; (B,D,F) represents the score plot in negative ion mode.
Figure 3. The PCA (A,B), PLS-DA(C,D), and OPLS-DA(E,F) score plots of five Panax species (RS, XYS, SQ, ZJS, and ZZS). (A,C,E) represents the score plots in positive ion mode; (B,D,F) represents the score plot in negative ion mode.
Molecules 28 02745 g003
Molecules 28 02745 g004 550
Figure 4. The VIP plots of five Panax species (RS, XYS, SQ, ZJS, and ZZS). (a) Represents the VIP plot in positive ion mode and (b) represents the VIP plot in negative ion mode.
Figure 4. The VIP plots of five Panax species (RS, XYS, SQ, ZJS, and ZZS). (a) Represents the VIP plot in positive ion mode and (b) represents the VIP plot in negative ion mode.
Molecules 28 02745 g004
Molecules 28 02745 g005 550
Figure 5. Heatmap visualization of metabolomics data with hierarchical clustering analysis (HCA). The red color represents the peak value that is relatively large; the blue color represents the peak value that is relatively small. The more similar the color, the more similar the peak value. The units in the abscissa axis represent the sample names and their groups; the panel on the right represents the different metabolites. The upper dendritic structure is clustered according to the degree of metabolite similarity across samples.
Figure 5. Heatmap visualization of metabolomics data with hierarchical clustering analysis (HCA). The red color represents the peak value that is relatively large; the blue color represents the peak value that is relatively small. The more similar the color, the more similar the peak value. The units in the abscissa axis represent the sample names and their groups; the panel on the right represents the different metabolites. The upper dendritic structure is clustered according to the degree of metabolite similarity across samples.
Molecules 28 02745 g005
Molecules 28 02745 g006 550
Figure 6. Bar charts illustrating the content differences in 12 important ginsenoside markers among RS, XYS, SQ, ZJS, and ZZS. (* p < 0.05, ** p < 0.01).
Figure 6. Bar charts illustrating the content differences in 12 important ginsenoside markers among RS, XYS, SQ, ZJS, and ZZS. (* p < 0.05, ** p < 0.01).
Molecules 28 02745 g006
Molecules 28 02745 g007 550
Figure 7. Appearance of P. ginseng (A), P. quinquefolius (B), P. notoginseng (C), P. japonicus (D), and P. japonicus var. major (E).
Figure 7. Appearance of P. ginseng (A), P. quinquefolius (B), P. notoginseng (C), P. japonicus (D), and P. japonicus var. major (E).
Molecules 28 02745 g007
Table
Table 1. The information of potential chemical marker diagnostics for differentiating among RS, XYS, ZJS, ZZS, and SQ.
Table 1. The information of potential chemical marker diagnostics for differentiating among RS, XYS, ZJS, ZZS, and SQ.
NO.Ion ModeAdductsVIPtR (min)Observed m/zError ppmFormulaIdentificationRSXYSZJSZZSSQ
1+M+H5.430.49175.11739C6H14N4O2L-Arginine++++++
2-M-H4.820.6387.10743C20H20O85-Hydroxy-3,3′,4′,7,8-pentamethoxyflavone++++-
3-M-H5.940.61341.10038C19H18O6Tetramethoxyluteolin++++++++
4+M+H-2H2O5.950.64360.14622C19H25NO8Ethyl 2-hydroxy-3-(3-indolyl)propanoate glucoside+++++++
5+M+NH45.022.41372.18572C14H26O10Isopropyl apiosylglucoside--++-
6+/9.692.81322.2084//Unknown 1++++++-
7+/8.562.94386.1982//Unknown 2--++--
8+/9.723.07363.169//Unknown 3+++++++++
9+/4.353.68400.2183//Unknown 4-+---
10+M+CH3OH+H5.673.82579.29933C26H42O12Canavalioside++++-
11+M+H4.424.42407.200714C22H30O71-(4-Hydroxy-3,5-dimethoxyphenyl)-7-(4-hydroxy-3-methoxyphenyl)-3,5-heptanediol+-++++
12+M+NH4-H2O6.084.65381.18347C15H27NO10Pantothenic acid 4′-O-b-D-glucoside----+
13+M+H-2H2O8.074.78365.18710C22H28N2O5Rotundifoline----++
14-/5.314.78262.1171//Unknown 5----++
15-M-H4.085.5961.53820C48H82O19Ginsenoside M6a----+
16-M-H7.735.63931.51929C47H80O18Notoginsenoside R1----+++
17+M+H8.575.73817.482515C42H72O15Majonoside R1--+++-
18-M-H4.915.83799.48015C42H71O14Ginsenoside Rg1----++
19-/7.135.88991.5361//Unknown 6++++++++-
20-M+FA-H4.916.05845.48912C43H75O15Ginsenoside Rg2---++
21+M+NH47.296.06904.516211C45H74O17Melongoside G+-+-+
22+M+ACN+Na6.576.31580.32552C30H44O7Ganoderic acid--++-
23+/8.56.8817.4781//Unknown 7--++++-
24-M+Hac-H8.036.81861.469518C41H70O15Vinaginsenoside R14--+++++-
25+M+2ACN+H4.0171259.681322C57H92O25Asterlingulatoside D+---+
26-M-H4.467.031239.64586C59H100O27Ginsenoside Ra3----+
27+M+2H5.67.16563.30384C54H92O24Koryoginsenoside Rg2--++-
28-M+FA-H8.197.16845.474219C43H75O15Ginsenoside Rf++++++++-
29-M-H4.487.3769.47095C41H70O13Ginsenoside F3----+
30-M+F4.537.351193.589925C56H90O27Agavasaponin D++++++-
31+M+NH44.087.441126.63551C54H92O23Ginsenoside Rb1-+--+
32-M-H4.067.65783.48763C42H72O13Ginsenoside Rg3----+
33-M+Na-2H5.317.82925.468310C45H76O18Torvoside--++-
34+M+NH45.118.131168.64861C56H94O24Quinquenoside R1-+---
35-/6.368.13793.4239//Unknown 8--++++-
36-M+Hac-H6.358.361031.534327C48H76O20Calendasaponin B++++--
37+M+NH45.028.43964.5822C48H82O18Ginsenoside B2-+--+
38-M-H7.228.45945.53369C48H81O18Ginsenoside Rd----+++
39+/4.698.661050.5853//Unknown 9+---+
40+M+NH45.4710.32330.26225C18H32O4octadec-9-enedioic acid++--+
41+M+NH45.3510.88316.281310C18H34O312-Hydroxyoctadec-9-enoic acid++-+-+
42-M-H+HCOONa4.5811.78295.22357C15H32OPentadecanol++---
43+/5.5112.1622.3771/C36H60O8Ginsenoside Rh4++--+
44+M+H4.0712.37496.33774C24H50NO7P1-Palmitoylphosphatidylcholine++++-
45+/6.6914.64378.2217//Unknown 10--++++-
46-M-H4.3315.07271.224512C16H32O3Hydroxyhexadecanoic acid+---+
47+M+H9.9515.32380.2367C23H29N3O24-Isopropylphenserine-++++++-
48-M-H+HCOONa5.2515.79279.22944C15H32Pentadecane++++++
49-/6.4616.82833.5103//Unknown 11--+++++-
+++: the content of the compound is very high, peak area ≥ 260; ++: the content of the compound is high, peak area ≥ 130 but <260; +: the content of the compound is low, peak area ≥ 10 but <130; - the content of the compound is very low or even none, peak area < 10.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ji, R.; Garran, T.A.; Luo, Y.; Cheng, M.; Ren, M.; Zhou, X. Untargeted Metabolomic Analysis and Chemometrics to Identify Potential Marker Compounds for the Chemical Differentiation of Panax ginseng, P. quinquefolius, P. notoginseng, P. japonicus, and P. japonicus var. major. Molecules 2023, 28, 2745. https://doi.org/10.3390/molecules28062745

AMA Style

Ji R, Garran TA, Luo Y, Cheng M, Ren M, Zhou X. Untargeted Metabolomic Analysis and Chemometrics to Identify Potential Marker Compounds for the Chemical Differentiation of Panax ginseng, P. quinquefolius, P. notoginseng, P. japonicus, and P. japonicus var. major. Molecules. 2023; 28(6):2745. https://doi.org/10.3390/molecules28062745

Chicago/Turabian Style

Ji, Ruifeng, Thomas Avery Garran, Yilu Luo, Meng Cheng, Mengyue Ren, and Xiuteng Zhou. 2023. "Untargeted Metabolomic Analysis and Chemometrics to Identify Potential Marker Compounds for the Chemical Differentiation of Panax ginseng, P. quinquefolius, P. notoginseng, P. japonicus, and P. japonicus var. major" Molecules 28, no. 6: 2745. https://doi.org/10.3390/molecules28062745

APA Style

Ji, R., Garran, T. A., Luo, Y., Cheng, M., Ren, M., & Zhou, X. (2023). Untargeted Metabolomic Analysis and Chemometrics to Identify Potential Marker Compounds for the Chemical Differentiation of Panax ginseng, P. quinquefolius, P. notoginseng, P. japonicus, and P. japonicus var. major. Molecules, 28(6), 2745. https://doi.org/10.3390/molecules28062745

Article Metrics

Back to TopTop