Metabonomics screening of serum identifies pyroglutamate as a diagnostic biomarker for nonalcoholic steatohepatitis

Metabonomics screening of serum identifies pyroglutamate as a diagnostic biomarker for nonalcoholic steatohepatitis

Accepted Manuscript Metabonomics screening of serum identifies pyroglutamate as a diagnostic biomarker for nonalcoholic steatohepatitis Suwen Qi, Dep...

831KB Sizes 4 Downloads 49 Views

Accepted Manuscript Metabonomics screening of serum identifies pyroglutamate as a diagnostic biomarker for nonalcoholic steatohepatitis

Suwen Qi, Depeng Xu, Qiaoliang Li, Ni Xie, Jun Xia, Qin Huo, Pu Li, Qiwen Chen, Si Huang PII: DOI: Reference:

S0009-8981(17)30324-8 doi: 10.1016/j.cca.2017.08.022 CCA 14842

To appear in:

Clinica Chimica Acta

Received date: Revised date: Accepted date:

30 January 2017 11 August 2017 20 August 2017

Please cite this article as: Suwen Qi, Depeng Xu, Qiaoliang Li, Ni Xie, Jun Xia, Qin Huo, Pu Li, Qiwen Chen, Si Huang , Metabonomics screening of serum identifies pyroglutamate as a diagnostic biomarker for nonalcoholic steatohepatitis, Clinica Chimica Acta (2017), doi: 10.1016/j.cca.2017.08.022

This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

ACCEPTED MANUSCRIPT Title page Article title Metabonomics screening of serum identifies pyroglutamate as a diagnostic biomarker for nonalcoholic steatohepatitis

PT

Running title

RI

Metabolic biomarker of NASH

SC

Authors

Suwen Qia, Depeng Xua, Qiaoliang Li*, Ni Xie*b, Jun Xia b, Qin Huo b, Pu Lic, Qiwen

NU

Chena, Si Huanga

a

MA

Affiliations and addresses

Department of Biomedical Engineering, School of Medicine, Shenzhen University,

Shenzhen Second Hospital, The First Affiliated hospital of Shenzhen University,

PT E

b

D

Shenzhen 518060, China

Shenzhen 518060, P.R.China

Department of Laboratory Medicine, The Second Hospital Affiliated to Chongqing

CE

c

AC

Medical University, 410006, Chongqing, China *Corresponding author Name: Qiaoliang Li, PhD Email: [email protected] Phone & fax: +86 755 2562 6750 Address: Shenzhen University, Shenzhen, Guangdong Province, 518020, P.R. China

1

ACCEPTED MANUSCRIPT Abstract Objective A key step in managing non-alcoholic fatty liver disease (NAFLD) is to differentiate nonalcoholic steatohepatitis (NASH) from simple steatosis (SS). Method Serum samples were collected from three groups: NASH patients (N = 21),

PT

SS patients (N = 38) and healthy controls (N = 31). High performance liquid

RI

chromatography-mass spectrometry (HPLC-MS) was used to analyse the metabolic

SC

profile of the serum samples. The acquired data were processed by multivariate principal component analysis (PCA) and orthogonal partial least-squares-discriminant

NU

analysis (OPLS-DA) to identify novel metabolites. The potential biomarkers were

MA

quantitatively determined and their diagnostic power was further validated. Results A total of 56 metabolites were capable of distinguishing NASH from SS

D

samples based on the OPLS-DA model. Pyroglutamate was found to be the most

PT E

promising factor in distinguishing the NASH from SS groups. With an optimal cut-off value of 4.82 mmol/L, the sensitivity and specificity of the diagnosis of NASH were

CE

72% and 85%, respectively. The area under the receiver operating characteristic

AC

(AUROC) of the pyroglutamate levels of NASH versus SS patients was more than those of tumor necrosis factor-α, adiponectin and interleukin-8. Conclusion These data suggest that pyroglutamate may be a new and useful biomarker for the diagnosis of NASH. Key words: Non-alcoholic fatty liver disease; nonalcoholic steatohepatitis; simple steatosis; metabolic profiling; pyroglutamate

2

ACCEPTED MANUSCRIPT 1. Instruction Non-alcoholic fatty liver disease (NAFLD) is a spectrum of diseases that covers simple steatosis (SS) and the more severe progressive form of NAFLD—non-alcoholic steatohepatitis (NASH) [1]. NASH is a hallmark of NAFLD,

PT

which shows typical pathologic characteristics, such as lobular inflammation,

RI

hepatocyte clarification or ballooning and fibrosis progression. Hence, early diagnosis

SC

of NASH is extremely important for the treatment and prognosis of NAFLD patients [1-3].

NU

The actual gold standard for diagnosing NASH is to conduct a histologic evaluation

MA

of an invasive liver biopsy [4]. However, it is not always logistically feasible due to the high risk of complications. Imaging techniques and serum biomarkers are two

D

non-invasive approaches for the diagnosis of NASH. Ultrasound, computed

