Abstract

Dendrobium denneaum paxt., which has been widely used for health prevention in traditional Chinese medicine (TCM), is one of the most popular tonic herbs in China. In order to analyze its flavonoids, characterization and antitumor activity of crude extract and flavonoids rich fractions from D. denneaum paxt. were investigated. Flavonoids extracted from D. denneaum paxt. were clearly enriched in fraction II after determining the total flavonoids content; there were 15 characteristic peaks which have been detected; ultra-high performance liquid chromatography-electrospray ionization/mass spectrometry (UHPLC-ESI-MS/MS) was applied for structural elucidation of compounds. 13 characteristic peaks including flavonoid-O-glycosides and flavonoid-C-glycosides were determined or preliminarily characterized through comparing retention times and UV and MS spectra with standard compounds or documented literature. The antitumor activity of fraction II on human liver cancer cells HepG2 was investigated. MTT assay method was used to test the antiproliferation activity and to confirm the appropriate treatment concentration as well as inducing time. The morphological changes of the apoptosis cells after being induced by fraction II were observed by a Hoechst reagent and the apoptosis rate was tested by flow cytometry. The results showed that fraction II can inhibit HepG2 cells from proliferation in a dose-dependent and time-dependent manner. The apoptosis experiments indicated that fraction II can significantly induce apoptosis in HepG2 cells in a concentration over 50 μg/mL for 48 h and the most effective level was 150 μg/mL for 48 h.

1. Introduction

Dendrobium genus are used as traditional herbal remedies called shihu in China for a long period, which has been deemed as a valuable folk medicine with diversity medicinal activities in traditional Chinese medicine (TCM) theory, including promoting stomach fluid, nourishing yin, and clearing heat [14]. Modern pharmacological studies manifest that Dendrobium genus have various functions on human body such as antitumor, immune-regulation, antioxidant, antiplatelet aggregation, dilation of blood vessels, and hypoglycemic activity [58]. In the Dendrobium genus, D. officinale and Dendrobium huoshanense are the most popular species in traditional Chinese market. Since the 1950s, these two species were on the edge of extinction due to the tremendous demand in health preservation of the public. And since then, the herbal medicinal markets were filled with other Dendrobium species such as D. denneaum paxt. which accounts for the vast majority named “Zipi Shihu.” Traditionally the stems of D. denneaum paxt. will be heated in adequate temperature to be softened and then twisted into a spiral like a spring and dried before selling in the market [9]. After these processes the products are commonly known as Zipi Fengdou. As time goes by, D. denneaum paxt. has become the major source of Dendrobium species in traditional Chinese medicine market which occupies half of the market share. The traditional herbal medicinal practitioners have confirmed that it has similar therapeutic effect with D. officinale. In recent decades, the plantation area of D. denneaum paxt. is expanding rapidly in Longling County of Yunnan Province, China [3].

The existence and activity of free radicals are the main cause of metabolic disorder which may lead to many diseases including angiocardiopathy, diabetes, and cataracts. It is feasible to introduce antioxidant into the body for the prevention or treatment of these diseases. Flavonoids are a widespread group of phytochemistry with diverse biological functions and the significant substances in plants, which not only play a key role in the pharmaceutical industry, but also serve as the excellent chemical markers to control the quality of medicinal plants [10, 11]. Plenty of studies showed that the flavonoid glycosides have obvious pharmaceutical function of antioxidant, hypoglycemic and cardiovascular prevention activity [1219]. Many researchers have studied the mechanism of the antioxidant and anticancer effect of flavonoids extract in plants. Thong et al. studied the mechanism of antioxidant activities of chalcone, flavones, and flavanone and found that the mechanism of antioxidant activities was activated through hydrogen atom transfer and single electron transfer mechanism [20]. Jiang et al. found that the lyophilized aqueous extracts which contained a high level of flavonoids exhibit significant antioxidant activities on mesenchymal stem cells. This kind of activities was through direct ROS-scavenging and indirect Fe2+- chelating [21]. Yuan et al. found that Isoorientin could enhance the expression level of GSH and increase phase II detoxifying enzyme activities in liver cancer cell line which may be the mechanism for Isoorientin anticancer activities [22]. Comparing to the chemical antioxidants, the natural products are healthier with fewer side effects. At present, there have been no reports on the pharmacological studies of the flavonoids from D. denneaum paxt. We speculated that the flavonoids in D. denneaum paxt. would exhibit similar activities as the flavonoids from other plant extract. Therefore, we decided to study the antitumor activity of flavonoids in D. denneaum paxt. by inducing HepG2 cell apoptosis.

