Abstract

Triterpenoids are the major secondary metabolites and active substances in Ganoderma, considered as the “marker compounds” for the chemical evaluation or standardization of Ganoderma. A response surface methodology was used to optimize the ultrasonic-assisted extraction of triterpenoids. The extraction rate was 7.338 ± 0.150 mg/g under the optimum conditions: 87% ethanol, ratio of solid to liquid (w : v) 1 : 28, and ultrasound extraction time 36 min. Based on the high sensitivity and selectivity of HPLC-LTQ-Orbitrap-MSn, 24 components of triterpenoids were tentatively identified in the negative mode. Then, the global chemical profiling consisting of HPLC and TLC fingerprints generated by ChemPattern™ software was developed for evaluation of Ganoderma species. For fingerprint analysis, 11 peaks of triterpenoids were selected as the characteristic peaks to evaluate the similarities of different samples. The correlation coefficients of similarity were greater than 0.830. The cluster analysis showed a clear separation of three groups, and 11 peaks played key roles in differentiating these samples. The developed global chemical profiling method could be applied for rapid evaluation, quality control, and authenticity identification of Ganoderma and other herbal medicines.

1. Introduction

Ganoderma, a popular edible and medicinal mushroom, is commonly used as dietary supplements and touted as a remedy to promote health and longevity [1]. So far, more than 200 species of Ganoderma have been found in the world, and Ganoderma lucidum (Leyss. ex Fr.) Karst. and Ganoderma sinense Zhao, Xu et Zhang are officially recorded in Chinese pharmacopoeia [2]. Previous studies have demonstrated that Ganoderma possesses various biological properties, such as antitumor, antiaging, antioxidant, hyperglycemic, and regulating immunity [35]. Owing to its satisfactory clinical effects, more and more Ganoderma products as health foods or medicines have appeared at the market. However, it is hard to say whether its quality is good or bad, and it is also difficult to identify whether the raw materials are authentic or adulterant.

For Ganoderma, triterpenoids are the major secondary metabolites and active substances [6]. The pioneering has isolated and identified more than 300 triterpenoids from the spores, fruiting bodies, and cultured mycelia of Ganoderma [79]. Triterpenoids could be considered as the “marker compounds” for the chemical evaluation or standardization of Ganoderma. Professor Guo and his team’s researches concluded that the content and composition of triterpenoids vary significantly due to difference in the strain, geographic origin, cultivation method, extraction process, and other factors [10, 11]. Besides, scholars have made a lot of exploratory work on the chromatographic fingerprint of Ganoderma products [12, 13]. However, the traditional fingerprints were always processed by the fingerprint similarity software (2004 or 2012 version), only limited to HPLC profiles without TLC. After that, researchers always apply other statistical software (such as SPSS and SAS) to process the data for cluster analysis or principal component analysis, which is relatively cumbersome and time-consuming. Therefore, it is essential to develop a global chemical profiling method for rapidly evaluating the quality of Ganoderma to ensure the efficacy.

Chemical profile reflects the totality of intrinsic chemical compounds of herbal medicines and emphasizes the integral characterization of a complex system [14, 15]. Here, the global chemical profiling method contained the establishment of HPLC and TLC fingerprints, the characterization of common peaks and statistical analysis. For processing the large scale of physical properties characteristic in the complex herbal extract, an advanced chemometric and chemical fingerprinting software was developed by Chemmind Technology Co., Ltd. ChemPattern is an advanced chemometric and chemical fingerprinting software, which endeavors to provide solutions for qualitative and quantitative quality evaluation and characteristic analysis. The software was employed to calculate the correlation coefficients between different chromatographic profiles, as well as to generate the representative standard fingerprint by mean simulation.

In the present work, triterpenoids were extracted from different Ganoderma samples under the ultrasonic-assisted condition optimized by response surface methodology (RSM) with the Box-Behnken design (BBD) [16]. The extract was analyzed by high-performance liquid chromatography coupled with linear ion trap-Orbitrap mass spectrometry (HPLC-LTQ-Orbitrap-MSn) for a comprehensive study of the multiple chemical constituents. In addition, the chromatographic data of HPLC and TLC were submitted into the professional ChemPattern software for establishing global chemical profiling and evaluating similarity. Furthermore, Ganoderma samples from different regions could be distinguished by clustering analysis and principal component analysis of fingerprint data.

2. Experimental

2.1. Materials and Reagents

