|Year : 2019 | Volume
| Issue : 63 | Page : 378-384
Research on mechanism of traditional hot processing on velvet antler based on 1H nuclear magnetic resonance and multivariate statistical analysis
Jingfeng Li1, Hui Sun2, Xuefeng Bian1, Nanxi Zhang1, Yaping Wang3, Ying Zong4, Hui Zhang1, Jiaming Sun1
1 Jilin Ginseng Academy, Changchun University of Traditional Chinese Medicine, Changchun, China
2 Eastern District of Clinical Coordination Department, Hangzhou Simo Co., LTD, Hangzhou, P. R. China
3 College of Pharmacy, Changchun University of Chinese Medicine, Changchun, China
4 College of Chinese Medicine Materials, Jilin Agricultural University, Changchun, China
|Date of Submission||23-Nov-2018|
|Date of Decision||02-Jan-2019|
|Date of Web Publication||16-May-2019|
Jilin Ginseng Academy, Changchun University of Traditional Chinese Medicine, Changchun
Jilin Ginseng Academy, Changchun University of Traditional Chinese Medicine, Changchun
Source of Support: None, Conflict of Interest: None
| Abstract|| |
Background: The processing of Chinese materia medica is an important part for the preparation of traditional Chinese medicine (TCM) which has been believed that can change the chemical composition and physiological function. Velvet antler is a kind of precious TCM and widely used after hot processing as folk medicines in China. However, no strategy has been presented to reveal the processing principle of velvet antler. Objective: The aim of this study is to develop a method for exploring the mechanism of hot processing on velvet antler. Materials and Methods: In this study, the chemical compositions of fresh velvet antler (FVA) and hot processed velvet antler (HPVA) were compared based on the 1H nuclear magnetic resonance (1H NMR) and multivariate statistical analysis methods to explore the mechanism of hot processing on velvet antler. Results: It showed that hot processing increased phosphorylcholine, taurine, alanine, uridine, phenylalanine, uracil, tyrosine, accompanied by decreased bile acid, choline, lipid 2 and succinic acid in comparison to the FVA, which changed significantly in 52 common metabolites of FVA and HPVA. In addition, a clear separation between FVA and HPVA was obtained by using principal component analysis and 11 “chemical markers” were found by using orthogonal partial least squares discriminant analysis (OPLS-DA) approaches. Particular attention was given to the “chemical markers” with potentially beneficial bioactivities. Conclusion: The study demonstrated that 1H NMR-based chemometric approach could be a promising tool for investigation of the hot processing of velvet antler in a holistic way of TCM.
Keywords: 1H nuclear magnetic resonance, hot processing, multivariate statistical analysis, traditional Chinese medicine, velvet antler
|How to cite this article:|
Li J, Sun H, Bian X, Zhang N, Wang Y, Zong Y, Zhang H, Sun J. Research on mechanism of traditional hot processing on velvet antler based on 1H nuclear magnetic resonance and multivariate statistical analysis. Phcog Mag 2019;15:378-84
|How to cite this URL:|
Li J, Sun H, Bian X, Zhang N, Wang Y, Zong Y, Zhang H, Sun J. Research on mechanism of traditional hot processing on velvet antler based on 1H nuclear magnetic resonance and multivariate statistical analysis. Phcog Mag [serial online] 2019 [cited 2020 Feb 22];15:378-84. Available from: http://www.phcog.com/text.asp?2019/15/63/378/258396
- The chemical compositions of fresh velvet antler (FVA) and hot processed velvet antler (HPVA) were analyzed by 1H nuclear magnetic resonance (1H NMR)-based metabolomics approach for the first time
- 52 compounds were identified and 11 compounds were considered as potential “chemical markers” concerning FVA and HPVA
- 1H NMR spectroscopy and chemometric tools are recommended as a holistic way of traditional Chinese medicine for studying the hot processing mechanism of velvet antler.
Abbreviations used: TCM: Traditional Chinese medicine; FVA: Fresh velvet antler; HPVA: Hot processed velvet antler; 1H NMR: 1H nuclear magnetic resonance; PCA: Principal component analysis; OPLS-DA: Orthogonal partial least squares discriminant analysis.
| Introduction|| |
Velvet antler, the immature unossified horn of Cervus nippon Temminck or Cervus elaphus Linnaeus, a traditional Chinese medicine (TCM), has been recorded in the 2015 edition of the Chinese Pharmacopoeia, which is believed to invigorate kidney yang, promote blood flow, and strengthen the bones and muscles. Modern pharmacological studies show that velvet antler exhibits wide functions, such as protection of myocardial injury, anti-aging, immune modulation, and facilitation of learning and memory. Velvet antler has been reported to contain amino acids, proteins, peptides, polyamines, vitamins, hormones, nucleic acids, bases, lipids, aromatic compounds, and inorganic elements.
Processing of medicinal materials is a traditional pharmaceutical technique used to fulfil the different requirements of storage, dispensing, and making preparations according to TCM theory. The aims of processing are to increase medical potency and reduce toxicity and side effects of crude drugs., The growth inhibition rate of hot proceed antler hepatoma carcinoma cell was higher than fresh antler. Processed velvet antler is widely used in the clinic in TCM for the treatment of invigorating kidney yang. In this work, the velvet antler metabolites of two differently processed products were compared.
As an animal medicine, the processing method of velvet antler has an important influence on its effective components and curative effect and directly relates to the quality of velvet antler. Conventionally, hot processed velvet antler (HPVA) was made by boiling fresh velvet antler (FVA) at 100°C for a short time and draining the residual blood, which have been repeated several times and then drying at 70°C. The purpose of hot processing is to make velvet antler dry and antiseptic, which is convenient for storage, transportation, and utilization. For analyzing various compositions of the velvet antler well-established techniques currently used include the Folin-phenol method, the sodium dodecyl sulphate–polyacrylamide gel electrophoresis, reversed-phase high-performance liquid chromatography, and liquid chromatography-mass spectrometry.,,, Beyond these traditional methods, which are mostly investigated for the changes of a certain type components between FVA and HPVA, the 1 H nuclear magnetic resonance (1 H NMR) technique can be used to characterize the overall variation of a variety of water-soluble small molecules.
1 H NMR spectroscopy, which has an informative, nondestructive and nontargeted nature, coupled with multivariate analyses such as principal component analysis (PCA) and orthogonal partial least squares discriminant analysis (OPLS-DA) have been applied for the discrimination of different donkey-hide gelatin samples produced by different manufacturers. Donkey-hide gelatin, like antler, is a kind of precious animal TCM. Thus, the metabolite profiling of FVA and HPVA could be investigated with the 1 H NMR method combined with multivariate analysis methods.
In the present study, a comparative analysis of the chemical constituents of the processed and unprocessed velvet antler was performed based on 1 H NMR spectra combined with PCA and OPLS-DA, to clarify the mechanism of the hot processing of velvet antler in accordance with the holism view of TCM. Furthermore, a reference point for a convenient method was provided for investigation of the processing of animal medicine in TCM.
| Materials and Methods|| |
Chemicals and materials
The FVA of Cervus nippon Temminck was purchased from Dongfeng in Jilin Province of China. Chromatographic grade acetonitrile, analytical grade Na2 HPO4·12H2O, and NaH2 PO4·2H2O were purchased from Beijing Chemical Works (Beijing, China). Sodium 3-trimethlysilyl [2, 2, 3, 3-d4] propionate (TSP) was acquired from Cambridge Isotope Laboratories Inc., (Andover, MA, USA).
Preparation of fresh velvet antler and hot processed velvet antler
The fresh sika deer antlers (two-branched) were burned the hair with alcohol lamp and cut into two portions in the longitudinal direction. One portion of FVA was wrapped with plastic wrap and boiled at 100°C for 1 min and then drained the residual blood, which have been repeated 5 times. Subsequently, that portion was dried at 70°C as HPVA sample, accompanied by using the other portion as FVA sample for the experiments. Then ten samples were gotten including five FVA samples and five HPVA samples after five sika deer antlers were proceed according to the above method.
Sample preparation for 1 H nuclear magnetic resonance analysis
Twenty grams of FVA sample or HPVA sample was weighed and 200 mL of deionized water was added for homogenizing, followed by extraction conducted at 4°C for 24 h, the supernatant was collected by centrifugation at 3600 r/min for 10 min and freeze-dried. Forty milligrams of lyophilized powder was weighed into a 10 mL centrifuge tube and added to 5 mL 80% acetonitrile (v/v). The solution was dissolved by ultrasound for 1 min and centrifuged at 4000 r/min for 30 min and then, the solvent was evaporated at 75°C in a water bath. Dried extract was added to 500 μL buffer (pH 7.4 with Na2 HPO4 and NaH2 PO4, 0.1 M) and 100 μL D2O (containing 0.05% TSP [m/v]). The solution (600 μL) was transferred into a 5-mm NMR tube for analysis. All of 10 samples, including five FVA samples and five HPVA samples, were treated with the above method for 1 H NMR analysis.
Spectra acquisition and processing
1 H NMR experiment was recorded at 25°C on a Bruker AVANCE III NMR spectrometer (Bruker, Karlsruhe, Germany) operating at 500 MHz.1 H NMR spectra used standard NOESYPRGP1D pulse sequence. Waiting time (3s) and mixing time (300 m) used lower powered continuous wave pulse for water peak suppression. With deuteroxide lock field, the internal standard is TSP. Each 1 H NMR spectrum consisted of 64 scans requiring 5-min acquisition time with the following parameters: 0.18 Hz/point, pulse width = 5498.53 Hz, at 90° (12.08 μs), sampling time 3.997s, and relaxation delay = 5.0 s.
The 1 H NMR spectrum was processed using Mest Re Nova (version 9.0.0, Mestrelab Research, Santiago de Compostella. Spain). Every spectrum was subsequently segmented at 0.04 ppm intervals across the chemical shift δ 0.00–9.00 ppm and the signal intensity in each region were integrated. The regions of δ 4.70–5.10 were excluded from the analysis because of the residual signal of HDO. SIMCA-P 11.0 (Umetrics, Umeå, Sweden) was used for PCA, projections to latent structures-discriminant analysis (PLS-DA) and orthogonal projections to latent structures-discriminant analysis (OPLS-DA) of the NMR data to find differential metabolites.,
| Results and Discussion|| |
Metabolites identification by 1 H nuclear magnetic resonance analysis
To illustrate the differences in metabolites between FVA and HPVA, 10 samples from five different velvet antlers were purchased and analyzed with the 1 H NMR spectroscopic technique. Identification of the metabolites in the two groups of velvet antler unveiled diversity in metabolite profiles.
Each compound is not only a chemical shift and couple constant but also most of the compounds are several sets of chemical shifts, split pattern, and couple constant that constitute a “three-dimensional data set” for the identification, which is based on the identification results obtained by referring to metabolomics literature.1 H-NMR data and assignments of the metabolites of FVA and HPVA have the same peak profils including chemical shift, split pattern and couple constant of the common peaks, so only the spectrum of FVA is used as a representative to illustrate the identification of all metabolits in FVA and HPVA [Figure 1]. The metabolites were confirmed by comparing with the chemical shifts, split pattern and couple constant of standard compounds from the Biological Magnetic Resonance Data Bank, Human Metabolome Database, as well as the literature [Table 1].,,,,,,,,,,, The main metabolites are amino acids, organic acids, choline, and nucleoside compounds, which are consistent with the literatures.,, The results of the identification are shown in [Table 1], and the comparison of FVA and HPVA is shown in [Figure 2].
|Table 1: 1H nuclear magnetic resonance data of major metabolites in fresh velvet antler and hot processed velvet antler|
Click here to view
|Figure 2: 1H nuclear magnetic resonance comparison chart of fresh velvet antler and hot processed velvet antler|
Click here to view
Multivariate statistical analysis of 1 H-NMR data
The 1 H NMR data were subjected to various multivariate data analyses to determine the chemical differences between FVA and HPVA.
Analysis of PCA is a commonly used unsupervised learning method in multi-dimensional statistical analysis. In search of the principal metabolites, which emphasize the chemical changes in processing, a comparison of FVA and HPVA through PCA was conducted for the metabolomic study. Its scatter plot can directly indicate the difference between FVA and HPVA and it was clear that processed and unprocessed could be divided into two groups, indicating that the variety and/or content of the chemical components in HPVA had significantly changed [Figure 3]. Five HPVA samples were proceed by a series of complex treatments including heating, bleeding, and drying from five different FVA, so there were more differences within the HPVA group.
|Figure 3: Principal component analysis score plots of fresh velvet antler and hot processed velvet antler, fresh velvet antler (), hot processed velvet antler ()|
Click here to view
PLS-DA is a supervised discriminant technique to process large data set and can reduce the differences of inter-group to the lowest level. In general, the fitness of model was indicated by R2 and the predictive ability was represented by Q2. In order to improve the separation between FVA and HPVA groups and interpret the potential discriminant metabolites of the two groups, OPLS-DA is performed and validated by the corresponding PLS-DA models. It showed that the model had excellent predictive abilities, which R2 = 0.637 and Q2 = −0.36 were obtained [Figure 4]a.
|Figure 4: Projections to latent structures-discriminant analysis model (a), score plot (b) and S-plot (c) of Orthogonal projections to latent structures-discriminant analysis results obtained from 1H nuclear magnetic resonance spectra of the different samples: Fresh velvet antler () and hot processed velvet antler ()|
Click here to view
A clear separation in the OPLS-DA score plot between FVA and HPVA groups was acquired, which revealed that they had different metabolic profiling [Figure 4]b. From the S-score plot constructed from OPLS-DA [Figure 4]c, metabolites which had significant differences between FVA and HPVA groups could be found. In this study, 11 significant differences metabolites, as “chemical markers” for discriminating between FVA and HPVA, were selected based on their distance from the origin in the S-plot, variable importance in projection values larger than 1.0 and P < 0.05 [Table 2].
Comparison of metabolite content between fresh velvet antler and hot processed velvet antler
Taking TSP as the internal standard, the 15 components of the identified compounds with obvious characteristic peaks without overlapped peaks with other substances were manually integrated [Table 2]. The differences in the content of 11 metabolites had the significance of statistics (P < 0.01) by using of t-test with SPSS version 16.0 (SPSS Inc., Chicago, IL). Levels of fatty acid, phosphoryl-choline, taurine, alanine, uridine, phenylalanine, uracil, and tyrosine in HPVA were significantly (P < 0.01) higher than in the samples from FVA, on the contrary, the contents of choline, lipid 2, and succinic acid were higher in FVA than in HPVA.
The results were basically consistent with the multivariate statistical analysis, which 11 extremely significant differences components are illustrated in the S-plots [Figure 4]c.
The most significant components to discriminate FVA and HPVA are amino acids, organic acids, and nucleosides. The contents of alanine, phenylalanine, and tyrosine were higher in HPVA, which were hydrophobic amino acids. This may be due to the breaking of peptide bonds of protein and polypeptide contained in FVA by heat treatment, which produced more free hydrophobic amino acids. The decrease of the content of hydrophilic threonine may be due to dissolving in hot water during boiling and discharging blood. In addition, the content of taurine, uridine, and uracil was increased in HPVA. Uridine, uracil, and succinic acid had a certain immune-enhancing activity., Moreover, taurine, high concentration of saturated fatty acid could significantly promote the secretion of testosterone in the rat Leydig cells in vitro,, which may be caused by taurine involved in the biosynthesis pathway of steroid hormone. Kidney-yang deficiency is related to the neuroendocrine-immune system in modern medicine including enhancing immune function and regulating sex hormone content. It is suggested that the change in the content of the above “chemical markers” is one of the reasons for the invigorating Kidney Yang activity of HPVA.
| Conclusion|| |
Hot processing of velvet antler is a traditional method, and it has great significance for the efficacy of antler. This is the first time that a metabolic profiling analysis, using 1 H NMR and multivariate data analysis, was applied to explain the processing mechanism of antler. The results presented here suggested that HPVA differed from FVA greatly on the chemical composition. The information obtained in this work is valuable for identifying fresh or processed products of velvet antler and establishing the analytical requirements for quality control in future. It is of practical significance for the modern research of velvet antler and to ensure the safety and efficacy of clinical medicine. In addition, further research of their respective bioactivity is also important and could be conducted to fully understand the essence of hot processing, which be linked to clinical practice.
The authors would like to thank Prof. Chunyu Wang, State Key Laboratory of Supramolecular Structure and Materials, Jilin University, for his kind help in determination of 1 H NMR.
Financial support and sponsorship
This work was supported by the National Natural Science Foundation of China (grant no. 81373936), the key Scientific and Technological project of Jilin province (grant no. 20160209006YY) and the Science and Technology Develop Project of Jilin Province (grant no. 20140203015YY).
Conflicts of interest
There are no conflicts of interest.
| References|| |
The Medicine Science and Technology Press of China. Chinese Pharmacopoeia Commission. Part. I. Beijing: The Medicine Science and Technology Press of China; 2015. p. 323-4.
Liu D, Gao JT, Sun JM, Zhang H. The pharmacological action research deer nearly five years. Jilin J Tradit Chin Med 2015;9:927-30.
Wu F, Li H, Jin L, Li X, Ma Y, You J, et al.
Deer antler base as a traditional Chinese medicine: A review of its traditional uses, chemistry and pharmacology. J Ethnopharmacol 2013;145:403-15.
Ke LJ, Lin DY, Huang XN, Huo YS, Rao PF, Ye XY. Comparison of protein composition and activities of pilose antler processed by different methods. Zhong Yao Cai 2008;31:11-4.
Zhao Z, Liang Z, Chan K, Lu G, Lee EL, Chen H, et al.
A unique issue in the standardization of Chinese materia medica: Processing. Planta Med 2010;76:1975-86.
Zhao J, Wang M, Avula B, Zhong L, Song Z, Xu Q, et al.
Effect of processing on the traditional Chinese herbal medicine Flos lonicerae
: An NMR-based chemometric approach. Planta Med 2015;81:754-64.
Zhang R, Li Y. Processing technology of pilose antler. Guangming J Chin Med 2012;27:822-3.
Wang Y. Study on the Quality Analysis Method of Pilose Antler Oral Liquid. Heilongjiang University; 2015.
Wang X, Zhao L, Wang CT. Extraction and in vitro
simulated gastrointestinal digestion of velvet antler water-soluble protein. Sci Technol Food Ind 2014;12:81-6.
Zhou J. Construction of two-dimensional liquid chromatography and its application in separation of velvet antler protein. Mod J Integr Tradit Chin Western Med 2016;25:437-9.
Liu Y, Liu C, Li CY, Tan E, Su YW, Fan G, et al
. Construction of two-dimensional liquid chromatography and its application in separation of velvet antler protein. China: The Eighth Chinese Medicine Association; 2015.
Li J, Wei YH, Qin XM, Li ZY. Chemical comparison on different parts of Angelica sinensis
Radix based on NMR metabolomics. Chin Tradit Herb Drugs 2017;48:1409-15.
Na LD, Chen JL, Qin XM, Tian JS. Research on analysis method for chemical variability of E'Jiao based on 1
H-NMR metabolomics. Chin Tradit Herb Drugs 2015;46:2573-9.
Bing GL. Processing technology of pilose antler. China Rural Sci Technol 1998;8:36-7.
Zhang ZZ, Zhi HJ, Qin XM, Li ZY. Chemical comparison of different Farfarae Flos by NMR-based metabolomic approaches. Yao Xue Xue Bao 2015;50:599-604.
Han XJ, Fan ML, Qin XM, Li ZY. Effect of scutellariae radix on metabolism in mice by metabolomics. Chin J Exp Tradit Med Formulae 2016;103:85-91.
Li ZY, Zhang ZZ, Du GH, Qin XM. Comparative analysis of Danggui
and European Danggui
using nuclear magnetic resonance-based metabolic fingerprinting. J Pharm Biomed Anal 2015;103:44-51.
Fan ML, Zhang ZZ, Zhang XQ, Li ZY, Qin XM. Comparation on effect between Phellodendri Chinensis Cortex and Coptis Rhizoma
on organism metabolism in mice Chin Tradit Herb Drugs 2015;46:2892-9.
Tian JS, Shi BY, Feng GM, Qin XM.1
H-NMR metabonomic study on fecal of chronic unpredicted mild stress model of depression in rats. Chin Tradit Herb Drugs 2013;44:3170-6.
Chang YF, Gong WX, Zheng YH, Li JW, Zhou YZ, Qin XM, et al.
Urinary metabolomics study of the effects of Scutellaria baicalensis
Georgi ethanol extract on D-galactose-induced rats. Yao Xue Xue Bao 2016;51:86-92.
Chen L, Liu H, Chen JL, Gao XX, Zhou YZ, Tian JS, et al.
Anti-depressive mechanism of Fufang Chaigui prescription based on neuroendocrine hormone and metabolomic correlation analysis. Zhongguo Zhong Yao Za Zhi 2015;40:4080-7.
Xing J, Yuan SC, Sun HM, Fan ML, Li ZY, Qin XM. Discrimination of patients with Xiao-Chaihu Tang syndrome using 1
H NMR metabonomics and partial least square analysis. Yao Xue Xue Bao 2015;50:1032-7.
Niu QY, Li ZY, Du GH, Qin XM. (1) H NMR based metabolomic profiling revealed doxorubicin-induced systematic alterations in a rat model. J Pharm Biomed Anal 2016;118:338-48.
Xue Y, Li XW, Li ZY, Zeng ZP, Ai FS, Li AP, et al
. UPLC/Q-TOF MS and NMR plant metabolomics approach in studying the effect of growth year on the quality of Polygala tenuifolia
. Acta Pharm Sin 2015;50:340-7.
Mi X, Li ZY, Qin XM, Zhang LZ. Monitoring of chemical components with different color traits of Tussilago farfara
using NMR-based metabolomics. Yao Xue Xue Bao 2013;48:1692-7.
Tian JS, Na LD, Xiang H, Chen JL, Ma XQ, Qin XM. Identification of donkey-hide and its counterfeits based on 1
H-NMR metabolomics. Chin Tradit Herb Drugs 2015;46:255-61.
Peng GJ, Shi BY, Tian JS, Gao S, Qin XM.1
HNMR based metabonomics study on the antidepressant effect of genipin in rat hippocampus. Acta Pharm Sin 2014;49:209-16.
Tian JS, Zuo YM, Sun HF, Qin XM. Analysis on endogenous metabolites in cecal tissue of CUMS-induced depression rats after Xiaoyao San intervention by GC-MS metabolomics. Chin Tradit Herb Drugs 2015;46:1931-6.
Guo XQ, Tian JS, Shi BY, Xing J, Jia JP, Zhang LZ, et al
H-NMR metabonomics on antidepression of Xiaoyao Powder composed of Bupleurum scorzonerifolium
or Bupleurum chinense
. Chin Tradit Herb Drugs 2012;43:2209-16.
Sui Z, Zhang L, Huo Y, Zhang Y. Bioactive components of velvet antlers and their pharmacological properties. J Pharm Biomed Anal 2014;87:229-40.
Wang YH, Sun YS, Zhang L, Liu C, Gong RZ. Research on chemical compositions in different processed products of Cervi Cornu Pantotrichum. Zhongguo Zhong Yao Za Zhi 2018;43:1145-55.
Zong Y, Wang Y, Li H, Li N, Zhang H, Sun J, et al.
Simultaneous quantification and splenocyte-proliferating activities of nucleosides and bases in Cervi Cornu Pantotrichum. Pharmacogn Mag 2014;10:391-7.
Yin RG, Wei S, Li H, Yu H. Introduction of unsupervised learning methods in deep learning. Comput Syst Appl 2016;25:1-7.
Chen JL, Shi BY, Xiang H, Hou WJ, Qin XM, Tian JS, et al.
(1)H NMR-based metabolic profiling of liver in chronic unpredictable mild stress rats with genipin treatment. J Pharm Biomed Anal 2015;115:150-8.
Zong Y, Zhang H, Niu XH, Li H, Sun JM. Alkaloids from Cervi Cornu Pantotrichum and its effect on murine splenocytes proliferation. Zhong Yao Cai 2014;37:752-5.
Duan Y, Wang Y, Zhang J, Sun Y, Wang J. Dietary effects of succinic acid on the growth, digestive enzymes, immune response and resistance to ammonia stress of Litopenaeus vannamei
. Fish Shellfish Immunol 2018;78:10-7.
Feng Y, Yang JC, Shi J, Yang QH, Jiang Y, HU JM. Effects of taurine on testosterone secretion and its mechanism. Acta Vet Zootech Sin 2006;37:1293-6.
Hu XC. Effect of Fatty Acids on Survival Rate of Sertoli Cells in Testis of Mice and Its Mechanisms [doctor]. China: Nanjing University; 2018.
Liu Q, Zhao HW, Zhang AH, Sun H, Zhao XH, Nan Y, et al.
Chinmedomics strategy to discover effective constituents and elucidate action mechanism of Nanshi capsule against kidney-yang deficiency syndrome. Zhongguo Zhong Yao Za Zhi 2016;41:2901-14.
[Figure 1], [Figure 2], [Figure 3], [Figure 4]
[Table 1], [Table 2]