There is a substantial content of flavonoid glycosides in the plant D. denneaum paxt.; it is significant for searching and screening out some strong antioxidant activities components to develop the new drug. No reports can be found in fingerprint chromatography establishment by UHPLC-MS as well as antitumor activity study of the flavonoid in D. denneaum paxt. so far. Therefore, it is instructive to build a method of fingerprint chromatography by UHPLC-MS and study the in vitro antitumor activity to evaluate the plant D. denneaum paxt. This study will become an experimental reference in developing a new health preserving and functional drug.

2. Materials and Methods

2.1. Chemicals and Reagents
2.1.1. Chemicals and Reagents in UHPLC-MS Analysis

Apigenin-6,8-di-C-β-D-glucoside and Apigenin-6-C-β-D-xyloside-8-C-β-D-glucoside were isolated from the leaves of D. officinal by preparative liquid chromatography as reference substances for the experiments. Rutin and Schaftoside were obtained from the National Institute of the Control of Pharmaceutical and Biological Products (Guangzhou, China). The purity of all compounds mentioned above is over 96% in HPLC level and their chemical structures were identified by comparing their UV, IR, ESI/MS, and NMR data with other published literatures.

Chromatographic grade methanol was purchased from Merck (Darmstadt, Germany). Preparation of pure water was done by using A Milli-Q water purification system (MA, USA). Methanol and ammonium acetate (CH3CO2NH4) with analytical reagent grade were obtained from Damao Chemical Corporation, Tianjin, China. All standard analytical solutions were stored at 4°C.

2.1.2. Materials and Chemicals in Cell Experiments

The human liver cancer HepG2 cell line was obtained from Guangzhou Jini Ou Biological Technology Co., Ltd. (Guangzhou, China). High glucose Dulbecco’s modified Eagle’s medium (H-DMEM), fetal bovine serum (FBS), penicillin-streptomycin, trypsin, trypsin (non-EDTA), and phosphate-buffered saline (PBS) were purchased from Gibco Chemical Co. (Rockville, MD, USA). MTT, 5-Fluorouracil, and dimethyl sulfoxide (DMSO) were acquired from Sigma-Aldrich (St. Louis, MO, USA). The Annexin V-FITC apoptosis assay kits and Hoechst 33258 assay kits were purchased by KeyGen Biotech Co., Ltd. (Jiangsu, China). The fraction II was a fraction from stem of D. denneaum Paxt. enriched with flavonoids.

2.2. Materials Sample Preparation

50 g of the stem of D. denneaum Paxt. was extracted by refluxing with 800 mL 70% ethanol for 1 h twice. The extraction was combined and concentrated under vacuum to 20 mL. Petroleum Ether was used to remove pigment and other impurities. Then the extraction was dried with water bath and dissolved into suspension with 30 mL of pure water. The BuOH-soluble fraction was chromatographed on AB-8 macroporous adsorption resin eluting with ethanol in different concentrations of 10%, 30%, 50%, and 90%, respectively. The four elution parts named fractions I, II, III, and IV were collected and weighed after being dried. The weights of fractions (I–IV) were 599.23 mg, 345.38 mg, 155.13 mg, and 153.64 mg. The four fractions were dissolved in HPLC grade methanol to test flavonoids.

2.3. UHPLC Instrumentation and Chromatographic Conditions

The UHPLC system consisting of a vacuum degasser, quaternary pump, autosampler, and ultraviolet detector (Thermo Separation Products Inc., Riviera Beach FL, USA) was used for acquiring chromatograms. The studies were carried out on a Hypersil GOLD C18 (100 × 2.1 mm ID, 1.9 μml, Thermo, USA). The mobile phase consisted of methanol (A) and 10 mM (v/v) ammonium acetate aqueous solution (B). Gradient elution at a flow rate of 200 μL·min−1 was applied to establish fingerprint chromatography. The suitable methanol gradient was performed at the following procedure: 0–5 min, 16–20% B; 5–13 B, 20–23% B; 13–19 B, 23-24% B; 19–25 min, 24–26% B; 26–31 min, 26–28% B; 31–45 min, 28–35% B; 45–55 min, 45–55% B. The column temperature was 30°C in the whole analysis process. Sample introducing amounts were 3 uL.