PT E

tomography (CT) and magnetic resonance imaging (MRI) are commonly employed imaging techniques for qualitative or quantitative assessments of hepatic steatosis [5,

CE

6]. Controlled attenuation parameter (CAP), which is a high accuracy ultrasound

AC

based technique, is the most promising method for NASH diagnosis, but it needs to be implemented with a probe in patients [7]. CT is inappropriate for assessing hepatic steatosis as it suffers from poor accuracy in detecting mild hepatic steatosis. MRI-proton density fat fraction (PDFF) is more accurate and sensitive, though the high costs associated with this technique limit its availability [8]. The biomarker-based liver fat score of NAFLD is the preferred method for quantification identification of steatosis. Serum biomarkers often are analyzed in 3

ACCEPTED MANUSCRIPT panels to establish live fat score models of NAFLD. SteatoTest TM (Biopredictive, Paris, France) is a score algorithm for the quantification of steatosis. The score ranges from 0 to 1 based on the values of twelve parameters that can be determined via biochemical assays, such as alanine transaminase [9], haptoglobin, total bilirubin,

PT

gamma-glutamyl transferase (GGT) and body mass index (BMI) [10]. Meta-analysis

RI

has shown that the SteatoTest TM has an area under receiver operating characteristic

SC

curve (AUROC) of 0.80 for diagnosing individuals with more than 33% steatosis [11]. Some biomarkers that correlate with NAFLD, such as tumor necrosis factor-α

NU

(TNF-α), adiponectin and Interleukin-8 (IL-8), had the same diagnostic score.

MA

However, no panel or single biomarker could effectively discriminate NASH from SS [12-14]. The aim of the proposed metabonomics study was to identify highly sensitive

PT E

diagnosis of NASH.

D

and specific biomarkers that are convenient and applicable for the differential

2. Materials and Methods

CE

2.1 Patients and controls

AC

A total of 59 participants (38 SS and 21 NASH) were enrolled from the Department of Gastroenterology and Infectious Diseases of the Second Hospital Affiliated to Chongqing Medical University from November 2015 to March 2016. NASH and SS were diagnosed based on the results of histological examinations or the combined results of clinical, laboratory and imaging examinations. Patients with known causes of chronic liver disease, including viral hepatitis (HBV and HCV), Wilson disease, haemochromatosis, or autoimmune hepatitis, were excluded. Patients with a history of 4

ACCEPTED MANUSCRIPT alcoholic consumption > 20 g/d for a period of > 5 years or > 50 g/d for a period of > 2 years were excluded [15]. Thirty-one healthy volunteers served as the controls. The baseline clinical characteristics of the NAFLD patients and controls are summarized in Table 1. Approximately 5 mL of fasted peripheral venous blood from patients and

PT

controls was collected, and then gently mixed by inverting the tube 8 to 10 times.

RI

After 30-40 min at room temperature for a clot to form, the whole blood were

SC

centrifuged 15 min at 3000 rpm at room temperature. One thousand microliters of serum were pipetted carefully into labeled tubes and stored in -80°C freezer. The

NU

procedure was conducted by the same personnel and the whole operation should be

MA

completed within 1 hr. If hemolysis or lipidemia was observed, the sample should be excluded.

D

This study was performed according to the guidelines of Ethics Committee of

PT E

Chongqing Medical University (Chongqing, China; approval No. CMU-2016-07) , which abided by the Declaration of Helsinki on ethical principles for medical

CE

research involving human subjects.

AC

2.2 HPLC-MS analysis

Approximately 500 µL of serum was mixed with 800 µL of cold methanol and precipitated at 12,000 RPM at 4°C for 15 min. Two-hundred microliters of supernatant was analyzed using high performance liquid chromatograph-mass spectrometer (HPLC-MS, Ultimate 3000 LC, Orbitrap Elite; Thermo Fisher Scientific, Waltham, MA, USA) analysis. The mass spectrum (MS) parameters in positive ion mode (ESI+ MODE) were as follows: Heater Temp, 300°C; Sheath Gas Flow rate, 45 5

ACCEPTED MANUSCRIPT arb; Aux Gas Flow Rate, 15 arb; Sweep Gas Flow Rate, 1arb; spray voltage, 3.0 KV; Capillary Temp, 350°C; and S-Lens RF Level, 30%. Standard controls for pyroglutamate were purchased from Sigma (St. Louis, MO, USA) and were used to prepare calibration curves.

PT

HPLC-MS data were preprocessed using SIEVE software (Thermo Fisher Scientific)

RI

and then normalized and edited in EXCEL 2010 software (Microsoft Corporation,

SC

Redmond, WA, USA). The results were visualized as a 2D data matrix, including metabolites’ retention time/mass-to-charge ratio (rt/mz), molecular weights

NU

(compMW) and peak intensities.

MA

2.3 Pattern recognition analysis

The metabolic profiles of the serum were first examined by unsupervised principal

D

