Xinxin Dai
Shulan Su*
Hongdie CaiDandan WeiHui YanTianyao ZhengZhenhua ZhuEr-xin ShangSheng GuoDawei Qian
Jin-ao Duan*- Jiangsu Collaborative Innovation Center of Chinese Medicinal Resources Industrialization, State Key Laboratory Cultivation Base for TCM Quality and Efficacy, Nanjing University of Chinese Medicine, Nanjing, China
Rehmannia glutinosa Libosch (RG), is officially listed in the Chinese Pharmacopoeia and is widely used in China. The leaves of RG (LR) is an important vegetative organ of the plant. At present, the total glycosides of RG (TLR) were extracted from RG, and developed a national second class of new drugs to the Dihuangye total glycoside capsule (DTG). Additionally, DTG has the effect of nourishing yin and tonifying kidney, promoting blood circulation and blood cooling, and applicable to chronic glomerulonephritis mild to Qi and Yin Deficiency. Moreover, diabetic nephropathy (DN) rats model was induced by intraperitoneal injection of a small dose of streptozotocin (45 mg/kg) and high-fat diet and plus 5% glucose drinking water. Over 15 days, after oral administration TLR and DTG in DN rats, samples from serum, urine and kidney were collected for biochemical indicators measurements, pathological analysis, western blotting and metabolomics. Therefore, the analytical results of biochemical indicators, histopathological observations and western blotting showed that TLR and DTG exhibited a significant effect in renal protection. And 27 endogenous metabolites (12 in serum and 15 in urine) could be tentatively identified in the process of DN in rats using metabolomics method. Those endogenous metabolites were chiefly involved in sphingolipid metabolism; pentose, glucuronate interconversion; terpenoid backbone biosynthesis; purine metabolism and retinol metabolism. After drug intervention, these endogenous metabolites turned back to normal level some extent (P < 0.05). Furthermore, TLR and DTG prevent high glucose-induced glomerular mesangial cells (GMCs) by inhibiting TGF-β1 and Wnt/β-catenin signaling pathway, providing a powerful supports to develop a new therapeutic agent for DN. This study paved the way for further exploration of the pathogenesis of DN, early diagnosis and the evaluation of curative effect.
Introduction
Traditional Chinese medicine (TCM) with multi-components and multi-targets has been used in clinical for over 4000 years in China due to its efficacy and low side-effects (Yokozawa et al., 1986). RG, belonging to Scrophulariaceae family, is widely distributed in China, and its root tubers are commonly called “Dihuang" in China. So far large quantity of iridoids and phenylethanoid glycosides have been isolated from the root and leaves of RG, together with other types of active principles. “Dihuang" has the significant effects on cardiovascular diseases (Liu et al., 2011), neuroprotection (Cai et al., 2011), diabetes and its complications (Zhang et al., 2014), osteoporosis (Zhang Z. et al., 2013), and hyperlipidemia (Zhang R. et al., 2013). Currently, with the deep exploration of RG, the leaves of RG (LR) is gaining increased attention around the world. In order to find the possibility to be used as medicine, people are beginning to pay attention to the development and utilization of LR. At present, LR is officially listed in the Beijing standard of Chinese herbal medicines, and it has the effect of heat-clearing, promoting blood circulation, supplementing Qi and nourishing yin and tonifying kidney (Beijing Municipal Health Bureau, 1988). The total glycosides of RG (TLR) were extracted from RG. The major active components of DTG is the phenylethanoid glycosides, which is extracted from LR, and phenylethanoid glycosides are a class of natural glycosides containing hydroxyl groups, methoxy substituted phenethyl or cinnamoyl groups, usually containing β-glucose as the parent nucleus, and are widely present in dicotyledonous plants (Jiménez and Riguera, 1994). DTG can increased glomerular permeability and reduced glomerular hyperfiltration, thus reducing proteinuria and the protection of renal function (Shen et al., 2010).
Diabetic nephropathy (DN) is one of the most common microangiopathy in diabetic patients, and renal fibrosis is the main cause of death in patients with DN (Lv et al., 2015). At present, the clinical diagnosis of DN mainly depending on the urinary albumin excretion rate detection, however, urinary microalbumin content may be affected by many factors such as obesity, insulin resistance and there are limitations (Weir, 2007). Therefore, it is necessary to find an early diagnostic marker with high sensitivity and specificity. To the best of our knowledge, there is no unified theory about the pathogenesis of DN. The current research focuses on the aspects of glucose metabolism, lipid metabolism, oxidative stress (Oberg et al., 2004; Chen and Quilley, 2008; Balakumar et al., 2009), and so on. Studies have shown that Wnt/β-catenin signaling plays a prominent role in cell differentiation, adhesion, survival, and apoptosis and is involved in glomerular cell proliferation and renal fibrosis, which affected the occurrence and development of DN (Xiao et al., 2013).
Metabolomics, based on the dynamic changes of endogenous metabolites in organisms, revealing the overall physiological status in responding to pathophysiological stimuli or genetic, environmental, or lifestyle factors (Pereira and Chang, 2004), and it has been widely used in the diagnosis of diseases, marker screening and pathogenesis research as a kind of technical means to study the dynamic metabolism of biological organism or tissue cells in recent years (Odunsi et al., 2005; Hodavance et al., 2007; Hwang et al., 2010). Metabolomics is regarded as “biochemical phenotypes” of the whole functional state of organism, and can be real-time, sensitive and real to express the response and regulation of the overall functional state of the organism under various external factors (Nicholson and Lindon, 2008).
Therefore, this study parallel simultaneous analysis of DN rats and normal rats of serum and urine, markers of different metabolites, and to study the effect of DTG and TLR on DN by analyzing the regulation of different metabolites. Moreover, the high glucose-induced GMCs model in vitro evaluate the efficacy of DTG and TLR on DN by determination of TGF-β1, Wnt4 and β-catenin proteins expression levels.
Materials and Methods
Chemicals and Instruments
UPLC-grade acetonitrile was purchased from Merck (Darmstadt, Germany), formic acid and STZ were purchased from Sigma-Aldrich (Sigma, St. Louis, MO, United States). HPLC-grade acteoside was purchased from the National Institutes for Food and Drug Control (Beijing, China) and the purity is above 98%. BUN reagent kit, LDL reagent kit, T-CHO reagent kit, TG reagent kit, Scr reagent kit, UP reagent kit, serum β2-microglobulin (β2-MG) reagent kit were bought from Nanjing Jiancheng Bioengineering Institute (Nanjing, China). Huangkui capsule was purchased from SZYY Group Pharmaceutical Limited; Irbesartan was purchased from Shenzhen Haibin Pharmaceutical Co., Ltd. DMEM and F12 were purchased from GIBCO, America; TGF-β1, Wnt4 and β-catenin first antibody were purchased from Abcam (Cambridge, United Kingdom). All other chemicals and reagents used in this study were of analytical grade and made in China.
Waters AcquityTM Ultra Performance LC system (Waters, United States) equipped with a Quattro Micro MS spectrometer and a Waters Xevo TM G2 QTof MS (Waters MS Technologies, Manchester, NH, United States). Deionized water was purified on a Milli-Q system (Millipore, Bedford, MA, United States). Mass Lynx v4.1 workstation was adopted to analyze the data, and Ultra-high speed centrifuge at low temperature (Thermo Scientific, United Kingdom); DMI3000M microscope (Leica, Germany) were used.
Preparation of TLR and DTG
The plant material, LR were purchased from Henan farmers market, and identified by the Prof. Jin-Ao Duan (Department of Nanjing University of Chinese Medicine). Fresh LR were vacuum-dried in 80°C and ground into powder. The vacuum-dried LR (500 g) was extracted with 6 L 80% alcohol three times and 2 h each time with reflux extraction. In rotary evaporation instrument reduced pressure concentration, as TLR low dose oral solution. Concentrated on dose in the twice as high dose oral solution.
Before the experiment, the chemical of TLR and DTG were previously established by the UPLC-TQ-MS. And the main components of TLR identified mainly as catalpol, ajugol, and acteoside, the contents of above components were 0.6326%, 0.4105%, and 0.6833%, respectively. DTG were bought from Sichuan Meidakang Pharmaceutical Co., Ltd., the main component is acteoside, and the content of the acteoside in the DTG was 13.61%. The MRM chromatogram and structure of catalpol, ajugol and acteoside were presented in Figure 1. Moreover, we made the content of acteoside in TLR oral solution and DTG oral solution is consistent for rats.
FIGURE 1. MRM chromatogram and structure of catalpol, ajugol, and acteoside.
Experimental Animals
All experiments were conducted with male Sprague–Dawley (SD) rats weighing 200–240 g purchased from Experimental animal center of Zhejiang Province (Zhejiang, China; Certificate no. SCXK 2014-0001). Animals were housed in cages with a constant humidity (ca. 60% ± 2%) and temperature (ca. 23 ± 2°C) and with a light/dark cycle of 12 h. The animals were underwent an adaptation period of 7 days, during which they were allowed unlimited access to chow and tap water.
After the 7-day acclimation period, animals were randomly divided into eight groups with 6 in each group: control group (normal saline, C), model group (normal saline, M), Huangkui capsule group (0.75 g/kg⋅d, HK), Irbesartan group (27 mg/kg⋅d, YX), TLR low dose group (4.3 g/kg⋅d, DHYL), TLR high dose group (7.2 g/kg⋅d, DHYH), DTG low dose group (216 mg/kg⋅d, JNL), DTG high dose group (360 mg/kg⋅d, JNH). Expect for control group rats, the rats in other groups were fed with high-fat diet (59% basic mice feed, 20% sugar, 18% lard, and 3% egg yolk) and 5% glucose water for 21 days, then, the rats were injected intraperitoneally with STZ (45 mg/kg) dissolved in 0.1 mol/L sodium citrate buffer (pH 4.5) to induce DN model. 72 h after the injection of STZ, rats with FBG levels of 16.7 mmol/L and above were considered to be in diabetic state and were used for the study. Moreover, in this study, Huangkui capsule and Irbesartan were selected as positive control drug. The animal dose of Huangkui capsule and Irbesartan were extrapolated from the human daily dose, using the body surface area normalization method. The formula for dose translation was as follows: human clinical equivalent dose of medicine × 70 × 0.018/0.2. The dose of JNL and JNH were three times and five times of the clinical equivalent dose, respectively. Moreover, TLR are about 20 times the corresponding dose of DTG. Then, the rats were given corresponding drug once a day by gastric gavage for 15 days, and the control group rats were given the same volume of saline at the same time. This study was carried out in accordance with the recommendations of the Institutional Animal Ethics Committee of Nanjing University of Chinese Medicine (Nanjing, China). The protocol was approved by the Institutional Animal Ethics Committee of Nanjing University of Chinese Medicine.
Sample Collection
During the experimental period, body weights were recorded at 3-days intervals, FBG levels were measured by a One Touch Ultra II blood glucose monitoring system (Life Scan, Milpitas, CA, United States) by tail vein at 3-days intervals (at 8 a.m). Prior to the experiments, rats were put individually in the metabolism cages for 24-h urinary collection, all urine samples were immediately centrifugated at 3,000 rpm for 10 min after collection, and the supernatants were separated and stored at -80°C until analysis. At the end of study, rats were sacrificed under anesthesia by intraperitoneal injection of chloral hydrate (330 mg/kg body weight), and blood was collected for concentration of biochemical parameters and metabolomics study. Then, the blood samples were centrifugated at 3,000 rpm for 10 min, and the serum samples were separated and stored at -80°C until analysis. The right kidney was removed with one piece fixed with 10% neutral-buffered formalin for histologic examination. Renal cortex was isolated from the other piece and frozen in liquid nitrogen for western blotting analysis.
Biochemical Indicators Measurements, Pathological Analysis and Western Blotting
The therapeutic efficacy of TLR and DTG on DN rats was evaluated for the levels of UP in urine and FBG, LDL, T-CHO, TG, Scr, BUN and β2-MG in serum, which were detected followed by the description supplied by kits manufacturer. Part of the renal cortex were performed for HE staining in order to observe pathological changes in renal tissue, degrees of fibrosis tissue hyperplasia, structures of glomeruli and tubules through electron microscope (200×), and another part of the renal tissues samples were analyzed by western blotting for detecting the expression levels of Wnt4, β-catenin and TGF-β1 protein.
GMC Cells Culture
Glomerular mesangial cells were provided by Nanjing KeyGen Biotechnology Development Co., Ltd. The complete medium was low glucose DMEM with 10% FBS and 1% penicillin-streptomycin solution. The GMCs were cultured in 37°C, 5% CO2 and saturated humidity. The GMCs were plated on 96-well plates, pretreated with DMEM low glucose (5 mmol⋅L-1), incubated at 37°C and 5% CO2 for 6 h, and then were synchronized. The experimental groups were divided as follows: control group (C), which without intervention factor; model group (M), which cells were cultured in high glucose (30 mmol⋅L-1) DMEM solution; Acteoside (MRHTG) group (acteoside is stable in 0.1% DMSO in 24 h by analyzing its contents), TLR group (DHY) and DTG group (JN) at five dosages (5, 10, 25, 50, and 100 μmol⋅L-1), respectively. Cellular count was used to detect the number of cells in each group after incubation for 48 h, and then determined the total protein content of each group by BCA kit. The ratio of total cell protein/cell number was calculated and expressed as the total amount of protein (μg) contained in 103 cells as an indicator of GMCs hypertrophy. Moreover, these cells were prepared for the western blotting to examine the expression levels of TGF-β1, Wnt4 and β-catenin protein.
Sample Preparation and UPLC-QTOF/MS Analysis
All serum and urine samples were thawed at room temperature before preparation. The 300 μL acetonitrile was added to 100 μL serum samples to precipitated protein, and vigorously mixed for 60 s. The 100 μL acetonitrile was added to 100 μL urine samples to remove proteins and vigorously mixed for 30 s. All the samples were centrifuged at 13,000 rpm for 10 min at 4°C. Finally, the supernatant was injected to UPLC-QTOF/MS analysis.
Metabolites separation was performed using a Waters AcquityTM Ultra Performance LC system (Waters, United States) equipped with a Waters XevoTM G2 Q/TOF-MS (Waters MS Technologies, Manchester, NH, United States). An aliquot of 2 μL of sample solution was injected on an Acquity UPLC BEH C18 (100 mm × 2.1 mm, 1.7 μm, Waters Corporation, Milford, CT, United States) at 35°C and the flowrate was 0.4 mL/min. The optimal mobile phase consisted of water (A) (containing 0.1% formic acid) and acetonitrile (B). The optimized UPLC elution conditions for serum were: 0 ∼ 3 min, 95% ∼ 55% A; 4 ∼ 13 min, 55% ∼ 5% A; 13 ∼ 14 min, 5% A. The optimized UPLC elution conditions for urine were: 0 ∼ 8 min, 95% ∼ 70% A; 8 ∼ 11 min, 70% ∼ 30% A; 11 ∼ 13 min, 30% ∼ 5% A; 13 ∼ 14 min, 5% A. The autosampler was maintained at 4°C.
Mass spectrometry was performed using a XevoTM G2 QTof (Waters MS Technologies, Manchester, NH, United States), a quadrupole and orthogonal acceleration time-of- flight tandem mass spectrometer. Leucine-enkephalin was used as the lock mass generating an [M+H]+ ion (m/z 556.2771) and [M-H]- ion (m/z 554.2615) in positive and negative modes, respectively. The concentration of Leucine-enkephalin 200 pg/mL and the infusion flow rate was 100 μL/min to ensure accuracy during the MS analysis via a syringe pump. Data were collected in centroid mode from 100 to 1000 m/z. For both positive and negative electrospray modes, the capillary and cone voltage were set at 3.0 kV and 30 V, respectively. The desolvation gas was set to 600 L/h at a temperature of 350°C, the cone gas was set to 50 L/h and the source temperature was set to 120°C. The data acquisition rate was set to 30 ms, with a 0.02 s interscan delay.
In addition, 10 serum (or urine) samples were randomly selected from each group and mixed together as the quality control (QC) samples, respectively. The QC sample was used to optimize the condition of UPLC-QTOF/MS, as it contained most information of whole samples. The QC samples were injected six times at the beginning of the running in order to condition or equilibrate the system and then every ten samples to further monitor the stability of the analysis. Every day, after the instrument was calibrated, the QC sample was firstly analyzed to test the stability of the instrument in order to ensure consistent performance of the system.
Metabolomics Data Processing and Analysis
All of the data acquisition and analyses of data were controlled by Waters MassLynx v4.1 software. The multivariate data matrix was analyzed by EZinfo software 2.0 (Waters Corp., Milford, CT, United States). The main parameters include: retention time range 0 ∼ 14 min; mass ratio m/z 100 ∼ 1000, mass tolerance range 0.01 Da, peak intensity threshold 50, quality window 0.05 Da, retention time windows 0.20 min, automatic detection of 5% peak height and noise. The intensity of each ion was normalized with respect to the total ion count to generate a data matrix that consisted of the retention time, m/z value, and the normalized peak area.
The resultant data matrices were introduced to EZinfo 2.0 software for principal component analysis (PCA), partial least-squares-discriminant analysis (PLS-DA) and orthogonal projection to latent structures (OPLS-DA) analysis. From the OPLS-DA, various metabolites could be identified as being responsible for the separation between control group and model group and were therefore viewed as potential markers. Potential markers of interest were extracted from S-plots constructed following analysis with OPLS-DA, and variables that had significant contributions to discrimination between groups were subjected to further identification of the molecular formula.
The variable importance (VIP) in the projection value is a weighted sum of squares of the PLS weights, and the variables with VIP >1 were considered to be influential for the separation of samples in the score plots generated from OPLS-DA analysis. In all experiments, confidence level was set at 95% to determine the significance of difference (P < 0.05). Those variables were eventually selected as potential biomarkers. The PLS-DA score plots were described by the cross-validation parameter R2Y and Q2, which represents the total explained variation for the X matrix and the predictability of the model, respectively. Excellent models are obtained when the cumulative values of R2Y and Q2 are above 0.8. The relative distances between administration groups and control group from PLS-DA score plot was calculated with the average value (x-axis and y-axis) of all samples of the control group as the referenced point (Li et al., 2010). This quantitative value was used as an indicator of pharmacodynamic evaluation of metabolomics, and it solves the problem of lack of accurate and quantitative evaluation methods for many pharmacodynamics.
Metabolites Identification and Construction of Metabolic Pathway
Potential metabolites of DN selected were identified according to the determination of the accurate m/z, retention time, and typical MS/MS fragment and pattern of the potential biomarkers above, which were obtained in the positive and negative ion modes. The construction of metabolic pathway was performed with Metabo Analyst, which is a web-based tool for visualization of metabolomics1 based on database source including the KEGG2 and the HMDB3 databases for searching. The retention time and typical MS/MS fragment and pattern had great avail to narrow the range of possible molecules.
Results
Determination of Biochemical Indicators, Histopathological Observations and Western Blotting
The analytical results of biochemical indicators were presented in Figure 2A. After 15 days of administration, the level of FBG, LDL, T-CHO, TG, Scr, BUN, β2-MG in serum, and UP in urine of DN rats increased significantly compared with that of control rats (P < 0.05 or P < 0.01 or P < 0.001). After treatment of HK (0.75 g/kg⋅d), YX (27 mg/kg⋅d), DHYL (4.3 g/kg⋅d), DHYH (7.2 g/kg⋅d), JNL (216 mg/kg⋅d), JNH (360 mg/kg⋅d), the levels of these indicators were nearly restored to normal. The results of the histological examination were consistent with the biochemical analysis (see Figure 2B). Compared with the control group, the DN model group showed renal tubular epithelial vacuolar degeneration, interstitial inflammatory cell infiltration, glomerular and tubular atrophy and interstitial fibrosis. Western blotting analysis (Figure 2C) showed the expression of Wnt4, β-catenin, TGF-β1 in renal tissue of control rats were significant from that in DN rats (P < 0.05). Compared with DN model group, the levels of creatinine, urea nitrogen, 24-h urinary protein and renal histopathology were significantly improved in the administration groups, and the expression levels of Wnt4, β-catenin and TGF-β1 protein was significantly decreased, especially TLR treatment.
FIGURE 2. Determination of biochemical indicators and analysis of renal pathological section and Western blotting among control group, model group and administered group. (A) Determination of FBG, LDL, T-CHO, TG, Scr, β2-MG in serum, and BUN, UP in urine. (B) HE staining of kidney biopsy (200×). (C) Western blotting analysis of renal tissues of control group, HK group, YX group, DHYL group, DHYH group, JNL group and JNH group. #P < 0.05, ##P < 0.01, ###P < 0.001: DN model group vs. control group; ∗P < 0.05, ∗∗P < 0.01, ∗∗∗P < 0.001: after treatment by HK, YX, DHYL, DHYH, JNL and JNH vs. DN model group. Data were presented as the mean ± SE.
Effect of TLR and DTG on GMCs Hypertrophy
As shown in Figure 3, high glucose increased the ratio of total protein to total cells in GMCs, resulting in cell hypertrophy, whereas TLR, DTG and the major active ingredient acteoside (5, 10, 25, 50, and 100 μmol⋅L-1), the ratio of total protein to total cell number decreased to a certain extent in a dose-dependent manner.
FIGURE 3. Effects of DHY, JN, and MRHTG on hypertrophy of high glucose-induced GMCs. ###P < 0.001: model group vs. control group; ∗∗P < 0.01, ∗∗∗P < 0.001: after treatment by DHY, JN, and MRHTG vs. model group. Data were presented as the mean ± SE.
TGF-β1 and Wnt/β-Catenin Signaling Pathway on GMCs Treated With TLR and DTG
As shown in Figure 4, the protein expression levels of Wnt4, β-catenin and TGF-β1 were significantly increased in high glucose-induced GMCs (P < 0.05). The expression of Wnt4, β-catenin and TGF-β1 returned to the normal level after intervention of TLR and DTG, but there was no significant effect of the MRHTG group on the expression levels of Wnt4, β-catenin and TGF-β1.
FIGURE 4. Effects of DHY and JN on the expression of Wnt-4, β-catenin and TGF-β1 in GMCs. #P < 0.05: model group vs. control group; ∗P < 0.05: after treatment by DHY, JN, and MRHTG vs. model group. Data were presented as the mean ± SE.
Metabolomics Insights
QC samples Analysis
The relatively clustering of QC samples (Figure 5A) and relative standard deviations (RSD%) of ion intensity (Table 1) expound and prove the quality of QC data. The trend plot showing the variation of t[1] over all observations (Figure 5B). The extracted ion chromatographic peaks of 10 ions were selected for method validation. The repeatability of method was evaluated by using six replicates of QC sample. This type of results demonstrated that the method had excellent repeatability and stability.
FIGURE 5. Assessment of QC samples (A) PCA score plot (PC1 vs. PC2) of test samples and QC samples. (B) Trend plot showing the variation of t[1] over all observations. QC samples were colored as red boxes and test samples were colored as black triangle. X-axis numbers represented sample number (54 injections). Y-axis was arbitrary (3 SD).
TABLE 1. Coefficient of variation of ion intensity of selected ions present in the QC samples covering the range of retention times.
Multivariate Data Analysis
All the data containing the retention time, peak intensity and exact mass were imported in the MassLynxTM software for multiple statistical analyses. The supervised OPLS-DA model, a pattern recognition approach, were established to separate serum or urine samples into two blocks between DN model group and control group. The supervised OPLS-DA with 100% sensitivity and no less than 95% specificity using a leave one out algorithm showed a better discrimination between the two groups (Figures 6A1–A4), which demonstrated that the DN model was built successfully. Based on these results, the OPLS-DA score plot (Figures 6B1–B4) were used to look for potential markers associated with DN progress. The R2Y and Q2 of PLS-DA model in positive and negative modes for serum and urine samples were suggested that the PLS-DA model was good to fitness and prediction.
FIGURE 6. OPLS-DA scores plots (A), S-plot of OPLS-DA (B) for serum (A1,A2,B1,B2) and urine (A3,A4,B3,B4) samples of DN model group (red) vs. control group (black) in positive (A1,B1,A3,B3) and negative (A2,B2,A4,B4) ion mode. (A1, R2Y = 0.9920, Q2 = 0.8994; A2, R2Y = 0.9882, Q2 = 0.9269; A3, R2Y = 0.9932, Q2 = 0.9363; A4, R2Y = 0.9933, Q2 = 0.9513).
The UPLC-QTOF/MS analysis platform provided the retention time and precise molecular mass within measurement errors (<5 ppm) as well as the fragments of corresponding production for the structural identification of metabolites. According to the precise molecular mass, predicted elemental composition was predicted and potential molecular formula could be searched out in Human Metabolome Database4. Twenty-seven endogenous metabolites (12 in serum and 15 in urine) were tentatively identified by comparing with authentic standards or based on the protocol detailed above method. The information about the detected endogenous metabolites was summarized in Table 2.
TABLE 2. Potential metabolites selected and identified between DN model group and control group.
Changes of Relative Intensity of Endogenous Metabolites
In order to study the efficacy and action mechanism of TLR and DTG for treating DN disease, PLS-DA model analysis was built to obtain the changes during the control group, DN model group and administration group rats (Figure 7). The variations of metabolic profiling in serum and urine for administration group rats had the tendency to restore back to the levels of controls, especially TLR and DTG groups (Table 3). Furthermore, the relative quantities of 27 endogenous metabolites (12 in serum and 15 in urine) were ultimately identified by comparing with authentic standards or based on the protocol detailed above method. Twenty-four endogenous metabolites (except for taurochenodesoxycholic acid, chenodeoxycholic acid glycine conjugate and L-gulonolactone) in serum and urine significantly affected by TLR and DTG were restored back to a control-like level. The detailed information was showed in Figure 8.
FIGURE 7. PLS-DA scores plots for serum (S1,S2) and urine (U1,U2) samples from DN model rats and controls and administration group rats in positive and negative ion mode.
TABLE 3. The relative distance between treatment groups and control group from the PLS-DA score plot (mean ± SE).
FIGURE 8. Changes in the relative intensity of target metabolites identified by UPLC-QTOF/MS. A two tailed, parametric t-test was used to determine the significance of the change in relative intensity for each metabolite. Bars represent the mean relative metabolite concentration and standard deviations. ##P < 0.01: model group vs. control group; ∗P < 0.05, ∗∗P < 0.01: after treatment by HK, YX, DHYL, DHYH, JNL, and JNH vs. model group. Data were presented as the mean ± SE.
Metabolic Pathway Analysis
The metabolic pathway was established by importing the potential metabolites into the web-based database MetPA. The pathway impact value calculated from pathway to topology analysis with MetPA above 0.1 was screened out as the potential target pathway. Here in Figure 9, for the five pathways, sphingolipid metabolism with the impact-value 0.31; pentose, glucuronate interconversion with the impact-value 0.27; terpenoid backbone biosynthesis with the impact-value 0.21; purine metabolism with the impact-value 0.15 and retinol metabolism with the impact-value 0.15 were filtered out as the most important metabolic pathways.
FIGURE 9. Summary of pathway analysis of serum and urine samples of rats. (a) sphingolipid metabolism; (b) pentose and glucuronate interconversions; (c) terpenoid backbone biosynthesis; (d) purine metabolism; (e) retinol metabolism.
Correlation Analysis
Correlation networks is a useful tool to elucidate the relationship between biochemical indicators and potential biomarkers, providing supports for clinical diagnosis, medical treatment and pathophysiology research. Pearson correlation coefficient r > 0.65 and r < -0.65 separately represented for significant positive correlation and negative correlation (Cai et al., 2018). Here in Figure 10, BUN, UP, Scr, β2-MG positively correlated with ceramide (d18:1/12:0) (r > 0.65); and negatively correlated with sphingosine 1-phosphate and phytosphingosine (r < -0.65). Ceramide (d18:1/12:0), phytosphingosine and sphingosine 1-phosphate have been found and used to explain the sphingolipid metabolism. Therefore, BUN, UP, Scr, β2-MG closely related to sphingolipid metabolism.
FIGURE 10. Correlation heat-map between biochemical indicators and potential biomarkers.
Discussion
The local and traditional uses of “Dihuang” are nourishing yin and tonifying kidney. The pharmacological studies stated that “Dihuang” possesses the significant effects on cardiovascular diseases, neuroprotection, diabetes and its complications, osteoporosis, and hyperlipidemia. Moreover, acteoside, as the main active component, was reportedly possesses widely pharmacological activities, such as kidney protection, anti-oxidation, neuroprotection, liver protection, anti-cancer, and so on. But the action mechanisms still unclear completely. In our this study, the results in vivo showed that TLR and DTG can decrease 24-h urinary protein, serum creatine and blood urea nitrogen, alleviating the degree of renal interstitial fibrosis, reducing the expression of Wnt4, β-catenin and TGF-β1, played a role in renal protection. Moreover, the effect of the low dose group was more significant. Therefore, the optimal dose should be determined by further studies. Studies have shown that, the in vivo model of GMCs induced by high glucose simulates the pathophysiological changes of DN (Abboud, 2012). These in vitro results reveal that inhibitory effects of TLR and DTG on the hypertrophy and fibrosis of GMCs, and consistent with the in vivo results in this study. It is confirmed that the TLR and DTG had a better inhibiting effect on TGF-β1 and Wnt/β-catenin signaling pathway.
Notably, there are several distinct metabolic pathways in various stages of progress of DN. Ceramide (d18:1/12:0), phytosphingosine and sphingosine 1-phosphate have been found and used to explain the sphingolipid metabolism. Ceramides are one of the hydrolysis byproducts of sphingomyelin by the enzyme sphingomyelinase (sphingomyelin phosphorylcholine phosphohydrolase E.C.3.1.4.12) which has been identified in the subcellular fractions of human epidermis and many other tissues (Summers, 2006). Besides, the accumulation of ceramide can cause many lipid metabolic diseases (Li et al., 2002), and some pathogen selectively alter metabolic pathways and promote the absorption of fat into ceramides. Studies have shown that high-fat diet for 3 weeks can significantly cause ceramide and sphingomyelin in the mouse liver and nucleus accumulation, and ultimately make the liver browning and lead to fatty liver (Chocian et al., 2010). Scientists have recently demonstrated that inhibition and control of sphingolipid synthase contribute to the treatment of atherosclerosis, insulin resistance, diabetes and cardiomyopathy (Li et al., 2014). In DN, the content of ceramide increased significantly compared to normal rats and recovered to normal after administration of TLG and DTG. In recent years, studies have found that sphingosine 1-phosphate (S1P) can participate in the regulation of GMC growth and differentiation process (Spiegel and Kolesnick, 2002). Phytosphingosine is involved in diverse cell processes, including cell-cell interaction, cell proliferation, differentiation, and apoptosis. The contents of S1P and phytosphingosine decreased in DN, which was consistent with the injured of GMCs induced by DN. In addition, sphingolipid metabolic disorders may lead to cell death and are closely related to insulin-related diseases (Pang et al., 2008).
Cholic acid glucuronide has been found and used to explain the pentose and glucuronate interconversions. Cholic acid glucuronide is the glucuronidated metabolite of cholic acid, one of the four main acids produced by the liver where it is synthesized from cholesterol. Upon formation, the glucuronide is rapidly and effectively cleared from the circulation and excreted via urine (Little et al., 1985). In DN, the decreased amount of cholic acid glucuronide in urine indicated that the renal excretion was perturbed. Guanosine triphosphate, guanosine diphosphate, inosinic acid and adenine have been found and used to explain the purine metabolism. Inosinic acid is a purine nucleotide which has hypoxanthine as the base and one phosphate group esterified to the sugar moiety. The hypoxanthine in the purine metabolism oxidize to xanthine, adenosine monophosphate can degrade into uric acid and accumulate in the urine into uremic toxins, eventually leading to renal failure (Sánchez-Lozada et al., 2005; Seaman et al., 2008). Results suggested that these target pathways showed the marked perturbations over the time-course of the treatment and could contribute to the development of DN.
Conclusion
This study showed that TLR and DTG could prevent its occurrence and development of DN by combining biochemical and pathological indicators with metabolomics technology, and by assessing comprehensive curative effect. Through the application of metabolomics technology, 27 endogenous metabolites (12 in serum and 15 in urine) could be identified in the process of DN. After drug intervention, these markers varied to some extent (P < 0.05). The results showed that TLR and DTG played a significant role in regulating the expression levle of TGF-β1 protein in renal tissue and mesangial cells. Therefore, to further study the mechanism of TLR on TGF-β/Smad signaling pathway, it is necessary to systematically study the protective mechanism of TLR on DN. The study paved the way for further exploration of the pathogenesis of DN, early diagnosis and the evaluation of curative effect.
Author Contributions
SS and J-aD designed the project. XD, HC, TZ, and HY performed the experiments. ZZ, E-xS, SG, DQ, and DW provided technology guidance. XD and SS wrote and modified the manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (Nos. 81373889 and 81673533) and supported by program for excellent talents in school of pharmacy of Nanjing University of Chinese Medicine (15ZYXET-2). This work was also financially supported by Jiangsu Collaborative Innovation Center of Chinese Medicinal Resources Industrialization (No. ZDXM-1-10), Construction Project for Jiangsu Key Laboratory for High Technology of TCM Formulae Research (BM2010576), and a project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (ysxk-2014).
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.
Abbreviations
β2-MG, serum β2-microglobulin; BUN, serum urea nitrogen; C, control group; DHYH, TLR high dose group; DHYL, TLR low dose group; DN, diabetic nephropathy; DTG, Dihuangye total glycoside capsule; FBG, Fasting blood glucose; GMC, glomerular mesangial cells; HE, hematoxylin eosin; HK, Huangkui capsule group; JNH, DTG high dose group; JNL, DTG low dose group; LDL, low density lipoprotein; M, model group; MRHTG, Acteoside; RG, Rehmannia glutinosa Libosch; Scr, serum creatinine; STZ, streptozotocin; T-CHO, total cholesterol; TCM, Traditional Chinese medicine; TG, triglyceride; TLR, total glycoside extracted from leaves of Rehmannia; UP, urine protein; YX, Irbesartan group.
Footnotes
- ^http://www.metaboanalyst.ca/faces/upload/PathUploadView.xhtmlpc
- ^http://www.genome.jp/kegg/
- ^http://www.hmdb.ca/
- ^http://www.hmdb.ca/
References
Abboud, H. E. (2012). Mesangial cell biology. Exp. Cell Res. 318, 979–985. doi: 10.1016/j.yexcr.2012.02.025
Balakumar, P., Arora, M. K., and Singh, M. (2009). Emerging role of PPAR ligands in the management of diabetic nephropathy. Pharmacol. Res. 60, 170–173. doi: 10.1016/j.phrs.2009.01.010
Beijing Municipal Health Bureau (1988). Beijing Standard of Chinese Herbal Medicines. Beijing: Capital Normal University press.
Cai, H. D., Su, S. L., Li, Y. H., Zeng, H. T., Zhu, Z. H., Guo, J. M., et al. (2018). Salvia miltiorrhiza protects adenine-induced chronic renal failure by regulating the metabolic profiling and modulating the NADPH oxidase/ROS/ERK and TGF-β/Smad signaling pathways. J. Ethnopharmacol. 212, 153–165. doi: 10.1016/j.jep.2017.09.021
Cai, Q. Y., Chen, X. S., Zhan, X. L., and Yao, Z. X. (2011). Protective effects of catalpol on oligodendrocyte death and myelin breakdown in a rat model of chronic cerebral hypoperfusion. Neurosci. Lett. 497, 22–26. doi: 10.1016/j.neulet.2011.04.013
Chen, Y. J., and Quilley, J. (2008). Fenofibrate treatment of diabetic rats reduces nitrosative stress, renal cyclooxygenase-2 expression, and enhanced renal prostaglandin release. J. Pharmacol. Exp. Ther. 324, 658–663. doi: 10.1124/jpet.107.129197
Chocian, G., Chabowski, A., Zendzian-Piotrowska, M., Harasim, E., Łukaszuk, B., and Górski, J. (2010). High fat diet induces ceramide and sphingomyelin formation in rat’s liver nuclei. Mol. Cell. Biochem. 340, 125–131. doi: 10.1007/s11010-010-0409-6
Hodavance, M. S., Ralston, S. L., and Pelczer, I. (2007). Beyond blood sugar: the potential of NMR-based metabonomics for type 2 human diabetes, and the horse as a possible model. Anal. Bioanal. Chem. 387, 533–537. doi: 10.1007/s00216-006-0979-z
Hwang, G. S., Yang, J. Y., Ryudo, H., and Kwon, T. H. (2010). Metabolic profiling of kidney and urine in rats with lithium-induced nephrogenic diabetes insipidus by (1) H-NMR-based metabonomics. Am. J. Physiol. Renal Physiol. 298,F461–F470. doi: 10.1152/ajprenal.00389.2009
Jiménez, C., and Riguera, R. (1994). Phenylethanoid glycosides in plants: structure and biological activity. Nat. Prod. Rep. 11, 591–606. doi: 10.1039/NP9941100591
Li, H., Junk, P., Huwiler, A., Burkhardt, C., Wallerath, T., Pfeilschifter, J., et al. (2002). Dual effect of ceramide on human endothelial cells: induction of oxidative stress and transcriptional upregulation of endothelial nitric oxide synthase. Circulation 106, 2250–2256. doi: 10.1161/01.CIR.0000035650.05921.50
Li, W. X., Tang, Y. P., Guo, J. M., Shang, E.-X., Qian, Y. F., Wang, L. Y., et al. (2014). Comparative metabolomics analysis on hematopoietic functions of herb pair Gui-Xiong by ultra-high-performance liquid chromatography coupled to quadrupole time-of-flight mass spectrometry and pattern recognition approach. J. Chromatogr. A 1346, 49–56. doi: 10.1016/j.chroma.2014.04.042
Li, X., Becker, K. A., and Zhang, Y. (2010). Ceramide in redox signaling and cardiovascular diseases. Cell Physiol. Biochem. 26, 41–48. doi: 10.1159/000315104
Little, J. M., Chari, M. V., and Lester, R. (1985). Excretion of cholate glucuronide. J. Lipid Res. 26, 583–592.
Liu, C. L., Cheng, L., Kwok, H. F., Ko, C. H., Lau, T. W., Koon, C. M., et al. (2011). Bioassay-guided isolation of norviburtinal from the root of Rehmannia glutinosa, exhibited angiogenesis effect in zebrafish embryo model. J. Ethnopharmacol. 137, 1323–1327. doi: 10.1016/j.jep.2011.07.060
Lv, M., Chen, Z., Hu, G., and Qian, B. L. (2015). Therapeutic strategies of diabetic nephropathy: recent progress and future perspectives. Drug Discov. Today 20, 332–346. doi: 10.1016/j.drudis.2014.10.007
Nicholson, J. K., and Lindon, J. C. (2008). Systems biology: metabonomics. Nature 7216:1054. doi: 10.1038/4551054a
Oberg, B. P., McMenamin, E., Lucas, F. L., McMonagle, E., Morrow, J., Ikizler, T. A., et al. (2004). Increased prevalence of oxidant stress and inflammation in patients with moderate to severe chronic kidney disease. Kidney Int. 65, 1009–1016. doi: 10.1111/j.1523-1755.2004.00465.x
Odunsi, K., Wollman, R. M., Ambrosone, C. B., Hutson, A., McCann, S. E., Tammela, J., et al. (2005). Detection of epithelial ovarian cancer using 1H-NMR-based metabonomics. Int. J. Cancer 113, 782–788. doi: 10.1002/ijc.20651
Pang, L. Q., Liang, Q. L., Wang, Y. M., Ping, L., and Luo, G. A. (2008). Simultaneous determination and quantification of seven major phospholipid classes in human blood using normal-phase liquid chromatography coupled with electrospray mass spectrometry and the application in diabetes nephropathy. J. Chromatogr. B 869, 118–125. doi: 10.1016/j.jchromb.2008.05.027
Pereira, T., and Chang, S. W. (2004). Semi-automated quantification of ivermectin in rat and human plasma using protein precipitation and filtration with liquid chromatography/tandem mass spectrometry. Rapid Commun. Mass Spectrom. 18, 1265–1276. doi: 10.1002/rcm.1485
Sánchez-Lozada, L. G., Tapia, E. J., Avila, C. C., Avila-Casado, C., Soto, V., Nepomuceno, T., et al. (2005). Mild hyperuricemia induces vasoconstriction and maintains glomerular hypertension in normal and remnant kidney rats. Kidney Int. 67, 237–247. doi: 10.1111/j.1523-1755.2005.00074.x
Seaman, C., Moritz, J., Falkenstein, E., Van Dyke, K., and Klandorf, H. (2008). Inosine ameliorates the effects of hemin-induced oxidative stress in broilers. Comp. Biochem. Physiol. A 151, 670–675. doi: 10.1016/j.cbpa.2008.08.014
Shen, X., Li, D. F., Zong, G. Z., Wu, Z. L., and He, W. (2010). Effects of total saponins extracted from leaves of rehmannia on accelerated nephrotoxic nephritis induced by rabbit IgG in rat. Chin. J. Exp. Tradit. Med. Form. 16, 179–181.
Spiegel, S., and Kolesnick, R. (2002). Sphingosine 1-phosphate as a therapeutic agent. Leukemia 16, 1596–1602. doi: 10.1038/sj.leu.2402611
Summers, S. A. (2006). Ceramides in insulin resistance and lipotoxicity. Prog. Lipid. Res. 45, 42–72. doi: 10.1016/j.plipres.2005.11.002
Weir, M. R. (2007). Microalbuminuria and cardiovascular disease. Clin. J. Am. Soc. Nephro. 2, 581–590. doi: 10.2215/CJN.03190906
Xiao, L., Wang, M., Yang, S., Liu, F. Y., and Sun, L. (2013). A glimpse of the pathogenetic mechanisms of Wnt/β-catenin signaling in diabetic nephropathy. Biomed. Res. Int. 2013:987064. doi: 10.1155/2013/987064
Yokozawa, T., Suzuki, N., Oura, H., Nonaka, G., and Nishioka, I. (1986). Effect of extracts obtained from rhubarb in rats with chronic renal failure. Chem. Pharm. Bull. 34, 4718–4723. doi: 10.1248/cpb.34.4718
Zhang, R., Zhou, J., Li, M., Ma, H., Qiu, J., Luo, X., et al. (2014). Ameliorating effect and potential mechanism of Rehmannia glutinosa oligosaccharides on the impaired glucose metabolism in chronic stress rats fed with high-fat diet. Phytomedicine 21, 607–614. doi: 10.1016/j.phymed.2013.11.008
Zhang, R., Zhao, Y., Sun, Y., Lu, X., and Yang, X. (2013). Isolation, characterization, and hepatoprotective effects of the raffinose family oligosaccharides from Rehmannia glutinosa Libosch. J. Agric. Food Chem. 61, 7786–7793. doi: 10.1021/jf4018492
Keywords: total glycoside extracted from leaves of Rehmannia, Dihuangye total glycoside capsule, metabolomics, diabetic nephropathy, TGF-β1, Wnt/β-catenin signaling pathway
Citation: Dai X, Su S, Cai H, Wei D, Yan H, Zheng T, Zhu Z, Shang E-x, Guo S, Qian D and Duan J-a (2018) Protective Effects of Total Glycoside From Rehmannia glutinosa Leaves on Diabetic Nephropathy Rats via Regulating the Metabolic Profiling and Modulating the TGF-β1 and Wnt/β-Catenin Signaling Pathway. Front. Pharmacol. 9:1012. doi: 10.3389/fphar.2018.01012
Received: 07 February 2018; Accepted: 20 August 2018;
Published: 11 September 2018.
Edited by:
Jiang Bo Li, Second People’s Hospital of Wuhu, ChinaReviewed by:
Jianxin Chen, Beijing University of Chinese Medicine, ChinaYanjun Zhang, Tianjin University of Traditional Chinese Medicine, China
Copyright © 2018 Dai, Su, Cai, Wei, Yan, Zheng, Zhu, Shang, Guo, Qian and Duan. 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(s) 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: Shulan Su, sushulan1974@163.com Jin-ao Duan, dja@njutcm.edu.cn