2.4. MS Instrumentation Conditions

The MS analysis was performed on LCQ ion trap instrument (Thermo Finnigan, San Jose, CA, USA) with a heated electrospray ionization (ESI) source coupled to the LC apparatus for mass spectrometric identification of characteristic peaks. The same conditions as described above were used. For optimum MS results, the negative ion mode for MS analyses was selected and worked under the following conditions: spray voltage 3.0 kV, capillary temperature 300°C, sheath gas flow rate at 35 (arbitrary units), auxiliary gas flow rate at 10 (arbitrary units), and split ratio 4 : 1 (v/v), respectively. The full-scan mass spectra were set in the range of 100–1,000. All data were manipulated by Finnigan Xcalibur 2.0 Advanced chromatography workstation (Thermo Quest Corporation, San Jose, CA, USA).

2.5. Total Flavonoids Content (TFC)

The TFC was measured following a described colorimetric method [23]. Rutin was used as standard and expressed as milligrams of Rutin equivalent per gram of sample (mg RE/g sample). The liner regression was   , where was the absorbance and was the concentration of Rutin equivalent (RE) mg/ml.

2.6. Analysis of Antitumor Activity
2.6.1. The Choice of Flavonoids Fraction and Conducting Concentration

In the cell experiment section, we chose the flavonoids fraction which contained most of the flavanoids. The conducting concentrations (10, 25, 50, 100, and 150 μg/mL) and duration time for cell viability were based on a routine practice of the antitumor in vitro experiments. And the conducting concentration and duration time for further experiments included cell morphological change and flow cytometry analysis was based on the cell viability results.

HepG2 cells were cultured at 37°C in a humidified atmosphere of 5% CO2 in H-DMEM medium supplemented with 10% FBS (v/v) and 1% penicillin-streptomycin. The control group were cultured without fraction II in serum-free medium while the positive drug group were cultured in serum-free medium with 5-Fluorouracil (5-Flu). The stock solution of flavonoids fraction II was prepared by dissolving flavonoids fraction II in DMSO with a concentration of 200 mg/mL in order to guarantee the DMSO content was less than 1% in the final fraction II solution. The cell experiments were done by diluting the stock solution with H-DMEM medium (without FBS) to different concentrations (10, 25, 50, 100, and 150 μg/mL).

2.6.2. Cell Viability Assay

The cell viability was estimated using MTT assays. Briefly, HepG2 cells were seeded in a 96-well plate (3000 cells/well) and cultured for 24 h using H-DMEM (with FBS). Then the culture medium was replaced by H-DMEM (without FBS) with fraction II in different concentrations (10, 25, 50, 100, and 150 μg/mL) for 24 and 48 hours, respectively. The control group and positive drug group were under the same conditions. After the treatment, 20 μl MTT (5 mg/mL) was added to each well and sequentially incubated for an additional 4 h at 37°C. The cell supernatant was removed and 100 μL DMSO was added to dissolve the dark blue formazan. The absorbance was measured at 490 nm using a microplate reader (Bio-Rad Laboratories, Inc., Hercules, CA, USA) after the plate was shaken for 5 minutes. The cell viability was calculated as follows:

2.6.3. Cell Morphological Changes

The morphological changes of HepG2 cells after being induced by fraction II were evaluated by staining with the chromatin dye Hoechst 33258. The HepG2 cells were cultured in 6-well plates in the concentration of 1 × 105 cells/mL for 24 h and then exposed to various concentrations of fraction II (0, 50, and 150 μg/mL) and 5-Flu (50 μg/mL) for 48 h under the same condition. Subsequently, the cells were immobilized with 4% paraformaldehyde for 10 min at 4°C and then stained with Hoechst 33258 for 10 min in the dark after being washed twice with PBS. Morphological changes and the images of the apoptotic cells were observed using a fluorescence microscope (Leica, Microsystems GmbH, Wetzlar, Germany).

2.6.4. Flow Cytometry Analysis of Cell Apoptosis