component analysis (PCA) to obtain the score plots of NASH, SS and healthy controls

PT E

[16]. After an initial overview of the PCA analysis, orthogonal partial least-squares-discriminant analysis (OPLS-DA) was performed to develop a more

CE

sophisticated model with the specific discriminant information between the different

AC

groups [17, 18]. Variable importance in the projection (VIP) values of the OPLS-DA model (threshold value > 1) combined with t-test p values (p < 0.05) were adopted to find differences in the expression of metabolites. The metabolites were assigned according to the mass-to-charge ratio found in an online database (https://metlin.scripps.edu). 2.4 Measurements of biomarkers The quantitative determination of pyroglutamate was tested by HPLC-MS. 6

ACCEPTED MANUSCRIPT Calibration curves for pyroglutamate were prepared at different concentrations from 0.1 mmol/L to 20 mmol/L. TNF-α and IL-8 were quantified using a Luminex 200 (Austin, Texas, USA). The adiponectin levels were measured using an ELISA assays kit from Phoenix Pharmaceuticals, Inc. (Belmont, CA, USA).

PT

2.5 Data analysis

RI

Statistical analysis was performed using SPSS11.5 software (SPSS Company,

SC

Chicago, Illinois, USA). The data were compared using the t-test, and the threshold P value was set at 0.05 throughout the study.

NU

3. Results

MA

The HPLC-MS signals of all common metabolites were identified, and the extracted data were subsequently analyzed using multivariate statistics, including PCA and

D

OPLS-DA. From all samples, we obtained 1761 features in ESI+ MODE. Mass

PT E

chromatograms of serum obtained from control, simple steatosis nonalcoholic steatohepatitis were displayed in Figure 1. The edited data matrix was imported into

CE

SIMCA-P software (v13.0) and subsequently analyzed using multivariate statistics.

AC

3.1 Principal component analysis PCA was performed on samples to determine the overall metabolic differences between groups. In the PCA score plot, a single point represented a serum metabolic profile. Figure 2A displays the two-dimensional PCA score plots of controls and SS. The first two principal components (PC1 and PC2) exhibit the most efficient features contained in the data set. Twelve principal components were obtained with a R2X of 0.734 and a Q2 Y of 0.235. In general, for the unsupervised model, an R2X value 7

ACCEPTED MANUSCRIPT greater than 0.4 means that it is reliable [19]. Similarly, Figure 2B shows the PCA score plots of SS and NASH with an R2X of 0.694 and a Q2Y of 0.258. 3.2 OPLS-DA OPLS-DA is a supervised pattern recognition algorithm that is more focused on

PT

discriminatory variations. Healthy control and SS samples were discriminated with an

RI

R2 X of 0.306, R2 Y of 0.844, and Q2 Y of 0.717(Figure 3A). SS and NASH were

SC

discriminated with an R2 X of 0.324, R2 Y of 0.972, and Q2 Y of 0.956 (Figure 3B). According to a threshold value of VIP greater than 1, a total of 108 metabolites were

NU

capable of distinguishing NAFLD from healthy controls and 56 metabolites were

MA

capable of distinguishing NASH from SS. The metabolic profiles showed that the main altered metabolites were glutamate, kynurenic acid, among others. To screen

D

remarkable metabolites, we applied a stricter criterion that stated that the VIP of

PT E

metabolites should be greater than 1.3. Five metabolites, including uracil, α-linolenic acid, glutamate, L-glutamine and pyroglutamate, were identified as prominent factors

CE

in discriminating NASH from SS. The fold changes of pyroglutamate between SS

AC

versus control and NASH versus SS were 1.56 and 2.26, respectively, which were higher than those of the other metabolites (Table 2). 3.3 Biomarkers levels in serum A standard curve of y= 3366.8x + 5826.7 (peak area vs concentration) for pyroglutamate was obtained with a correlation coefficient of 0.9779. Two different levels of the sera (3.57 mmol/L, 9.61 mmol/L) were measured twenty times to evaluate the imprecision. The variable coefficient (CV) were 4.53% and 3.27%, 8

ACCEPTED MANUSCRIPT respectively. A known amount of pyroglutamate (15.0 mmol/L) was added to two human sera at different levels as spike-and-recovery controls. The average recoveries were 97.56% and 102.35%. We selected 25 control, 25 SS and 25 NASH samples to quantitatively determine the

PT

levels of pyroglutamate, TNF-α, IL-8 and adiponectin in serum. The results showed

RI

that the level of pyroglutamate was highest in the NASH, followed by the SS and

SC

control samples.

The median pyroglutamate level in the NAFLD (including NASH and SS) group was

NU

5.19 mmol/L (95% CI 3.89–8.19 mmol/L), which was dramatically higher than that of

MA

the control group, 1.48 mmol/L (95% CI 0.48–3.97 mmol/L) (Figure 4A). ROC analysis of the pyroglutamate levels of NAFLD samples versus controls yielded

D

AUROC of 0.880 (Figure 4B). With an optimal cut-off value of 3.20 mmol/L, the

PT E

sensitivity and specificity for the diagnosis of NAFLD were 73% and 87%, respectively.

CE

The pyroglutamate levels in the NASH group were 7.89 mmol/L (95% CI 5.35–10.41

AC

mmol/L), which was higher than that of the SS group at 3.57 mmol/L (95% CI 0.98– 7.34 mmol/L). With a cut-off value of 4.82 mmol/L, ROC analysis yielded an AUROC of 0.846; the sensitivity and specificity for the differential diagnosis of NASH from SS were 72% and 85%, respectively. The level of adiponectin, IL-8 and TNF-α in patients with SS and NASH are shown in Table 3. Groupwise comparisons showed decreased adiponectin (Se 64%, Sep 56%, AUROC 0.723), decreased IL-8 (Se 65%, Sep 68%, AUROC 0.757) and increased 9

ACCEPTED MANUSCRIPT TNF-α (Se 72%, Sep 76%, AUROC 0.810), which were associated with NASH. The ROC curves of pyroglutamate, adiponectin, IL-8 and TNF-α of NASH versus SS samples are shown in Figure 5. In summary, pyroglutamate performed better than the others as a single diagnostic biomarker of NASH. A combination of 2-3 biomarkers is

PT

usually used in clinical lab to improve diagnosis sensitivity. The diagnostic

RI

performances of pyroglutamate combination with adiponectin, IL-8, TNF-α are

SC

summarized in Table 4. In addition, the combination of pyroglutamate with standard laboratory parameters such as GGT and triglycerides also are displayed in Table 4.

NU

4. Discussion

MA

Distinguishing NASH from SS is a key challenge in the management of NAFLD patients. Serum biomarkers are the most available noninvasive metrics for the

D

diagnosis of NASH [20]. In this study, we studied the serum metabolic profiles of

PT E

NAFLD based on HPLC-MS to identify NASH biomarkers. PCA score plots clearly distinguished the control, SS and NASH groups. As a highly sensitive analytical

CE

method, OPLS-DA, identified the altered metabolites in the NASH group compared

AC

with the SS group. Five metabolites, including uracil, α-Linolenic acid, glutamate, L-glutamine and pyroglutamate, were identified as the most prominent biomarkers. The fold changes of pyroglutamate between SS and control and NASH and SS were 1.56 and 2.26, respectively, which were higher than those of the other biomarkers. Pyroglutamate is also known as 5-oxoproline. The metabolism of pyroglutamate is related to the gamma-glutamyl cycle, which is the primary pathway for glutathione synthesis and degradation [21, 22]. The level of pyroglutamate is often used to screen 10

ACCEPTED MANUSCRIPT for glutathione synthetase deficiency in pediatric patients [23]. Children with deficiencies of either glutathione synthase or 5-oxoprolinase will present with metabolic acidosis or other abnormalities due to the accumulation of pyroglutamic acid [24-26]. When glycine is insufficiently available to support glutathione recovery

PT

via the renal gamma-glutamyl cycle, the rising level of pyroglutamate is regarded as

RI

an index of glycine insufficiency. In addition, the components of gamma-glutamyl

SC

cycle such as glutathione (CompMW: 307.08, retention time: 0.94) and cysteine (CompMW: 283.28, retention time: 7.23) were identified in the LC-MS data.

NU

Although both were capable of distinguishing between SS or NASH group and

MA

controls, glutathione and cysteine did not show significant changing between NASH and SS.

D

Furthermore, pyroglutamate has been proposed as a potential biomarker for enhanced

PT E

oxidative stress in liver cells [27, 28]. Studies have shown enhanced levels of pyroglutamate in various biofluids or tissues following administration of hepatic

CE

toxicants, such as acetaminophen [29] and bromobenzene [30]. A growing number of

AC

reports have documented increased pyroglutamate in the blood, urine, or both in cases of metabolic acidosis complicating paracetamol poisoning [31]. However, there is no study focusing on the role of pyroglutamate in the progress of NAFLD. The metabolic profiling conducted in this study showed that pyroglutamate is a remarkable metabolite that is correlated with NAFLD. Its sensitivity and specificity are equal to those of SteatoTestTM, which is a steatosis score algorithm based on 12 clinical parameters. The AUROC value (0.880) for the diagnosis of NAFLD was more 11

ACCEPTED MANUSCRIPT than 0.80 of that of SteatoTestTM [11]. Pyroglutamate profiling has the capacity to identify NAFLD from healthy controls. Nevertheless, more should be done to discriminate patients with SS and NASH because both are very similar. The performance of pyroglutamate in diagnosing

PT

NASH was carefully compared with that of adiponectin, IL-8 and TNF-α, which are

RI

correlated with NASH. Low adiponectin is a characteristic of steatosis progression

SC

and leads to an increase of intra-cellular fatty acids [12, 32]. In this study, the level of adiponectin in the SS group was higher than that in controls, whereas there was no

NU

significant difference between the NASH and SS groups. IL-8 is an independent

MA

factor that is associated with NASH and was noticeably decreased in NASH compared to SS samples. The levels of TNF-α increased gradually from healthy

D

controls to SS, then to NASH. Although adiponectin, IL-8 and TNF-α showed some

PT E

diagnostic significance, the results showed that their sensitivity and specificity were inferior to those of pyroglutamate.

CE

In addition to using single serum biomarkers related to NASH, some non-invasive

AC

models based on several biomarkers have been used to diagnosis NASH. The fatty liver index (FLI) is a routinely used tool based on four components, including BMI, waist circumference, triglycerides, and GGT [33, 34]. However, the unreliability of the components limited the interpretation of FLI [35]. Munkhzul et al. [33] established an index of NASH (ION) model based on anthropometric, clinical, ultrasound, laboratory, and mortality data, and its sensitivity (81%) was better, but its specificity (82%) was less than the pyroglutamate biomarker. In consideration of the 12

ACCEPTED MANUSCRIPT cost and operating performance, pyroglutamate still has an advantage in diagnosing NASH. The panel of pyroglutamate with other biomarkers contributed to enhance the diagnostic performance. The sensitivity of combination of pyroglutamate, adiponectin and TNF-α reached 92%.

PT

The difference in expression of pyroglutamate between the SS and NASH groups may

RI

be a reflection of their role in disease pathogenesis. It was reported that increased

SC

pyroglutamate excretion is an indication of abnormal response to oxidative stress [36]. NAFLD is accompanied by broad effects on mitochondria, which is the main location

NU

for fatty acid oxidation. When SS progresses to NASH, mitochondria have impaired

MA

respiratory coupling, which may increase the reductive load on respiration and form reactive oxygen species (ROS);ROS damage hepatocytes, trigger inflammation, and

D

elevate pyroglutamate levels [37, 38] .In conclusion, we identified pyroglutamate to

PT E

be a novel biomarker for the diagnosis of NASH. Future studies with larger sample sizes in clinical settings are worth pursuing.

CE

Acknowledgements

AC

This work was supported by the Project of the National Science Foundation of China (61401285

and

81401750),

Natural

Science

Foundation

of

Guangdong

(2016A030313029), Scientific Plan Program of Guangdong (2014A020212038), Innovation Program of Shenzhen (JCYJ20150330102720122), and International Cooperation Foundation of Shenzhen (GJHZ20160301163138685), Shenzhen Science Plan

(JCYJ2016030711492524,

JCYJ20150330102720122,

CXZZ201418182638764) and The Fund of Shenzhen University (2015006). 13

and

ACCEPTED MANUSCRIPT References [1]

V. Ajmera, E.R. Perito, N.M. Bass, N.A. Terrault, K.P. Yates, R. Gill, et al.,Novel plasma biomarkers associated with liver disease severity in adults with nonalcoholic fatty liver disease, Hepatology. 65(2017) 65-77. E. Bugianesi, C. Rosso, and H. Cortez-Pinto,How to diagnose NAFLD in

PT

[2]

M.E. Rinella and A.J. Sanyal,Management of NAFLD: a stage-based approach,

SC

[3]

RI

2016, J Hepatol. 65(2016) 643-4.

Nat Rev Gastroenterol Hepatol. 13(2016) 196-205.

C. Stern and L. Castera,Non-invasive diagnosis of hepatic steatosis, Hepatol

NU

[4]

[5]

MA

Int. 2016)

