Scholarly article on topic 'Disease-syndrome combination modeling: metabolomic strategy for the pathogenesis of chronic kidney disease'

Disease-syndrome combination modeling: metabolomic strategy for the pathogenesis of chronic kidney disease Academic research paper on "Veterinary science"

CC BY
0
0
Share paper
Academic journal
Sci. Rep.
OECD Field of science
Keywords
{""}

Academic research paper on topic "Disease-syndrome combination modeling: metabolomic strategy for the pathogenesis of chronic kidney disease"

SCIENTIFIC REPORTS

Received: 17 May 2017 Accepted: 25 July 2017 Published online: 18 August 2017

Disease-syndrome combination modeling: metabolomic strategy for the pathogenesis of chronic kidney disease

Shasha Li 1, Peng Xu1, Ling Han1, Wei Mao\Yiming Wang1'2, Guoan Luo1'2 & Nizhi Yang1

Conventional disease animal models have limitations on the conformity to the actual clinical situation. Disease-syndrome combination (DS) modeling may provide a more efficient strategy for biomedicine research. Disease model and DS model of renal fibrosis in chronic kidney disease were established by ligating the left ureter and by ligating unilateral ureteral combined with exhaustive swimming, respectively. Serum metabolomics was conducted to evaluate disease model and DS model by using ultra performance liquid chromatography coupled with quadrupole time-of-flight mass spectrometry. Potential endogenous biomarkers were identified by multivariate statistical analysis. There are no differences between two models regarding their clinical biochemistry and kidney histopathology, while metabolomics highlights their difference. It is found that abnormal sphingolipid metabolism is a common characteristic of both models, while arachidonic acid metabolism, linolenic acid metabolism and glycerophospholipid metabolism are highlighted in DS model. Metabolomics is a promising approach to evaluate experiment animal models. DS model are comparatively in more coincidence with clinical settings, and is superior to single disease model for the biomedicine research.

As an absolutely essential part to most of the modern medicine research, animal models were widely used to investigate disease or injury states especially in ways which would be inaccessible in patients. Most of the animal models serving in research at present are disease models, which may have an existing, inbred or induced disease similar to a human condition, and thus may meet most of the requirement of western medicine researches.

However, people come to realize the limitation of an animal model1, that the duplicated disease in animal is often different from the disease in human being. Human diseases, in most cases, are multifactorial and thus often accompanied by complications or a series of other symptoms. For instance, most patients of renal fibrosis are also characterized with weakness, fatigue, dizziness, deficient sweating, stuffy chest, short breath and blood stasis, which are, according to the theory of traditional Chinese medicine (TCM), typical observations of syndrome of qi-deficiency and blood-stasis (QDBS). The most common rat model of renal fibrosis is established by ligation of the left ureter2, and validated with blood biochemistry and kidney histopathology indicators, without too much attention on other symptoms that should be present. As a result, this kind of simple disease animal model cannot reflect the real condition of human being, which might be, as least in part, responsible to the possible inconsistency between animal experiments and clinical trials3.

In theory and practice, TCM believed that disease does not exist in isolation, often accompanied by some specific "syndrome" (more accurately, maybe should be transliterated as "Zheng")4. Disease-Syndrome (DS) combination modeling5, taking both disease and syndrome into consideration, is expected to be an efficient solution deal with the limitation of disease model. Then, how to evaluate the validity of an established DS model? A disease model can be evaluated by determination on some addressed biochemical indicators, which is, nevertheless, may not suffice to a DS model. Syndrome is a holistic concept with global pathological or physiological abnormalities occurred in patients or model animal, and can hardly be described by limited biochemical index6. Metabolomics

Guangdong Provincial Hospital of Chinese Medicine, No. 111 Dade Road, Guangzhou, Guangdong, 510120, China. department of Chemistry, Tsinghua University, No. 30 Shuangqing Road in Haidian Distric, Beijing, 100084, China. Shasha Li and Peng Xu contributed equally to this work. Correspondence and requests for materials should be addressed to G.L. (email: luoga@mail.tsinghua.edu.cn) or N.Y. (email: yangnizhi_gzucm@163.com)

is a holistic analysis strategy widely applied in the field of biomedicine research and in clinical application, and may be used as a powerful tool for the evaluation of TCM and DS models7' 8.

Chronic kidney disease (CKD) is a serious disease that jeopardizes human health worldwide9. The common pathological alteration in progressive CKD was tubulointerstitial fibrosis10. Tubulointerstitial fibrosis was complex, involving various independent and overlapping cellular and molecular pathways11. It was reported that metabolomics has been used to study patients with CKD for the last several years12-16. Compared with clinical aplication, metabolomics has been applied widely to the study of different animal models of CKD such as 5/6 nephrectomized CKD, adenine-induced chronic renal failure and other drug-induced CKD17-23. However, most publications were just based on the studies on animal disease models24. Statistic study of clinical cases in our hospital found that renal fibrosis patients had mostly a syndrome of qi-deficiency and blood-stasis (QDBS), and this conclusion has also been supported by literature in year 201025. Ultra-high performance liquid chromatography-quadrupole time-of-flight mass spectrometry (UPLC-QTOF/MS) with rapidly separation and high reproducibility is regarded as a suitable for untargeted metabolomics in disease biomarker discovery and drug research26-29. In this present study, a DS model was established in rats as a combination model of QDBS and renal fibrosis in CKD, and metabolomic approach was used to comparatively evaluate this model. It is expected to discover novel diagnostic biomarkers and to highlight pathophysiological mechanisms of renal fibrosis in CKD.

Materials and Methods

Animals and Treatment. Thirty-two specific pathogen-free Sprague Dawley (SD) rats (male, weight, 180220 g) were provided by the Experimental Animal Center of Guangdong Medical (certificate number, 0079361). All experiments were performed in accordance with the internationally accepted standard guidelines for the use of animals. This study was conducted in accordance with the Chinese national legislation and local guidelines, and the care and handling of rats were also approved by the Ethical Committee of Guangdong Provincial Hospital of Chinese Medicine. The rats were housed under standard environmental conditions (23 ± 2 °C, 55% ± 5% humidity and 12 h/12 h light/dark cycle) and were allowed free access to water as well as standard laboratory rat food.

After 3 days of acclimatization, the rats were randomly divided into three groups of control group (n = 10), disease model group (n = 10) and DS model group (n = 12), respectively. Rats in disease and DS model group underwent ligation of the left ureter to establish a model of renal fibrosis, and rats in DS group were further forced to swim exhaustively every day to establish a DS model; they swam in a pool of water with a depth of 60 cm at room temperature and were picked up when 50% of the rats were exhaustively submerged in water for 10 s. This experiment lasted for 28 days.