HepG2 was seeded in 6-well plates (1 × 105 cells/mL) and cultured for 24 h and then treated with fraction II in different concentration (0, 50, and 150 μg/mL) for 48 h. The experiments was performed according to the manufacturer’s instructions of kit; HepG2 cells were digested with trypsin (non-EDTA) and collected after being washed by PBS twice. Afterwards, the cells were suspended in 500 μl binding buffer then stained with 5 μl Annexin V-FITC and 5 μl propidium iodide (PI) for 15 min at room temperature in the dark. The samples were immediately detected using BD FACS Canto II flow cytometer.

2.7. Data Analysis

Data were expressed as mean as the mean ± SEM and the results from each group were obtained using triplicate samples independently. One-way analysis of variance (ANOVA) and Fisher’s least significant difference (LSD) test were used to determine significant differences for different samples. A probability value of was considered to be statistically significant while was considered to be highly statistically significant and are indicated in the figures. All statistical analysis were conducted using SPSS 20.0 and Prism 5 software package (GraphPad Software Inc., La Jolla, CA, USA).

3. Results

Take adequate amounts of the four fractions for HPLC test. Most of the flavones were in fraction II while very few of them can be detected in fraction III. No flavones can be spotted in fractions I and IV. And, as a result, the fraction II was applied for further experiments including qualitative analysis of the flavonoids and cell experiments. From Figure 1 we can see that, at the wavelength of 340 nm, 15 chromatographic peaks can be found in fraction II. In order to clarify the structure of these components, UHPLC-ESI-MS/MS was applied.

3.1. Identification of Flavonoid Glycosides

Comparing the RT (retention time) of UV and ESI-MSn spectra with literature data, we speculated and identified the potential chemical structures of 13 main essential ingredients (expect peaks (1) and (2) which were not detected), including flavonoid-O-glycosides and flavonoid-C-glycosides. The RT, MS, and MS2 spectral information and the identification results of the flavonoids were listed in Table 1.