S. Kinner, S.B. Reeder, and T. Yokoo,Quantitative Imaging Biomarkers of

J. Pirhonen, J. Arola, S. Sadevirta, P. Luukkonen, S.M. Karppinen, T.

PT E

[6]

D

NAFLD, Dig Dis Sci. 61(2016) 1337-47.

Pihlajaniemi, et al.,Continuous Grading of Early Fibrosis in NAFLD Using

CE

Label-Free Imaging: A Proof-of-Concept Study, PLoS One. 11(2016)

[7]

AC

e0147804.

R.P. Myers, A. Pollett, R. Kirsch, G. Pomier-Layrargues, M. Beaton, M. Levstik, et al.,Controlled Attenuation Parameter (CAP): a noninvasive method for the detection of hepatic steatosis based on transient elastography, Liver Int. 32(2012) 902-10.

[8]

E. Hashimoto, M. Taniai, and K. Tokushige,Characteristics and diagnosis of NAFLD/NASH, J Gastroenterol Hepatol. 28 Suppl 4(2013) 64-70. 14

ACCEPTED MANUSCRIPT [9]

J. Vulevic, A. Juric, G.E. Walton, S.P. Claus, G. Tzortzis, R.E. Toward, et al.,Influence of galacto-oligosaccharide mixture (B-GOS) on gut microbiota, immune parameters and metabonomics in elderly persons, Br J Nutr. 114(2015) 586-95. T. Poynard, V. Ratziu, S. Naveau, D. Thabut, F. Charlotte, D. Messous, et

