Uncovering Active Constituents Responsible for Different Activities of Raw and Steamed Panax notoginseng Roots

Although Panax notoginseng (PN) roots in raw and steamed forms were historically supposed to be different in the efficacies, the raw materials and steamed ones were often undifferentiated in the use and market circulation, which might bring unstable curative effects or even adverse reactions. To uncover chemical constituents responsible to different activities of raw and steamed PN, chemometrics analyses including partial least squares regression (PLSR) and multi-linear regression analysis (MLRA) were used to establish the relationships between the chromatographic fingerprints and activities of PN samples. Chemical fingerprints of PN were determined by HPLC. Anticoagulant and antioxidant activities were evaluated by the thromboplastin inhibiting test and hydroxyl radical scavenging assay, respectively. Results showed that there was a significant difference in the chemical composition between raw and steamed PN, which could be discriminated by principle component analysis according to different steaming temperatures. Compared with the steamed PN, raw PN exhibited stronger anticoagulation and weaker antioxidation. By chemometrics analyses, notoginsenoside R1, ginsenosides Rg1, Re, Rb1, and Rd were found to be the major active constituents of raw PN, whereas ginsenosides Rh1, Rk3, Rh4, and 20(R)-Rg3 had the key role in the activities of steamed PN, which could be used as new markers for the quality control (QC) of steamed PN.


INTRODUCTION
The processing of herbal medicines, including special crafts of steaming, baking, cooking, and other methods with liquid or solid supplementary materials, plays an important role in the application of traditional Chinese medicine (TCM). The purposes of processing include transforming the properties of medicines, strengthening the curative efficacy, generating new effects and reducing the toxicity or side-effects Li et al., 2013). In recent decades, it has been found that the main mechanism underlying the property changes of herbs is mainly related to the alteration in the chemical composition and/or bioactivity of constituents in herbs (Cao et al., 2011;Li S. L. et al., 2012). For example, diester diterpene alkaloids responsible for the toxicity of Aconitum Radix could be decomposed into less or non-toxic derivatives by boiling the raw materials at 100 • C for 8 h before drying it (Sun et al., 2012). Thus, the uses of raw and processed medicines cannot be mixed up for different curative or toxic effects.
Panax notoginseng (PN) Burk., a plant in genus Panax (Araliaceae), is a well-known medicinal herb used to treat blood disorders for more than 400 years (Ge et al., 2015;Li et al., 2016). Based on the US Dietary Supplement Health and Education Act of 1944, PN and the relative products were also classified as dietary supplements (103rd Congress, 1994). Traditionally, there are two forms of PN, namely raw PN and steamed PN, of which the former one is used as a hemostatic, traumatic, and cardiovascular medicine removing blood stasis, while the steamed one is used as a tonic for nourishing the body and improving the health (Ge et al., 2015;Gu et al., 2015). Pharmacologic studies have also shown that the effects of PN changed when steamed. Lau et al. (2009) reported that the treatment of raw PN extract resulted in shorter bleeding time compared with rats treated with steamed PN, which is consistent with the traditional use of raw PN as a hemostatic. While for steamed PN, it could significantly increase the levels of hemoglobin and white blood cells, as well as the organ index of mice with blood deficiency caused by cyclophosphamide, which were inapparent when treated with raw PN (Zhou et al., 2014). Triterpenoid saponins including ginsenosides and notoginsenosides are considered to be the major bioactive constituents of PN, among which ginsenosides Rg 1 , Rb 1 , Re, and notoginsenoside R 1 show higher levels than others in raw PN (Kim, 2012).
Despite the difference in traditional uses and pharmacological effects, trade in raw and steamed PN has still been complicated by the mixed use. Until now, only raw PN material or powder is recorded in the pharmacopeias of different countries or regions, in which the standard of processed PN is not included yet (British Pharmacopoeia Commission, 2014;Chinese Pharmacopoeia Commission, 2015; European Pharmacopoeia 8.0 1 ; U.S.A Herbal Medicines Compendium 1.0 2 ), indicating that the difference between the raw and steamed PN has not been authorized by codex standards. Even for some institution standards involving steamed PN, the marker constituents for its quality control (QC) are consistent with those of raw PN, despite the significant change in the chemical composition of PN during the steaming process reported in several studies Ge et al., 2015). For example, in the China Food and Drug Administration (2012), the marker constituents for the QC of steamed PN powder are ginsenoside Rb 1 , ginsenoside Rg 1 , and notoginsenoside R 1 , which are same as the markers of raw PN involved in Chinese Pharmacopoeia of 2015 edition. Due to the lack of respective quality standards to differentiate raw and steamed PN, there would be risks of unstable curative effects or even adverse reactions for customers. Therefore, besides the investigation on changes in the chemical composition and pharmacological effects of PN during the steaming process, constituents responsible to different activities should also be uncovered for the individualized QC of raw and steamed PN, which have not been reported yet.
Since the traditional efficacy of removing blood stasis and promoting blood circulation of raw PN could be related to the anticoagulant effect , and the body tonifying function of herbal medicines is partly attributed to their antioxidant and immunomodulatory effects by modern pharmacological researches (Yim and Ko, 2002), the anticoagulation (obtained by thromboplastin inhibiting test) and antioxidation (obtained by hydroxyl radical scavenging assay) activities of PN during the steaming process were studied in this research. Meanwhile, we developed the HPLC chromatographic fingerprints of PN under different steaming conditions, and investigated the correlation between the activities and fingerprints of PN samples by using multivariate regression techniques including principal component analysis (PCA), partial least squares regression (PLSR), and multi-linear regression analysis (MLRA). Constituents/peaks predicted to be responsible for different activities of PN were then identified, of which the activities were finally verified by pharmacologic tests.