Ganoderma collected from different regions is shown in Figure 1. Samples were pulverized into powder and kept in a vacuum dryer. The standard sample of Ganoderma was identified by National Institutes for Food and Drug Control. Cellulase and oleanolic acid were obtained from Aladdin. Ganoderic acid A was purchased from Chengdu Must Bio-technology Co., Ltd. Acetonitrile and methanol were of HPLC grade and obtained from Anhui Tedia High Purity Solvents Co., Ltd. Ethanol, vanillin, sulfuric acid, acetic acid, perchloric acid, phosphoric acid, formic acid, petroleum ether, ethyl acetate, and other chemicals were of analytical grade.

2.2. Heating Reflux Extraction

Sample powder (2.0 g) and 95% ethanol (30 mL) were extracted under reflux at 75°C for 30 min. Then, the extraction solution was filtered through a filter paper and evaporated to dryness at 60°C.

2.3. Ultrasonic Extraction of Total Triterpenoids

Sample powder (2.0 g) and 95% ethanol (30 mL) were placed in an ultrasonic bath (300 W) at 75°C for 30 min. The suspension was cooled to room temperature and filtered. The filtrate was vacuum-dried at 60°C.

2.4. Determination of the Total Content of Triterpenoids

The total content of triterpenoids was determined according to the method of Hou with some modifications [17]. The oleanolic acid (2.0 mg) was dissolved in methanol (10 mL) to produce a standard solution. 0 μL, 100 μL, 200 μL, 300 μL, 400 μL, 500 μL, and 600 μL standard solutions were added into a test tube, respectively, and evaporated in a water bath. 5% vanillin-acetic acid reagent (400 μL) and perchloric acid (1000 μL) were added, and the tube was placed in a water bath for 30 min at 65°C. When the reaction solution was cooled, acetic acid (5 mL) was added. The absorbance was determined at 546 nm using a microplate reader (Infinite®200 Pro NanoQuant, Tecan, Switzerland). The sample (100 μL) was determined following the abovementioned method. Then, the weight of triterpenoids was calculated according to the standard curve. The extraction yield of triterpenoids was calculated as follows: yield (mg/g, w/w) = weight of triterpenoids/weight of raw materials. All determinations were performed in triplicates. The standard curve was y = 223.2x + 4.473, R2 = 0.9950. The oleanolic acid weight in 20–120 μg range showed a good linear relationship.

2.5. Experimental Design of RSM for Ultrasonic Extraction

Three independent factors of ultrasonic extraction were investigated using the RSM of BBD, including liquid-solid ratio (A: 15, 20, and 35 mL/g), extraction time (B: 35, 45, and 60 min), and ethanol concentration (C: 75, 85, and 95%), as shown in Table S1. The three levels were designated as −1, 0, and +1 for low, intermediate, and high values.

In order to predict the conditions of ultrasound extraction, experimental data were analyzed using the software Design-Expert version 8.06 and explained using the following nonlinear computer-generated quadratic model [18]:where β0 is the constant coefficient, βi, βii, and βij are the coefficients for the linear, quadratic, and interaction effect, xi and xj are the independent variables, and ε is the error.

The adequacy of the model was tested through analysis of variance (ANOVA). The coefficients of determination R2 and adj R2 expressed the quality of fit of the resultant polynomial model, and the statistical significance was checked by F-value and lack of fit [19].

2.6. LC-LTQ-Orbitrap-MSn Conditions

For accurate mass measurements, an Agilent 1290 HPLC instrument was coupled with a LTQ Orbitrap Velos mass spectrometer (Thermo Scientific, Hemel Hempstead, UK) equipped with an ESI source. An Eclipse plus C18 (50 mm × 4.6 mm, 1.8 μm) column was used for chromatographic separation. The column temperature was kept at 35°C. The mobile phase consisted of acetonitrile (A) and 0.03% phosphoric acid solution (B). The gradient elution was as follows: 75–68% B at 0–15 min, 68–60% B at 15–20 min, 60–40% B at 20–25 min, 40–0% B at 25–40 min, and 0% B at 40–125 min. The DAD was set at 254 nm. The injection volume was 10 μL, and the flow rate was 0.6 mL/min.