PT

[10]

RI

al.,The diagnostic value of biomarkers (SteatoTest) for the prediction of liver

[11]

SC

steatosis, Comp Hepatol. 4(2005) 10.

T. Poynard, G. Lassailly, E. Diaz, K. Clement, R. Caiazzo, J. Tordjman, et

NU

al.,Performance of biomarkers FibroTest, ActiTest, SteatoTest, and NashTest

MA

in patients with severe obesity: meta analysis of individual patient data, PLoS One. 7(2012) e30325.

L. Abenavoli and V. Peta,Role of adipokines and cytokines in non-alcoholic

D

[12]

[13]

PT E

fatty liver disease, Rev Recent Clin Trials. 9(2014) 134-40. M.H. Jarrar, A. Baranova, R. Collantes, B. Ranard, M. Stepanova, C. Bennett,

CE

et al.,Adipokines and cytokines in non-alcoholic fatty liver disease, Aliment

[14]

AC

Pharmacol Ther. 27(2008) 412-21. S. Abiru, K. Migita, Y. Maeda, M. Daikoku, M. Ito, K. Ohata, et al.,Serum cytokine and soluble cytokine receptor levels in patients with non-alcoholic steatohepatitis, Liver Int. 26(2006) 39-45. [15]

V.R. Mas, R. Fassnacht, K.J. Archer, and D. Maluf,Molecular mechanisms involved in the interaction effects of alcohol and hepatitis C virus in liver cirrhosis, Mol Med. 16(2010) 287-97. 15