In the present study, 5 flavonoid-c-glycosides have been identified in the Dendrobium aurantiacum var. denneaum including peaks (3), (4), (5), (6), and (8). A series of characteristic fragment patterns of flavonoid-c-glycosides arise from the fracture of the glucosyl as follows: the major fragmentation pathways were parent fragment ions loss of [(M-H)-60], [(M-H)-90] and [(M-H)-120], [(M-H)-90-120], [M-H-2 × 120], [(M-H)-120-CO], [(M-H)-2 × 120-2CO], and so on in MS or MS/MS. The glycosides of flavonoids-C-glycosides in D. denneaum paxt. are usually replaced in the C-6 and C-8 position. Furthermore, the glycosides on the C-6 position revealed more commonly ion fragmentation than C-8 position. [M+H]+ and [M+Na]+ ions were the primary positive ESI/MS in D. denneaum paxt., in which further fragment was [M+H-18]+ resulting from continuous losses of one molecule of H2O. The characteristic fragment indicated that the positive negative ion mode could not provide much help to structure determination. Peak (3) showed [M-H] at 593, the fragment ion peak showed in Figure 2(a), including 473 [(M-H)-120] and 353 [(M-H)-120-90]; it was identified as Apigenin-6,8-di-C-β-D-glucoside by comparing with the standard. Peaks (4) and (5) had a molecular ion [M-H] at 563, which produced similar MS2 major fragments at 473 [(M-H)-90] and MS3 at 353 [(M-H0-90-120], as shown in Figure 2(b); upon comparison with standards, peaks (4) and (5) were established as Apigenin-6-C-β-D-xyloside-8-C-β-D-glucoside and Isoschaftoside. Peaks (6) and (8) (in Figure 2(c)) also showed the same molecular ion [M-H] at 563, yielding produced fragments at 443 [(M-H)-120] and 353 [(M-H)-120-90]. Meanwhile, accompanied with ions at 503 [(M-H)-60], 383 [(M-H)-120-60] and 353 [(M-H)-120-90] have also been found. Moreover, these findings are consistent with the literature data. Compounds 6 and 8 were unambiguously identified as Schaftoside and Apigenin-6,8-di-C-α-L-arabinoside by comparing with the standard.

Peaks (7), (11), (12), (13), (14), and (15) showed similar fragmentation behaviours with the loss of 308, 176, 162, and 146 mass units showing the potential existence of flavonoid-o-glycosides, while 308 Da was proved as the typical fragment of Rutinose. The glycosides ordinarily took place on the position of C-3 and C-7; the relative intensity of [(M-2H)-gly] higher than [(M-H)-gly] indicated that the glycosides superseded on the C-3 position while the C-7 position is completely the opposite. For instance, Rutin, [M-H] ion, at 300 possessed higher intensity than 301. Peaks (7), (12), and (13) were assigned as 3-O-rutinosides, whereas peaks (11), (14), and (15) were assigned as 7-O-Rutinosides. Peak (7) showed [M-H] ion at 771, and further fragmentation at 609 ([M-H]-162) (Figure 2(d)), 300 ([M-H]-162-308), and 301 ([M-2H]-162-308), suggesting the presence of Rutinose and Quercetin in the structure; it is tentatively identified as Quercetin-3-O-rutinoside-7-O-glucoside. The parent ions at 609.34 [M-H] and 610.82 [M+H]+ of peak (11) were similar to that of peak (12). It means that molecular mass of peaks (11) and (12) was 610, which produced fragments at 301 [(M-H)-308], 300 [(M-2H)-308]. The 303 that appeared in positive ion mode of peaks (11) and (12) was the proof that Quercetin existed. Peak (11) was tentatively determined as Quercetin-7-O-rutinoside. Peak (12) (showed in Figure 2(e)) was unambiguously identified as Rutin through comparison with the standard compound. Peaks (13) and (14) shared ions of [M-H] at 593 and exhibited fragments at 285, 284, 163, and 151; 285 and 284 indicated that these compounds contain Kaempferol. Based on the chromatographic column behaviour, peak (13) was assigned as Kaempferol-3-O-rutinoside; compound (14) was identified as Kaempferol-7-O-rutinoside (Figure 2(f)). Peak (15) had a [M-H] ion at 577, which produced second-level fragment at 269.06 (confirmed as Apigenin). Thus, compound (15) was deduced as Apigenin-7-O-rutinoside (Figure 2(g)); it has been reported in reference articles.

Peaks (9) and (10) showed the same [M-H] ion at 563.13, which generated the peaks at 413 [M-H-150], 443 [M-H-120], and 293 [M-H-120-150]; this study led to the identification of C-glucosyl-2′′-O-xyloside. Compound (9) was assigned as Apigenin-6-C-glucosyl-2′′-O-xyloside and compound (10) was identified as Apigenin-8-C-glucosyl-2′′-Oxyloside. Both of them were initially identified through comparing their mass spectrum and UV spectral values, which have been identified in Citrus fruit juices. The results were shown in Figure 2(h).

3.2. Total Flavonoids Content

The total content of flavonoids components enriched fraction from D. denneaum paxt. was 1792.11 ug RE/g.

3.3. Antitumor Activity of Fraction II
3.3.1. Antiproliferative Effect of Fraction II in HepG2 Cell

The antitumor effect of fraction II on human liver cancer cell HepG2 was measured by MTT assay. Different concentrations of fraction II were used to investigate the effects on the cell viability of HepG2. As it is shown in Figure 3, fraction II inhibited the growth of HepG2 cells in a dose-dependent and time-dependent manner after being induced for 24 h and 48 h. When comparing fraction II treatment group with the control group and positive drug group (5-Flu), the appropriate duration time of fraction II treatment was 48 h. After being treated for 48 h, the cell viability will decrease to 50.91%. Interestingly, even though the treatment concentration level reaches over 50 μg/mL, the cell viabilities will not decrease markedly and are only maintained at a level of about 50% in 48 h. According to the statistical analysis, there was no statistical significance between the cell viability of the treatment groups in three different concentration levels at 50, 100, and 150 μg/mL after 48 h. Therefore, we would like to use the concentration levels at 50 μg/mL and 150 μg/mL to proceed with further experiments to find out whether the morphological changes and apoptosis rate of the HepG2 cells would be the same like the MTT assay results.

3.3.2. Anticancer Effect of Fraction II on the Morphological Changes of HepG2 Cell

The apoptotic morphology and morphological changes of HepG2 cells were observed by fluorescence microscopy. As it is shown in Figure 4, after the treatment with fraction (50 and 150 μg/mL) for 48 h, the number of cells decreased and typical characteristics of apoptosis phenomenon such as membrane blebbing, cytoplasm shrinkage, nuclear chromatin condensation, and nucleoli fragmentation of cells can be observed. However, there were no significant differences between the cells treated with 50 and 150 μg/mL which is similar to MTT assay. That is to say fraction II can effectively inhibit the growth of human liver cancer HepG2 cells by inducing apoptosis.

3.3.3. Anticancer Effect of Fraction II on the Apoptosis Rate of HepG2 Cells

In order to explore the apoptosis rate of HepG2, we applied the Annexin V and PI double-staining assay kit to study using flow cytometry. As it is shown in Figure 5(b), after incubation with fraction II at the concentration of 50 and 150 μg/mL for 48 h, the apoptosis rates were 25.4% and 27.0%, respectively. Although there were no significant differences between the cells induced with 50 and 150 μg/mL fraction II, they were both remarkably higher than the control group, and almost reaching the apoptosis rate of the positive drug group, which exhibited a significant increase in the cell apoptosis rate. This result was also similar to MTT assay and demonstrates once again that fraction II has an effectively apoptotic effect on HepG2 cell when the concentration was over 50 μg/mL for 48 h.

4. Discussion

Flavonoid is a group of phenols with low molecular weight which widely exists in plants. Some of them have been proved to exhibit certain antitumor activities, and the mechanism was clarified by researchers [24]. Dendrobium genus is the largest in genus in Orchid and they have been used for thousand years in the mainland of China, as it was recorded that they exhibit certain pharmacology function such as nourishing the stomach and enhancing production of body fluids [25]. In the modern pharmacology study, researchers started to recognize the high value plants from Dendrobium genus that exhibit high efficacy on antitumor activities, such as Dendrobium officinale. Fan and Luo studied the antitumor activity of the polysaccharides from Dendrobium denneaum and found that the purified polysaccharides fractions DDP1-1 exhibit strong antitumor activities [26]. Haizhen et al. studied the antitumor activities of Dendrobium officinale polysaccharides and found that polysaccharides fraction DOPA can inhibit the growth of MCF-7 cell line. And the mechanism may be related to the PI3 K/AKT pathway [27]. However, most studies of the antitumor activities of Dendrobium genus mainly focus on the polysaccharides or crude extract, not the purified flavonoids fractions. So, in this study, we investigated the antitumor activities of the flavonoids from D. denneaum paxt. and found that the flavonoids’ fraction II exhibits strong antitumor activities on human liver cancer cells HepG2. Furthermore, we preliminarily confirmed that the antitumor activities effects were through inducing apoptosis of HepG2 cells and the effective concentration and duration time was at a concentration above 50 μg/mL for 48 h of the flavonoids from D. denneaum paxt.

We applied separation technology to obtain flavonoids fractions from D. denneaum paxt. and tested their antitumor activities. AB-8 macroporous resin was used to enrich and purify the crude extract. After being eluted with ethanol in different concentration, four fractions were collected and most of flavones were in fraction II. The results of UHPLC-MS showed that there were 15 chromatographic peaks in the flavones of the D. denneaum paxt. And the structure of 13 peaks was elucidated. When comparing the characterized peaks to the former studies published previously by Ye et al. [28], some compounds such as two Apigenin glucosides, Isoschaftoside, Schaftoside, and Rutin identified in this study were also identified by the previous paper. In this study, four Apigenin glucosides, two Kaempferol glucosides, and two Quercetin were identified from D. denneaum paxt. for the first time. Apart from the qualitative analysis of flavonoids, the total flavonoids content was also tested. There is a substantial content of flavones constituents in D. denneaum paxt.; the content of the Rutin especially is exceptionally high, which possesses many potential bioactive properties.

5. Conclusion

The study of natural plants products is becoming more and more popular all around the world. As the human society is developing faster and faster, the tendency of population to age would raise more attention to the medical care issue. In this study we proved that the flavonoids in D. denneaum paxt. exhibit certain antitumor and proapoptotic activities. So this plant has the potential to be developed into a pharmaceutical and nutraceutical products.

Abbreviations

D. denneaum paxt.: Dendrobium denneaum paxt.
TCM:Traditional Chinese medicine
UHPLC-ESI-MS/MS:Ultra-high performance liquid chromatography-electrospray ionization/mass spectrometry
RT:Retention times.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Chunhua Zhou and Yingyi Luo contributed equally to this work.

Acknowledgments

This study was funded by the Special Foundation of 2017 High Level University Construction (no. A1-AFD018171Z11024) and key project at central government level-the ability establishment of sustainable use for valuable Chinese medicine resources (no. 2060302).