The operation parameters of mass spectrometry were as follows: source voltage, 4.0 kV; sheath gas, 20 (arbitrary units); auxiliary gas, 12 (arbitrary units); sweep gas, 2 (arbitrary units); and capillary temperature, 350°C. Default values were used for most other acquisition parameters: Fourier transformation (FT) automatic gain control (AGC) target 5 × 105 for the MS mode and 5 × 104 for the MSn mode. Perfusion samples were analyzed in the data-dependent scan mode at a resolving power of 60,000 at m/z 400. The most intense ions were selected, and parent ions were fragmented by high-energy C-trap dissociation (HCD) with a normalized collision energy of 45% and an activation time of 100 ms. The maximum injection time was set to 100 ms with two microscans for the MS mode and to 1000 ms with one microscan for the MSn mode. The mass range was from m/z 100 to 1500. Each sample was analyzed both in negative and positive modes. Data were analyzed using Xcalibur software version 2.2 (Thermo Fisher Scientific).

2.7. HPLC Chromatographic Fingerprint Analysis Conditions

HPLC chromatographic fingerprint analysis was conducted on a liquid chromatography system (1260, Agilent, America) equipped with a quaternary solvent deliver system, an autosampler, and a DAD (Agilent Technologies). The mobile phase consisted of acetonitrile (A) and 0.03% phosphoric acid solution (B) using a gradient elution of 75–68% B at 0–40 min, 68–60% B at 40–60 min, 60–40% B at 80–120 min, 40–0% B at 80–120 min, and 0% B at 120–125 min. Chromatographic separation was carried out at an Agilent Zorbax Extend-C18 column (4.6 mm × 250 mm, 5 μm) with a solvent flow rate of 1.0 mL/min at a temperature of 35°C. The wavelength was set at 254 nm. The injection volume was 10 μL.

2.8. TLC Chromatographic Conditions

Ganoderic acid A solution and samples a to k (5 μL) were spotted using a microinjector on a 20 × 20 cm silica gel plate (GF254, Qingdao, China). The silica gel was activated in an oven at 80°C for 30 min before use. The mobile phase consisting of petroleum ether : ethyl acetate : formic acid (1 : 1 : 0.02, v/v/v) was added into a twin-trough chamber and saturated for 10 min. The plate in the chamber was developed upward over a path of 15 cm and sprayed with 1% vanillin-sulfuric acid solution. The plate was placed in an oven at 80°C for 10 min until the color of the triterpenoid spots was distinct. The image of TLC was reverse-phase processed, and the information of TLC spots was turned into the gray curve by software.

2.9. Statistical Analysis

All experiments were performed at least in triplicate. The values were expressed as means ± standard deviation (SD).

3. Results and Discussion

3.1. Comparison of Ultrasonic-Assisted Extraction and Reflux Extraction

Comparing heating reflux extraction (6.404 mg/g) with ultrasonic-assisted extraction (6.869 mg/g), the extraction yield of triterpenoids was increased significantly by ultrasonic-assisted extraction in the same extraction time with simple operation. The ultrasonic wave produced a strong cavitation, mechanical crushing, and thermal effect, which dissolved the active ingredients into the solvent more adequately and saved more energy. So, the ultrasonic-assisted extraction was selected for further optimization.

3.2. RSM for Optimization of Ultrasonic-Assisted Extraction

The main factors which affected the ultrasonic extraction yield of triterpenoids were investigated by single-factor experiments, including the temperature, liquid-solid ratio, extraction time, and ethanol concentration. When the procedures were conducted at 30 min with a concentration of 95% ethanol, the maximum yield of triterpenoids was 0.6906% at the liquid-solid ratio 25 mL/g (Figure 2(a)). The liquid-solid ratio and concentration of ethanol were fixed at 25 mL/g and 95%. The result showed that the extraction efficiency increased to the maximum amount of 0.7097% at 45 min (Figure 2(b)). The yield of triterpenoids was significantly increased with the ethanol concentration varying from 15% to 95%, and the optimal ethanol concentration was from 75% to 95% (Figure 2(c)). Therefore, the liquid-solid ratio 25 mL/g, extraction time 45 min, and 85% ethanol were selected as the center points for each factor in the RSM experiments. The single-factor experiment of temperature showed that the extraction rate of triterpenoids was almost unchanged with increasing temperature after 50°C. So, the extraction temperature was not considered as the experimental factor of BBD and set at 50°C. An experimental program for optimizing the extraction of triterpenoids using RSM with BBD is shown in Table S2. The predicted values were obtained from the model fitting technique using the software Design-Expert version 8.06.