ACCEPTED MANUSCRIPT [16]

B. Worley and R. Powers,PCA as a practical indicator of OPLS-DA model reliability, Curr Metabolomics. 4(2016) 97-103.

[17]

S. Wagner, K. Scholz, M. Donegan, L. Burton, J. Wingate, and W. Volkel,Metabonomics and biomarker discovery: LC-MS metabolic profiling

PT

and constant neutral loss scanning combined with multivariate data analysis

Y. Wang, H. Man, J. Gao, X. Liu, X. Ren, J. Chen, et al.,Plasma metabonomics

SC

[18]

RI

for mercapturic acid analysis, Anal Chem. 78(2006) 1296-305.

study on toxicity biomarker in rats treated with Euphorbia fischeriana based

S. Qi, X. Ouyang, L. Wang, W. Peng, J. Wen, and Y. Dai,A pilot metabolic

MA

[19]

NU

on LC-MS, Biomed Chromatogr. 30(2016) 1386-96.

profiling study in serum of patients with chronic kidney disease based on (1)

H. Yki-Jarvinen,Diagnosis of non-alcoholic fatty liver disease (NAFLD),

PT E

[20]

D

H-NMR-spectroscopy, Clin Transl Sci. 5(2012) 379-85.

Diabetologia. 59(2016) 1104-11. D.B. Liss, M.S. Paden, E.S. Schwarz, and M.E. Mullins,What is the clinical

CE

[21]

AC

significance of 5-oxoproline (pyroglutamic acid) in high anion gap metabolic acidosis following paracetamol (acetaminophen) exposure?, Clin Toxicol (Phila). 51(2013) 817-27. [22]

Y.M. Yu, C.M. Ryan, Z.W. Fei, X.M. Lu, L. Castillo, J.T. Schultz, et al.,Plasma L-5-oxoproline kinetics and whole blood glutathione synthesis rates in severely burned adult humans, Am J Physiol Endocrinol Metab. 282(2002) E247-58. 16

ACCEPTED MANUSCRIPT [23]

C.L. Foot, J.F. Fraser, and D.V. Mullany,Pyroglutamic acidosis in a renal transplant patient, Nephrol Dial Transplant. 20(2005) 2836-8.

[24]

I.A. Almaghlouth, J.Y. Mohamed, M. Al-Amoudi, L. Al-Ahaidib, A. Al-Odaib, and F.S. Alkuraya,5-Oxoprolinase deficiency: report of the first human

L. Zand, A. Muriithi, E. Nelsen, P.M. Franco, E.L. Greene, Q. Qian, et

RI

[25]

PT

OPLAH mutation, Clin Genet. 82(2012) 193-6.

SC

al.,Severe anion gap metabolic acidosis from acetaminophen use secondary to 5-oxoproline (pyroglutamic acid) accumulation, Am J Med Sci. 344(2012)

P. Tailor, T. Raman, C.L. Garganta, R. Njalsson, K. Carlsson, E. Ristoff, et

MA

[26]

NU

501-4.

al.,Recurrent high anion gap metabolic acidosis secondary to 5-oxoproline

E.J. Rolleman, E.J. Hoorn, P. Didden, and R. Zietse,Guilty as charged:

PT E

[27]

D

(pyroglutamic acid), Am J Kidney Dis. 46(2005) e4-10.

unmeasured urinary anions in a case of pyroglutamic acidosis, Neth J Med.

