Dan-Dan Wei
Jun-Song Wang
Jin-Ao Duan
Ling-Yi Kong- 1Jiangsu Collaborative Innovation Center of Chinese Medicinal Resources Industrialization, National and Local Collaborative Engineering Center of Chinese Medicinal Resources Industrialization and Formulae Innovative Medicine, Nanjing University of Chinese Medicine, Nanjing, China
- 2State Key Laboratory of Natural Medicines, Department of Natural Medicinal Chemistry, China Pharmaceutical University, Nanjing, China
- 3Center for Molecular Metabolism, Nanjing University of Science and Technology, Nanjing, China
Huang-Lian-Jie-Du-Decoction, a traditional Chinese formula, has been reported to protect liver from various injuries. Two cholestasis models of rats induced by thioacetamide and by bile duct ligation were established and treated with Huang-Lian-Jie-Du-Decoction. Nuclear Magnetic Resonance-based urinary metabolic profiles were analyzed by orthogonal partial least squares discriminant analysis and univariate analysis to excavate differential metabolites associated with the injuries of the two models and the treatment effects of Huang-Lian-Jie-Du-Decoction. The two cholestatic models shared common metabolic features of excessive fatty acid oxidation, insufficient glutathione regeneration and disturbed gut flora, with specific characteristics of inhibited urea cycle and DNA damage in thioacetamide-intoxicated model, and perturbed Kreb's cycle and inhibited branched chain amino acid oxidation in bile duct ligation model. With good treatment effects, Huang-Lian-Jie-Du-Decoction could regain the balance of the disturbed metabolic status common in the two cholestasis injuries, e.g., unbalanced redox system and disturbed gut flora; and perturbed urea cycle in thioacetamide-intoxicated model and energy crisis (disturbed Kreb's cycle and oxidation of branched chain amino acid) in bile duct ligation model, respectively.
Introduction
Traditional Chinese medicine (TCM) has been used for healthcare for thousands of years in China (Cheung, 2011). As different from western medicine in both methodology and philosophy, TCM concentrates on the overall functional state of the patient and the adjustment of its balance, and has gained ever-increasing interest worldwide, especially for the treatment of complex diseases (Jiang, 2005). The major challenge to the adoption and popularity of TCM is the mechanistic study of its efficacy due to its complicated principles (Wang R. et al., 2012). Similar to TCM, systems biology also treats the body as a whole and investigates the life holistically (Van der Greef, 2011), which was deemed as a technique that has the potential to make a real breakthrough in TCM (Qiu, 2007). Metabolomics, an important component of systems biology, is a sensitive tool to detect various stimuli by snapshotting metabolic status of the body (Nicholson et al., 1999; Lindon et al., 2003; Lucio, 2009). Metabolomics has been widely used for the assessment of the toxicity of TCM (Zira et al., 2013), disease diagnosis (Manna et al., 2013; Zeng et al., 2014), and phenotyping (An et al., 2012). The most common analytical techniques for metabolic profiling are nuclear magnetic resonance (NMR) spectroscopy, gas chromatography-mass spectrometry (GC-MS) and liquid chromatography mass spectrometry (LC-MS). Despite of a lower sensitivity than MS, NMR has been widely applied in metabolic profiling for its associated advantages such as non-destructive, non-biased, separation free, short analysis time, and rich in structural information (Barba et al., 2008; Wishart, 2008; Wheelock et al., 2013).
Liver is critical for bile formation, amino acid utilization and ammonia detoxification, and is also the tissue where glycolysis, gluconeogenesis, and the synthesis of certain plasma proteins happen (Beyoǧlu and Idle, 2013). The metabolic rate and energy production and consumption in liver are extremely fast, which makes liver vulnerable to drugs and toxins and sensitive to many pathological insults (Whitfield et al., 2005; Vinaixa et al., 2010). Flowing from liver to duodenum through hepatobiliary transport systems, bile plays an important role in the digestion of fats. Any obstruction in bile duct systems could result in an accumulation of bile acids and other harmful compounds, leading to cholestatic injury. Cholestatic injuries could be formed by various means. The hepatoxin thioacetamide could inhibit bile secretion, producing intrahepatic cholestatic injuries (Shukla et al., 1992); ligation of the common bile duct could lead to a mechanical blockage in the duct system, inducing obstructive extrahepatic cholestatic injuries.
Huang-Lian-Jie-Du-Decoction (HLJDD; also known as oren-gedoku-to or Hwangryun-Hae-Dok-Decotion in Japan), representative ancient antipyretic and detoxifying traditional Chinese medicine recipe, is comprised of Rhizoma coptidis, Radix scutellariae, Cortex phellodendri, and Fructus gardeniae in a weight ratio of 3:2:2:3, functioning as clearing “heat” and removing “poison.” It has been used to treat in clinical for nearly two thousand years, was mentioned firstly in Wang Tao's “Wai Tai Mi Yao” two thousand years ago, it has been widely used to treat inflammatory-related diseases, such as liver injury (Wei et al., 2015; Lv et al., 2017), ischemic stroke (Wang et al., 2014; Zhang et al., 2017a), arthritis (Hu et al., 2013), and dermatitis (Chen et al., 2016). Its component herb, Rhizoma coptidis and its main alkaloid berberine exhibited hepatoprotective activity on acute hepatotoxicity and liver fibrosis induced by various drugs and chemicals, e.g., acetaminophen (Janbaz and Gilani, 2000), carbon tetrachloride (Ye et al., 2009; Feng et al., 2010), diethylnitrosamine and phenobarbital (Zhao et al., 2008). Their protection on liver injuries was mainly mediated by inhibition of various inflammatory factors: tumor necrosis factor-α, cyclooxygenase-2, and inducible nitric oxide synthase expression (Zhao et al., 2008; Domitrovi et al., 2011). Baicalin, a major flavone from RS, has also showed protective effect on liver injuries induced by acetaminophen (Jang et al., 2003), lipopolysaccharide/d-galactosamine (Liu et al., 2008; Wan et al., 2008), ischemia/reperfusion (Kim et al., 2010), and iron overloading (Zhao et al., 2005).
Cholestatic injury was associated with inflammation, and our previous study has demonstrated anti-inflammatory effects of HLJDD both in vitro and in vivo (Lu et al., 2011). HLJDD, its main herb and representative components have been indicated to possess therapeutic potential in the treatment of stroke (Wang et al., 2013, 2014; Zhang et al., 2015, 2016, 2017a,b; Zhu et al., 2018), chronic cholestasis injury (Wei et al., 2015) and multiple organ failure during sepsis (Liao et al., 2016a; Li et al., 2017; Lv et al., 2017; Xu et al., 2017). Herein the aim of the present work was to comparatively study the two cholestatic injuries induced by thioacetamide and by bile duct ligation, and their treatment by HLJDD using an NMR-based metabolomic approach. Multivariate analysis technique orthogonal partial least squares-discriminant analysis (OPLS-DA), combined with univariate analysis, was applied to excavate metabolic features associated with metabolic and mechanical cholestatic injuries of the two models and the treatment effects of HLJDD on them.
Materials and Methods
Chemicals and Reagents
Rhizoma Coptidis, Radix Scutellariae, Cortex Phellodendri Chinensis, and Fructus Gardeniae were obtained from Jiangsu medicine company (Nanjing, China) and were authenticated by Professor Jin-ao Duan, department of Resources and Authentication of Chinese Medicine, school of Pharmacy, Nanjing University of Chinese Medicine (Nanjing, China) as Coptis teeta Wall (Ranunculaceae), Scutellaria baicalensis Georgi (Lamiaceae), Phellodendron chinense Schneid. (Rutaceae) and Gardenia jasminoides Ellis (Rubiaceae), respectively. Voucher specimens were deposited at the herbarium of Jiangsu Collaborative Innovation Center of Chinese Medicinal Resources Industrialization, Nanjing University of Chinese Medicine, with the voucher number of RC201335, RS201323, CP201329, and FG201319.
Geniposidic acid, geniposide, chlorogenic acid, phellodendrine, magnoflorine, tetrahydropalmatine, coptisine, jatrorrhizine, baicalein, palmatine, berberine, baicalin and wogonin were bought from Tianjin Chemical Industrial Co. (Tianjin, China). Formic acid, methanol and acetonitrile were of HPLC grade, and were obtained from Merck KGaA (Darmstadt, Germany). Sodium 3-trimethylsilyl-1-(2,2,3,3-2H4) propionate was bought from Sigma Chemical Co. (St. Louis, MO, USA). Deuterium oxide (99.9%) was purchased from Sea Sky Bio Technology Co. Ltd (Beijing, China). Thioacetamide were bought from Shanghai Lingfeng Chemical Reagent Co., Ltd (Shanghai, China). All reagents were of analytical grade with purity over 98%. Ultra-pure distilled water was prepared from a Milli-Q purification system.
Preparation of HLJDD Extract Solution and the Reference Solution
Rhizoma Coptidis, Radix Scutellariae, Cortex Phellodendri Chinensis, and Fructus Gardeniae in a weight ratio of 3:2:2:3, reaching a total weight of 5.0 kg, was dried, powdered, and extracted with 70% ethanol (1:10, 1:8 and 1:5, w/v) under reflux three times for 1 h each as our previous report (Wang et al., 2013). The extract solution was combined, filtered through 5-layer gauze to give a filtrate, which was concentrated on a rotary vacuum evaporator, and freeze-dried in a vacuum to collect HLJDD extract (1,525.432 g, yield: 30.51%). 1.00 mg HLJDD extract was accurately weighed and dissolved in 5.00 ml methanol to obtain a sample solution at the concentration of 200 μg/ml, while references of approximately 2 mg were accurately weighed and placed in a 10 ml volumetric flask to prepare a mixture of stock solution, with concentrations of 209 μg/ml (geniposidic acid), 308 μg/ml (chlorogenic acid), 106 μg/ml (phellodendrine), 135 μg/ml (magnoflorine), 263 μg/ml (tetrahydropalmatine), 143 μg/ml (coptisine), 120 μg/ml (jatrorrhizine), 268 μg/ml (baicalein), 134 μg/ml (palmatine), 106 μg/ml (berberine), 83 μg/ml (baicalin), and 243 μg/ml (wogonin).
Qualitative and Quantitative Analysis of HLJDD Extract Solution by Ultra-High Pressure Liquid Chromatography With a Linear Ion Trap-High Resolution Orbitrap Mass Spectrometry System (UPLC-LTQ-Orbitrap MS)
The UPLC analyses were performed on a Dionex Ultimate 3000 UPLC system (Thermo Fisher, Hopkinson, MA, USA), hyphenated with a linear ion trap orbitrap XL mass spectrometry system, which was equipped with an ESI source. Chromatographic separations were carried out using a Waters Acquity UPLC column (BEH C18 100 × 2.1 mm, 1.7 μm) at 30°C. The mobile phase was composition of 0.1% formic acid-water (v/v, phase A) and acetonitrile (phase B), in a gradient program: 0–1 min, 3% B; 1–10 min, 3–20% B; 10–25 min, 20–50% B; 25–27 min, 50–95% B. The flow rate was set at 0.4 ml/min and the injection volume was 2 μL.
The Orbitrap instrument was calibrated using a solution containing caffeine, methionine-arginine-phenylalanine-alanine, and Ultramark according to the manufacturer's instructions before use, and signals from each of the orbitrap outer electrodes were amplified by a differential amplifier and transformed into a frequency spectrum by fast Fourier transformation. Mass spectrometer was operated with an ESI probe in positive and negative ion modes, the mass range was between 100 and 1,000, and the resolution was 30,000. All data were acquired in centroid mode, and operation of the entire LC/MS instrumentation was controlled using LCQ Xcalibur software (Thermo Fisher Scientific, Waltham, MA, USA). Parameters of the ion source were as follows: capillary temperature and voltage at 275°C and 40 V, respectively, source voltage at 5 kV and tube lens voltage at 80 V. Nitrogen was used as the sheath gas with a flow rate of 8 arbitrary units.
Animal Treatment
Seventy-two male Sprague-Dawley rats (with an average age of 8 weeks and weight of 220–240 g), were bought from Laboratory Animal Research Center, Nanjing University (Nanjing, China). Rats were housed in a well-ventilated room at constant ambient temperature (25 ± 2°C) and air humidity (50 ± 10%) with a light/dark cycle of 12 h, and were allowed ad libitum access to water throughout the study. Animal experiments were reviewed, approved and conducted in accordance with guidelines set by National Institutes of Health and the Institutional Animal Care and Use Committee at Nanjing University of Chinese Medicine (approbation number: NJUCM2013068). Standard biosecurity procedures were followed throughout the study in a Good Laboratory Practice-compliant experimental trials platform in Nanjing University of Chinese medicine, which maintained a strict rodent control policy. All the rooms and housing environment (ie walls, floor, bars, bottle and roof of the cage) were completely disinfected every 3 days. All food, water and bedding were packed in a bulk sterilizing chamber and vacuum sterilized with a vacuum autoclave before use. All the operation of animal tissues, organs, body parts, carcasses waste, fluid, discarding urine, sharps and disposal of experimental refuse were in accordance with the Centers for Disease Control's recommendations and the Environmental Protection Agency's recommendations to prevent introduction of pathogens.
After acclimatization for a week, rats were randomly divided into six groups (n = 12): two model groups TAA (intoxicated with thioacetamide) and BDL (bile duct ligation), two HLJDD treatment groups THC (on TAA model) and BHD (on BDL model), and two control groups NC and HLD. The HLD, THC and BHD groups were intragastrically (i.g.) administrated with HLJDD extract suspended in 0.5% carboxymethyl cellulose sodium salt (1.47 g/kg/day), while the NC, TAA and BDL groups were administrated with an equivalent amount of vehicle for 1 week.
After an overnight fasting, the TAA and THC rats were injected intraperitoneally with thioacetamide (dissolved in saline, 200 mg/kg), the HLD and NC rats were intraperitoneally administrated with the same amount of saline. BDL and BHD rats were anesthetized by intraperitoneal injection with chloral hydrate (300 mg/kg) and kept in anesthetization throughout the experiment. After a midline incision in the stomach under sterile conditions, gastroduodenal ligament was exposed to allow the isolation of the common bile duct, which was ligated singly with 4-0 nylon silk sutures, followed by careful suturing of the peritoneum and muscle layers as well as the skin wound.
Urine Samples Collection
All rats were ad libitum accessed to standard rat pellet diet and tap water. HLD, THC, and BHD rats were i.g. administrated with HLJDD solution once a day after the operation, while NC, TAA, and BDL rats were treated with an equivalent amount of vehicle. Animals were housed individually in metabolic cages to collect urine samples at the following time periods: 0–8 h (T1), 8–24 h (T2), 24–32 h (T3), 32–48 h (T4), 48–72 h (T5), 72–96 h (T6), 96–120 h (T7) and 120–144 h (T8) after treatment. All urine samples were collected over dry ice with the addition of sodium azide (0.1% w/v urine solution) to the collection vessel to minimize bacterial contamination. At the end of the collection, urine was centrifuged at 12,000 g, 10 min, at 4°C to remove food. All the urine samples for NMR spectroscopic analysis were stored in multiple aliquots at −80°C before use.
Pathological Assessment
Overnight fasted rats were euthanized at 48 and 144 h after operation, and hearts, livers, spleens, lungs, kidneys and brains tissues were collected, and then preserved in 10% buffered formalin. A series of adjacent 3-μm-thick sections were cut and stained with hematoxylin & eosin for independent histological assessment by an experienced pathologist from Jiangsu Provincial Institute of Materia Medica (Nanjing, China).
Acquisition of 1H NMR Spectra of Urine
The deep-frozen urine samples were thawed at 4°C overnight, and then were vortexed to dissolve any precipitates. An aliquot of 400 μL urine was mixed with 200 μL phosphate buffered saline (0.2 M Na2HPO4 and 0.2 M NaH2PO4, pH 7.0), allowing to equilibrate for 10 min prior to centrifugation (12,000 g for 10 min, at 4°C). The supernatant was added with 50 μL Sodium 3-trimethylsilyl-1-(2,2,3,3-2H4) propionate solution (dissolved in deuterium oxide, 1 mg/ml).
All NMR spectra were recorded using a Bruker AV 500 MHz spectrometer, with a nuclear overhauser enhancement spectroscopy presaturation pulse sequence (relaxation delay-90°-t1-90°-tm-90°-acquire-FID). Typically, 64 free induction decay was collected into 32 K data points, using a spectral width of 10 kHz, an acquisition time per scan of 2.54 s, recycle delay of 2 s and a mixing time of 100 ms. With TopSpin software (version 3.0, Bruker Biospin, Germany), all the spectra were automatically phased and baseline corrected, and calibrated to Sodium 3-trimethylsilyl-1-(2,2,3,3-2H4) propionate at 0.00 ppm.
Processing of NMR Data
With the removal of signals of water and urea (4.2–5.9 ppm), the NMR spectra were aligned in R, a freely available, open-source software (R Development Core Team, http://cran.r-project.org/). Regions between 0.5 and 9.8 ppm were binned, probabilistic quotient normalized, mean-centered and pareto scaled in R prior to multivariate statistical analysis.
Identification of Metabolites
Aided by Statistical total correlation spectroscopy, metabolites were identified by comparison with the standard references in the Human Metabolome Database (www.hmdb.ca) under the help of Chenomx NMR Suite, version 4.0 (Chenomx Inc., Edmonton, Canada). These assignments were finally confirmed by 2D NMR techniques, e.g., 1H-1H total correlation spectroscopy and 1H-13C heteronuclear single-quantum correlation.
Multivariate Analysis
Supervised orthogonal partial least squares-discriminant analysis (OPLS-DA) was performed using scripts written in R. A repeated 2-fold cross-validation and permutation test (n = 2,000) were then performed to validate the OPLS-DA models. Goodness-of-fit (R2) and goodness of prediction (Q2) of established models was calculated: typically a model with a Q2 value of greater than 0.40 is considered good, and with Q2 values over 0.70 is robust (Lourenço et al., 2011). Receiver operating characteristic (ROC) plots were delineated to evaluate classification performance using R-package ROCR (http://rocr.bioinf.mpi-sb.mpg.de). The area under the ROC (AUROC) is a measure of the overall accuracy, with a value of 1.000 and 0.500 representing perfect discrimination and random classification, respectively.
Univariate Analysis
Levels of metabolites of each group throughout the experiment were calculated corresponding to the normalized average integration of area (Zira et al., 2013). To deduct the metabolic alterations of healthy rats among different time points, the variations of metabolites in HLD, TAA, BDL, THC, and BHD groups were normalized by NC group. Considering the potential effect of HLJDD on metabolic status of healthy rats, the variations of metabolites in THC and BHD groups were normalized by HLD group. According to the conformity to normality, the statistical significance was calculated using either parametric (t-test) or non-parametric statistical test (Wilcoxon signed rank test) for metabolites between groups using R, with p < 0.05 as significant.
Results
Histopathological Inspection
It revealed no pathological alterations in hearts, spleens, lungs and brains for all groups and also no noticeable changes in the livers and kidneys (Figures 1a–h) of NC and HLD rats. At 48 h, the livers of TAA rats exhibited inflammatory foci, characteristics of monocytes in mesenchymal (Figure 1i), and kidneys of TAA rats were slightly swollen and lobulated with active glomerulogenesis (Figure 1k). The pathological alterations of livers and kidneys in TAA rats were both recovered at 144 h (Figures 1j,l). On the treatment of HLJDD, THC rats showed attenuated inflammatory cell infiltration in hepatocytes at 48 h (Figures 1m,n), in consistent with the previously reported anti-inflammatory effect of HLJDD (Lu et al., 2011). Slight necrosis was observed in kidneys of THC rats at 48 h (Figures 1o,p).
Figure 1. Treatment effect of Huang-lian-jie-du decoction (HLJDD) on the histology of livers and kidneys of thioacetamide-intoxicated (TAA) and bile duct ligation (BDL) rats sacrificed at 48 and 144 h. (a–d) Normal control rats (NC). (e–h) NC rats administrated with HLJDD (HLD): showing no signs of pathological change. (i–l) TAA rats: inflammatory foci in liver (i) and slightly swollen kidneys (k). (m–p) TAA rats treated with HLJDD (THC): slight inflammatory cell infiltration in hepatocytes (m,n) and partial necrosis in kidneys (o,p). (q–t) BDL rats: intense focal necrosis (q), proliferation of small bile duct and fibrogensis (r), collagenic exudation (s) and necrosis of epithelial cell (t). (u–x) BDL rats treated with HLJDD (BHD): slight necrosis (u), partial proliferation of small bile duct and fibrogensis (v), mild collagenic exudation (w), and no signs of pathological changes (x).
Livers of BDL rats exhibited severe hepatocellular injury with intense focal necrosis and degeneration at 48 h (Figure 1q), and striking cholangiocyte proliferation, stellate cell activation and fibrogensis at 144 h after operation (Figure 1r), in consistent with those previously reported (Chen et al., 2012). Renal slices of BDL rats presented significant intense collagenic exudation, around by severe necrosis of epithelial cell at 48 h (Figure 1s), and decreased volume of glomeruli at 144 h (Figure 1t). On the treatment of HLJDD, the collagenic exudation was strikingly attenuated in livers of BHD rats at 48 h (Figure 1u) and recovered at 144 h (Figure 1v), and renal slices exhibited slight necrosis and moderate cholangiocyte proliferation at 48 h (Figure 1w), which was recovered to normal at 144 h (Figure 1x). The pathological changes induced by hioacetamide and bile duct ligation were mainly in livers and kidneys. Thioacetamide produced a mild injury as compared with a relatively severer damage induced by bile duct ligation, both of which peaked at 48 h and attenuated at 144 h, suggesting a self-regulating and self-regenerating capacity of liver and kidney. HLJDD could ameliorate the lesions produced in both models.
UPLC-LTQ-Orbitrap MS Analysis of HLJDD Extract Solution
Thirty-two main base peaks in the chromatograms of HLJDD extract solution were successfully identified according to the retention time, accurate mono-molecular weight and by comparison with the standard references, including 2 phthalate of bis(2-ethylhexyl) phthalate and dioctyl phthalate, 2 phenols of chlorogenic acid and sinapyl aldehyde, 3 iridoid glycosides of geniposide, geniposidic acid and obacunone, 11 flavonoids of baicalin, baicalein-7-O-β-D-glucopyranoside, apigenin 7-O-β-D-glucuronide, baicalein, norwogonin, luteolin/kaempferol, wogonin, 3′,4′,5′-O-trimethyltricetin, prunetin, puararin, dihydrowogonin, and 14 alkaloids of phellodendrine, magnoflorine, thalifendine, berberastine, 8-oxypalmatine, tetrahydropalmatine, epiberberine, coptisine, jatrorrhizine, columbamine, berberrubine, worenine, palmatine and berberine (Figures 2A,C, Table 1). Thirteen standard reference chemcials were employed as external standard solutions, and calibration curves were generated for each reference by linear regression analysis to determine the contents in HLJDD extract (Table 2). Thirteen major components were successfully identified and determined: Geniposidic acid (1), geniposide (2), chlorogenic acid (3), phellodendrine (4), magnoflorine (6), tetrahydropalmatine (11), coptisine (13), jatrorrhizine (14), baicalein (17), palmatine (21), berberine (22), baicalin (25), and wogonin (28) (Figures 2B,D).
Figure 2. (A,C) UPLC chromatogram of the HLJDD extract with its secondary metabolites assigned. (B,D) The typical UPLC chromatograms of a mixture of thriteen standards (A,B in ESI+ mode, C,D in ESI- mode): Geniposidic acid (1, 3.82 min), geniposide (2, 6.13 min), chlorogenic acid (3, 6.25 min), phellodendrine (4, 7.71 min), magnoflorine (6, 8.30 min), tetrahydropalmatine (11, 12.65 min), coptisine (13, 12.67 min), jatrorrhizine (14, 12.83 min), baicalein (17, 13.79 min), palmatine (21, 14.28 min), berberine (22, 14.36 min), baicalin (25, 17.91 min), and wogonin (28, 20.67 min).
Table 1. Metabolites putatively identified by UPLC-LTQ-orbitrap-MS analysis in the extract of HLJDD.
Table 2. The quantification of 13 reference compounds in the extract of HLJDD (n = 6).
Urinary Metabolites Identification
A combination of statistical total correlation spectroscopy, 1H-1H total correlation spectroscopy and 1H-13C heteronuclear single-quantum correlation was used to assign metabolites in urine, in a manner the same as before (Wei et al., 2014a,b). As a result of the joined effort of these techniques, 58 metabolites were identified and assigned (Table S1 and Figure 3).
Figure 3. Typical 500 MHz 1H NMR spectra for urine of two cholestatic injuries model groups BDL (A, bile duct ligation) and TAA (B, intoxicated with thioacetamide), two HLJDD treatment groups BHD (C, on BDL model) and THC (D, on TAA model), and two control groups HLD (E) and NC (F) with metabolites assigned. MMA, methylamine; DMA, dimethyl amine; TMA, trimethylamine; TMAO, trimethylamine N-oxide; DMG, dimethyl glycine; Cr, creatine; PCr, phosphocreatine.
Comparison of Two Cholestatic Injuries and the Effects of HLJDD
Direct comparison of the two models was made by OPLS-DA analysis of urinary 1H NMR data of TAA, BDL, and NC groups. The BDL group dominated the variations, being far away from both the NC and TAA groups. NC and TAA groups, on the other hand, could only be marginally differentiated along the vertical axis. In consistent with the pathological observations, the score and the trajectory plots (Figures S1a,b) suggested that bile duct ligation induced mechanical cholestatic injuries was severer than thioacetamide induced metabolic cholestatic injuries, which might due to the persistent obstruction of bile duct in BDL model versus administration of thioacetamide only once in TAA model. Color-coded loading plot (Figure S1c) provided a direct visualization of the contribution of variables to group separation. Compared with TAA and NC groups, BDL group exhibited significantly increased levels of branched-chain amino acids (BCAAs: leucine, isoleucine and valine), lysine/arginine, phosphocholine, glycine and N-phenylacetyl-glycine, and strikingly decreased levels of citrate, 2-oxoglutarate (2-OG), fumarate and hippurate.
The two cholestatic injuries and the treatment effects of HLJDD were compared according to the metabolomics profiles of 1H NMR data for HLD, NC, TAA, and THC groups (TAA model) and HLD, NC, BDL and BHD groups (BDL model). In the OPLS-DA score (Figure 4A) and trajectory plots (Figure 4B) for TAA model, the four groups formed three distinct clusters along the horizontal axis (T1 score): NC and HLD to the right, TAA to the left, and THC in the middle. TAA and THC groups showed a zigzag trajectory: being furthest away from the NC and HLD groups at T1, then moving back gradually till T4; from T5 and T6, moving toward the NC and HLD groups. The metabolomic profiles of groups in the BDL model resembled those in the TAA model (Figure 4C). However, the trajectories for the BDL and BHD groups (Figure 4D) showed trends different from those for the TAA and THC groups: being the furthest away from the NC and HLD groups at T1 and then moving steadily toward them. The NC and HLD rats could be differentiated along the vertical axis of corresponding score plots in both models. Therefore, the horizontal axis represented pathological changes in the two models; and the vertical axis reflected variations irrelevant to the effects of the two models, e.g., the effects of HLJDD on healthy rats.
Figure 4. OPLS-DA analysis of 1H NMR data in urine for normal control rats (NC), healthy rats administrated with HLJDD (HLD), cholestatic injury rats induced by thioacetamide (TAA) and treated with HLJDD (THC) (A,B), and NC, HLD, cholestatic injury rats induced by bile duct ligation (BDL) and treated with HLJDD (BHD) (C,D). (A,B) scores plots; (C,D) mean trajectory plots.
Influence of HLJDD on Healthy Rats
An OPLS-DA analysis was performed on HLD and NC groups of all time-points to explore the metabolic alterations of HLJDD on healthy rats. The score (Figure 5A) and the mean trajectory plots (Figure 5B) showed a good separation between the two groups. In the loading plot (Figure 5C) and S-plot (Figure 5D), HLD group showed significantly increased levels of propionate, 2-hydroxybutyrate (2-HB), lysine/arginine, glutamate, succinate, choline and guanidoacetate, and noticeably decreased levels of citrate, 2-OG, pyridoxine, phenylalanine and hippurate. The OPLS-DA model constructed for HLD and NC groups demonstrated a good prediction with R2 at 0.620 and Q2 at 0.500 (Figure 5E) and a favorable discrimination with AUROC of 0.928 (Figure 5F). The variations of metabolites in HLD group were shown in lattice plots (Figure 6), and their levels relative to the NC group were presented in color-coded table (Table 3). Compared with NC, levels of 2-HB, glutamate and glycine were noticeably increased in the urine of HLD group, while concentration of pyridoxine, benzoate, hippurate, and 3-indoxylsulfate were significantly decreased.
Figure 5. OPLS-DA analysis of urinary metabolomic profiles of normal control rats (NC) and healthy rats administrated with HLJDD (HLD). (A) scores plots; (B) mean trajectory plots; (C) loading plots color-coded according to the absolute value of correlation coefficients. 2-OG, 2-oxoglutarate; Cr, creatine; PCr, phosphocreatine; (D) S-plots; (E) scatter plots of statistical validation obtained by permutation test, with R2 and Q2 values in the vertical axis, the correlation coefficients in the horizontal axis, and the ordinary least squares (OLS) line for the regression of R2 and Q2 on the correlation coefficients, intercepts: R2 = (0.0, 0.0921), Q2 = (0.0, −0.230); (F) receiver operating characteristic (ROC) curves of classifier performance of OPLS-DA models, the x-axis denoting the false positive rate, the y-axis the true positive rate.
Figure 6. Lattice plots of metabolites levels in urine for two model groups TAA (intoxicated with thioacetamide) and BDL (bile duct ligation), two HLJDD treatment groups THC (on TAA model) and BHD (on BDL model), and two control groups NC and HLD. The horizontal axis denoted the time points and the vertical axis denoted the integrated areas of the metabolites.
Table 3. Levels of integrated areas of metabolites in HLJDD administrated healthy rats (HLD) relative to normal control rats (NC) as depicted by z-scores (z-scores are colored for easy interpretation).
Treatment of HLJDD on Thioacetamide Induced Cholestatic Rats
NMR data of TAA vs. NC, and THC vs. TAA groups throughout the experiment were subjected to OPLS-DA analysis to unravel metabolic outcomes associated with cholestatic injuries induced by thioacetamide, and to find out metabolites that were directly related to the treatment effects of HLJDD on TAA model. TAA and NC groups showed good separation in score (Figure 7A) and trajectory plots (Figure 7B), with favorable R2 of 0.510 and Q2 of 0.450 (Figure S2a), and AUROC of 0.947 (Figure S3a), demonstrating a good distinction of TAA from NC group. THC was partly separated from TAA in the score (Figure 7D) and trajectory plots (Figure 7E), together with unfavorable values of Q2 (0.035, Figure S2b) and AUROC (0.763, Figure S3b), implying an unobvious metabolite change to TAA rats on the treatment of HLJDD. The loading plot (Figure 7C) revealed markedly increased levels of lactate/threonine, lysine/arginine and glutamate, and noticeably reduced levels of N-acetylglutamate, N-phenylacetyl-glycine, pyridoxine, phenylalanine and hippurate in TAA group, as compared with NC rats. Compared with TAA rats, THC rats exhibited increased levels of lysine/arginine, glutamate, phosphocholine, taurine, guanidoacetate, N-phenylacetyl-glycine, and Cr/PCr, and decreased levels of glutarate, 2-OG, citrate, benzoate, hippurate, 3-indoxylsulfate and trigonelline (Figure 7F).
Figure 7. OPLS-DA analysis of 1H NMR data in urine for normal control rats (NC) vs. cholestatic injury rats induced by thioacetamide (TAA) (A–C), and TAA rats vs. HLJDD treated TAA rats (THC) (D–F). (A,D) scores plots; (B,E) mean trajectory plots; (C,F) loading plots color-coded according to the absolute value of correlation coefficients. BCAAs: branched-chain amino acids, leucine, isoleucine and valine; 2-OG: 2-oxoglutarate; PC: Phosphocholine; TMAO: trimethylamine N-oxide; Cr: creatine; PCr: phosphocreatine.
However, most of these alterations couldn't support the beneficial effects of HLJDD on TAA model; neither did the univariate analysis of metabolites in THC group relative to NC group. To deduct the effect of HLJDD on metabolic status of healthy rats, metabolites of THC group were normalized by HLD group. Levels of metabolites in THC group relative to HLD group showed great difference to those in TAA group relative to NC group, especially those of glutamate, glycine, 2-HB, benzoate, pyridine and N-acetylglutamate. Variations of metabolites in TAA group were visualized relative to NC rats (Table S2). Throughout the whole experiment, significant high levels of allantoin (P < 0.01) were observed in TAA rats. Glutamate and 2-HB were observed increased in THC group relative to NC group (Table S3); however, their levels in THC group were decreased from T5 to T8 relative to HLD group (Table S4). N-acetylglutamate and hippurate were decreased in TAA rats from T2 to T7, which were significantly increased in THC group relative to HLD group, suggestive the protective effect of HLJDD on the damaged liver mitochondria.
Treatment of HLJDD on Bile Duct Ligated Cholestatic Rats
NMR data of BDL vs. NC, and BHD vs. BDL groups from T1 to T8 were subjected to OPLS-DA analysis to explore the metabolic feature of cholestatic injuries induced by bile duct ligation and the treatment effects of HLJDD on it. Clear separation was achieved between the two groups in both constructed OPLS-DA models in score (Figures 8A,D) and trajectory plots (Figures 8B,E). The favorable model parameters showcased a good performance and prediction ability of the constructed model (R2 = 0.640, Q2 = 0.590 in Figure S2c, and AUROC = 0.942 in Figure S3c for BDL vs. NC groups; R2 = 0.470, Q2 = 0.420 in Figure S2d, and AUROC = 0.986 in Figure S3d for BHD vs. BDL groups). Compared with NC rats, BDL rats exhibited significantly increased levels of BCAAs, fucose, lysine/arginine, glutamate, glycine, N-phenylacetyl-glycine and creatine/phosphocreatine (Cr/PCr), and decreased levels of glutarate, 2-OG, fumarate, hippurate and trigonelline (Figure 8C). BHD rats showed increased levels of lactate/threonine, lysine/arginine, acetate and fumarate, and decreased levels of 2-OG, citrate, 3-indoxylsulfate, phenylalanine, pyridoxine and hippurate as compared with the NC group (Figure 8F).
Figure 8. OPLS-DA analysis of 1H NMR data in urine for normal control rats (NC) vs. cholestic injury rats induced by bile duct ligated (BDL) (A–C), and BDL vs. HLJDD treated BDL rats (BHD) (D–F). (A,D) scores plots; (B,E) mean trajectory plots; (C,F) loading plots color-coded according to the absolute value of correlation coefficients. BCAAs, branched-chain amino acids, leucine, isoleucine and valine; 3-HB, 3-hydroxybutyrate; 2-OG, 2-oxoglutarate; PC, Phosphocholine; TMA, trimethylamine; Cr, creatine; PCr, phosphocreatine.
Color-coded table (Table S5) was used to visualize important variations of metabolites in BDL group relative to NC group. Glutamate and 2-HB showed noticeable sustained increase from T1 to T8 in BDL rats. Their levels were still increased in BHD group relative to NC group (Table S6), and were decreased relative to HLD group (Table S7). Glycine showed higher levels in BDL rats than that in NC rats throughout the experiment but gradually declined from T1 to T8. The high levels of glycine in BDL group from T4 to T5 were decreased in BHD group relative to HLD group (Figure 9). Trigonelline showed decreased level in BDL group throughout the whole experiments, which exhibited increased level in BHD group at T5. Similar to trigonelline, 2-OG, hippurate and fumarate were also decreased in BDL rats. Increased hippurate and decreased formate and N-phenylacetyl-glycine were found in BHD group relative to HLD group.
Figure 9. Noticeably altered metabolites in two cholestatic injury model groups TAA (intoxicated with thioacetamide) and BDL (bile duct ligation), two HLJDD treatment groups THC (on TAA model) and BHD (on BDL model), and HLJDD administrated healthy rats (HLD). NC, normal control group; FA, fatty acid; 2-HB, 2-hydroxybutyrate; GSH, glutathione; ROS, reactive oxygen species; BCAAs, branch-chained amino acids, leucine, isoleucine and valine; ATP, triphosphate adenosine; SAM, S-adenosyl-methionine; Kreb's cycle, tricarboxylic acid cycle; 2-OG, 2-oxoglutarate; NAG, N-acetylglutamate; NAPG, N-phenylacetyl-glycine. ↑↑/↓↓ and ↑/↓ stands for significantly and partly altered (increased or decreased), respectively.
Discussion
The extract was analyzed by UPLC-LTQ-Orbitrap-MS to successfully identify the major chemical component in HLJDD as phthalate, iridoid and its glycosides, flavonoids and alkaloids. Among the quantified compounds in HLJDD extract, baicalin, baicalein and berberine were exhibited to be the most abundant (Table S2), which was in accordance with previously reported (Xu et al., 2017; Zhang et al., 2017a). As the represent hepatic protective chemical components of Radix Scutellariae and HLJDD, baicalin was reported to ameliorate experimental cholestatic liver fibrosis in mice by modulation of oxidative stress, inflammation, and NRF2 transcription factor (Shen et al., 2017), attenuate non-alcoholic steatohepatitis by suppressing key regulators of lipid metabolism, inflammation and fibrosis in mice (Zhang et al., 2018), attenuate the liver injury induced by alcohol through modulation of hepatic oxidative stress, inflammation and sonic hedgehog pathway (Wang H. et al., 2016), and prevent against acetaminophen-induced liver injury via ERK signaling pathway (Liao et al., 2017). Moreover, it showed that baicalein and baicalin could alleviate acetaminophen-induced liver injury by activating Nrf2 antioxidative pathway (Shi et al., 2018). As the major and represent alkaloid of both Rhizoma Coptidis, Cortex Phellodendri Chinensis, and HLJDD, berberine could inhibit the liver fibrogenesis both in vivo and in vitro (Wang N. et al., 2016), alert hepatoprotective effect against methotrexate induced liver toxicity in rats (Mehrzadi et al., 2018), and promote male-specific expression of a bile acid uptake transporter to ameliorate liver cholestasis by inactivation of signal transducer and activator of transcription 5 signaling (Bu et al., 2017). HLJDD and its major representative components of baicalin, baicalein and berberine could protect against various hepatic injuries by ameliorating metabolic disturbances and improving hepatic function.
Excessive Oxidative Stress
The two acute cholestasis model shared common excessive fatty acid oxidation, as revealed by the notably elevated levels of lysine and 2-HB both in the urine of TAA and BDL rats. Lysine and 2-HB were indicative of oxidative stress induced by lipid peroxide (O'Doherty et al., 2014) and by excessive fatty acid oxidation (Gall et al., 2010), respectively. Reactive oxygen species (ROS) facilitated excessive fatty acid oxidation, as evidenced by the increased levels of 2-HB in TAA rats (Gall et al., 2010). As one of the representativ phenols and iridoid glycosides of HLJDD, chlorogenic acid and geniposide were reported to decrease the serum concentration of 2-HB (Ruan et al., 2014) and to confront glutamate-induced obese and exhibit antiobesity effects (Shi et al., 2011). The findings were consistent with the results of THC and BHD rats after administration with HLJDD in our study.
As a biomarker for oxidative stress, allantoin was originated from antioxidant uric acid in case of ROS (Yardim-Akaydin et al., 2004). Its elevated levels in TAA rats thus disclosed a rise in the production of endogenous ROS after thioacetamide intoxication, in consistent with previous report (Stanková et al., 2010). The levels of allantoin in the serum and urine were both increased in collagen-induced arthritis, which were decreased after supplementation with HLJDD and its combination formula with chemical constituents including 14% geniposide, 0.6% coptisine, 0.9% phellodendrine, 3% jatrorrhizine, 3.2% magnoflorine, 5.1% palmatine, 32% berberine, 26.1% baicalin, 2.6% chlorogenic acid, 1% crocin, 5.9% wogonoside, 3.6% baicalein, and 2.2% wogonin (Zhang et al., 2014).
Glutathione Depletion
Excessive oxidative stress in the two acute cholestasis model resulted in depletion of glutathione, as revealed by the notably elevated level of glutamate and disturbed contents of glycine, precursors of the most import natural antioxidant glutathione (GSH), implying the depletion of GSH (Choucha Snouber et al., 2013). Chlorogenic acid was reported to increase serum concentrations of glycine and hepatic concentrations of glutathione (Ruan et al., 2014). Therefore, comes naturally for the decrease of glutamate in BHD and THC rats after administration of HLJDD since that HLJDD is rich in chlorogenic acid.
Trigonelline is the methylation product of nicotinic acid in the demethylation process from S-adenosyl-methionine (SAM) to S-adenosyl-homocysteine (SAH) (Sun et al., 2008). S-adenosyl-homocysteine after hydrolysis afford homocysteine, which could be converted eventually into glutathione. Therefore, the increase of trigonelline in BHD group suggested an increased conversion from S-adenosyl-methionine to S-adenosyl-homocysteine, thus a regeneration of glutathione to regain the balance of redox system. Trigonelline was reported to protect the liver-kidney functions of diabetic rats efficiently (Hamden et al., 2013), thus increased trigonelline in BHD group suggested regained function of livers and kidneys of BDL rats, in consistent with the histopathological observation.
DNA Damage
ROS could also led to mitochondrial injury and even cell death of hepatocytes (Stanková et al., 2010). Bilirubin in dead hepatocytes entered blood, exhibiting symptoms of jaundice in TAA rats, as evidenced by the obviously consistent high levels of deoxycytidine, biomarker for jaundice syndrome (Wang X. et al., 2012), in TAA rats. Deoxycytidine is a composition of DNA, thus its continuous increase also suggested DNA damage in TAA rats, which has also been observed in the toxicity of thioacetamide S-oxides, metabolites of thioacetamide (Lozano et al., 2014).
Disturbed Urea Cycle
ROS induced injury of mitochondria in livers, led to the obstacle of urea cycle as evidenced by decreased level of N-acetylglutamate from T2 to T7 in TAA rats (Ishihara et al., 2006), which was increased after treatment with HLJDD. N-acetylglutamate is the most important cofactor in the control of ureagenesis in the mitochondria of liver (Tuchman et al., 2008), thus its increase in THC group relative to HLD group demonstrated that HLJDD could reverse the disturbed urea cycle in cholestasis TAA model. Synthesized from glutamate, the decrease of N-acetylglutamate has a direct association with the increased level of glutamate.
Disturbance of Gut Flora
As an antioxidant produced by the intestinal bacteria (Cornelli et al., 2010), the level of Pyridoxine (ca. VB6) was decreased in HLD group compared with NC rats, which suggested a dysfunction of gut flora that also evidenced by significantly decreased levels of benzoate, hippurate, and 3-indoxylsulfate (mammalian-microbial co-metabolites). Their decrease was consistent with oral administration of antibiotics of penicillin and streptomycin sulfate (Swann et al., 2011), therefore, comes naturally for their decrease in HLD rats since that HLJDD is rich in antibacterial berberine and other principles (Lewis and Ausubel, 2006). Therefore, it is strongly recommended to avoid taking any natural or artifact bacteriostatic agents including TCM such as HLJDD unless it is necessary.
Hippurate was generated in rat liver mitochondria by the conjugation of glycine with benzoate that was synthesized by bacterial action in gut. The decreased hippurate in BDL groups thus indicated disturbed gut flora, also supported by increased levels of two other gut microbial-host co-metabolites such as formate and N-phenylacetyl-glycine (Dong et al., 2012). The bile in cholestatic BDL rats couldn't flow to duodenum, resulting in alteration of microenvironment of gut and dysfunction of the gut flora. Gut microbiota was able to metabolize polyphenols, such as chlorogenic acid, into more absorbable compounds such as benzoic acid which was then detoxied to form hippurate (Claus et al., 2008). After administration of HLJDD, increased hippurate and decreased formate and N-phenylacetyl-glycine in BHD group relative to HLD group were observed, in accordance with previously reported, suggesting that HLJDD and its main principle chlorogenic acid could regain the unbalanced gut flora.
Dysfunction of Energy Metabolism
As the intermediates of Kreb's cycle, their significant decrease and slight decrease of citrate and succinate in BDL rats demonstrated a limited supply of substrates to the Kreb's cycle and thus an insufficient energy production. The energy crisis was even aggravated as evidenced by the increased level of BCAAs in BDL rats from T1 to T8. BCAAs play essential role in energy supply, especially for muscles and other organs by their oxidation (Choudry et al., 2006). The accumulation of BCAAs in BDL rats thus suggested an inhibited oxidation to compensate energy insufficiency. The increase of BCAAs in BDL group indicated renal injury, in consistent with those of nephrotoxicity study induced by D-serine (Williams and Lock, 2005), as evidenced by the observed severe pathological alterations in kidneys. The increased levels of fumarate and decreased levels of BCAAs in BHD group relative to HLD group suggested that HLJDD could regain the energy balance by regulating disturbed Kreb's cycle and inhibited BCAAs oxidation, in consistent with the protection of baicalin on liposaccharide induced damage (Liao et al., 2016b).
Conclusion
The two cholestatic injury models induced by thioacetamide and bile duct ligation, and the treatment effect of HLJDD were investigated and compared using an NMR-based metabolomic approach, complemented with histopathological inspection. The common metabolic features in both cholestatic models were excessive fatty acid oxidation, insufficient glutathione regeneration and disturbed gut flora. Model specific metabolic characteristics were inhibited urea cycle and DNA damage in TAA model, and perturbed Kreb's cycle and inhibited BCAAs oxidation in BDL model. Rich in various bioactive polyphenols, flavonoids, iridoid glycosides and alkaloids, such as chlorogenic acid, baicalin, geniposide and berberine, HLJDD could regain the balance of the disturbed metabolic status common in the two cholestasis injuries with good treatment effects, e.g. unbalanced redox system and disturbed gut flora; and perturbed urea cycle in TAA model, and disturbed Kreb's cycle and oxidation of BCAAs) in BDL model, respectively (Figure 9).
Author Contributions
D-DW made substantial contributions to the conception, design of experiment, data acquisition and analysis, interpretation of data, and agreement to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. J-SW made substantial contributions to the analysis of the NMR-based metabolomics data, drafting the work, revising it critically for important intellectual content, and agreement to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. J-AD and L-YK made substantial contributions to the approval of the version to be published, and agreement to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
Funding
This work was funded by the Key Project of National Natural Science Foundation of China (No. 81430092), the National Natural Science Foundation of China (No.81603273) and the Key program of Jiangsu Collaborative Innovation Center of Chinese Medicinal Resources Industrialization (ZDXM-3-15).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer SC and handling Editor declared their shared affiliation.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2018.00458/full#supplementary-material
Abbreviations
2-OG, 2-oxoglutarate; 2-HB, 2-hydroxybutyrate; AUROC, area under the receiver operating characteristic plot; BDL, bile duct ligation; BHD, bile duct ligation animals treated with Huang-Lian-Jie-Du-Decoction; BCAAs, branched-chain amino acids; Cr, creatine; GC-MS, gas chromatography-mass spectrometry; R2, goodness-of-fit; Q2, goodness of prediction; HE, hematoxylin & eosin; HLJDD, Huang-Lian-Jie-Du-Decoction; HMDB, Human Metabolome Database; LC-MS, liquid chromatography mass spectrometry; HLD, normal animals administrated with Huang-Lian-Jie-Du-Decoction; NC, normal control; NMR, nuclear magnetic resonance; OPLS-DA, orthogonal partial least squares-discriminant analysis; PCr, phosphocreatine; ROC, receiver operating characteristic; TAA, thioacetamide-intoxicated; THC, thioacetamide-intoxicated animals treated with Huang-Lian-Jie-Du-Decoction; TCM, Traditional Chinese medicine; UPLC, ultra-high pressure liquid chromatography.
References
An, Y. J., Xu, W. J., Jin, X., Wen, H., Kim, H., Lee, J., et al. (2012). Metabotyping of the C. elegans sir-2.1 mutant using in vivo Labeling and 13C-Heteronuclear multidimensional NMR Metabolomics. ACS Chem. Biol. 7, 2012–2018. doi: 10.1021/cb3004226
Barba, I., Fernandez-Montesinos, R., Garcia-Dorado, D., and Pozo, D. (2008). Alzheimer's disease beyond the genomic era: nuclear magnetic resonance (NMR) spectroscopy-based metabolomics. J. Cell. Mol. Med. 12, 1477–1485. doi: 10.1111/j.1582-4934.2008.00385.x
Beyoǧlu, D., and Idle, J. R. (2013). The metabolomic window into hepatobiliary disease. J. Hepatol. 59, 842–858. doi: 10.1016/j.jhep.2013.05.030
Bu, P., Le, Y., Zhang, Y., Zhang, Y., and Cheng, X. (2017). Berberine-induced Inactivation of signal transducer and activator of transcription 5 signaling promotes male-specific expression of a bile acid uptake transporter. J. Biol. Chem. 292, 4602–4613. doi: 10.1074/jbc.M116.757567
Chen, W.-Y., Lin, S.-Y., Pan, H.-C., Liao, S.-L., Chuang, Y.-H., Yen, Y.-J., et al. (2012). Beneficial effect of docosahexaenoic acid on cholestatic liver injury in rats. J. Nutr. Biochem. 23, 252–264. doi: 10.1016/j.jnutbio.2010.11.022
Chen, Y., Xian, Y., Lai, Z., Loo, S., Chan, W. Y., and Lin, Z.-X. (2016). Anti-inflammatory and anti-allergic effects and underlying mechanisms of Huang-Lian-Jie-Du extract: implication for atopic dermatitis treatment. J. Ethnopharmacol. 185, 41–52. doi: 10.1016/j.jep.2016.03.028
Choucha Snouber, L., Bunescu, A., Legallais, C., Brochot, C., Dumas, M. E., Elena-Herrmann, B., et al. (2013). Metabolomics-on-a-chip of hepatotoxicity induced by anticancer drug flutamide and its active metabolite hydroxyflutamide using HepG2/C3a microfluidic biochips. Toxicol. Sci. 132, 8–20. doi: 10.1093/toxsci/kfs230
Choudry, H. A., Pan, M., Karinch, A. M., and Souba, W. W. (2006). Branched-chain amino acid-enriched nutritional support in surgical and cancer patients. J. Nutr. 136, 314S–318S. doi: 10.1093/jn/136.1.314S
Claus, S. P., Tsang, T. M., Wang, Y., Cloarec, O., Skordi, E., Martin, F. P., et al. (2008). Systemic multicompartmental effects of the gut microbiome on mouse metabolic phenotypes. Mol. Syst. Biol. 4, 219. doi: 10.1038/msb.2008.56
Cornelli, U., Belcaro, G., Nardi, G. M., Cesarone, M. R., Dugall, M., Hosoi, M., et al. (2010). Action of an antioxidant complex on the antioxidant power of saliva. Panmin. Med. 52(2 Suppl. 1), 69–73.
Domitrović, R., Jakovac, H., and Blagojević, G. (2011). Hepatoprotective activity of berberine is mediated by inhibition of TNF-α, COX-2, and iNOS expression in CCl4-intoxicated mice. Toxicology 280, 33–43. doi: 10.1016/j.tox.2010.11.005
Dong, F., Wang, B., Zhang, L., Tang, H., Li, J., and Wang, Y. (2012). Metabolic response to klebsiella pneumoniae infection in an experimental rat model. PLoS ONE 7:51060. doi: 10.1371/journal.pone.0051060
Feng, Y., Siu, K.-Y., Ye, X., Wang, N., Yuen, M.-F., Leung, C.-H., et al. (2010). Hepatoprotective effects of berberine on carbon tetrachloride-induced acute hepatotoxicity in rats. Chin. Med. 5, 1–6. doi: 10.1186/1749-8546-5-33
Gall, W. E., Beebe, K., Lawton, K. A., Adam, K.-P., Mitchell, M. W., Nakhle, P. J., et al. (2010). α-Hydroxybutyrate is an early biomarker of insulin resistance and glucose intolerance in a nondiabetic population. PLoS ONE 5:10883. doi: 10.1371/journal.pone.0010883
Hamden, K., Mnafgui, K., Amri, Z., Aloulou, A., and Elfeki, A. (2013). Inhibition of key digestive enzymes related to diabetes and hyperlipidemia and protection of liver-kidney functions by trigonelline in diabetic rats. Sci. Pharm. 81, 233–246. doi: 10.3797/scipharm.1211-14
Hu, Y., Hu, Z., Wang, S., Dong, X., Xiao, C., Jiang, M., et al. (2013). Protective effects of Huang-Lian-Jie-Du-Tang and its component group on collagen-induced arthritis in rats. J. Ethnopharmacol. 150, 1137–1144. doi: 10.1016/j.jep.2013.10.038
Ishihara, K., Katsutani, N., and Aoki, T. (2006). A metabonomics study of the hepatotoxicants galactosamine, methylene dianiline and clofibrate in rats*. Basic Clin. Pharmacol. Toxicol. 99, 251–260. doi: 10.1111/j.1742-7843.2006.pto_455.x
Janbaz, K. H., and Gilani, A. H. (2000). Studies on preventive and curative effects of berberine on chemical-induced hepatotoxicity in rodents. Fitoterapia 71, 25–33. doi: 10.1016/S0367-326X(99)00098-2
Jang, S. I., Kim, H. J., Hwang, K. M., Jekal, S. J., Pae, H. O., Choi, B. M., et al. (2003). Hepatoprotective effect of baicalin, a major flavone from scutellaria radix, on acetaminophen-induced liver injury in mice. Immunopharmacol. Immunotoxicol. 25, 585–594. doi: 10.1081/IPH-120026443
Jiang, W.-Y. (2005). Therapeutic wisdom in traditional Chinese medicine: a perspective from modern science. Trends Pharmacol. Sci. 26, 558–563. doi: 10.1016/j.tips.2005.09.006
Kim, S.-J., Moon, Y.-J., and Lee, S.-M. (2010). Protective effects of baicalin against ischemia/reperfusion injury in rat liver. J. Nat. Prod. 73, 2003–2008. doi: 10.1021/np100389z
Lewis, K., and Ausubel, F. M. (2006). Prospects for plant-derived antibacterials. Nat. Biotech. 24, 1504–1507. doi: 10.1038/nbt1206-1504
Li, P., Liao, S.-T., Wang, J.-S., Zhang, Q., Xu, D.-Q., Lv, Y., et al. (2017). Protection by Huang-Lian-Jie-Du decoction and its constituent herbs of lipopolysaccharide-induced acute kidney injury. FEBS Open Bio 7, 221–236. doi: 10.1002/2211-5463.12178
Liao, C.-C., Day, Y.-J., Lee, H.-C., Liou, J.-T., Chou, A.-H., and Liu, F.-C. (2017). ERK signaling pathway plays a key role in baicalin protection against acetaminophen-induced liver injury. Am. J. Chin. Med. 45, 105–121. doi: 10.1142/S0192415X17500082
Liao, S., Li, P., Wang, J., Zhang, Q., Xu, D., Lv, Y., et al. (2016a). Huang-Lian-Jie-Du decoction treated sepsis via regulating ERK and SRC/STAT3 pathways and ameliorating metabolic status. RSC Adv. 6, 89855–89866. doi: 10.1039/c6ra17380b
Liao, S., Li, P., Wang, J., Zhang, Q., Xu, D., Yang, M., et al. (2016b). Protection of baicalin against lipopolysaccharide induced liver and kidney injuries based on 1 H NMR metabolomic profiling. Toxicology Research 5, 1148–1159. doi: 10.1039/C6TX00082G
Lindon, J. C., Holmes, E., and Nicholson, J. K. (2003). Peer reviewed: so what's the deal with metabonomics? Anal. Chem. 75, 384 A-391A. doi: 10.1021/ac031386+
Liu, L.-L., Gong, L.-L., Wang, H., Xiao, Y., Wu, X.-F., Zhang, Y.-H., et al. (2008). Baicalin inhibits macrophage activation by lipopolysaccharide and protects mice from endotoxin shock. Biochem. Pharmacol. 75, 914–922. doi: 10.1016/j.bcp.2007.10.009
Lourenço, A., Ascenso, J., and Sá-Correia, I. (2011). Metabolic insights into the yeast response to propionic acid based on high resolution 1H NMR spectroscopy. Metabolomics 7, 457–468. doi: 10.1007/s11306-010-0264-1
Lozano, E., Sanchez-Vicente, L., Monte, M. J., Herraez, E., Briz, O., Banales, J. M., et al. (2014). Cocarcinogenic effects of intrahepatic bile acid accumulation in cholangiocarcinoma development. Mol. Cancer Res. 12, 91–100. doi: 10.1158/1541-7786.mcr-13-0503
Lu, J., Wang, J.-S., and Kong, L.-Y. (2011). Anti-inflammatory effects of Huang-Lian-Jie-Du decoction, its two fractions and four typical compounds. J. Ethnopharmacol. 134, 911–918. doi: 10.1016/j.jep.2011.01.049
Lucio, M. (2009). Datamining Metabolomics: The Convergence Point of Non-target Approach and Statistical Investigation. Doctoral dissertation, Technische Universität München.
Lv, Y., Wang, J., Xu, D., Liao, S., Li, P., Zhang, Q., et al. (2017). Comparative study of single/combination use of Huang-Lian-Jie-Du decoction and berberine on their protection on sepsis induced acute liver injury by NMR metabolic profiling. J. Pharm. Biomed. Anal. 145, 794–804. doi: 10.1016/j.jpba.2017.07.062
Manna, S. K., Krausz, K. W., Bonzo, J. A., Idle, J. R., and Gonzalez, F. J. (2013). Metabolomics reveals aging-associated attenuation of noninvasive radiation biomarkers in mice: potential role of polyamine catabolism and incoherent DNA damage-repair. J. Proteome Res. 12, 2269–2281. doi: 10.1021/pr400161k
Mehrzadi, S., Fatemi, I., Esmaeilizadeh, M., Ghaznavi, H., Kalantar, H., and Goudarzi, M. (2018). Hepatoprotective effect of berberine against methotrexate induced liver toxicity in rats. Biomed. Pharmacother. 97, 233–239. doi: 10.1016/j.biopha.2017.10.113
Nicholson, J. K., Lindon, J. C., and Holmes, E. (1999). ‘Metabonomics’: understanding the metabolic responses of living systems to pathophysiological stimuli via multivariate statistical analysis of biological NMR spectroscopic data. Xenobiotica 29, 1181–1189.
O'Doherty, P. J., Lyons, V., Tun, N. M., Rogers, P. J., Bailey, T. D., and Wu, M. J. (2014). Transcriptomic and biochemical evidence for the role of lysine biosynthesis against linoleic acid hydroperoxide-induced stress in Saccharomyces cerevisiae. Free Radic. Res. 48, 1454–1461. doi: 10.3109/10715762.2014.961448
Qiu, J. (2007). Traditional medicine: a culture in the balance. Nature 448, 126–128. doi: 10.1038/448126a
Ruan, Z., Yang, Y., Zhou, Y., Wen, Y., Ding, S., Liu, G., et al. (2014). Metabolomic analysis of amino acid and energy metabolism in rats supplemented with chlorogenic acid. Amino Acids 46, 2219–2229. doi: 10.1007/s00726-014-1762-7
Shen, K., Feng, X., Pan, H., Zhang, F., Xie, H., and Zheng, S. (2017). Baicalin ameliorates experimental liver cholestasis in mice by modulation of oxidative stress, inflammation, and NRF2 transcription factor. Oxid. Med. Cell. Longev. 2017:6169128. doi: 10.1155/2017/6169128
Shi, L., Hao, Z., Zhang, S., Wei, M., Lu, B., Wang, Z., et al. (2018). Baicalein and baicalin alleviate acetaminophen-induced liver injury by activating Nrf2 antioxidative pathway: the involvement of ERK1/2 and PKC. Biochem. Pharmacol. 150, 9–23. doi: 10.1016/j.bcp.2018.01.026
Shi, L., Sun, J., and Li, L. (2011). Antiobesity effects of geniposide on monosodium glutamate-induced obese mice. China Pharmacy 22, 3671–3673.
Shukla, B., Visen, P. K. S., Patnaik, G. K., and Dhawan, B. N. (1992). Reversal of thioacetamide induced cholestasis by picroliv in rodents. Phytother. Res. 6, 53–55. doi: 10.1002/ptr.2650060114
Stanková, P., Kučera, O., Lotková, H., Roušar, T., Endlicher, R., and Cervinková, Z. (2010). The toxic effect of thioacetamide on rat liver in vitro. Toxicol. In vitro 24, 2097–2103. doi: 10.1016/j.tiv.2010.06.011
Sun, J., Schnackenberg, L. K., Holland, R. D., Schmitt, T. C., Cantor, G. H., Dragan, Y. R., et al. (2008). Metabonomics evaluation of urine from rats given acute and chronic doses of acetaminophen using NMR and UPLC/MS. J. Chromatogr. B Analyt. Technol. Biomed. Life Sci. 871, 328–340. doi: 10.1016/j.jchromb.2008.04.008
Swann, J. R., Tuohy, K. M., Lindfors, P., Brown, D. T., Gibson, G. R., Wilson, I. D., et al. (2011). Variation in antibiotic-induced microbial recolonization impacts on the host metabolic phenotypes of rats. J. Proteome Res. 10, 3590–3603. doi: 10.1021/pr200243t
Tuchman, M., Caldovic, L., Daikhin, Y., Horyn, O., Nissim, I., Nissim, I., et al. (2008). N-carbamylglutamate markedly enhances ureagenesis in N-acetylglutamate deficiency and propionic acidemia as measured by isotopic incorporation and blood biomarkers. Pediatr. Res. 64, 213–217. doi: 10.1203/PDR.0b013e318179454b
Vinaixa, M., Ángel Rodrí-guez, M., Rull, A., Beltrán, R., Bladé, C., Brezmes, J., et al. (2010). Metabolomic assessment of the effect of dietary cholesterol in the progressive development of fatty liver disease. J. Proteome Res. 9, 2527–2538. doi: 10.1021/pr901203w
Wan, J.-Y., Gong, X., Zhang, L., Li, H.-Z., Zhou, Y.-F., and Zhou, Q.-X. (2008). Protective effect of baicalin against Lipopolysaccharide/d-galactosamine-induced liver injury in mice by up-regulation of Heme oxygenase-1. Eur. J. Pharmacol. 587, 302–308. doi: 10.1016/j.ejphar.2008.02.081
Wang, H., Zhang, Y., Bai, R., Wang, M., and Du, S. (2016). Baicalin attenuates alcoholic liver injury through modulation of hepatic oxidative stress, inflammation and sonic hedgehog pathway in rats. Cell. Physiol. Biochem. 39, 1129–1140. doi: 10.1159/000447820
Wang, N., Xu, Q., Tan, H. Y., Hong, M., Li, S., Yuen, M.-F., et al. (2016). Berberine inhibition of fibrogenesis in a rat model of liver fibrosis and in hepatic stellate cells. Evid. Based Complement. Alternat. Med. 2016:8762345. doi: 10.1155/2016/8762345
Wang, P.-R., Wang, J.-S., Yang, M.-H., and Kong, L.-Y. (2014). Neuroprotective effects of Huang-Lian-Jie-Du-Decoction on ischemic stroke rats revealed by H-1 NMR metabolomics approach. J. Pharm. Biomed. Anal. 88, 106–116. doi: 10.1016/j.jpba.2013.08.025
Wang, P.-R., Wang, J.-S., Zhang, C., Song, X.-F., Tian, N., and Kong, L.-Y. (2013). Huang-Lian-Jie-Du-Decotion induced protective autophagy against the injury of cerebral ischemia/reperfusion via MAPK-mTOR signaling pathway. J. Ethnopharmacol. 149, 270–280. doi: 10.1016/j.jep.2013.06.035
Wang, R., Xiong, A.-Z., Teng, Z.-Q., Yang, Q.-W., Shi, Y.-H., and Yang, L. (2012). Radix paeoniae rubra and radix paeoniae alba attenuate CCl4-induced acute liver injury: an ultra-performance liquid chromatography-mass spectrometry (UPLC-MS) based metabolomic approach for the pharmacodynamic study of Traditional Chinese Medicines (TCMs). Int. J. Mol. Sci. 13, 14634–14647. doi: 10.3390/ijms131114634
Wang, X., Zhang, A., Han, Y., Wang, P., Sun, H., Song, G., et al. (2012). Urine metabolomics analysis for biomarker discovery and detection of jaundice syndrome in patients with liver disease. Mol. Cell. Proteom. 11, 370–380. doi: 10.1074/mcp.M111.016006
Wei, D., Liao, S., Wang, J., Yang, M., and Kong, L. (2015). Cholestatic liver injury model of bile duct ligation and the protection of Huang-Lian-Jie-Du decoction by NMR metabolomic profiling. RSC Adv. 5, 66200–66211. doi: 10.1039/c5ra12224d
Wei, D.-D., Wang, J.-S., Li, M.-H., Guo, P.-P., Dong, G., Yang, M.-H., et al. (2014a). A pilot study of the onset of hepatic encephalopathy (OHE) in mice induced by thioacetamide and the protective effect of taurine by holistic metabolic characterization. Metabolomics 11, 559–570. doi: 10.1007/s11306-014-0715-1
Wei, D.-D., Wang, J.-S., Wang, P.-R., Li, M.-H., Yang, M.-H., and Kong, L.-Y. (2014b). Toxic effects of chronic low-dose exposure of thioacetamide on rats based on NMR metabolic profiling. J. Pharm. Biomed. Anal. 98, 334–338. doi: 10.1016/j.jpba.2014.05.035
Wheelock, C. E., Goss, V. M., Balgoma, D., Nicholas, B., Brandsma, J., Skipp, P. J., et al. (2013). Application of ‘omics technologies to biomarker discovery in inflammatory lung diseases. Eur. Respir. J. 42, 802–825. doi: 10.1183/09031936.00078812
Whitfield, P. D., Noble, P.-J. M., Major, H., Beynon, R. J., Burrow, R., Freerman, A. I., et al. (2005). Metabolomics as a diagnostic tool for hepatology: validation in a naturally occurring canine model. Metabolomics 1, 215–225. doi: 10.1007/s11306-005-0001-3
Williams, R., and Lock, E. (2005). Sodium benzoate attenuates D-serine induced nephrotoxicity in the rat. Toxicology 207, 35–48. doi: 10.1016/j.tox.2004.08.008
Wishart, D. S. (2008). Quantitative metabolomics using NMR. TrAC Trends Anal. Chem. 27, 228–237. doi: 10.1016/j.trac.2007.12.001
Xu, D., Lv, Y., Wang, J., Yang, M., and Kong, L. (2017). Deciphering the mechanism of Huang-Lian-Jie-Du-Decoction on the treatment of sepsis by formula decomposition and metabolomics: enhancement of cholinergic pathways and inhibition of HMGB-1/TLR4/NF-kappa B signaling. Pharmacol. Res. 121, 94–113. doi: 10.1016/j.phrs.2017.04.016
Yardim-Akaydin, S., Sepici, A., Özkan, Y., Torun, M., Simsek, B., and Sepici, V. (2004). Oxidation of uric acid in rheumatoid arthritis: is allantoin a marker of oxidative stress? Free Radic. Res. 38, 623–628. doi: 10.1080/10715760410001694044
Ye, X., Feng, Y., Tong, Y., Ng, K.-M., Tsao, S., Lau, G. K. K., et al. (2009). Hepatoprotective effects of Coptidis rhizoma aqueous extract on carbon tetrachloride-induced acute liver hepatotoxicity in rats. J. Ethnopharmacol. 124, 130–136. doi: 10.1016/j.jep.2009.04.003
Zeng, Z., Liu, X., Dai, W., Yin, P., Zhou, L., Huang, Q., et al. (2014). Ion fusion of high-resolution LC MS-based metabolomics data to discover more reliable biomarkers. Anal. Chem. 86, 3793–3800. doi: 10.1021/ac500878x
Zhang, H., Fu, P., Ke, B., Wang, S., Li, M., Han, L., et al. (2014). Metabolomic analysis of biochemical changes in the plasma and urine of collagen-induced arthritis in rats after treatment with Huang-Lian-Jie-Du-Tang. J. Ethnopharmacol. 154, 55–64. doi: 10.1016/j.jep.2014.03.007
Zhang, J., Zhang, H., Deng, X., Zhang, N., Liu, B., Xin, S., et al. (2018). Baicalin attenuates non-alcoholic steatohepatitis by suppressing key regulators of lipid metabolism, inflammation and fibrosis in mice. Life Sci. 192, 46–54. doi: 10.1016/j.lfs.2017.11.027
Zhang, Q., Fu, X., Wang, J., Yang, M., and Kong, L. (2017a). Treatment effects of ischemic stroke by berberine, baicalin, and jasminoidin from Huang-Lian-Jie-Du-Decoction (HLJDD) explored by an integrated metabolomics approach. Oxid. Med. Cell. Longev. 2017:9848594. doi: 10.1155/2017/9848594
Zhang, Q., Guo, P., Wang, J., Yang, M., and Kong, L. (2015). Gender-specific metabolic responses in focal cerebral ischemia of rats and Huang-Lian-Jie-Du decoction treatment. RSC Adv. 5, 95558–95575. doi: 10.1039/c5ra19934d
Zhang, Q., Wang, J., Liao, S., Li, P., Xu, D., Lv, Y., et al. (2017b). Optimization of Huang-Lian-Jie-Du-Decoction for ischemic stroke treatment and mechanistic study by metabolomic profiling and network analysis. Front. Pharmacol. 8:165. doi: 10.3389/fphar.2017.00165
Zhang, Q., Wang, J., Zhang, C., Liao, S., Li, P., Xu, D., et al. (2016). The components of Huang-Lian-Jie-Du-Decoction act synergistically to exert protective effects in a rat ischemic stroke model. Oncotarget 7, 80872–80887. doi: 10.18632/oncotarget.12645
Zhao, X., Zhang, J.-J., Wang, X., Bu, X.-Y., Lou, Y.-Q., and Zhang, G.-L. (2008). Effect of berberine on hepatocyte proliferation, inducible nitric oxide synthase expression, cytochrome P450 2E1 and 1A2 activities in diethylnitrosamine- and phenobarbital-treated rats. Biomed. Pharmacother. 62, 567–572. doi: 10.1016/j.biopha.2007.02.009
Zhao, Y., Li, H., Gao, Z., and Xu, H. (2005). Effects of dietary baicalin supplementation on iron overload-induced mouse liver oxidative injury. Eur. J. Pharmacol. 509, 195–200. doi: 10.1016/j.ejphar.2004.11.060
Zhu, B., Cao, H., Sun, L., Li, B., Guo, L., Duan, J., et al. (2018). Metabolomics-based mechanisms exploration of Huang-Lian Jie-Du decoction on cerebral ischemia via UPLC-Q-TOF/MS analysis on rat serum. J. Ethnopharmacol. 216, 147–156. doi: 10.1016/j.jep.2018.01.015
Keywords: metabolomics, Huang-Lian-Jie-Du-Decoction, acute cholestatic injuries, thioacetamide, bile duct ligation, orthogonal partial least squares-discriminant analysis
Citation: Wei D-D, Wang J-S, Duan J-A and Kong L-Y (2018) Metabolomic Assessment of Acute Cholestatic Injuries Induced by Thioacetamide and by Bile Duct Ligation, and the Protective Effects of Huang-Lian-Jie-Du-Decoction. Front. Pharmacol. 9:458. doi: 10.3389/fphar.2018.00458
Received: 31 January 2018; Accepted: 18 April 2018;
Published: 08 May 2018.
Edited by:
Adolfo Andrade-Cetto, Universidad Nacional Autónoma de México, MexicoReviewed by:
Sol Cristians, Universidad Nacional Autónoma de México, MexicoJianxin Chen, Beijing University of Chinese Medicine, China
Copyright © 2018 Wei, Wang, Duan and Kong. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jin-Ao Duan, dja@njutcm.edu.cn
Ling-Yi Kong, cpu_lykong@126.com
†These authors have contributed equally to this work.