Sample Preparation
Samples were obtained from a single batch of PN root in Yunnan, China. Steamed PN samples were prepared by steaming the crushed raw PN in an autoclave (Shanghai, China) for 2, 4, 6, 8, and 10 h at 105, 110, and 120 • C, respectively. The steamed powder was then dried in a heating-air drying oven at about 45 • C to constant weight, then powdered and sieved through a 40-mesh sieve.

Animals
Kunming mice, male and female, weighing 18-22 g, were purchased from Tianqin Biotechnology Co. Ltd., Changsha, Hunan (SCXK (Xiang) 2014-0011). Before the experiments, the mice were given 1-week acclimation period in a laboratory at room temperature (20-25 • C) and constant humidity (40-70%), and fed with standard rodent chow and tap water freely. Animal experimental procedures in the study were strictly conformed to the Guide for the Care and Use of Laboratory Animals and related ethics regulations of Kunming University of Science and Technology. The protocol was approved by the Experimental Animal Welfare and Ethics Committee, Kunming University of Science and Technology.

HPLC Analyses
The sample solutions were prepared according to the method described in Chinese Pharmacopoeia Commission (2015). HPLC analyses were done on an Agilent 1260 series system (Agilent Technologies, Santa Clara, CA, USA) consisting of a G1311B pump, a G4212B DAD detector and a G1329B autosampler. A Vision HT C 18 column (250 × 4.6 mm, 5 µm) was adopted for the analyses. The mobile phase consisted of A (ultra pure water) and B (acetonitrile). The gradient mode was as follows: 0-20 min, 80% A; 20-45 min, 54% A; 45-55 min, 45% A; 55-60 min, 45% A; 60-65 min, 100% B; 65-70 min, 80% A; 70-90 min, 80% A. The flow rate was set at 1.0 mL/min. The detection wavelength was set at 203 nm. The column temperature was set at 30 • C and sample volume was set at 10 µL.