ANOVA was applied to optimize the extraction conditions of ultrasonic-assisted extraction for the triterpenoid yield and evaluate the relationship between response and variables. ANOVA for the response surface quadratic regression model showed that the F-value of model was 30.56 and the value of model was smaller than 0.0001 (Table S3), suggesting the model was significant [16]. The coefficient of determination (R2) of the model was 0.9752, and the adjusted determination coefficient (adj R2) was 0.9433, which indicated good agreement between the experimental parameters and the predicted values of triterpenoids. The sequence of three factors influencing the triterpenoid yield was the ethanol concentration (C), extraction time (B), and liquid-solid ratio (A).

By statistically processing, the multiple second-order equation for the extraction yield of triterpenoids was obtained as follows:where is the extraction yield of triterpenoids, , , and are the coded values of the liquid-solid ratio, extraction time, and ethanol concentration, respectively.

The extraction yield and the interaction of different variables could be predicted from the three-dimensional (3D) response surface (Figures 2(d)2(f)). The steeper slope represented that the factor had a more significant effect on the extraction yield. The steeper slope of ethanol concentration implied that it had greater effect on the extraction yield of triterpenoids than on the extraction time and liquid-solid ratio.

According to the regression equation, the optimal parameters were the liquid-solid ratio 28.32 mL/g, extraction time 35.64 min, and ethanol concentration 87.25%. The theoretical highest yield of triterpenoids was 7.3118 mg/g predicted by the model. In order to validate the adequacy of the model, verification experiments were carried out by slightly modified conditions: liquid-solid ratio 28 mL/g, extraction time 36 min, and ethanol concentration 87%. The yield of triterpenoids of 7.338 ± 0.150 mg/g could be attained, which was 6.83% higher than the previous ultrasonic extraction method.

3.3. Characterization and Identification of Triterpenoids by LC-LTQ-Orbitrap-MSn

Qualitative analysis of triterpenoids was performed on the HPLC-LTQ-Orbitrap-MSn system. ESI-MS spectra in both negative and positive modes were examined in this study. Negative-mode ESI was found to be sensitive for triterpenoids. All triterpenoids gave [M – H] ions in their negative ion mass spectra. The total ion chromatograms (TICs) of triterpenoids in the negative ion mode by LC-LTQ-Orbitrap-MSn are shown in Figure 3. The fragmentation pathway of triterpenoids is summarized by using ganoderic acid A as the standard compound. The mass spectrum of ganoderic acid A and its major fragmentation pathways are given in Figure 4. In the negative ion mode, the prominent fragmentation pathways begin with the prominent losses of H2O or CO2; then, the cleavages took place on the A, C, and D rings. The [M – H] ion at m/z 515.30 () of ganoderic acid A produced a prominent ion at m/z 497.34 by eliminating a molecule of H2O (18 Da). The m/z 497.34 ion was further subjected to produce signals at m/z 479.32 or 453.36 by the sequential losses of H2O or CO2 (44 Da). The [M – H] ion at m/z 515.30 produced an ion at m/z 417.32 () by direct cleaving on ring A. The ion at m/z 355.29 could be obtained by the process of losing H2O and CO2, followed by cleavage of ring A. Ganoderic acid A was also cleaved on ring C to give the product ion at m/z 249.12 (). The cleavage of ring D could be observed in ganoderic acid A, besides the cleavage of rings A and C. The [M – H] ion at m/z 515.30 produced an ion at m/z 301.25 (). The m/z ion 301.25 then underwent losses of CH3 (15 Da) to generate an ion at m/z 285.31.

According to the nontarget compound identification strategy based on the accurate mass measurement (<5 ppm), MS/MS fragmentation patterns, diagnostic product ions, and different chromatographic behaviors, 24 compounds were unambiguously identified from triterpenoids [9, 10]. Table 1 summarizes the retention times (tR), molecular formula, [M – H] and CAS number of each compound, and MS/MS ions. The structures of 24 compounds are shown in Figure S1. These results provided the critical information for constructing chemical fingerprints of triterpenoids.

3.4. Validation of HPLC and TLC Methods

Prior to the establishment of the HPLC fingerprint, the precision, repeatability, and stability were chosen to validate the reliability of HPLC, which were expressed by the relative standard deviations (RSDs) of the retention time (tR) and peak area (Pa). For the precision test, the working solutions were analyzed in triplicate, and RSD values of tR and Pa were lower than 0.2% and 4% (Table 2). To confirm the repeatability, five different working solutions prepared from the same batch of the sample were analyzed. The repeatability test (Table 2) demonstrated that the developed assay was reproducible (RSD < 5%). Stability of 11 analytes was detected at 2, 4, 8, 12, and 24 h, respectively. RSD values of Pa for stability tests were 1.46–4.81% (Table 2), which indicated that the sample was stable in 24 h.

