The plasma profile of subjects with non-alcoholic fatty liver disease (NAFLD), steatosis and steatohepatitis (NASH), was examined using an untargeted global metabolomic analysis in order to identify specific disease-related pattern/s and to identify potential non-invasive biomarkers. Plasma samples were obtained after an overnight fast from histologically confirmed non-diabetic subjects with hepatic steatosis (N=11) or NASH (N=24), and compared with healthy, age and sexmatched controls (n=25). Subjects with NAFLD were obese, were insulin resistant and had higher plasma concentration of homocysteine and total cysteine and lower plasma concentrations of total glutathione. Metabolomic analysis showed markedly higher levels of glycocholate, taurocholate and glycochenodeoxycholate in subjects with NAFLD. Plasma concentrations of long chain fatty acids were lower and concentrations of free carnitine, butyrylcarnitine and methylbutyryl carnitine were higher in NASH. Several glutamyl dipeptides were higher, while cysteine-glutathione levels were lower in NASH and steatosis. Other changes included higher branched chain amino acids, phosphocholine, carbohydrates (glucose, mannose), lactate, pyruvate, and several unknown metabolites. Random forest analysis and recursive partitioning of the metabolomic data could separate healthy subjects from NAFLD with an error rate of ~8%, and NASH from healthy controls with an error rate of 4%. Hepatic steatosis and steatohepatitis could not be separated using the metabolomic profile.Conclusion-Plasma metabolomic analysis revealed marked changes in bile salts and in biochemicals related to glutathione in subjects with non-alcoholic fatty liver disease. Statistical analysis identified a panel of biomarkers that could effectively separate healthy controls from NAFLD and healthy controls from NASH. These biomarkers can potentially be used to follow response to therapeutic interventions.