Anticoagulation Test in Vitro
Blood was collected from healthy mice and directly transferred into citrated tubes (0.109 mol citrate, 9:1). The supernatant platelet-poor plasma (PPP) was obtained by centrifuging the blood samples above at 3,000 rpm for 10-15 min. The mixture of PPP and thrombokinase of various concentrations at the proportion of 2:1 (v/v) of total 50 µL was added into the test cup and incubated for 3 min at 37 • C in a Blood Coagulation Instrument (XN06 series, Diagnostic Technology Ltd of Wuhan Jingchuan, China). 10 U/mL thrombokinase of 100 µL dissolved in 0.1 mol/L Tris-HCl buffer solution (pH 7.4) was subsequently added and incubated at the same condition. The prothrombin time (PT) was determined in accordance with the manufacturer's recommended protocols. The prolongation rate of PT was calculated according to the following equation: where PT 0 was the PT of control (blank, the normal saline replaced of thrombokinase), PT was the PT in the presence of thrombokinase.
The standard curve was drawn with the concentration of thrombokinase (U i ) as the X axis and the lg [prolongation rate of PT (%)] as the Y axis. PN samples of 5 g, in the powdered form, were extracted with pure water (50.0 mL) by refluxing twice for 2 h at 80 • C. The combined solution was filtered and concentrated under reduced pressure to the extract containing 0.1 g/mL of PN. The extract was then diluted with the normal saline to different concentrations. The prothrombin time of the mixed plasma sample containing PPP and PN extract (PT') of different concentrations was determined. The prolongation rate of PT' was calculated according to the following equation: Prolongation rate of PT'(%) = (PT' − PT 0 ')/PT 0 '×100% (2) where PT 0 ' was the prothrombin time of control (blank, the normal saline replaced of extracts), PT' was the prothrombin time in the presence of extracts.
The corresponding concentration of thrombokinase (U i ) was determined according to the standard curve. And the thromboplastin inhibition rate (%) was calculated according to the following equation: Thromboplastin inhibition rate (%) = (U i −10)/10×100% (3) where U i was the concentration of thromboplastin determined by the standard curve.

Antioxidation Test in Vitro
The extracts prepared in "Anticoagulation Test in vitro" were diluted with normal saline to 0.5, 1, 1.5, 2, 2.5, 3, and 3.5 mg/mL, respectively. The scavenging activity for hydroxyl radicals was measured according to the procedure described by Zhao et al. (2006). Reaction mixture contained 60 µL of 1.0 mmol FeCl 2 , 90 µL of 1 mmol 1,10-phenanthroline, 2.4 mL of 0.2 mol phosphate buffer (pH7.8), 150 µL of 0.17 mol H 2 O 2 , and 1.5 mL of extracts prepared. The reaction was started by adding H 2 O 2 . After incubation at room temperature for 5 min, the absorbance of the mixture at 560 nm was measured with a spectrophotometer. The scavenging activity for hydroxyl radicals was calculated according to the following equations: Where A 0 was the absorbance of the control, A 1 was the absorbance in the presence of the extract and A 2 was the absorbance without 1,10-phenanthroline.

Pharmacologic Verification
Constituents/peaks predicted to be responsible for different activities of PN were then identified by reference standards, of which the activities were finally verified by pharmacologic evaluation using the methods described in section "Anticoagulation Test in vitro" and "Antioxidation Test in vitro."

Statistical Analyses
All data were expressed as means ± SD. SPSS 21.0 software (Statistical Program for Social Sciences, SPSS inc, Chicago) was applied to carry out the two-tailed unpaired t-test and PCA. DPS 9.50 software (Data Processing System, China) was used for MLRA. Umetrics SIMCA-P 11.5 software (Sartorius Stedim Biotech, Sweden) was applied for PLS analysis. A value of P < 0.05 was considered to be significant difference. A value of P < 0.01 was considered to be highly significant difference. EC 50 value was fitted by Probit regression with Origin 7.5 for windows (OriginLab Corporation, USA) software.

THEORY PCA
PCA is applied for data compression and visualization. PCA produces so-called latent variables, here called principal components (PCs), which are linear combinations of the original manifest variables. The orthogonal PCs are constructed in such way that they maximize the description of the data variance in the n × p data matrix X. The projections of the objects (PN samples) on PC i are the scores on PC i , and the projections of the variables (HPLC fingerprint signals here) on PC i are the loadings on PC i . Thus, score plots give information related to the (dis)similarity of the objects, while on loading plots information about the contribution of the original variables to a given PC i can be found (Nguyen Hoai et al., 2009).
MLRA MLRA attempts to model the relationship between two or more variables and a response by fitting a linear equation to observed data (Noori et al., 2010;Placca et al., 2010). The general purpose of MLRA is to learn about the relationship between several independent variables and a dependent variable. MLRA can be generally represented in the following form: where Y is the estimated value and represents the dependent variable. X 1 , X 2 , X 3 ,..., X n are measures of not correlated variables that may help in estimating Y. For example, X 1 is the known score of the first independent variable, X 2 is the known score of the second independent variable, etc. The coefficient b 0 is the estimated constant, and b 1 , b 2 , b 3 ..., b n are called the regression coefficients (Hair et al., 1999).

PLSR
PLSR is used to find the inner relationship between independent variables (X) and dependent variables (Y), which are simultaneously modeled by taking into account not only X variance, but the covariance between X and Y (Martens and Naes, 1989). In our study, the X matrix is composed of the enhanced fingerprints and the Y vector is constructed with the reference values of anticoagulation and antioxidation activities (EC 50 ) obtained by the thromboplastin inhibition rate and hydroxyl radicals assay, respectively. Then, X and Y are decomposed in a product of another two matrices of scores and loadings; as described by the following equations: where TP T approximates to the chromatographic data and UQ T to the true Y values; notice that the relationship between T and U scores is a summary of the relationship between X and Y. The terms E and F from the equations are error matrices. Hence, the PLS algorithm attempts to find factors (called Latent Variables) that maximize the amount of variation explained in X that is relevant for predicting Y; i.e., capture variance and achieve correlation (Brereton, 2007).

HPLC Fingerprints
The results of methodology validation showed that the relative standard deviation values for precision, reproducibility and storage stability were less than 3.0, 4.0, and 3.0%, respectively. All the results indicated that the method of HPLC for the fingerprint analyses was valid and satisfactory. The optimized conditions for the 90-min HPLC fingerprints were described in sections "HPLC Analyses." The chromatograms were generated for all batches of PN samples (Figure 1A), and for a typical raw PN sample and a steamed PN sample ( Figure 1B). Peaks with good segregation, which also occupied large areas from consecutive peaks, were determined as the common peaks of PN samples. Therefore, fifteen peaks were selected by comparing their ultraviolet spectra and HPLC retention time. The method used to identify common peaks refers to reports in similar researches (Zheng et al., 2014;Shi et al., 2016). Along with the duration of steaming time and rise of temperature, the area and height of major peaks (peaks 1-3, and 6-8) in the raw PN were decreased gradually, while other peaks (peaks 5, 9-15) were increased or formed ( Figure 1B). The areas of 15 peaks in 18 batches of PN samples were listed in Table 1. The peak area was defined as 0 for  peaks lacked in chromatograms. The coefficients of variance for almost all common peaks were higher than 46.6%. This is due to the diversity in the levels of constituents contained in samples under different process conditions. The areas of 15 common peaks were used for the following analysis.

Clustering Results by PCA
PCA is a classical technique to reduce the dimensionality of the data set by transforming to a new set of variables, named PCs to summarize the features of the data. Since PCs are uncorrelated and ordered, the first few PCs, which contain most of the variations in the data, are usually used in cluster analysis (Yeung and Ruzzo, 2001). As shown in Figure 2, the fingerprints of PN samples were separated into four clusters according to the peak area. Samples 1-3 in cluster one were raw PN. Samples 4-8 in cluster two were PN steamed at 105 • C. Samples 9-13 in cluster three were PN steamed at 110 • C. Samples 14-18 in cluster four were PN steamed at 120 • C. It indicated that PN samples steamed at the same temperature had similar chemical fingerprints. And PCA could initially separated raw and steamed PN samples at different temperatures from the chemical level.
The results also suggested that compared with the steaming time, the steaming temperature had a more important role in the change of chemical composition of PN. Based on the results, a combination of PCA and HPLC methods could, at least roughly, discriminate raw and PN samples under different process conditions.

Anticoagulation Test
PT is used to evaluate the overall efficiency of extrinsic clotting pathway. A prolonged PT indicates a deficiency in coagulation factors V, VII, and X (Chan et al., 2007). In the study, the EC 50 determined by the logarithm to base 10 of PT prolongation rate was applied to evaluate the anticoagulant effect of PN. The standard curve between the concentration of thrombokinase and logarithm of PT prolongation rate showed a good linearity (R = −0.9991). As shown in Figure 3A, raw PN samples (S1-S3) exhibited lower EC 50 values of anticoagulation, suggesting that the anticoagulant effect of raw PN was stronger compared with steamed PN samples. The EC 50 values of PN steamed at the same temperature were generally increased along with the increase of steaming time. For PN steamed for the same time, the higher the steaming temperature, the higher was the EC 50 values of samples. Among these samples, S18 steamed for the longest time of 10 h at the highest temperature of 120 • C showed the highest EC 50 value, suggesting that its anticoagulant activity was the weakest compared with other ones.

Antioxidation Test
PN roots are mostly consumed as popular food tonic in the soup form by people in the southern region of China. Various studies have suggested that the tonifying functions of Chinese herbal medicines could be due to, at least partially, the protective effects against oxidation (Yim and Ko, 2002). Hydroxyl radical is very reactive which can be generated in biological cells through the Fenton reaction. Meanwhile, hydroxyl radical scavenging assay is commonly used for the determination of antioxidant activities of plant extracts. And PN showed higher sensitivity of scavenging hydroxyl radicals than other ones like 1,1-diphenyl-2-picrylhydrazyl free radicals in our previous work ( Figure S1). Therefore, the method was applied to investigate the antioxidant effect of raw and steamed PN roots, with results shown in Figure 3B, where EC 50 was the concentration of PN scavenging 50% hydroxyl radicals. According to the results, raw PN samples (S1-S3) exhibited much higher EC 50 values of antioxidation, suggesting that the antioxidant effect of raw PN was much weaker compared with steamed PN samples. PN steamed at 120 • C showed general lower EC 50 values than samples steamed at lower temperatures, suggesting that higher steaming temperature was related to stronger antioxidation activity of PN.

MLRA
The relationship between the fifteen independent variables X 1 , X 2 , X 3 ,..., X 15 (the values of normalized peak areas) and the dependent variable Y of each activity was established by fitting a linear equation to observed data with multiple linear regression model. The regression equations of anticoagulation and antioxidation were shown as follows:$$ where Y anticoagulation was the EC 50 of lg [prolongation rate of PT (%)], Y antioxidation was the EC 50 of hydroxyl radicals scavenging activity, X 1 -X 15 were the normalized peak areas of peaks 1-15 (Figure 1), respectively. The F-values for the two equations were 8.82 and 13.16, respectively. And the corresponding Pvalues were <0.01 (R 2 = 0.9996), and <0.05 (R 2 = 0.9999), respectively, showing that the established MLRA models were satisfied. According to the equations, EC 50 values of a new PN sample could be obtained by inputting the corresponding peak areas to preliminarily evaluate the anticoagulant and antioxidant activities of this sample. For the anticoagulant activity, values of Pr > |T| of X 2 , X 3 , X 5 , X 6 , X 7 , X 10 , and X 11 were all <0.05, suggesting that constituents corresponding to peaks 2, 3, 5, 6, 7, 10, and 11 had more important role in the anticoagulation of PN. From Table 1 and Figure 1B, peaks 2, 3, 5, 6, and 7 were observed in the chromatographic fingerprints of raw PN. Among them, only peak 5 showed an increase trend along with the increase of steaming temperature and time. For steamed PN, peaks 10 and 11 were exclusively existed in the fingerprints. Thus, constituents corresponding to peaks 2, 3, 6, and 7 might play the major role in the anticoagulation of raw PN, whereas constituents corresponding to peaks 5, 10, and 11 could be the major active ones for the anticoagulation of steamed PN.
For the antioxidant activity, values of Pr > |T| of X 2 , X 3 , X 5 , X 7 , X 10 , and X 13 were all <0.05, indicating that constituents corresponding to peaks 2, 3, 5, 7, 10, and 13 had more significant influence on the scavenging activity of hydroxyl radicals. From Table 1 and Figure 1B, peaks 2, 3, 5, and 7 were observed in the chromatographic fingerprints of raw PN. And peaks 10 and 13 are exclusively existed in the fingerprints of steamed PN. The area of peak 5 was increased along with the increase of steaming temperature and time. Therefore, constituents corresponding to peaks 2, 3, and 7 had the major role in the antioxidation of raw PN, whereas constituents corresponding to peaks 5, 10, and 13 were the major active ones for the antioxidation of steamed PN. The variations in the contents and contribution degrees of above constituents to the activities of PN may lead to the difference in the anticoagulant and antioxidant effects of raw and steamed PN samples.

PLSR
The PLSR models to correlate chromatographic data and the activities of PN were constructed with the 18 batches of PN samples. Since the total number of samples (18) was small and since the prediction for new samples was not our first concern, no division was made into a calibration set to build a PLSR model and a test set to validate the predictive properties. Our main concern was to focus on the indication of anticoagulant and antioxidant peaks from the modeling results. PLSR models were built from the normalized data matrix X containing the 18 PN fingerprints and the response matrix Y, i.e., either the EC 50 of lg [prolongation rate of PT (%)] or the EC 50 of hydroxyl radicals scavenging activity. For the anticoagulation model, two principle components were achieved, accounting for an explained variance of 89.9% for X variable, 84.3% for Y variable, and a predictive ability (Q 2 ) of 85.3% (Table S1), indicating the obtained model was excellent. As shown in the regression coefficients plot (Figure 4A), peaks 4, 5, and 9-15 were positively correlated with the EC 50 of lg [prolongation rate of PT (%)], whereas peaks 1-3 and 6-8 were negatively correlated with the EC 50 value. It should be noted that the predicted EC 50 values could not be defined if these variables increased or decreased, because a negative coefficient did not necessarily mean that the relevant variable has the opposite effect on the anticoagulant activity. Besides, the importance of the X-variables for the model could be summarized by variable importance for the projection (VIP) values (usually with a threshold >1.0). Thus, constituents corresponding to peaks 1, 2, 5, 10, and 11, of which the VIP values were >1.0 (Table S2) with high absolute values of coefficients were considered to be highly related to the anticoagulant activity of PN samples. For the antioxidant model, two principle components were achieved, accounting for an explained variance of 89.9% for X variable, 83.1% for Y variable, and a predictive ability (Q 2 ) of 65.8% (Table S1), indicating the obtained model was excellent. As shown in the regression coefficients plot (Figure 4B), peaks 1, 2, 4, 6, 8, and 11 were positively correlated with the EC 50 of lg [prolongation rate of PT (%)], whereas peaks 3, 5, 7, 9, 10, and 12-15 were negatively correlated with the EC 50 value. Besides, the VIP value of each peak was shown in Table S2. Therefore, constituents corresponding to peaks 3, 5, 10, and 13, of which the VIP values were >1.0 with high absolute values of coefficients were considered to be highly related to the antioxidant activity of PN samples.
Next, the contents of those constituents in raw and steamed PN were determined, as shown in Figure 6. The contents of all the constituents were significantly different between the raw PN and PN samples steamed at 120 • C, and could be markers for the QC of raw and steamed PN.

Pharmacological Verification
In order to verify the predicted active constituents and determine their contributions to each activity, the anticoagulation and antioxidation of the nine identified constituents were tested. As shown in Figure 7A, the sequence of the anticoagulant activity of the constituents was ginsenoside Rd > ginsenoside Rg 1 > ginsenoside Re > ginsenoside Rb 1 > ginsenoside Rh 1 > ginsenoside Rh 4 > notoginsenoside R 1 > ginsenoside Rk 3 > 20 (R)-Rg 3 . Meanwhile, Figure 7B showed that the antioxidant activity in descending order was ginsenoside Rg 1 > ginsenoside Rd > ginsenoside Rk 3 > ginsenoside Rh 1 > ginsenoside Re > ginsenoside 20 (R)-Rg 3 > notoginsenoside R 1 > ginsenoside Rh 4 > ginsenoside Rb 1 . According to the results, ginsenosides Rd, Rg 1 , Re , and Rb 1 , and notoginsenoside R 1 , with stronger anticoagulant activities than other constituents, and higher levels in raw PN than steamed ones, were the major active constituents for the anticoagulation of raw PN, which was consistent with the predicted result of chemometrics analyses. Among the five constituents, Rg 1 , Rd, and Re also showed certain degrees of antioxidant activity, which should be the major antioxidant constituents of raw PN. Conversely, the levels of the three constituents were decreased in PN along with the increase of steaming temperature and duration of time. And other constituents of ginsenosides Rk 3 , Rh 1 , 20 (R)-Rg 3 , and Rh 4 with higher levels or exclusively existed in steamed PN should have more important role in the activities of steamed PN.

DISCUSSION
PN is widely used as a herbal medicine or food tonic in the global market. In order to control the quality of PN, several constituents including notoginsenoside R 1 , and some or all of ginsenosides Rg 1 , Re, Rb 1 , and Rd are determined as markers in the quality standards of different countries or regions (British Pharmacopoeia Commission, 2014;Chinese Pharmacopoeia Commission, 2015; European Pharmacopoeia 8.0; U.S.A Herbal Medicines Compendium 1.0). PN in raw and steamed forms are considered to be different in the medicinal qualities by practitioners of Oriental medicine: the raw materials eliminate and the steamed ones tonify. The so-called "eliminate" means raw PN can move stagnant blood, promote blood circulation, stopping bleeding, and resolving swelling. And the "tonify" means steamed PN can tonify the blood, enhance the immunity and anti-aging (Ge et al., 2015;Gu et al., 2015). Besides the differentiated use of raw and steamed PN by traditional medicine practitioners, the differences in the chemical composition and pharmacologic effects between raw PN and steamed PN have also been verified by modern researches (Lau et al., 2009;Wang et al., 2012). In our research, along with the duration of steaming, some peaks in the chromatograms of raw PN were decreased, whereas other ones were increased or formed. That transformation might be due to the hydrolyzation or dehydration of constituents induced by high temperature. Besides, raw PN was found to be much better in the anticoagulation than the antioxidation compared with the steamed PN, suggesting that raw PN was more suitable to treat coagulation disorders. With the increase of steaming time and temperature, the anticoagulant activity of PN weakened and the antioxidant effect strengthened, which was consistent with the traditional description of medicinal properties of raw and steamed PN. The difference in the pharmacologic effects between raw and steamed PN could be attributed to the change in the chemical composition of PN during the steaming process. However, such differences have still not been acknowledged by national statutes, in terms of the fact that QC markers for raw and steamed PN were undifferentiated at present. This may be due to that the relationships between efficacies and specific constituents of raw and steamed PN remain ambiguous.
The typical approach to select bioactive chemical markers in herbal medicines usually involves two sequential procedures: phytochemical isolation and purification to obtain pure compounds, in vitro and/or in vivo bioactivity evaluation of the individual isolates. There are some shortfalls of this reductionist methodology, such as ignorance of the possible synergism of multiple constituents, missing of the minor components, and too much time-and labor-consumption (Jiang et al., 2010). To overcome these shortcomings, the analysis of fingerprinteffect relationship has been applied to screen characteristic constituents with bioactivity related to the efficacy of a herbal FIGURE 6 | The contents of notoginsenoside R 1 , ginsenosides Rg 1 , Re, Rh 1 , Rb 1 , Rd, Rk 3 , Rh 4 , and 20 (R)-Rg 3 in raw and steamed Panax notoginseng (PN) roots at different steaming temperatures. The points on the scatter plots represent the quantities of the nine constituents in 18 batches of PN samples. Error bars represent mean ± s.e.m. From a two-tailed unpaired Student's t-test, *P < 0.05, **P < 0.01. Frontiers in Pharmacology | www.frontiersin.org medicine. Multivariate regression techniques such as PLSR, PCR, and MLRA are often used to correlate the bioactivity of a herbal medicine and its fingerprint. Regression is a generic term for all methods intended to adjust a model to the observed data, with the purpose of quantifying the relationship between two groups of variables. The adjusted model can then be used either to describe the relationship between the two variables or to predict new variables. Among these techniques, MLRA and PLSR are frequently used to specify a linear relationship between a set of dependent variables from a large set of independent variables, especially when the sample size is small relative to the dimension of these variables (Garza-Juárez et al., 2011;Wu et al., 2015). In this study, the two chemometrics modeling methods, MLRA and PLSR, were preliminarily applied to predict bioactive constituents of raw and steamed PN. Combined with the verification of pharmacologic tests, constituents responsible to different bioeffects of PN were rapidly uncovered. According to the results, nine constituents differently distributed in raw and steamed PN were predicted to be active ones related to different activities of raw and steamed PN. Notoginsenoside R 1 , ginsenosides Rg 1 , Re, Rb 1 , and Rd were predicted by MLRA and PLSR to be the major constituents related to the anticoagulation of raw PN, which was consistent with the measured EC 50 values of anticoagulation. Among them, ginsenosides Rg 1 , Re, and Rd also showed a certain degree of antioxidation, giving evidence for determining these constituents as QC markers of raw PN. Whereas, for steamed PN samples, we found that ginsenosides Rh 1 , Rk 3 , 20 (R)-Rg 3 , and Rh 4 with higher levels or exclusively existed in them could be the major constituents contributing to the activities of steamed PN. Conversely, notoginsenoside R 1 , ginsenosides Rg 1 , Re, Rb 1 , and Rd, as QC markers of raw PN, showed little or no contents in PN steamed at higher temperatures (110 and 120 • C) or for a longer time, which were samples with strong antioxidation (Table 1, Figure 3). Therefore, active constituents as markers for the QC of steamed PN should be different from raw PN, because QC markers of a herbal medicine should be correlated with its safety and efficacy (Capasso et al., 2000).
Ensuring the safety and efficacy of drugs involves multiple considerations and the quality of the drug must be fundamentally guaranteed. An important part of drug QC is to ensure consistent medical and biological effects are delivered by the same drug dosage (Busse, 2000). A differentiated QC standard of steamed PN from raw PN is necessary to ensure the accuracy and safety in clinic use. Meanwhile, the anticoagulation and antioxidation effects are just parts of the major medicinal properties of PN. For better uncovering bioactive constituents of raw and steamed PN, the relationships between chemical information and other efficacies such as anti-inflammation, hemostasis, blood-tonifying, and immunoregulation, need to be further studied.

CONCLUSIONS
In the research, there were divergences in the chemical composition and bioactivities between raw and steamed PN based on the fingerprints and pharmacologic results. Notoginsenoside R 1 , ginsenosides Rg 1 , Re, Rb 1 , and Rd with higher levels in raw PN were verified to be its active constituents, whereas ginsenosides Rh 1 , Rk 3 , 20 (R)-Rg 3 , and Rh 4 with higher levels or exclusively existed in steamed PN were found to be its active constituents. Ginsenosides Rh 1 , Rk 3 , 20 (R)-Rg 3 , and Rh 4 could be used in the future as new markers for the QC of steamed PN. Future research is needed to uncover bioactive constituents related to other efficacies of raw and steamed PN.

AUTHOR CONTRIBUTIONS
YX did the writing of paper and a part of statistical analysis; LC did the anticoagulation test and pharmacological verification as well as the statistical analysis; YH did the HPLC experiments and antioxidation test; YX and XC supervised the project. All the authors read and approve the final manuscript.