T. Soga, R. Baran, M. Suematsu, Y. Ueno, S. Ikeda, T. Sakurakawa, et

AC

[28]

CE

66(2008) 351-3.

al.,Differential metabolomics reveals ophthalmic acid as an oxidative stress biomarker indicating hepatic glutathione consumption, J Biol Chem. 281(2006) 16768-76. [29]

D.T.

Lawrence,

L.K.

Bechtel,

N.P.

Charlton,

and

C.P.

Holstege,5-oxoproline-induced anion gap metabolic acidosis after an acute acetaminophen overdose, J Am Osteopath Assoc. 110(2010) 545-51. 17

ACCEPTED MANUSCRIPT [30]

N.J. Waters, C.J. Waterfield, R.D. Farrant, E. Holmes, and J.K. Nicholson,Integrated

metabonomic

analysis

of

bromobenzene-induced

hepatotoxicity: novel induction of 5-oxoprolinosis, J Proteome Res. 5(2006) 1448-59. S. Geenen, C. Guallar-Hoyas, F. Michopoulos, J.G. Kenna, K.L. Kolaja, H.V.

PT

[31]

RI

Westerhoff, et al.,HPLC-MS/MS methods for the quantitative analysis of

SC

5-oxoproline (pyroglutamate) in rat plasma and hepatic cell line culture medium, J Pharm Biomed Anal. 56(2011) 655-63.

R. Kumar, S. Prakash, S. Chhabra, V. Singla, K. Madan, S.D. Gupta, et

NU

[32]

MA

al.,Association of pro-inflammatory cytokines, adipokines & oxidative stress with insulin resistance & non-alcoholic fatty liver disease, Indian J Med Res.

M. Otgonsuren, M.J. Estep, N. Hossain, E. Younossi, S. Frost, L. Henry, et

PT E

[33]

D

136(2012) 229-36.

al.,Single non-invasive model to diagnose non-alcoholic fatty liver disease

CE

(NAFLD) and non-alcoholic steatohepatitis (NASH), J Gastroenterol Hepatol.

[34]

AC

29(2014) 2006-13. G. Calori, G. Lattuada, F. Ragogna, M.P. Garancini, P. Crosignani, M. Villa, et al.,Fatty liver index and mortality: the Cremona study in the 15th year of follow-up, Hepatology. 54(2011) 145-52. [35]

M.A. Borman, F. Ladak, P. Crotty, A. Pollett, R. Kirsch, G. Pomier-Layrargues, et al.,The Fatty Liver Index has limited utility for the detection and quantification of hepatic steatosis in obese patients, Hepatol Int. 7(2013) 18

ACCEPTED MANUSCRIPT 592-9. [36]

G. Brooker, J. Jeffery, T. Nataraj, M. Sair, and R. Ayling,High anion gap metabolic acidosis secondary to pyroglutamic aciduria (5-oxoprolinuria): association with prescription drugs and malnutrition, Ann Clin Biochem.

W.B. Ashworth, N.A. Davies, and I.D. Bogle,A Computational Model of

RI

[37]

PT

44(2007) 406-9.

SC

Hepatic Energy Metabolism: Understanding Zonated Damage and Steatosis in NAFLD, PLoS Comput Biol. 12(2016) e1005105.

NU

M. Ono, N. Okamoto, and T. Saibara,The latest idea in NAFLD/NASH

CE

PT E

D

MA

pathogenesis, Clin J Gastroenterol. 3(2010) 263-70.

AC

[38]

19

ACCEPTED MANUSCRIPT Table 1. Clinical and biochemical characteristics of NAFLD and controls Controls (n=31) AGE (yrs)

(n=38)

52.8±14.6 25:7

31:9

(n=21) 56.3±12.9 24:8

Body mass index

47.3±6.8

47.5±8.1

Alanine aminotransferase (U/L)

35.6±7.2

65.7±20.8*

Aspartate aminotransferase (U/L)

21.3±6.7

38.7±12.8

50.7±23.5*

gamma -glutamyl transferase (U/L)

36.8±8.4

78.8±14.3

152±15.9*†

Triglyceride (mg/dL)

136±66.4

150±84.1

198±95.2*†

180±30.7

197±31.8

192±32.5

104±20.6

110±25.9

134±40.5*

D

Glucose (mg/dL)

41.4±23.5

SC

NU

MA

Cholesterol (mg/dL)

PT E

*p < 0.05, compared to controls. †

58.7±14.5

NASH

RI

Gender (M:F)

SS

PT

Parameter

p < 0.05, compared to SS.

AC

CE

Data refer to the number of patients or the mean ± SD.

20

47.5±8.4

ACCEPTED MANUSCRIPT Table 2. Fold changes of identified metabolites between the SS versus control groups or NASH versus SS groups based on the OPLS-DA. SS vs Control

NASH vs SS

Uracil

1.34

1.63

α-Linolenic Acid

1.33

1.52

Glutamate

1.56

1.80