In contrast with HPLC, TLC fingerprints were more cost-effective and provided a vivid colorful image for parallel comparison [20]. At present, there are few literatures about TLC fingerprints applied to Ganoderma. In this study, 10 batches of Ganoderma and the standard herb were further investigated by TLC fingerprints. In the preliminary experiment, developing solvent, sample concentration, chromogenic reagent, chromogenic temperature, and chromogenic time were optimized to achieve the optimum effect of separation and coloration. Owing to the large polarity of triterpenoids, a small amount of formic acid was added in the developing agent. Then, 1% vanillin-sulfuric acid solution was chosen as a chromogenic agent to color the compounds of triterpenoids. In order to validate the reliability of the TLC method, precision, repeatability, and stability were determined. The operation methods of working solutions were the same as HPLC. As shown in Table 3, RSDs of the flow rate value (Rf) ranged from 0.33% to 2.48%, which proved that the TLC experiments were reliable. However, RSDs of the peak area value (Pa) of the main characteristic peaks were from 2.32% to 7.79% (Table 3), inferring that a certain amount of error was produced by manual spotting.

3.5. Establishment of Global Chemical Profile by ChemPattern Software

Both common patterns of HPLC and TLC fingerprints for triterpenoids were generated to represent the characteristic peaks of this authenticated herbal medicine by ChemPattern software (Chemmind Technologies, Beijing, China). The similarity and clustering analyses on different kinds of chromatographic fingerprint data comprehensively evaluated the types and quantities of triterpenoids by ChemPattern software.

The HPLC common pattern of 10 batches of samples is shown in Figures 5(a) and 5(b). Sample c from America was not suitable to establish the common model of the HPLC fingerprint because of its significant differences in chemical composition. 11 peaks existing in 10 batches of samples were found as common peaks. According to the compound database of triterpenoids obtained from HPLC-LTQ-Orbitrap-MSn, peaks 1 to 11 were identified as lucidenic acid LM1, ganoderic acid G, ganoderic acid B, lucidenic acid E, ganoderenic acid A, ganoderic acid A, lucidenic acid A, ganoderenic acid D, ganoderic acid D, lucidenic acid D, and ganoderic acid F, respectively. The relative retention time ( = retention time of the characteristic peak/retention time of the marker peak) and relative peak area (RPA = peak area of the characteristic peak/peak area of the marker peak) of the common peaks are shown in Tables S4 and S5. Peak 6, identified as ganoderic acid A, was selected as the reference peak. The results indicated that of 11 common peaks (between 0.13% and 0.57%) was invariable between samples, and was a valid parameter for constituent identification. However, the RPA (from 29.19% to 135.20%) showed significant differences between 10 batches of samples, indicating that the content of triterpenoids from various sources was different.

The similarity and clustering analyses of HPLC fingerprints were also analyzed by ChemPattern software. The advanced chemical pattern recognition module of ChemPattern could forecast the classification of complex samples. The similarities of samples a, b, d, e, f, g, h, and i were greater than 0.830 (Figure 5(c)). However, the profiles of samples j and k were different from the common pattern, suggesting that these two samples belonged to Ganoderma sinense. This method could distinguish Ganoderma sinense and Ganoderma lucidum quickly. When the clustering distance was extended to 1.0, the samples from different regions could be divided into three groups. As shown in Figure 5(d), samples f, d, h, b, and a belonged to class I. These five samples were collected from Anhui Province. Samples k, j, g, and e were classified into class II, for the reason that the triterpenoid content and species of these samples were similar. Sample i was classified into class III individually, due to the different categories of triterpenoids. These results illustrated that the origins of Ganoderma could be distinguished via clustering analysis. The HPLC fingerprints could not only reflect the chemistry information of Ganoderma but also distinguish the Ganoderma species from different geographical origins.

The image of TLC taken by a digital camera was reverse-phase processed before imported into ChemPattern software (Figure S2). Twelve tracks were set manually, and the information of TLC spots was turned into the gray curve. TLC fingerprints of different triterpenoid extracts and the common model were obtained as shown in Figures 6(a) and 6(b). The fingerprints exhibited that nine common peaks were likely to represent the major constituents of triterpenoids. The similarity and clustering analyses of TLC by ChemPattern software are shown in Figures 6(c) and 6(d). The similarities of the samples a, b, d, e, f, g, h, i, j, and k were all greater than 0.900, while samples j and k were relatively low, which was consistent with the result of the HPLC fingerprint (Figure 6(c)). The clustering analysis showed that the samples could be divided into three clusters if the Euclidean distance was equal to 25 (Figure 6(d)). The samples e, d, b, f, i, h, and a belonged to class I. Among these samples, only samples e and h were not collected from Anhui. Sample k was classified into class II individually, while samples j and g were classified into class III. The cluster analysis of 10 batches of Ganoderma showed a clear separation of the three groups, and common peaks played key roles in differentiating these samples. Both HPLC and TLC fingerprint classifications could provide a simple reference standard for quality identification of Ganoderma.