Sample Preparation. Three groups of rats were sacrificed after being anesthetized by intraperitoneal injection of pentobarbital (50 mg/kg of body weight) (GBCBIO Technologies, Guangzhou, China) and their blood and kidney were collected. Blood samples were drawn into tubes, allowed to stand for 30 min, and were centrifuged to obtain serum. Serum samples were snap frozen in liquid nitrogen for 5 min and then stored in a refrigerator at -80 °C. Kidney tissues were fixed with 10% formalin (Sigma, America) in PBS for 12 h after dehydration, cleared in xylene (Damao chemical reagent company, Tianjin, China), penetrated in wax and embedded in paraffin (Leica, Germany).

Thawed serum (400 |iL) was mixed with methanol (1200 |iL), then swirled for 2 min and centrifuged at 13,000 rpm for 15 min at 4 °C. The clear supernatant was transferred and evaporated to dryness by N2, followed by dilution with 300 (L of methanol-H2O (80:20, v/v) (Merck, Germany) and centrifugation at 13,000 rpm for 15 min at 4 °C. The supernatant was transferred to a new EP tube for UPLC-QTOF/MS analysis.

Chromatographic Condition. Chromatography was performed on a Waters Acquity UPLC BEH C18 column (2.1 mm x 100 mm, 1.7 |im) at 30 °C. The mobile phase consisted of acetonitrile (Merck, Germany) as solvent A and 0.1% formic acid (Merck, Germany) in water as solvent B. The gradient programs was as follows: 0-5 min, 2-50% A; 5-7 min, 50-60% A; 7-14 min, 60-65% A; 14-18 min, 65-80% A; 18-24 min, 80-90% A; 24-26 min, 90-95% A; 26-28 min, 95-2% A; and 28-30 min, 2% A. The injection volume of the test sample was 5 (L and each sample was injected three times. Each wash cycle consisted of 200 (L of strong wash solvent (80% CH3CN-H2O, 8:2, v/v) and 600 (L ofweak wash solvent (10% CH3CN-H2O, 1:9, v/v).

Mass spectrometry condition. MS was performed on a Waters Q-TOF Premier with an electrospray ionization system (Waters MS technologies, Manchester, UK). A preliminary experiment was conducted to optimize the experimental conditions. To achieve the desired detection results, the flow rate and column temperature for chromatography, as well as capillary voltage, flow and temperature of the desolvation gas for the mass spectrometry detector were optimized carefully. As a result, the optimal parameters were fixed as followed. The capillary voltages were set at 3.0 and 2.5 kV, and the sample cone voltages were 30 and 45 V for the positive and negative ion modes, respectively. The desolvation gas flow was set to 600 L/h at a desolvation temperature of 350 °C, the cone gas set to 50 L/h and the source temperature was 110 °C. The Q-TOF Premier acquisition rate was 0.1 s with a 0.02 s inter-scan delay. The instrument was operated with the first resolving quadrupole in a wide-pass mode (50-1500 Da). Argon was used as the collision gas. ESI-MS spectra were acquired in both ionization modes of positive and negative with full-scan detection.

The precision and stability of the UPLC-MS method was determined by repeated analysis of six injections of the same QC samples, and the repeatability of sample preparation was accessed by preparing six parallel samples using the same protocol. All assays were performed using a lock spray to ensure accuracy and reproducibility, with Leucine enkephalin amide acetate as the lock mass at 200 pg/|iL concentration and 20 |iL/min flow rate.

0 5 10 15 20 25 30 0 5 10 15 20 25 30

Times (day) Times(day)

Figure 1. The change of weight in all three groups and swimming time in DS model group. The weight of all rats was measured weekly and the swimming time of the rats for the DS model group was measured every 2 days from the start to the end of the experiment. Rats in DS model group were weighed before exhaustive swimming. Data were reported as the mean ± standard deviation in each group. *P < 0.05, compared with the control group. The swimming time was calculated from the point when rats were placed in water to the point when 50% of the rats were exhaustively submerged in water for 10 s.

Data processing. Data pre-processing was performed with MarkerLynx applications manager version 4.1 (Waters, MA, USA). The parameters included retention time from 1 min to 26 min, mass range from 50 Da to 1500 Da and mass tolerance was 0.01 Da. The minimum intensity was set at 10% of base peak intensity, the maximum mass per retention time was set at 6 and the retention time tolerance was 0.05 min. MarkerLynx incorporates an Apex-Track-peak detection package and the data were combined into a single matrix by aligning peaks with the exact mass/retention time pair together. The ion intensities for each peak were normalized to a constant integrated intensity to partially compensate for the concentration bias of each sample.

The processed data was then exported and processed by principle component analysis (PCA) and partial least square-discriminant analysis (PLS-DA) in the software package SIMCA-P 11.5 version (Umetrics AB, Umea, Sweden). PCA was used first to determine the general interrelation between the three groups, and PLS-DA was subsequently performed to maximize the difference in metabolic profiling. The value of Q2 and R2 were used to estimate the accuracy of the model and a typical 7-round cross-validation was used to validate the model against over fitting.

Statistical analyses were performed using SPSS software (Version 18.0, USA). Assumptions of normality and homogeneity of variance were first checked. Data were presented as the mean ± standard deviation for continuous variables with a normal distribution, as counts and percentages for categorical variables. The independent samples i-test or one-way ANOVA were used to analyze the differences among groups for continuous measures. Differences with P values less than 0.05 were considered to have statistical differences, and P values less than 0.01 were considered to have significant differences. All probability values were two-sided.

Data availability. All data generated during this study are included in this published article and the raw data are available from the corresponding author on reasonable request.

Results

Change of weight and swimming time of rats. The weight of all rats was recorded every week, and the swimming time of the rats for the DS model group was measured every 2 days. Figure 1a shows that the weight of all rats increased with time. Starting from the second week, the weight growth rate of the control group and disease model group was faster than that of DS model group, and this phenomenon persisted until the end of the experiment. The Dunnett's test showed that the weights of the rats in each group were equal at the beginning of the experiment; however, the final weights of the DS model group were significantly lower than that of the control group and the disease group (P = 2.92 x 10-6 and 1.05 x 10-3, respectively). No significant difference was observed between the control and disease model group. Figure 1b shows the exhaustive swimming time of the DS model group during the adaptive swimming period. After stress reduction in the second day, the swimming time of rats shortened gradually in followed two weeks. Then, the time remained relatively stable until the end of the experiment.

Serum biochemistry. Thirty-two rats were included in this study and the results of serum biochemistry are listed in Table 1. Compared with control rats, rats of both model groups show significant changes in levels of cholesterol, lactate dehydrogenase, creatine kinase, blood urea nitrogen and creatinine. However, excepting triglyceride, these biochemical parameters in rats showed no statistical differences between disease group and DS group.

Immunohistochemical staining. Paraffin embedded kidney tissue were cut into 4 |im sections. Hematoxylin & Eosin (H&E) and Masson staining were performed on sections for observation of kidney pathohistological changes. As illustrated in Fig. 2, compared with control group, those in disease and DS groups

Control group Disease model group DS model group

No. of subjects 10 10 12

Cholesterol, mmol/L 1.18 ± 0.20 2.56 ± 0.30** 2.66 ± 0.29**

Triglyceride, mmol/L 0.37 ± 0.10 0.57 ± 0.08** 0.37 ± 0.12aa

Lactate dehydrogenase, U/L 429.20 ± 89.17 1131.11 ± 122.93** 1111.86 ± 91.61**

Creatine kinase, U/L 225.80 ± 38.88 416.92 ± 95.67** 443.38 ± 80.69**

Blood urea nitrogen, mmol/L 7.40 ± 1.55 9.90 ± 1.06** 9.32 ± 1.13*

Creatinine, ^mol/L 14.75 ± 2.75 54.60 ± 8.40** 55.29 ± 13.46**

Table 1. Comparison of serum biochemistry analysis among groups"- b. aThe serum levels of biochemical indices of rats are presented as the mean ± standard deviation. bThe serum concentrations of biochemical indices were statistically tested by Dunnett's i-test when the variance in different groups was homogeneous and by the Games-Howell test when the variance was non-homogeneous. *P< 0.05, **P< 0.01 compared with the control group; AP < 0.05, AAP < 0.01 compared with the disease model group.

H&E staining MASSON staining

Tubular Glomerulus Tubular Glomerulus

Figure 2. H& E and Masson staining on kidney tissue of rats in control, disease and DS groups. H& E staining: rats in disease and DS groups were observed with renal cortex and medulla atrophied, glomerular sclerosis and atrophy, renal tubular expansion, structure diffusion, visible hyalinization, protein tube and a large number of cavity; Glomerular and capillary congestion and clear inflammatory cell infiltration was also observed. Masson staining: glomerular basement membrane and tubule interstitial fibers proliferated significantly in disease and DS groups of rats.

were observed with renal cortex and medulla atrophied, and microscopic examination showed the development of glomerular sclerosis and atrophy by H&E staining, together with renal tubular expansion, structure diffusion, visible hyalinization, protein tube and a large number of cavity; in addition, glomerular and capillary congestion and clear inflammatory cell infiltration was also observed in disease and DS groups of rats. Masson staining showed dramatically tubulointerstitial fibrosis. Immunohistochemical staining indicates both disease and DS models successfully established renal fibrosis models, while it cannot show any significant variations between them.

Metabolic Fingerprinting. Based on the optimal conditions, the UPLC-MS base peak intensity chromato-grams obtained in the ESI positive and negative ion modes are shown in Fig. 3. In this study, abundant chromatographic peak information was observed both in the positive and negative ion modes. To acquire the maximum amount of detectable metabolites, the detection data of two ionization patterns were all considered in further study. Quasi-molecular ion [M + H]+ was the most usually form with a higher signal/noise in the positive ion mode and [M — H] — or [M + HCOO] — in the negative ion mode.

Figure 3. Representative base peak intensity chromatogram of rat serum obtained in the ESI positive (a) and negative ion modes (b) based on UPLC-QTOF/MS. Con, control group; DM, disease model group; DSM, DS model group.

The relative standard deviations of the peak retention time and area value in methodological investigation were all less than 5.0%, indicating that the precision and repeatability of the proposed method were satisfactory for metabolomic analysis.

PCA and PLS-DA. In this study, PCA was used first to determine the general interrelation between the three groups. For a clearer display of clustering, PCA 3D score plots were provided in this study. As shown in Fig. 4, the model groups and control group were separated clearly, indicating that the rats of model groups have significantly altered endogenous metabolism both in the positive and negative ion modes. Furthermore, the assembly of DS model group was far from control group compared to the disease model group both in positive and negative ion modes, indicating a more obvious change in metabolic state of DS model group. This indicates that metabolic fingerprints could reflect the alternation between disease and normal rats. In addition, metabolic fingerprints could exhibit the difference between disease model group and DS model group that cannot be observed in detection of serum biochemistry and immunohistochemical staining.

Subsequently, PLS-DA was performed to maximize the difference in metabolic profiling (Fig. 5). The R2Y and Q2Y values calculated from the results in SIMCA-P package were 0.990 and 0.944 in the positive ion mode, meaning 99% of data fit the model and 94.4% of data could be predicted by this model. Besides, the R2Y and Q2Y values in the negative ion mode were 0.989 and 0.919 respectively. Both the Q2Y and R2Y close to 1 indicate an excellent model which is good to fitness and prediction. In order to maximize the difference, Data of control model and disease model or DS model groups were performed PLS-DA, respectively.

Biomarkers. After PLS-DA processing, metabolites that significantly contributed to the clustering and discrimination were identified when the variable importance in the projection values was more than 1 and p-value of i-test between groups was less than 0.05. These endogenous metabolites were selected for further structural identification. In next process, the calculated mass, mass deviation (mDa and ppm), double bond equivalent,

.30 .20 -10 0 10 20 30 -30 -20 -10 0 10 20 30 40

Figure 5. PLS-DA score plots of control group and model groups. (a) Comparison of groups in positive ion mode (R2X(cum) = 0.349, R2Y(cum) = 0.990, Q2Y(cum) = 0.944); (b) Comparison of groups in negative ion mode (R2x(cum) = 0.279, R2Y(cum) = 0.989, Q2Y(cum) = 0.919). Con, control group; DM, disease model group; DSM, DS model group.

formula, and isotopic pattern of the selected ion (¿-fit value) were used to evaluate the accuracy of possible formulas30' 31. A smaller mass deviation and a lower ¿-fit value indicate a more accurate elemental composition. The structural information of the metabolites selected was subsequently obtained by searching the freely accessible databases HMDB (http://www.hmdb.ca) and KEGG (http://www.genome.jp). Finally, twenty significantly differential metabolites were identified (Table 2). Moreover, the relative levels of biomarkers were also analyzed and were shown in Table 2, and the heatmap of these 20 potential biomarkers is shown in Fig. 6.

Eleven and fifteen potential biomarkers were identified in disease model and DS model, respectively. Leukotriene E3(LTE3), 15(S)-Hydroxyeicosatrienoic acid (15(S)-HETrE), LysoPE (20:1), Tyrosyl-histidine, Glutaminyl-glutamate and Sphingosine 1-phosphate(S1P) were observed as potential biomarkers in both models and these biomarkers show the same variation trends. No. 7-11 and No. 12-20 metabolites were particular biomarkers obtained by PLS-DA analysis in disease model and DS model, respectively (Table 2 and Fig. 7). The relative level of all twenty potential biomarkers was listed in Table 2 and a one-way ANOVA was conducted to evaluate the significance of difference between groups. The _p-value in Table 2 showed that the relative level of five among No.12-20 had significant difference in disease model compared to control group, and No.7-11 metabolites had significant difference in DS model. This result illustrated again that the DS model could exhibit the abnormal endogenous metabolic status comprehensively.

Metabolic pathway analysis revealed key pathways involved in disease model and DS model32, 33. According to Fig. 8, abnormal sphingolipid metabolism is a common characteristic of both models, while there are some other important pathways of arachidonic acid (AA) metabolism, linoleic acid metabolism and glycerophospholipid metabolism are involved highlight in DS model. The biosynthetic pathways for parts of potential biomarkers related to arachidonic acid metabolism, linoleic acid metabolism, glycerophospholipid metabolism and sphingo-lipid metabolism were shown in Fig. 9.

Discussion

In clinic, serum creatinine, blood urea nitrogen, cholesterol and triglyceride levels are the main biochemical biomarkers that are preferred to detect patients with nephropathy. Among these biomarkers, serum creatinine and blood urea nitrogen levels are most reliable for the detection of renal function, whereas cholesterol and triglyceride levels are related to lipid metabolism34, 35. In addition, creatine kinase and lactate dehydrogenase levels were used to access the physical condition of rats in the present experiment. This study applied two approaches to establish models for the same disease, including a common disease model and a disease model combined syndrome under the guidance of the TCM theory. Table 1 indicates that the two types of modeling methods can induce injury of the kidneys, which presented as damaged renal function, reduced physical condition and disordered the lipid metabolism. Except triglyceride, all biochemical indices in two models show significant variance compared with control rats, which indicate the validity of both models. The difference in triglyceride levels observed between the disease and DS model groups may be due to daily swimming, which requires a great amount of energy.

Disease model and DS model show no difference in biochemical biomarkers except triglyceride and suggests common approach is not sufficient to evaluate DS modeling. Metabolomics focuses on the quantitative and qualitative analysis of low-molecular-weight metabolites and is an efficient tool for globally understanding of diseases pathogeneses. Disease model and DS model have almost the same biochemical parameters, while metabolomics may highlight their difference. Besides 6 biomarkers in common, five metabolites were discovered only in disease model and control groups; while investigation on DS model and control groups identified nine particular

No. Identified potential biomarker Related metabolism Concentration Fc Sig.c p valued

Control Disease model0 DS model' Con. VS DM Con. VS DSM

1 Leukotriene E3 Arachidonic acid metabolism 34.54 ± 5.67 14.81 ± 6.44| 7.94 ± 2.274 81.62 1.23 x 10-12 2.52 x 10-8 2.33 x 10-8

2 15(S)-Hydroxyeicosatrienoic acid Y-Linolenic acid metabolism 124.22 ± 14.73 101.64 ± 29.31| 44.29 ± 7.874 53.00 2.07 x 10-10 2.27 x 10-2 2.37 x 10-8

3 LysoPE (20:1) Glycerophospholipid metabolism 14.93 ± 4.49 1.78 ± 0.82| 5.80 ± 1.404 63.35 2.61 x 10-11 2.33 x 10-8 4.16 x 10-8

4 Tyrosyl-histidine Tyrosine metabolism 16.88 ± 1.46 21.01 ± 2.26t 29.40 ± 3.921 56.02 1.10 x 10-10 5.56 x 10-3 2.34 x 10-8

5 Glutaminyl-glutamate Glutamate metabolism 17.29 ± 1.19 21.12 ± 3.58t 27.65 ± 3.60t 32.61 3.79 x 10-8 1.67 x 10-2 4.25 x 10-8

6 Sphingosine 1-phosphate Sphingolipid metabolism 34.68 ± 6.52 43.28 ± 5.38t 52.80 ± 12.55t 11.01 2.78 x 10-4 1.45 x 10-2 1.22 x 10-4

7 Docosahexaenoic acid a-Linolenic acid metabolism 215.29 ± 21.48 288.42 ± 37.48t 283.26 ± 32.33t 17.50 1.04 x 10-5 2.59 x 10-5 3.96 x 10-5

8 8,11,14-eicosatrienoyl ethanolamide a-Linolenic acid metabolism 26.96 ± 6.69 41.47 ± 9.72t 46.69 ± 9.00t 15.01 3.35 x 10-5 1.42 x 10-3 1.84 x 10-5

9 11b-PGF2a Arachidonic acid metabolism 26.17 ± 7.08 43.35 ± 11.741 41.89 ± 8.64t 10.71 3.29 x 10-4 5.52 x 10-4 9.03 x 10-4

10 Palmitoleic acid Fatty Acid Biosynthesis 57.01 ± 11.35 94.07 ± 18.911 85.33 ± 23.24t 10.66 3.38 x 10-4 2.65 x 10-4 2.87 x 10-3

11 Cholic acid Bile Acid Biosynthesis 39.07 ± 12.93 144.96 ± 31.70t 166.96 ± 26.461 78.00 2.14 x 10-12 2.38 x 10-8 2.33 x 10-8

12 Arachidonic acid Arachidonic acid metabolism 35.62 ± 3.81 30.06 ± 8.904 16.72 ± 5.574 25.73 3.74 x 10-7 1.09 x 10-1 2.96 x 10-7

13 Eicosapentaenoic acid A-Linolenic acid metabolism 245.94 ± 56.43 153.10 ± 39.184 122.45 ± 39.354 21.35 1.99 x 10-6 1.58 x 10-4 1.18 x 10-6

14 9-HODE Linoleic acid metabolism 35.18 ± 1.40 35.76 ± 7.15 21.41 ± 4.154 32.27 4.22 x 10-8 9.47 x 10-1 4.88 x 10-7

15 Phytosphingosine Sphingolipid metabolism 461.18 ± 55.03 403.14 ± 35.194 180.93 ± 66.804 81.25 1.30 x 10-12 4.57 x 10-1 2.33 x 10-8

16 LysoPC (16:0) Glycerophospholipid metabolism 13.97 ± 5.31 16.45 ± 2.99t 22.59 ± 9.59t 4.73 1.67 x 10-2 6.29 x 10-1 1.16 x 10-2

17 LysoPE (18:1) Glycerophospholipid metabolism 10.57 ± 1.24 4.47 ± 1.974 4.01 ± 1.274 60.62 4.38 x 10-11 2.47 x 10-8 2.34 x 10-8

18 LysoPC (15:0) Glycerophospholipid metabolism 42.71 ± 5.02 29.30 ± 8.414 17.36 ± 10.694 23.96 7.20 x 10-7 2.88 x 10-3 2.84 x 10-7

19 Tri-N-acetyl chitotriose Amino sugar metabolism 52.54 ± 10.38 38.94 ± 7.024 23.61 ± 4.994 39.54 5.20 x 10-9 7.98 x 10-4 2.52 x 10-8

20 17ß-Estradiol glucuronide Steroid hormone biosynthesis 41.38 ± 11.24 30.78 ± 6.304 17.92 ± 3.174 27.42 2.06 x 10-7 6.56 x 10-3 1.04 x 10-7

Table 2. Identification of significantly different potential endogenous metabolites in the sera of individuals in the control and model groups". "Metabolites were confirmed by tR and m/z with authentic chemicals. 'Compared with control group. "f " indicates a higher level of metabolites, whereas " j" represents a lower level of metabolites. All data are represented the intensity values of metabolites. cThe value of F and Sig. was obtained from a one-way ANOVA. The F value is the ratio of the variance between the groups and the variance in the group. F value greater than 1 and Sig. less than 0.05 means the differences between groups were statistically significant. dp values were calculated from the Dunnett test in multivariate statistical analysis.

biomarkers. MataboAnalyst 3.0 is efficient tool for metabolites related metabolic pathway analysis, and it indicates that abnormal sphingolipid metabolism is a common characteristic of both model, while AA metabolism, linoleic acid metabolism and glycerophospholipid metabolism are highlighted in DS model.

It has been reported that CKD was closely related to the abnormal sphingolipid metabolism36, 37, and kidney injury may be caused by excessive accumulation of sphingolipids38-40. In this study, two important metabolites phytosphingosine and sphingosine 1-phosphate (S1P) involved in sphingolipid pathway could be considered as important biomarker in both models. Sphingolipids can regulate cell growth, differentiation, death and many important signal transduction pathways41. It was reported that phytosphingosine inhibited fibrosis by regulating the expression and activity of peroxisome proliferator-activated receptor y (PPARy)42. PPARy has various biological functions, such as anti-inflammatory, anti-fibrotic, regulation of lipid metabolism by regulating expression of transforming growth factor [3 (TGF-p) and nuclear factor kappa B (NF-kB), inhibition of the proliferation of mesangial cells and reduction in the accumulation of extracellular matrix and type I collagen43. In addition, S1P plays an important role in the progression of lung fibrosis. Ikeda's study found significant changes in the levels of S1P and its receptor in the liver tissues of patients with fibrosis44. Another report demonstrated that the role of S1P in fibrosis was correlated with the site of action and extracellular S1P promoted fibrotic processes in a S1P receptor-dependent manner, whereas intracellular S1P had an opposite effect45.

The biomarkers indicate abnormal sphingolipid metabolism in both models34. Compared with disease model, DS model includes additional TCM syndrome (qi-deficiency and blood-stasis, QDBS) on the base of abnormal serum biochemistry and pathological characteristics. Metabolic pathway analysis suggests DS model display more distinct abnormality in AA metabolism, linoleic acid metabolism and glycerophospholipid metabolism.

AA metabolic pathway includes AA together with its metabolites, LTE3 and 11b-PGF2a. AA mediates inflammation and the function of several organs and systems, either directly or upon its conversion into eicosanoids46 47. The metabolites that are involved in the regulation of renal epithelial transport and vascular function, such as EETs and 20-HETE can regulate transport function in the proximal tubules, which is important to maintain renal function in patients with kidney disease by affecting Na+-K+-ATPase activity48. LTE3 is a by-product of the metabolism of LTC3, a cysteinyl leukotriene. Cysteinyl LTs (leukotrienes C4, D4 and E4) possess bronchoc-onstrictive and inflammatory effects in the lung via their action on specific leukotriene receptors. Montelukast, a potent cysteinylleukotriene1 (CysLT1) receptor antagonist, is therapeutically effective for the alleviation of further progression of lung fibrosis via the inhibition of a-smooth muscle actin (a-SMA). It was found that LT receptor blockers prevented necroinflammatory liver injury and fibrogenesis by modulating the hepatic

Figure 6. Heatmap of 20 potential biomarkers between control and model groups. The color of each section is proportional to the significance of change of metabolites (red, upregulated; green, downregulated)as the numbers listed in Table 2. Rows: metabolites; Columns: samples (C: control group; DM: disease model group; DSM: DS model group).

Figure 7. Identification of potential metabolites related to disease model rats and DS model rats. The color represent different classes of biomarkers: (a) arachidonic acid metabolism, linoleic acid metabolism; (b) fatty acid biosynthesis; (c) sphingolipid metabolism; (d) amino acid metabolism; (e) glycerophospholipid metabolism. DHA: docosahexaenoic acid; LTE3: leukotriene E3; 15S-HETrE: 15(S)-hydroxyeicosatrienoic acid; DHEA: 8,11,14-eicosatrienoyl ethanolamide; EPA: eicosapentaenoic acid; AA: arachidonic acid; S1P: Sphingosine 1-phosphate.

expression of TGF-p, NF-kB, tissue inhibitor of metalloproteinase 1 and matrix metalloproteinases49. 11p-PGF2a is irreversibly produced from prostaglandin D2 via the enzyme prostaglandin-F synthase. Recent study demonstrated that PGD2 was probably involved in multiple aspects of inflammation through its dual receptor systems, DP and CRTH250.

Linoleic acid pathway involves five metabolites of 15S-HETrE, docosahexaenoic acid (DHA), eicosapentaenoic acid (EPA), 9-HODE and 8,11,14-eicosatrienoyl ethanolamide. Among the five metabolites, 15S-HETrE and 9-HODE play a role in regulating PPARy activity because 15S-HETrE can promote PPARy expression51 and 9-HODE is an endogenous agonist of PPARy52. 8,11,14-eicosatrienoyl ethanolamide is a DHA-derived N-acylethanolamine(NAE) metabolite, which modulates inflammation by reducing MCP-1 and NO production

Figure 8. Pathway analysis on biomarkers of disease model (a) and DS model (b). All matched pathways according to p-values from pathway enrichment analysis and pathway impact values from pathway topology analysis, using pathway library of Rattus norvegicus (rat).

Figure 9. Biosynthetic pathways for parts of potential biomarkers related to arachidonic acid metabolism, linoleic acid metabolism, glycerophospholipid metabolism and sphingolipid metabolism. The colored dots showed metabolites identified in the current study. DGLA: dihomo-^-linolenic acid; 15S-HETrE: 15(S)-hydroxyeicosatrienoic acid; AA: arachidonic acid; LTE3: leukotriene E3; EPA: eicosapentaenoic acid; DHA: docosahexaenoic acid; DHEA: 8, 11, 14-eicosatrienoyl ethanolamide; S1P: Sphingosine 1-phosphate.

and expression53. DHA and EPA both belong to the long-chain w-3 unsaturated fatty acid family, which maintains the normal physiological functions of cells and reduces blood cholesterol and triglycerides to improve blood microcirculation31, 54, 55. In addition, docosahexaenoic acid can suppress oxidative stress and fibrotic reaction in mice with non-alcohol fatty liver disease56. Moreover, EPA can regulate the ratio of urinary albumin/creatinine, as well as reduce the extracellular matrix accumulation of mesangial cells and renal tubulointerstitial fibrosis by adjusting monocyte chemoattractant protein-1 (MCP-1) levels57.

Lysophosphatidylcholines (LysoPCs) play a key role in process of glycerophospholipid metabolism, which have a pro-fibrotic effect because it can promote epithelial cell apoptosis, increase vascular permeability, promote fibroblast migration and anti-apoptosis, activate TGF-p activity and increase the secretion of the pro-fibrotic cytokines platelet-derived growth factor [3 and connective tissue growth factor in proximal tubule cells58-62. Moreover, LysoPCs can promote oxidative stress and enable cell structural and functional abnormalities that ultimately induce hemodynamic disorders by inhibiting Na+-K+-ATPase activity35, 63, 64.

Moreover, compared with disease model, some other particular biomarkers were also identifed in DS model. Phytosphingosine could regulate PPARy expression and activity42, and its reduced level was closely related to the progress of fibrosis. tri-N-acetyl chitotriose, a component of oligosaccharides, could protect mucosal and maintain normal digestive action, and the reduced level of tri-N-acetyl chitotriose in model rats may be correlated

with gastrointestinal dysfunction65. Estradiol is closely related to oxidative stress in patients with CKD as an anti-fibrosis component that inhibits the generation of reactive oxygen species and loss of anti-oxidant enzymes via the suppression of NADH/NADPH oxidase activity, attenuation of hydrogen peroxide induced TGF-p1 expression as well as proliferation and transformation of hepatic stellate cells66, 67.

The metabolomic study identified the similarities and differences of the two models. Common potential bio-markers related to abnormal sphingolipid metabolism were identified in disease model and DS model, which are associated with the formation of fibrotic cells or tissue and their cytokines, such as PPARy, TGF-p and MCP-1. Furthermore, compared with the disease model, DS model shows its differences on AA metabolism, linoleic acid metabolism and glycerophospholipid metabolism, which were mainly involved in the maintenance of cell structure and function, regulation of cytokines related to fibrosis and other aspects and moderation of oxidative stress and inflammation. Accumulated evidences have demonstrated the presence of inflammation, oxidative stress, lipid metabolism disorders and other symptoms in patients with CKD. Therefore, DS model is in more coincidence with clinical settings.

Conclusion

Disease model and DS model show no difference in clinical biochemistry, while metabolomics may highlight their difference. Abnormal sphingolipid metabolism is a common characteristic of both models, while AA metabolism, linoleic acid metabolism and glycerophospholipid metabolism are highlighted in DS model.

DS model not only involves the histopathological changes or abnormal biochemistry but can reveal potential biomarkers related to renal fibrosis, oxidative stress, inflammation and digestive disorders that are presented patients with CKD. This result is consistent with actual clinical situation. This consistence is beneficial for the translational medicine research and presents a greater advantage in the evaluation of disease and in the pharmacological study of drugs.

References

1. Jucker, M. The benefits and limitations of animal models for translational research in neurodegenerative diseases. Nature medicine 16, 1210-1214, doi:10.1038/nm.2224 (2010).

2. Jiang, G. T., Chen, X., Li, D., An, H. X. & Jiao, J. D. Ulinastatin attenuates renal interstitial inflammation and inhibits fibrosis progression in rats under unilateral ureteral obstruction. Molecular medicine reports 10, 1501-1508, doi:10.3892/mmr.2014.2323 (2014).

3. Perel, P. et al. Comparison of treatment effects between animal experiments and clinical trials: systematic review. Brit Med J 334, 197-200, doi:10.1136/bmj.39048.407928.BE (2007).

4. Wang, H. Z., Liu, X., Lv, B., Yang, F. & Hong, Y. Z. Reliable Multi-Label Learning via Conformal Predictor and Random Forest for Syndrome Differentiation of Chronic Fatigue in Traditional Chinese Medicine. Plos One 9, doi:10.1371/journal.pone.0099565 (2014).

5. Wang, J. & Xiong, X. J. Current Situation and Perspectives of Clinical Study in Integrative Medicine in China. Evid-Based ComplAlt, doi: 10.1155/2012/268542 (2012).

6. Liu, L. P. Differences in biochemical indexes of blood and relevant factors of obesity in hypertensive patients with different syndromes of traditional Chinese medicine. Chinese Journal of Clinical Rehabilitation 10, 12-13 (2006).

7. Chen, D. Q. et al. Metabolomic application in toxicity evaluation and toxicological biomarker identification of natural product. Chem Biol Interact 252, 114-130, doi:10.1016/j.cbi.2016.03.028 (2016).

8. Wang, M. et al. Metabolomics highlights pharmacological bioactivity and biochemical mechanism of traditional Chinese medicine. Chem Biol Interact 273, 133-141, doi:10.1016/j.cbi.2017.06.011 (2017).

9. Bates, J., Chitani, A. & Dreyer, G. Caring for patients with end-stage kidney disease. Lancet 386, 854-855 (2015).

10. Chen, L. et al. Role of RAS/Wnt/beta-catenin axis activation in the pathogenesis of podocyte injury and tubulo-interstitial nephropathy. Chem Biol Interact 273, 56-72, doi:10.1016/j.cbi.2017.05.025 (2017).

11. Lawson, J., Elliott, J., Wheeler-Jones, C., Syme, H. & Jepson, R. Renal fibrosis in feline chronic kidney disease: Known mediators and mechanisms of injury. Vet J 203, 18-26, doi:10.1016/j.tvjl.2014.10.009 (2015).

12. Goek, O. N. et al. Metabolites associate with kidney function decline and incident chronic kidney disease in the general population. Nephrol Dial Transpl 28, 2131-2138, doi:10.1093/ndt/gft217 (2013).

13. Rhee, E. P. et al. A Combined Epidemiologic and Metabolomic Approach Improves CKD Prediction. J Am Soc Nephrol 24, 1330-1338, doi:10.1681/Asn.2012101006 (2013).

14. Zhao, Y. Y. Metabolomics in chronic kidney disease. Clin Chim Acta422, 59-69, doi:10.1016/j.cca.2013.03.033 (2013).

15. Zhang, Z. H. et al. Metabolomic Signatures of Chronic Kidney Disease of Diverse Etiologies in the Rats and Humans. JProteome Res 15, 3802-3812, doi:10.1021/acs.jproteome.6b00583 (2016).

16. Hocher, B. & Adamski, J. Metabolomics for clinical use and research in chronic kidney disease. Nat Rev Nephrol 13, 269-284, doi:10.1038/nrneph.2017.30 (2017).

17. Zhao, Y. Y., Cheng, X. L., Wei, F., Bai, X. & Lin, R. C. Application of faecal metabonomics on an experimental model of tubulointerstitial fibrosis by ultra performance liquid chromatography/high-sensitivity mass spectrometry with MSE data collection technique. Biomarkers: biochemical indicators of exposure, response, and susceptibility to chemicals 17, 721-729, doi:10.3109/135475 0X.2012.724450 (2012).

18. Zhao, Y. Y. et al. Urinary metabonomics study on the protective effects of ergosta-4, 6, 8(14), 22-tetraen-3-one on chronic renal failure in rats using UPLC Q-TOF/MS and a novel MSE data collection technique. Process Biochem 47, 1980-1987, doi:10.1016/j. procbio.2012.07.008 (2012).

19. Zhao, Y. Y., Lei, P., Chen, D. Q., Feng, Y. L. & Bai, X. Renal metabolic profiling of early renal injury and renoprotective effects of Poria cocos epidermis using UPLC Q-TOF/HSMS/MSE. JPharmaceut Biomed 81-82, 202-209, doi:10.1016/j.jpba.2013.03.028 (2013).

20. Fang, J. W. et al. A urine metabonomics study of chronic renal failure and intervention effects of total aglycone extracts of Scutellaria baicalensis in 5/6 nephrectomy rats. Rsc Adv 5, 75612-75621, doi:10.1039/c5ra12710f (2015).

21. Zhang, Z. H. et al. Metabolomics insights into chronic kidney disease and modulatory effect of rhubarb against tubulointerstitial fibrosis. Sci Rep-Uk 5, doi: 10.1038/Srep14472 (2015).

22. Zhao, Y. Y. et al. Metabolomics analysis reveals the association between lipid abnormalities and oxidative stress, inflammation, fibrosis, and Nrf2 dysfunction in aristolochic acid-induced nephropathy. Sci Rep-Uk 5, doi: 10.1038/Srep12936 (2015).

23. Li, Y. H., Tan, Y. F., Cai, H. D. & Zhang, J. Q. Metabonomic study of the fruits of Alpinia oxyphylla as an effective treatment for chronic renal injury in rats. J Pharmaceut Biomed 124, 236-245, doi:10.1016/j.jpba.2016.02.035 (2016).

24. Zhao, Y. Y., Cheng, X. L., Lin, R. C. & Wei, F. Lipidomics applications for disease biomarker discovery in mammal models. Biomarkers in medicine 9, 153-168, doi:10.2217/bmm.14.81 (2015).

25. Qu, K. et al. [Common syndrome factors of traditional Chinese medicine in chronic renal failure based on the questionnaire investigation among experts]. Zhongxiyi jie he xue bao = Journal of Chinese integrative medicine 8, 1147-1152 (2010).

26. Zhao, Y. Y. & Lin, R. C. UPLC-MS(E) application in disease biomarker discovery: the discoveries in proteomics to metabolomics. Chem Biol Interact 215, 7-16, doi:10.1016/j.cbi.2014.02.014 (2014).

27. Miao, H. et al. Plasma lipidomics reveal profound perturbation of glycerophospholipids, fatty acids, and sphingolipids in diet-induced hyperlipidemia. Chem-Biol Interact 228, 79-87, doi:10.1016/j.cbi.2015.01.023 (2015).

28. Miao, H. et al. The antihyperlipidemic effect of Fu-Ling-Pi is associated with abnormal fatty acid metabolism as assessed by UPLC-HDMS-based lipidomics. Rsc Adv 5, 64208-64219, doi:10.1039/c5ra09766e (2015).

29. Miao, H. et al. Lipidomics Biomarkers of Diet-Induced Hyperlipidemia and Its Treatment with Poria cocos. Journal of agricultural andfood chemistry 64, 969-979, doi:10.1021/acs.jafc.5b05350 (2016).

30. Zhao, Y. Y., Liu, J., Cheng, X. L., Bai, X. & Lin, R. C. Urinary metabonomics study on biochemical changes in an experimental model of chronic renal failure by adenine based on UPLC Q-TOF/MS. Clin Chim Acta 413, 642-649, doi:10.1016/j.cca.2011.12.014 (2012).

31. Zhao, Y. Y. et al. Intrarenal Metabolomic Investigation of Chronic Kidney Disease and its TGF-beta 1 Mechanism in Induced-adenine Rats using UPLC Q-TOF/HSMS/MSE. JProteome Res 12, 692-703, doi:10.1021/pr3007792 (2013).

32. Xia, J., Sinelnikov, I. V., Han, B. & Wishart, D. S. MetaboAnalyst 3.0-making metabolomics more meaningful. Nucleic acids research 43, W251-257, doi:10.1093/nar/gkv380 (2015).

33. Xia, J. & Wishart, D. S. Using MetaboAnalyst 3.0 for Comprehensive Metabolomics Data Analysis. Current protocols in bioinformatics 55, 14 10 11-14 10 91, doi:10.1002/cpbi.11 (2016).

34. Chen, D. Q. et al. The link between phenotype and fatty acid metabolism in advanced chronic kidney disease. Nephrology, dialysis, transplantation: official publication of the European Dialysis and Transplant Association - European Renal Association, doi:10.1093/ ndt/gfw415 (2017).

35. Chen, D. Q. et al. Gene and protein expressions and metabolomics exhibit activated redox signaling and wnt/beta-catenin pathway are associated with metabolite dysfunction in patients with chronic kidney disease. Redox Biol 12, 505-521, doi:10.1016/j. redox.2017.03.017 (2017).

36. Xiang, Z., Sun, H., Cai, X. & Chen, D. The study on serum and urine of renal interstitial fibrosis rats induced by unilateral ureteral obstruction based on metabonomics and network analysis methods. Anal Bioanal Chem 408, 2607-2619, doi:10.1007/s00216-016-9368-4 (2016).

37. Zhao, Y. Y., Cheng, X. L., Vaziri, N. D., Liu, S. & Lin, R. C. UPLC-based metabonomic applications for discovering biomarkers of diseases in clinical chemistry. Clinical biochemistry 47, 16-26, doi:10.1016/j.clinbiochem.2014.07.019 (2014).

38. Chen, H. et al. Combined Clinical Phenotype and Lipidomic Analysis Reveals the Impact of Chronic Kidney Disease on Lipid Metabolism. J Proteome Res 16, 1566-1578, doi:10.1021/acs.jproteome.6b00956 (2017).

39. Zhao, Y. Y., Wu, S. P., Liu, S. M., Zhang, Y. M. & Lin, R. C. Ultra-performance liquid chromatography-mass spectrometry as a sensitive and powerful technology in lipidomic applications. Chem-Biol Interact 220, 181-192, doi:10.1016/j.cbi.2014.06.029 (2014).

40. Zhao, Y. Y. et al. A Pharmaco-Metabonomic Study on Chronic Kidney Disease and Therapeutic Effect of Ergone by UPLC-QTOF/ HDMS. Plos One 9, doi: 10.1371/journal.pone.0115467 (2014).

41. Kurek, K. et al. Metabolism, physiological role, and clinical implications of sphingolipids in gastrointestinal tract. BioMed research international 2013, 908907, doi:10.1155/2013/908907 (2013).

42. Kim, S. et al. Phytosphingosine stimulates the differentiation of human keratinocytes and inhibits TPA-induced inflammatory epidermal hyperplasia in hairless mouse skin. Mol Med 12, 17-24, doi:10.2119/2006-00001.Kim (2006).

43. Ruan, X., Zheng, F. & Guan, Y. PPARs and the kidney in metabolic syndrome. American journal of physiology. Renal physiology 294, F1032-1047, doi:10.1152/ajprenal.00152.2007 (2008).

44. Ikeda, H. et al. Plasma concentration of bioactive lipid mediator sphingosine 1-phosphate is reduced in patients with chronic hepatitis C. Clin Chim Acta 411, 765-770, doi:10.1016/j.cca.2010.02.063 (2010).

45. Schwalm, S., Pfeilschifter, J. & Huwiler, A. Sphingosine-1-phosphate: a Janus-faced mediator of fibrotic diseases. Biochimica et biophysica acta 1831, 239-250, doi:10.1016/j.bbalip.2012.07.022 (2013).

46. Harizi, H., Corcuff, J. B. & Gualde, N. Arachidonic-acid-derived eicosanoids: roles in biology and immunopathology. Trends Mol Med 14, 461-469, doi:10.1016/j.molmed.2008.08.005 (2008).

47. Pompeia, C., Lima, T. & Curi, R. Arachidonic acid cytotoxicity: can arachidonic acid be a physiological mediator of cell death? Cell Biochem Funct 21, 97-104, doi:10.1002/cbf.1012 (2003).

48. Zhao, X. & Imig, J. D. Kidney CYP450 enzymes: Biological actions beyond drug metabolism. Curr DrugMetab 4, 73-84, doi:10.2174/1389200033336892 (2003).

49. El-Swefy, S. & Hassanen, S. I. Improvement of hepatic fibrosis by leukotriene inhibition in cholestatic rats. Annals of hepatology 8, 41-49 (2009).

50. Hirai, H. [Prostaglandin D2 in allergy: PGD2 has dual receptor systems]. Nihonyakurigaku zasshi. Foliapharmacologica Japonica 123, 15-22 (2004).

51. Pham, H., Banerjee, T., Nalbandian, G. M. & Ziboh, V. A. Activation of peroxisome proliferator-activated receptor (PPAR)-gamma by 15S-hydroxyeicosatrienoic acid parallels growth suppression of androgen-dependent prostatic adenocarcinoma cells. Cancer Lett 189, 17-25, doi:10.1016/S0304-3835(02)00498-6 (2003).

52. Hampel, J. K. A. et al. Differential modulation of cell cycle, apoptosis and PPAR gamma 2 gene expression by PPAR gamma agonists ciglitazone and 9-hydroxyoctadecadienoic acid in monocytic cells. Prostag Leukotr Ess 74, 283-293, doi:10.1016/j.plefa.2006.03.002 (2006).

53. Meijerink, J. et al. The ethanolamide metabolite of DHA, docosahexaenoylethanolamine, shows immunomodulating effects in mouse peritoneal and RAW264.7 macrophages: evidence for a new link between fish oil and inflammation. The British journal of nutrition 105, 1798-1807, doi:10.1017/S0007114510005635 (2011).

54. Riediger, N. D., Othman, R. A., Suh, M. & Moghadasian, M. H. A Systemic Review of the Roles of n-3 Fatty Acids in Health and Disease. J Am Diet Assoc 109, 668-679, doi:10.1016/j.jada.2008.12.022 (2009).

55. Ruxton, C. H. S., Reed, S. C., Simpson, M. J. A. & Millington, K. J. The health benefits of omega-3 polyunsaturated fatty acids: a review of the evidence. J Hum Nutr Diet 20, 275-285 (2007).

56. Depner, C. M., Philbrick, K. A. & Jump, D. B. Docosahexaenoic Acid Attenuates Hepatic Inflammation, Oxidative Stress, and Fibrosis without Decreasing Hepatosteatosis in a Ldlr(-/-) Mouse Model of Western Diet-Induced Nonalcoholic Steatohepatitis. J Nutr 143, 315-323, doi:10.3945/jn.112.171322 (2013).

57. Furukawa, M., Gohda, T., Tanimoto, M. & Tomino, Y. Pathogenesis and Novel Treatment from the Mouse Model of Type 2 Diabetic Nephropathy. Sci World J, doi: 10.1155/2013/928197 (2013).

58. Zhang, Z. H. et al. An integrated lipidomics and metabolomics reveal nephroprotective effect and biochemical mechanism of Rheum officinale in chronic renal failure. Sci Rep-Uk 6, doi:10.1038/Srep22151 (2016).

59. Zhao, Y. Y. et al. Ultra performance liquid chromatography-based metabonomic study of therapeutic effect of the surface layer of Poria cocos on adenine-induced chronic kidney disease provides new insight into anti-fibrosis mechanism. Plos One 8, e59617, doi:10.1371/journal.pone.0059617 (2013).

60. Zhao, Y. Y. et al. Serum metabonomics study of adenine-induced chronic renal failure in rats by ultra performance liquid chromatography coupled with quadrupole time-of-flight mass spectrometry. Biomarkers: biochemical indicators of exposure, response- and susceptibility to chemicals 17, 48-55, doi:10.3109/1354750X.2011.637180 (2012).

61. Zhao, Y. Y. et al. Effect of ergosta-4, 6, 8(14), 22-tetraen-3-one (ergone) on adenine-induced chronic renal failure rat: a serum metabonomic study based on ultra performance liquid chromatography/high-sensitivity mass spectrometry coupled with MassLynx i-FIT algorithm. Clin Chim Acta 413, 1438-1445, doi:10.1016/j.cca.2012.06.005 (2012).

62. Shea, B. S. & Tager, A. M. Role of the lysophospholipid mediators lysophosphatidic acid and sphingosine 1-phosphate in lung fibrosis. Proceedings of the American Thoracic Society 9, 102-110, doi:10.1513/pats.201201-005AW (2012).

63. Chen, H. et al. Metabolomics insights into activated redox signaling and lipid metabolism dysfunction in chronic kidney disease progression. Redox Biol 10, 168-178, doi:10.1016/j.redox.2016.09.014 (2016).

64. Yasunari, K., Maeda, K., Minami, M. & Yoshikawa, J. HMG-CoA reductase inhibitors prevent migration of human coronary smooth muscle cells through suppression of increase in oxidative stress. Arterioscl Throm Vas 21, 937-942 (2001).

65. Silano, M. & De Vincenzi, M. Bioactive antinutritional peptides derived from cereal prolamins: A Review. Nahrung 43, 175-184, doi:10.1002/(Sici)1521-3803(19990601)43:3<175::Aid-Food175>3.0.Co;2-Z (1999).

66. Itagaki, T. et al. Opposing effects of oestradiol and progesterone on intracellular pathways and activation processes in the oxidative stress induced activation of cultured rat hepatic stellate cells. Gut 54, 1782-1789, doi:10.1136/gut.2005.053278 (2005).

67. Shimizu, I. et al. Inhibitory effect of oestradiol on activation of rat hepatic stellate cells in vivo and in vitro. Gut 44, 127-136 (1999).

Acknowledgements

We are very grateful to Dr. Min Huang in Tsinghua University and Dr. Hong-zhi Zhao, Dr. Jing Wang, M.D. Sheng-you Song in Nankai University for their comments and contributions. This work was supported by the National Natural Science Foundation of China (no. 81072784). The authors are also grateful for this financial support.

Author Contributions

N.Z.Y., G.A.L. and Y.M.W. conceived the idea, provided advice and interpretation of the analysis. S.S.L., P.X., L.H. and W.M. performed the survey and analysis of clinical biochemistry and kidney histopathology. S.S.L. carried out the metabonomic experiments and wrote the manuscript. All authors reviewed the manuscript.

Additional Information

Competing Interests: The authors declare that they have no competing interests.

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

l/JJv ® I Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2017