L-Glutamine

1.56

Pyroglutamate

1.56

RI

PT

Metabolite

AC

CE

PT E

D

MA

NU

SC

1.81

21

2.26

ACCEPTED MANUSCRIPT Table 3. Groupwise comparisons of serum biomarkers for healthy controls and patients with SS and NASH Control (n = 25)

SS (n = 25)

NASH (n = 25)

(mmol/L)

1.48±1.21

3.18±1.53*

7.85±3.22†

Adiponectin (ug/mL)

13.25±8.36

11.69±6.25

7.08±4.40†

TNF-α (pg/mL)

3.89±3.21

4.26±2.96

IL-8 (pg/ml)

25.15±13.87

26.02±10.69



p < 0.05, compared to SS.

RI

SC

NU MA

*p < 0.05, compared to controls.

PT

Pyroglutamate

AC

CE

PT E

D

Data refer to the number of patients or the mean ± SD.

22

8.99±4.63† 17.73±8.70†

ACCEPTED MANUSCRIPT Table 4. The diagnostic performance of biomarkers for NASH

Biomarker

Sensitivity

Specificity

Positive

Negative

predictive value

predictive value

Accuracy

64%

56%

59%

60%

60%

TNF-α

72%

76%

75%

73%

74%

IL-8

65%

68%

86%

triglyceride

54%

75%

68%

GGT

57%

80%

pyroglutamate

72%

85%

pyroglutamate +adiponectin

82%

80%

pyroglutamate+TNF-α

83%

77%

pyroglutamate+IL-8

85%

82%

66%

62%

65%

74%

65%

69%

62%

90%

82%

58%

93%

80%

70%

86%

79%

91%

69%

82%

86%

81%

78%

88%

85%

86%

93%

89%

91%

87%

83%

84%

86%

pyroglutamate+ triglyceride

83%

86%

85%

83%

85%

pyroglutamate +GGT

85%

80%

82%

88%

85%

91%

76%

79%

90%

84%

+adiponectin+IL-8

90%

AC

pyroglutamate+

CE

+TNF-α

92%

PT E

+adiponectin+TNF-α pyroglutamate+IL-8

triglyceride+GGT

NU

D

pyroglutamate

MA

pyroglutamate

RI

40%

SC

PT

adiponectin

23

ACCEPTED MANUSCRIPT Figure legends: Figure 1. Mass chromatograms of serum obtained from control (A), simple steatosis (B) nonalcoholic steatohepatitis (C).

PT

Figure 2. (A) The PCA score plots of healthy controls (blue square) and SS patients

RI

(red triangle). (B) The PCA score plots of the SS (red triangle) and NASH

SC

(inverse yellow triangle) groups.

NU

Figure 3. (A) The OPLS-DA score plots of healthy controls (blue square) and SS

D

(inverse yellow triangle) groups.

MA

patients (red triangle). (B) The OPLS-DA score plots of SS (red triangle) and NASH

PT E

Figure 4. (A) The serum pyroglutamate levels in controls and patients with SS and NASH. The short transverse line represents the median of diff erent groups. (B) The

CE

ROC analysis of pyroglutamate in serum of NAFLD versus control with AUROC of

AC

0.880 and SS versus controls with AUROC of 0.846, respectively.

Figure 5. The ROC analysis of pyroglutamate, adiponectin, TNF-α and IL-8 in serum of NASH versus SS.

24

ACCEPTED MANUSCRIPT

AC

CE

PT E

D

MA

NU

SC

RI

PT

Fig. 1

25

ACCEPTED MANUSCRIPT

AC

CE

PT E

D

MA

NU

SC

RI

PT

Fig. 2

26

ACCEPTED MANUSCRIPT

AC

CE

PT E

D

MA

NU

SC

RI

PT

Fig. 3

27

ACCEPTED MANUSCRIPT

AC

CE

PT E

D

MA

NU

SC

RI

PT

Fig. 4

28

ACCEPTED MANUSCRIPT

AC

CE

PT E

D

MA

NU

SC

RI

PT

Fig. 5

29

ACCEPTED MANUSCRIPT Highlights 1. High performance liquid chromatography-mass spectrometry (HPLC-MS) was

used to analyse the metabolic profile of the serum samples from patients with nonalcoholic steatohepatitis (NASH) from simple steatosis (SS), and controls.

PT

2. Five metabolites, including uracil, α-linolenic acid, glutamate, L-glutamine and

RI

pyroglutamate, were identified as prominent factors in discriminating NASH from

SC

SS.

3. Pyroglutamate was found to be the most promising factor in distinguishing the

NU

NASH from SS groups with an optimal cut-off value of 4.82 mmol/L, the

MA

sensitivity and specificity of the diagnosis of NASH were 72% and 85%, respectively.

D

4. The potential diagnostic power of pyroglutamate was superior to other biomarkers

AC

CE

PT E

such as factor-α (TNF-α), adiponectin and interleukin-8 (IL-8).

30