3.6. Discussion

In this article, RSM with the BBD method was successfully applied to optimize the factors for the ultrasonic extraction of triterpenoids. The suitable conditions were as follows: liquid-solid ratio 28 mL/g, ethanol concentration 87%, and extraction time 36 min at 50°C. The ultrasonic extraction technology has the advantage of accelerating the extraction time, lowering the temperature, causing less damage to the structure of plant materials, and increasing the extraction yield. It is more suitable for the extraction of triterpenoids. RSM is an effective statistical method useful for optimizing a complex process and evaluating the interaction between multiple parameters [21]. The yield of triterpenoids of 7.338 ± 0.150 mg/g could be attained, which was consistent with the theoretical predicted value (7.3118 mg/g). Therefore, the model was considered to be reliable, and RSM could be used for predicting the ultrasonic extraction yield of triterpenoids.

HPLC-LTQ-Orbitrap-MSn technique was applied to identify the chemical structure of triterpenoids in Ganoderma with higher sensitivity. It provided high resolution and abundant structural information for not only the pseudomolecular ions but also the fragment species. The application of HPLC-LTQ-Orbitrap-MSn could provide a large amount of information related to the chemical structure and make up the deficiency of UV and DAD detectors [8]. Based on the accurate mass measurement, MS/MS fragmentation patterns, and diagnostic product ions provided by HPLC-LTQ-Orbitrap-MSn and literatures, 24 compounds were identified, which was one of the major tools for the study of the chemical substance of Ganoderma.

To investigate triterpenoids further, the chromatographic fingerprint method was applied to analyze the complex composition of Ganoderma. The fingerprint analysis technology is different from traditional analysis methods because it analyzes objects from the perspective of whole component information. The chromatographic fingerprint method was established through importing HPLC and TLC profiles into ChemPattern software for the comprehensive quality control of Ganoderma. The common pattern of fingerprints provided the critical information for constructing chemical fingerprints of triterpenoids. In addition, the fingerprint method combined with the stoichiometry could be applied to identify the species and origins of Ganoderma. The results obtained by clustering analysis of HPLC and TLC fingerprints were slightly different due to the fact that HPLC possesses high resolution and reproducibility. The correct classification percentages of the pattern recognition method of HPLC fingerprint identification were higher than those of TLC. However, in contrast with the HPLC fingerprint, the TLC fingerprint was more cost-effective and provided a vivid colorful image for parallel comparison. In summary, the chemical fingerprint method could be adopted as a reliable tool for the authentication and quality control of Ganoderma.

4. Conclusion

In the present study, the chemical composition of triterpenoids was clarified by HPLC-LTQ-Orbitrap-MSn, and the global chemical profile consisting of HPLC and TLC fingerprints was established. Eleven triterpenoid peaks which differed significantly in all the analyzed samples were used as markers for origin identification and authenticity establishment of Ganoderma. This work suggested that the developed global chemical profiling method could provide a convenient approach, which might be applied for rapid evaluation, quality control, and authenticity establishment of Ganoderma products. In the future, the chemical constituents and pharmacological activities of Ganoderma would be explored in depth through the multidimensional fingerprints combined with chemometric methods, molecular biology, and pattern recognition techniques.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (Grant no. 81603249) and Scientific Research Project of Science Technology Department of Zhejiang Province (no. 2016C33076).

Supplementary Materials

Table S1: Factors and levels in extraction experiments. Table S2: Box-Behnken design of extraction parameters. Table S3: ANOVA for the response surface quadratic model of extraction parameters. Table S4: The a of the common peak in different Ganoderma species. Table S5: The RPAb of the common peak in different Ganoderma species. Figure S1: The structures of compounds identified in triterpenoids from Ganoderma. Figure S2: TLC image of triterpenoids (a) and the reverse-phase processing image (b). Twelve tracks from left to right were ganoderic acid A and samples a, b, c, d, e, f, g, h, i, j, and k, respectively. (Supplementary Materials)