Metabolic profiling reveals therapeutic effects of Herba Cistanches in an animal model of hydrocortisone-induced 'kidney-deficiency syndrome'
Chinese Medicine volume 3, Article number: 3 (2008)
Herba Cistanches (Roucongrong) is effective in treating Shenxu Zheng ('kidney-deficiency syndrome'). However, the mechanisms and systemic metabolic responses to the herbal intervention are unclear.
Using GC-MS-based metabolic profiling, we investigated the metabolic responses to Herba Cistanches intervention in a rat model of the hydrocortisone-induced 'kidney-deficiency syndrome'.
The metabolic profiles of the rats after hydrocortisone injection deviated from the pre-dose metabolic state at different time points, ranging from day 1 to day 10, whereas the metabolic profiles of the rats treated with both hydrocortisone and water extract of Herba Cistanches returned to the pre-dose state on day 10.
The intervention of Herba Cistanches caused a systemic recovery from the hydrocortisone-induced metabolic perturbation in rats. This study also demonstrates that metabolic profiling is useful in studying therapeutic mechanisms of herbal medicines.
Through modulation of biochemical reactions, control mechanisms, and enzyme activities, many drugs or chemicals cause fluctuation of metabolites present in single cells, tissues or body fluids . Metabolic profiling, i.e. probing metabolites of low molecular weight (MW < 1000 Da) by means of advanced analytical instrument coupled with multivariate statistics, can show systemic responses of living systems to xenobiotics. It is also technically feasible to catalogue all the multifactorial heritable and environmentally influenced metabolic profiles of an organism, including the physiopathological consequences of toxin and/or disease-induced disturbances or disequilibria in metabolic regulatory network on a systemic level. To date, metabolic profiling has been established in screening, diagnosis, prognosis of diseases [2–4], and safety evaluation of certain drugs and chemicals [5–11].
Nuclear magnetic resonance (NMR)  and mass spectrometry (MS) , applied alone or in combination, have been used in profiling and characterizing metabolic consequences of toxin and/or disease-induced disturbances. NMR, which does not require tedious sample preprocessing, is a fast and simple method to obtain intrinsic information from complex and intact biological samples. On the other hand, the wide application of hyphenated MS in metabolic profiling is due to its high sensitivity and availability [14, 15]. Notably, GC-MS-based metabolic profiling has been used in discovering mechanisms of drugs and herbicides in vivo, biomarkers of diseases  and effects of altered gene expression on metabolism, and monitoring performance of organisms in biotechnological applications [17–20].
Herba Cistanches (Roucongrong), a common Chinese tonic herb which grows in desert, exhibits marked activities for improving memory  and/or sexual potency , free radical scavenging, anti-aging [23–26] and neuroprotection [27, 28]. For centuries, Herba Cistanches has been effectively used in treating Shenxu Zheng ('kidney-deficiency syndrome') . Recently Herba Cistanches was demonstrated to have ameliorated hydrocortisone-induced kidney disorders ; however, its metabolic consequences are not clear. Our previous study  found that metabolic profiles of rats exposed to hydrocortisone at a high dosage (i.e. an animal model for the 'kidney-deficiency syndrome') , showed a unique biochemical pattern of endogenous metabolites in urine. These results inspired us to study the mechanisms of consistent biochemical changes following hydrocortisone modification, using GC-MS-based metabolic profiling to investigate whether Herba Cistanches could reverse or counteract aberrant metabolic effects of hydrocortisone.
Materials and instrument
Herba Cistanches was purchased from Shanghai Leiyunshang Pharmaceutical Co Ltd (China) and identified as Cistanche deserticola Y C Ma by Dr Mengyue Wang (Laboratory of Pharmacognostics, School of Pharmacy, Shanghai Jiao Tong University) according to a standard protocol . Hydrocortisone solution for injection (0.5%) was purchased from Shanghai Xinyi Pharmaceutical Co (China). The derivatizing reagents were N-methyl-N-trimethylsilyltrifluoroacetamide (MSTFA) (Sigma-Aldrich Inc, USA) and Trimethyliodosilane (TMSI) (Sigma-Aldrich Inc, USA) mixed at a ratio of 1000:1. All reagents used in the experiment were of analytical grade. Ultra-pure water was prepared with Millipore purification system (18.2 MΩ, USA). Metabolic cages were purchased from Suzhou Fengshi Laboratory Animal Experiment Co Ltd (China).
Preparation of Herba Cistanches extract
Five hundred grams of coarsely pulverized plant material was refluxed with 2 L of ultra-pure water for 2 hours. After filtration, the extract was evaporated to about one tenth of the original volume on a Buchi rotary evaporator and was diluted to 250 mL in a volumetric flask with ultra-pure water. The final concentration of crude Herba Cistanches extract was 2 g/mL.
Dosage and sampling
The handling of all animals in this study conformed to the national guidelines and was performed at the Center for Laboratory Animals, Shanghai University of Traditional Chinese Medicine, Shanghai, China. A total of 19 nine-week-old male Wistar rats were purchased from Shanghai Laboratory Animal Co Ltd (China). All animals were kept at a barrier system with regulated temperature (20–22°C) and humidity (60 ± 10%), and on a 12 hour dark/light cycle with lights on at 8:00 am. The rats were given ad libitum food and water. After two weeks' acclimation, the animals were transferred to individual metabolic cages and randomly divided into three groups: (1) treatment group (n = 7) in which hydrocortisone (5%) was injected i.p. at 1.5 mg/100 g of body weight followed by oral administration of Herba Cistanches extract for 10 days; (2) model group (n = 7) in which hydrocortisone (5%) was injected i.p. at 1.5 mg/100 g once daily for 10 days; and (3) control group (n = 5) in which the vehicle was injected i.p. at about 0.6 mL for 10 days . Herba Cistanches was administered to the treatment group at a dose of 20 g/kg following the recommendation by Shen et al. . Twenty-four hour urine samples were collected at specific time intervals: pre-dose (-24 – 0 h), day 1 (0 – 24 h), day 3, day 7, and day 10. All the urine samples were centrifuged (6383 × g, LG 16-W, Beijing Jingli Centrifuge Co Ltd, China) for 10 minutes to remove suspended debris and immediately stored at -80°C for subsequent GC-MS analysis.
Sample preparation and GC-MS
GC-MS was performed according to our previous study with minor modifications . Briefly, each 0.5 μL aliquot of trimethylsilyl (TMS) derivatized analyte was injected into a fused-silica capillary column (17 m × 220 μm inside diameter, 0.11 μm film thickness; HP Ultra-1, Agilent J&W Scientific, USA). GC-MS was conducted on a hyphenated PerkinElmer gas chromatograph and TurboMass-Autosystem XL mass spectrometer (PerkinElmer Inc, USA).
Data processing and multivariate analysis
The GC-MS data were converted into NetCDF format through DataBridge (PerkinElmer Inc, USA). Custom scripts were run in MATLAB 7.0 (The MathWorks Inc, USA) to perform baseline correction, peak deconvolution and alignment, internal standard exclusion, and normalization to the total sum of the chromatogram. The resultant 3-dimensional matrix encompassing arbitrary peak index (paired retention time-m/z), samples (observations), and normalized peak areas (variables) were imported into the SIMCA-P 11.0 software package (Umetrics, Sweden) for multivariate analysis.
Mean-centering was performed column-wise to remove the offsets. All the measured metabolites were treated on an equal level with auto-scaling (scaled to unit variance) prior to multivariate analysis. Principal component analysis (PCA) was performed using the SIMCA-P 11.0 software to reveal the general clustering, grouping, and trends among the subjects without prior knowledge. The first principal component (PC1) represents the most variance in the data. The second principal component (PC2) is orthogonal to PC1, and represents maximum amount of variance not explained by PC1. The remaining principal components were constructed in the same manner. Meanwhile, mean trajectories of PCA scores were used to provide a dynamic indication for the onset, progression, and/or recovery of the syndrome through time. Correlation coefficients from the partial least squares – discriminant analysis (PLS-DA) were used to rank the importance of each variable to further capture the differentially expressed metabolites accountable for the separation between groups. PLS-DA is derived from the partial least squares (PLS) method which is a generalized multiple regression method dealing with multiple collinear predictor and response variables . The PLS-DA was conducted using the SIMCA-P 11.0 software . A typical 7-round cross-validation was performed. One seventh of the samples were excluded from the model in each round so as to validate the model. This procedure was repeated in an iterative manner for cross-validation until each sample had been excluded once.
The differentially expressed metabolites identified from multivariate analysis were also verified in the MATLAB 7.0 software (The MathWorks Inc, USA) by nonparametric Kruskal-Wallis test with a significance level at P < 0.05.
Results and Discussion
Interpretation of GC-MS spectra
Typical GC-MS total ion current (TIC) chromatograms of rat urine on day 10 from the treatment group, model group and control group are shown in Figure 1. Using our optimized GC-MS analysis protocol in association with a software-based peak deconvolution procedure, a total of 117 individual metabolites were consistently detected in at least 90% of the urine samples. Compound identification of peaks of interest was performed by comparing the mass spectral fragment against the NIST (National Institute of Standards and Technology) reference libraries, Wiley libraries and reference standards. We were able to verify 23 out of the 117 metabolites (20%), the majority of which were amino acids, polyamines, fatty acids, purines and adrenal hormones that are mainly involved in energy metabolism, lipid metabolism and amino acid metabolism.
Time-dependent changes in urine samples
The mean trajectories of the PCA scores derived from the model group and treatment group were illustrated in Figure 2. Transient shift in the trajectory plot revealed the dynamic progress of the 'kidney-deficiency syndrome' induced by hydrocortisone alone or in combination with Herba Cistanches treatment. In the model group, metabolic patterns on day 1 and day 3 were different from those on day 7 and day 10, suggesting that the metabolic regulatory network on day 1 and day 3 might have undergone a transient period with high fluctuations, and that the perturbed network might have been restored on day 7 and day 10 which ultimately led to a stable pattern close to the pre-dose state. Analogously, the fact that metabolic pattern of day 1 and day 3 obviously deviated from that of pre-dose in the treatment group indicated the dominant 'kidney-deficiency syndrome' state. In this period, the effects of hydrocortisone were probably dominant over those of the Herba Cistanches extract. These findings were consistent with the general observation that rats from both groups showed less activity on day 1 and day 3. Interestingly, the metabolic patterns on day 7 and day 10 gradually and significantly approached the pre-dose state, suggesting that Herba Cistanches had some counteracting or therapeutic effects on the rats exposed to hydrocortisone. These results support the clinical findings that Herba Cistanches is effective in treating the 'kidney-deficiency syndrome'. In general, both trajectories provide a visual, overall and dynamic picture of the onset, progression and recovery of the 'kidney-deficiency syndrome'.
Comparative metabolic analysis of urine samples
To better understand the metabolic effects of hydrocortisone, we compared the metabolic profiles obtained from the control, model and treatment groups. General clustering of the three groups can be readily observed at various time points, i.e. pre-dose, day 3 and day 10 (Figure 3). While there is no separation tendency in pre-dose urine profiles, the metabolic profiles deviated from those of the control group on day 3 after hydrocortisone exposure. The metabolic perturbation by hydrocortisone appeared in both the model and treatment groups. However, after a consecutive 7-day treatment with Herba Cistanches, metabolic profiles of the treatment group became comparable to those of the control group again, indicating that Herba Cistanches did effectively restore the perturbed metabolism.
Differential identification of metabolic profiles
A cross-validated PLS-DA model was used to identify the key metabolites in different metabolic profiles for easier differentiation between the control group rats and model group (i.e. hydrocortisone-induced) rats with or without Herba Cistanches treatment on day 3 (Table 1). Fold changes in relative concentration of each key metabolite between groups were determined and the corresponding visualization of the changes between groups at pre-dose, day 3 and day 10 was produced (Figure 4). As shown in Figure 4 and Table 1, while most of the endogenous metabolites significantly increased or decreased in the model group, those in the treatment group underwent a transient period as observed on day 1 and day 3 and gradually approached the control (normal) level. For example, as compared to the slight variation (1.1–1.5 times) of the metabolites in the treatment group on day 10, the greatly enhanced levels (1.7–3.0 times) of urinary tyrosine, tyramine, dopamine, and noradrenaline were observed in the model group throughout the experiment. Our previous study showed that the enhanced catecholamine metabolism induced by glucocorticoids resulted in over-consumption of immune functions, thereby leading to the 'kidney-deficiency syndrome' . Herba Cistanches, a tonic herb which improves the immune system , may counteract against some effects of hydrocortisone. Herba Cistanches may also be able to restore a normal metabolic regulatory network. Further experiments with different approaches such as molecular biology, cell biology and plant chemistry are required to delineate the actions of Herba Cistanches (and its constituents) in the 'kidney-deficiency syndrome'.
The present metabolic profiling study using GC-MS showed that Herba Cistanches caused a systemic recovery from the hydrocortisone-induced metabolic perturbation in rats, an animal model for the 'kidney-deficiency syndrome'. This study also demonstrates that metabolic profiling is a useful method to study therapeutic effects of herbal medicines.
gas chromatography-mass spectrometry
nuclear magnetic resonance
principal component analysis
partial least squares – discriminant analysis
Nicholson JK, Lindon JC, Holmes E: 'Metabonomics': Understanding the metabolic responses of living systems to pathophysiological stimuli via multivariate statistical analysis of biological NMR spectroscopic data. Xenobiotica. 1999, 29 (11): 1181-1189. 10.1080/004982599238047.
Brindle JT, Antti H, Holmes E, Tranter G, Nicholson JK, Bethell HWL, Clarke S, Schofield PM, McKilligin E, Mosedale DE, Grainger DJ: Rapid and noninvasive diagnosis of the presence and severity of coronary heart disease using 1H-NMR-based metabonomics. Nat Med. 2002, 8 (12): 1439-1444. 10.1038/nm802.
Constantinou MA, Papakonstantinou E, Benaki D, Spraul M, Shulpis K, Koupparis MA, Mikros E: Application of nuclear magnetic resonance spectroscopy combined with principal component analysis in detecting inborn errors of metabolism using blood spots: A metabonomic approach. Anal Chim Acta. 2004, 511 (2): 303-312. 10.1016/j.aca.2004.02.012.
Beckonert O, Monnerjahn J, Bonk U, Leibfritz D: Visualizing metabolic changes in breast-cancer tissue using 1H-NMR spectroscopy and self-organizing maps. NMR Biomed. 2003, 16 (1): 1-11. 10.1002/nbm.797.
Mortishire-Smith RJ, Skiles GL, Lawrence JW, Spence S, Nicholls AW, Johnson BA, Nicholson JK: Use of metabonomics to identify impaired fatty acid metabolism as the mechanism of a drug-induced toxicity. Chem Res Toxicol. 2004, 17 (2): 165-173. 10.1021/tx034123j.
Waters NJ, Holmes E, Williams A, Waterfield CJ, Duncan Farrant R, Nicholson JK: NMR and pattern recognition studies on the time-related metabolic effects of α-naphthylisothiocyanate on liver, urine, and plasma in the rat: An integrative metabonomic approach. Chem Res Toxicol. 2001, 14 (10): 1401-1412. 10.1021/tx010067f.
Coen M, Lenz EM, Nicholson JK, Wilson ID, Pognan F, Lindon JC: An integrated metabonomic investigation of acetaminophen toxicity in the mouse using NMR spectroscopy. Chem Res Toxicol. 2003, 16 (3): 295-303. 10.1021/tx0256127.
Small-Howard A, Turner H: Exposure to tobacco-derived materials induces overproduction of secreted proteinases in mast cells. Toxicol Appl Pharmacol. 2005, 204 (2): 152-163. 10.1016/j.taap.2004.09.003.
Waters NJ, Waterfield CJ, Farrant RD, Holmes E, Nicholson JK: Metabonomic deconvolution of embedded toxicity: Application to thioacetamide hepato- and nephrotoxicity. Chem Res Toxicol. 2005, 18 (4): 639-654. 10.1021/tx049869b.
Robertson DG: Metabonomics in toxicology: A review. Toxicoll Sci. 2005, 85 (2): 809-822. 10.1093/toxsci/kfi102.
Robertson DG, Bulera SJ: High-throughput toxicology: Practical considerations. Curr Opin Drug Discovery Dev. 2000, 3 (1): 42-47.
Nicholson JK, Connelly J, Lindon JC, Holmes E: Metabonomics: A platform for studying drug toxicity and gene function. Nat Rev Drug Discov. 2002, 1 (2): 153-161. 10.1038/nrd728.
Taylor J, King RD, Altmann T, Fiehn O: Application of metabolomics to plant genotype discrimination using statistics and machine learning. Bioinformatics. 2002, 18 (SUPPL 2): S241-S248.
Wilson ID, Nicholson JK, Castro-Perez J, Granger JH, Johnson KA, Smith BW, Plumb RS: High resolution "ultra performance" liquid chromatography coupled to oa-TOF mass spectrometry as a tool for differential metabolic pathway profiling in functional genomic studies. J Proteome Res. 2005, 4 (2): 591-598. 10.1021/pr049769r.
Jonsson P, Gullberg J, Nordstrom A, Kusano M, Kowalczyk M, Sjostrom M, Moritz T: A strategy for identifying differences in large series of metabolomic samples analyzed by GC/MS. Anal Chem. 2004, 76 (6): 1738-1745. 10.1021/ac0352427.
Ohdoi C, Nyhan WL, Kuhara T: Chemical diagnosis of Lesch-Nyhan syndrome using gas chromatography-mass spectrometry detection. J Chromatogr, B: Anal Technol Biomed Life Sci. 2003, 792 (1): 123-130. 10.1016/S1570-0232(03)00277-0.
Fiehn O, Kopka J, Dormann P, Altmann T, Trethewey RN, Willmitzer L: Metabolite profiling for plant functional genomics. Nat Biotechnol. 2000, 18 (11): 1157-1161. 10.1038/81137.
Lafaye A, Junot C, Pereira Y, Lagniel G, Tabet JC, Ezan E, Labarre J: Combined proteome and metabolite-profiling analyses reveal surprising insights into yeast sulfur metabolism. J Biol Chem. 2005, 280 (26): 24723-24730. 10.1074/jbc.M502285200.
Schauer N, Steinhauser D, Strelkov S, Schomburg D, Allison G, Moritz T, Lundgren K, Roessner-Tunali U, Forbes MG, Willmitzer L, Fernie AR, Kopka J: GC-MS libraries for the rapid identification of metabolites in complex biological samples. FEBS Lett. 2005, 579 (6): 1332-1337. 10.1016/j.febslet.2005.01.029.
Willse A, Belcher AM, Preti G, Wahl JH, Thresher M, Yang P, Yamazaki K, Beauchamp GK: Identification of major histocompatibility complex-regulated body odorants by statistical analysis of a comparative gas chromatography/mass spectrometry experiment. Anal Chem. 2005, 77 (8): 2348-2361. 10.1021/ac048711t.
Wang XW, Wang XF, Wu LY: Improvement of rememberance of mice of phenylethanoid glycosides of Cistanche deserticola. Rep Chin Pharm. 2002, 19: 41-42.
Xie JH, Wu CF: Effect of ethanolic extract of Cistanche deserticola on the contents of monoamine neurotransmitters in rat brain. Zhongcaoyao. 1993, 24: 417-419.
Li LL, Wang XW, Wang XF: Antilipid peroxidation and antiradiative action of glycosides in herba Cistanches. Chin J Chin Mater Med. 1997, 22 (6): 364-367.
Shahat AA, Nazif NM, Abousetta LM, Ibrahim NA, Cos P, Van Miert S, Pieters L, Vlietinck AJ: Phytochemical investigation and antioxidant activity of Duranta repens. Phytother Res. 2005, 19 (12): 1071-1073. 10.1002/ptr.1766.
Gao J, Igarashi K, Nukina M: Three new phenylethanoid glycosides from Caryopteris incana and their antioxidative activity. Chem Pharm Bull. 2000, 48 (7): 1075-1078.
Kyriakopoulou I, Magiatis P, Skaltsounis AL, Aligiannis N, Harvala C: Samioside, a new phenylethanoid glycoside with free-radical scavenging and antimicrobial activities from Phlomis samia. J Nat Prod. 2001, 64 (8): 1095-1097. 10.1021/np010128+.
Deng M, Zhao JY, Ju XD, Tu PF, Jiang Y, Li ZB: Protective effect of tubuloside B on TNF alpha-induced apoptosis in neuronal cells. Acta Pharmacol Sin. 2004, 25: 1276-1284.
Geng XC, Song LW, Pu XP, Tu PF: Neuroprotective effects of phenylethanoid glycosides from Cistanches salsa against 1-methyl-4-phenyl-1, 2, 3, 6-tetrahydropyridine (MPTP) induced dopaminergic toxicity in C57 mice. Biol Pharm Bull. 2004, 27: 797-801. 10.1248/bpb.27.797.
He W, Shu X, Zong G, Shi M, Xiong Y, Chen M: Kidney reinforcing and yang supporting action of cistanche deserticola Y. C. Ma before and after preparation. Zhongguo Zhongyao Zazhi. 1996, 21 (9): 534-537. 575
Shen LZ, Zhong XY, Wang SX: Effect of Cistanche deserticola on normal and deficiency Shen-yang rats. Zhongyao Yaoli Yu Linchuang. 2001, 17 (1): 17-18.
Chen M, Zhao L, Jia W: Metabonomic study on the biochemical profiles of a hydrocortisone-induced animal model. J Proteome Res. 2005, 4 (6): 2391-2396. 10.1021/pr050158o.
Chen Q, Yi NY: The animal models and drugs for Yin-deficiency and Yang-deficiency. Experimental methodology of pharmacological research in traditional Chinese medicine. Edited by: Chen Q. 1993, Beijing: People's Health Publishing House, 982-984.
National Committee of Pharmacopoeia: Pharmacopoeia of the People's Republic of China. 2005, Beijing: Press of Chemical Industry, 1: 90-
Multi- and Megavariate Data Analysis Part I: Basic Principles and Applications, Second revised and enlarged edition. [http://www.umetrics.com/default.asp/pagename/training_literature/c/5]
SIMCA-P and SIMCA-P+ 11 User Guide. [http://www.umetrics.com/default.asp/pagename/downloads_userguide/c/3]
Chin HL, Su YC: Study on the pharmacological effects of Cistanche deserticola Ma. Zhongguo Zhongyao Zazhi. 1993, 19: 143-146.
The authors wish to thank Yumin Liu, Jingjing Gu, Yuan He and Congying Guo for their assistance with the GC-MS analysis and animal experiments. This study was financially supported by the Young Faculty Grant from Shanghai Jiao Tong University (No. A2816A), Shanghai Leading Academic Discipline Project (No. T0301) and the National Natural Science Fund (No. 30572286).
The author(s) declare that they have no competing interests.
YQ and MC performed the animal and instrumental experiments and drafted the manuscript. MS and GX analyzed the data and helped draft the manuscript. XL helped with animal experiments. MZ and AZ prepared Herba Cistanches extract. JJ provided suggestions for the experiment design and the manuscript. WJ supervised all the experiments and revised the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Qiu, Y., Chen, M., Su, M. et al. Metabolic profiling reveals therapeutic effects of Herba Cistanches in an animal model of hydrocortisone-induced 'kidney-deficiency syndrome'. Chin Med 3, 3 (2008). https://doi.org/10.1186/1749-8546-3-3