Figures
Abstract
Early detection is the most effective way to improve the clinical outcome of biliary atresia (BA). Emerging metabolomics provides a powerful platform for discovering novel biomarkers and biochemical pathways to improve early diagnosis. The aim of this study is to find the potential biomarkers to distinguish BA from neonatal hepatitis syndrome (NHS) by using a metabolomics method. We comprehensively analyzed the serum metabolites in a total of 124 blood samples from patients with BA or neonatal hepatitis syndrome (NHS) and from normal individuals using advanced metabolomic approaches, and found that the levels of glutarylcarnitine (C5DC) significantly increased in the BA group while the levels of threonine (Thr) significantly rose in the NHS group comparing with the other groups. The levels of glutamic acid (Glu) in the BA group were significantly elevated compared to those in the NHS group, but still lower than the hyperbilirubinemia and normal controls. The levels of propionyl carnitine (C3), isovaleryl carnitine (C5) and glutamine (Gln) were reduced in the BA group compared to those in the NHS group, but still higher than the hyperbilirubinemia and normal controls. This study demonstrates the possibility of metabolomics as non-invasive biomarkers for the early detection of BA and also provides new insight into pathophysiologic mechanisms for BA.
Citation: Zhao D, Han L, He Z, Zhang J, Zhang Y (2014) Identification of the Plasma Metabolomics as Early Diagnostic Markers between Biliary Atresia and Neonatal Hepatitis Syndrome. PLoS ONE 9(1): e85694. https://doi.org/10.1371/journal.pone.0085694
Editor: Jean-Marc A. Lobaccaro, Clermont Université, France
Received: October 1, 2013; Accepted: December 6, 2013; Published: January 8, 2014
Copyright: © 2014 Zhao et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by the Shanghai Natural Science Foundation (No. 11ZR1423800), the Program of “New One Hundred Talented People” of the Shanghai Jiaotong University School of Medicine, Science Foundation for The Excellent Youth Scholars of Xinhua Hospital (to Y. Zhang). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Biliary atresia (BA) is one of the most serious hepatobiliary diseases during infantile period, accounting for about 75% of liver transplantations in children younger than 2 years [1]. The reported incidence of BA varies from 5/100,000 to 32/100,000 live births [2]. The etiology and pathogenesis of BA are largely unknown. Virus infection and autoimmune disorders are the two major theoretic mechanisms [3]–[6]. Since the Kasai procedure and liver transplantation have developed over the years, BA can be treated or cured. However, the better outcomes for infants depend on early diagnosis and timely Kasai portoenterostomy. It is generally known that children operated before day 60 achieve biliary drainage and stabilization in up to 80% of cases, while in patients operated only after day 90, the success rate decreases to 20% [7]. Therefore it is extremely important to differentiate BA from other causes of neonatal cholestatic jaundice promptly and accurately. However, besides BA, neonatal hepatitis syndrome (NHS) is another major cause of cholestatic jaundice in infants. The clinical manifestations and liver function test between the two diseases are too similar to distinguish. Up to date, there is still no accuracy investigation to diagnose BA. For example, the abdominal ultrasound of BA has a high specificity but low sensitivity, while the hepatobiliary scintigraphy (HBS) has a high sensitivity but low specificity [8], [9]. The “gold standard” of diagnosing BA is liver biopsy, which is invasive. Therefore, more sensitive and specific biomarkers for the early diagnosis of BA are needed.
Metabonomics is concerned with the measurement of global sets of low molecular weight metabolites, which can be regarded as important indicators of various states of disease. Through analyzing and verifying the specific biomarkers of a disease, metabolomics enables us to better understand pathological processes, and substance metabolic pathways. Compared with traditional diagnostic methods, even small changes of metabolites can help to detect early pathologic changes more sensitively [10], [11]. Its utility has been demonstrated by the identification of new biomarkers for Hepatocellular Carcinoma, Kidney Injury, and Alzheimer’s disease [11]–[13].
Today, the improved efficiency and accuracy of metabolomic biomarkers and early diagnostic technologies have been increasingly used in a clinic setting. By applying advanced analytical and statistical tools, the “metabolome” is mined for biomarkers that are associated with the state of BA. It may help to understand the mechanism of BA occurrence and progression on the metabolic level and provide information for the identification of early and differential marker metabolites for BA. Thus, the aim of this study was to evaluate the usefulness of the metabolic alterations as a rapid metabolomics screening technique for BA. As such, the identification of serum biomarkers may provide useful information for early diagnosis of BA, and effectively distinguish BA from other diseases such as NHS and hyperbilirubinemia.
Methods
Ethics and Clinical Populations
This study was a prospectively randomized controlled trial, whose design has been approved by the Ethical Committees of the Xinhua Hospital affiliated to Shanghai Jiaotong University School of Medicine. All parents or legal guardians signed an agreement after being disclosed to necessary information on this study.
From January 2008 to December 2011, a total of 61 cholestasis infants were included in the study and 32 hyperbilirubinemia infants and 31 normal infants as control. For all patients, serum alanine transaminase (ALT), aspartate transaminase (AST), total bilirubin (TB), direct bilirubin (DB), γ-glutamyl transpeptidase (γ-GT), total bile acid (TBA) and alkaline phosphatase (ALP) were measured. Four dried blood spots samples were obtained from infants by heel prick and spotted on filter paper [S&S 903# (CFH)] at the patients’ first presentation to hospital and then stored at −20°C.
Cholestasis infants were defined as conjugated bilirubin more than 20% of the total bilirubin when the total bilirubin is 85 µmol/L or more and 17 µmol/L or more when the total bilirubin is 85 µmol/L or less [14]. First, history taking and physical examination of all patients were conducted carefully. And then they underwent differential diagnosis and etiologic work-up of cholestasis, including a complete blood count, virus biomarkers, serum-concentration of α1-antitrypsin, blood tandem mass spectrometry and urine gas chromatography inspection as well as magnetic resonance cholangiopancreatography, abdominal ultrasound and HBS. In addition, specific biochemical tests were done as needed. For all, there were 45 cases suspicious for BA underwent intraoperative cholangiography and liver biopsy. And all of them had a 3 to 6 month clinical follow-up. Exclusion criteria: 1. Cholestasis infants with any family history such as progressive familial intrahepatic cholestasis, (PFIC), α1-antitrypsin deficiency; 2. Cholestasis infants with other obvious accompany manifestation like endocrine disorders, chromosomal disorders and Alagille’s syndrome. 3. Cholestasis infants with any amino, organic, and fatty acid disorders such as neonatal cholestasis induced by citrin deficiency (NICCD), tyrosinemia. Diagnostic criteria: 1.Hyperbilirubinemia group: Infants had a serum bilirubin up to 221 µmol/L or clinical jaundice appearing greater than 14 days of life with a conjugated bilirubin less than 20% of the total bilirubin. 2. Biliary atresia group: the confirmed diagnosis of BA was based on laparoscopy with surgical cholangiography and liver biopsy. The gallbladders or the remnants of gallbladders were found, but extrahepatic biliary duct could not be detected during the operation. And liver biopsy showed the typical histopathologic features, including bile duct proliferation, bile plugs, giant cell transformation, canalicular and cellular bile stasis and periportal fibrosis [15]. At last, 30 patients were identified as BA. 3. Neonatal hepatitis syndrome group: There were 15 infants who had undergone laparoscopy with surgical cholangiography and liver biopsy. During the operation all of them showed biliary tree and the giant cell transformation, inflammatory cell infiltration, and lobular disorganization on liver biopsy. They were finally diagnosed as NHS. Additionally, 16 patients recovered from jaundice and had the normalization of laboratory values during the clinical follow-up for about 3–6 months. So it was unnecessary to perform an intraoperative cholangiography. We confirmed these patients as idiopathic neonatal hepatitis. Finally, 30 cases of BA and 31 cases of NHS were enrolled in this study.
Samples Preparation
At the time of LC-MS/MS analysis, the frozen blood spots samples were thawed in normal atmospheric temperature. According to the previous report [16], three millimeter discs were punched and incubated in the methanolic internal standard solution (100 µL) for 20 min. The internal standards of amino acids and acylcarnitine were purchased from Cambridge Isotope Labs, USA. The concentrations of the deuterium-labeled internal standards per liter of methanol were as follows: [15N,2-13C1]-glycine ([N1]-[C1]-Gly), 12.5 µmol/L; [2H4]- alanine([D4]-Ala), 2.5 µmol/L; [2H8]- valine ([D8]-Val), 2.5 µmol/L; [2H3]-leucine ([D3]-Leu), 2.5 µmol/L; [2H3]-methionine ([D3]-Met), 2.5 µmol/L; [2H5]-phenylalanine([D5]-Phe),2.5 µmol/L; [13C6]-tyrosine([C6]-Tyr), 2.5 µmol/L; [2H3]-aspartic acid ([D3]-Asp), 2.5 µmol/L; [2H3]-glutamic acid ([D3]-Glu), 2.5 µmol/L; [2H2]-ornithine ([D2]-Orn), 2.5 µmol/L; [2H2]-citrulline ([D2]-Cit), 2.5 µmol/L; [2H4-13C1]-arginine ([D4]-[C1]- Arg), 2.5 µmol/L; [2H9]-carnitine([D9]- C0), 0.76 µmol/L; [2H3]-acetylcarnitine([D3]-C2), 0.19 µmol/L; [2H3]-propionylcarnitine ([D3]-C3), 0.04 µmol/L; [2H3]-butyrylcarnitine([D3]-C4), 0.04 µmol/L; [2H9]-isovalerylcarnitine([D9]-C5), 0.04 µmol/L; [2H3]-octanoylcarnitine ([D3]-C8), 0.19 µmol/L; [2H9]-myristoylcarnitine([D9]-C14), 0.04 µmol/L; [2H3]-palmitoylcarnitine ([D3]-C16), 0.08 µmol/L. The diluted sample was then manually transferred to a 96 well polypropylene microtiter plate (membrane pore size 0.2/0.45 µm, Millipore Corporation, USA) and dried by a hot air blower (Dri-Block DB-30, Techne, Britain) at 55°C. Each diluted sample was placed in 65°C forced air oven for 15 min by the addition of 60 µl of n-butanol hydrochloride (3 mol/L), covered with a thin Teflon sheet. After the plate was removed from the oven, the hot air blower removed the butanol hydrochloride. The butanol-derivatized samples were reconstituted with 100 µL of acetonitrile and water (80∶20 by volume), and each plate was covered with aluminum foil. The samples were then ready for MS/MS analysis.
LC-MS/MS Measurements
An API Mass spectrometer 2000 (Applied Biosystems, USA) coupled with a HPLC system of Agilent 1000 (Agilent Technologies, USA) was used to detect amino acids and acylcarnitine in the samples. Twenty microliters of prepared sample and solvent were introduced by autosampler into MS1 of the API Mass spectrometer 2000. Mobile phase used 80% acetonitrile: water and quaternary pump were used to set flow speed: 140 µl/min, 0.2 min→30 µl/min,1 min→300 µl/min, 0.2 min. The obtained spectra of amino acids were analyzed with neutral loss scan (neutral loss fragment butyl formate, m/z102) and multiple reaction monitoring. Neutral loss scan range of m/z 140 to m/z 280 and multiple reaction monitoring were test for Gly, Orn, Arg, Cit and its internal standard. The obtained spectrum of acylcarnitine in every sample was using precursor ion scanning mode. The precursor ion scan range from m/z 210 to m/z 502 and the product ion fragment was m/z 85. Hexanoylearnitine (C6) has no internal standards. As the mass to charge ratio of butyl-esterified C6 was 316.3, we detected it with D9-C5 quantitatively. The scan time of all the amino acids and acylcarnitine was set to 1.4 minutes and the each test sample would be detected with four minutes. The mass to charge ratio of all the amino acid and acylcarnitine was shown in Table 1 and the spectra were seen in Fig. 1.
Blood samples were prepared and mass spectrometric analysis of individual sample was performed as described in “Methods” in details. Panel A shows a mass spectrum was acquired in the neutral loss scan mode. The scan was range from m/z 140 to m/z 280. It can detect most of the amino acids. Panel B shows a mass spectrum was acquired in the multiple reaction monitor mode. It was mainly used to detect Gly, Orn, Arg, Cit and its internal standard. Panel C shows a mass spectrum was acquired in precursor ion scanning mode. It used to detect the acylcarnitine in the sample.
Data Processing
The data were analyzed by Chem View b5, (U.S. Applied Biosystems), which was in-house analytical software in API 2000 MS. Based on a variety of amino acid ester and acylcarnitine and its isotopic ion peak intensity of the internal standard, this software enables peak identification and quantification using an in-house metabolite library. The markers were expressed as “concentration” or as a “ratio” and were exported in EXCEL format. The overall resultant data included 17 amino acids (Alanine [Ala], Aspartic acid [Asp], Glutamic acid [Glu], Methionine [Met], Phenylalanine [Phe], Tyrosine [Tyr], Leucine [Leu], Tryptophane [Trp], Valine [Val], Argnine [Arg], Citrulline [Cit], Glycine [Gly], Ornithine [Orn], Glutamine [Gln], Histidine [His], Serine [Ser], Threonine [Thr]) and 30 acylcarnitine concentrations (Free carnitine [C0], Acetyl carnitine [C2], Propionyl carnitine [C3], Malonyl carnitine [C3DC], Butyrylcholinesterase carnitine [C4], 3-hydroxybutyrate acyl carnitine [C4-OH], Succinic acid carnitine [C4DC], Isovaleryl carnitine [C5], Senecioyloxy carnitine [C5:1], 3-hydroxy-isovaleryl-carnitine [C5-OH], Glutaryl carnitine [C5DC], Hexanoyl carnitine [C6], Hexene acyl carnitine [C6:1], Adipic acid carnitine [C6DC], Octanoyl carnitine [C8], Octene carnitine [C8:1], Sim two acyl carnitine [C8DC], Kwai acyl carnitine [C10], Aoi enoyl-carnitine [C10:1], Lauroylcarnitine [C12], Myrcene carnitine [C12:1], Myristoyl carnitine [C14], Nutmeg enoyl-carnitine [C14:1], 3-hydroxy-myristoyl-carnitine [C14-OH], Palmitoyl carnitine [C16], Palm enoyl-carnitine [C16:1], 3-hydroxy-palmitoyl carnitine [C16-OH], Octadecadienoic acid carnitine [C18], Carnitine octadecene [C18:1], 3-hydroxy-octadecadienoic acid carnitine [C18-OH]) as well as 7 amino acid ratios (Arg/Orn, Cit/Arg, Orn/Cit, Met/Phe, Leu/Phe, Phe/Tyr, Gly/Phe) and 8 acylcarnitine ratios (C3/C0, C3/C2, C4/C2, C5/C2, C5-OH/C2, C5DC/C8, C8/C2, C14:1/C8:1).
Statistical Analysis
Orthogonal partial least squares discriminant analysis (OPLS-DA) was applied to analyze LC-MS spectral data by SIMCA-P+ (v12.0 version, Umetrics, Umeå, Sweden). The number of components for OPLS-DA models was auto-fitted using a cross validation. OPLS-DA, a supervised pattern recognition (PR) method, was conducted to maximize the difference between the BA and other groups. The variable importance for projection (VIP) score for each component was derived from the OPLS-DA analysis and served as an indication of contribution to the OPLS model. The quality of the models was assessed by a 10-fold cross-validation method. The obtained parameters R2 stands for the total explained variation of the model and describes how well the derived model fit the data, and Q2 represents the predictability of the model and provides a measure of the model quality. According to a greater-than-one threshold of VIP in the OPLS-DA model, a number of variables are responsible for the differences in the metabolic profiling between BA, INH, hyperbilirubinemia and normal groups. Variables that were normally distributed and had homogeneous variance were then analyzed by one-way ANOVA with a LSD test from SPSS Statistics 12.0 (SPSS, Chicago, IL, USA). We examined if different biomarker candidates obtained from the OPLS-DA models were statistically significant among the four groups at the univariate level. A p-value <0.05 (confidence level 95%) was considered statistically significant. Besides, other distributed data were also analyzed using one-way ANOVA followed by a LSD test by SPSS.
Results
Demographic Characteristics
There were a total of 124 patients including 67 boys and 57 girls with a mean age of 39.8±26.9 days and a mean weight of 4.3±1.1 kg. These infants were divided into four groups as described above. Table 2 shows clinical data including the patients’ demographic information and serum liver function values. There was no apparent difference in the age of presentation, sex distribution and the liver function parameters including TB, DB, ALT, AST, TBA and ALP between NHS and BA groups. However, the γ-glutamyl transpeptidase (γ-GT) value was extremely higher in BA group than other groups.
Discovery and Identification of Metabolic Biomarkers
The LC-MS variables obtained from all subjects were used to construct an OPLS-DA model. In those four groups, we analyzed two sets of pair-wise data to identify the metabolic biomarkers to discriminate the BA and NHS groups. Fig. 2 show a clear separation between BA and normal group, NHS and normal group, BA and NHS group, hyperbilirubinemia and normal group. Table 3 indicates the performance parameters of the four models, suggesting that these models were reliable and highly predictive.
A. OPLS-DA score plots of BA and normal; B. OPLS-DA score plots of with NHS and normal; C. OPLS-DA score plots of BA and NHS; D. OPLS-DA score plots of hyperbilirubinemia and normal. The fig. showing that the two populations are well separated respectively.
The VIP score for each compared model was generated for each OPLS-DA analysis. A VIP score >1.0 implied that the metabolite contributed significantly towards the differentiation of the two compared groups, shown in Fig. 3.
The Variables were aligned in descending order of their VIP values. From model A and model B, we can obtain the identified metabolites of BA and NHS. From model C, we can verify these metabolites obtained above can distinguish BA from NHS. From model D, we can exclude the interference of metabolites in hyperbilirubinemia infants.
As shown in Fig. 3, from model A and model B, we obtained the specific metabolites of BA and NHS group, compared with the normal group respectively. From model C, we verified whether these metabolites obtained above could distinguish BA from NHS. And from model D, we excluded the interference of metabolites in hyperbilirubinemia infants. Finally, according to the criterion for VIP statistics (VIP>1), we obtained a total of nine metabolites for their most contribution in discriminating profiles between BA and NHS cohorts. Their VIPs were shown in Table 4. Then One-way analysis of variance (ANOVA) test by LSD was used to examine the statistical significance of these metabolites among the four groups. At last, six of which were detected with significant differences. (Table 5). The results suggest that the levels of C5DC significantly increased in the BA group, while in the NHS group, the levels of Thr significantly rose. The levels of Glu were significantly elevated in the BA group compared to those in the NHS group, but still lower than the hyperbilirubinemia and normal controls. The levels of C3, C5 and Gln in BA were significantly reduced compared to those in the NHS group, but still higher than those in the hyperbilirubinemia and normal controls.
Discussion
Since the Kasai hepatoportoenterostomy has been evolving over the years, more and more studies have demonstrated that the earlier age at Kasai was a predictor of improved native liver survival and overall survival of BA [17], [18]. But it is still difficult to differentiate BA from other causes of neonatal cholestatic jaundice especially the NHS. Because of limitations in conventional approaches, several noninvasive tests have been developed to facilitate such a diagnosis. With recent advances in metabonomics, serum metabolites could be used as a novel diagnostic indicator for disease [19].
Accruing evidence suggests that metabolic profiles are useful for identification of diagnostic biomarkers in BA. Zhou and Mushtaq have suggested that bile acids elevated in biliary atresia by comparing with other forms of cholestatic hepatobiliary diseases or neonatal jaundice group basing on tandem mass spectrometry [20], [21]. In the future, we also intend to further investigate the metabolomics of amino acids and acylcarnitine profiles of BA. In this study, we attempted to identify sensitive and specific metabolites to distinguish BA form NHS based on LC/MS, which does not require the time-consuming derivatization procedure and is suitable for analysis of the metabolites compounds [22]. Because the symptoms of BA initially are indistinguishable from neonatal hyperbilirubinemia and later other symptoms like clay colored stools, dark urine, and large hardened liver are indistinguishable from NHS. We set both hyperbilirubinemia and NHS as the study groups.
In the previous studies, there were different amino acid and acylcarnitine metabolite profiles corresponding to the different forms of liver disease [11], [23]–[27]. However there were only a few studies on BA using LC-MS-based metabonomic methods to investigate the acylcarnitine profiles. As generally known, carnitine plays an important role in fatty acid oxidation and energy production. It helps long-chain fatty acids to enter into the mitochondria for the purpose of β-oxidation and regulating the proportion and coenzyme A and acylcoenzyme A. The carnitine bound with acyl-CoA to form acylcarnitine, which help the elimination of organic acid. Animal experiments indicate that biliary export was a major route for acylcarnitine clearance and long-chain acyl-carnitine was hardly excreted through bile duct [28], [29]. So short-chain carnitine are the mainly excretion of carnitine in the bile duct. Selimoglu MA et al reported that plasma carnitine levels were significantly increased in children with biliary atresia [23]. In this study, we also found that there was a set of short-chain carnitine metabolites elevating in the BA patients. Among them, the C5DC concentrations in BA were higher than the other three groups, and C2, C5 concentrations in BA are statistically lower than NHS group but still higher than hyperbilirubinemia and normal controls. The elevated glutarylcarnitine (C5DC) suggested the deficient or nonfunctional of glutaryl-CoA dehydrogenase. And this was usually seen in the Glutaric aciduria type I disorder, acting as an important neurological dysfunction resembling an episode of encephalitis with movement disorders such as dystonia and/or dyskinesia [30]. However, in our study, the BA patients did not show such clinical onset, so the further mechanisms that lead to the pathophysiological changes are to be studied.
As the liver plays a central role in amino acid metabolism, it is of great importance to investigate the metabolite profiles in BA and NHS. From the liver function tests, we found that there was no significantly different elevation of serum ALT and AST in BA and NHS patients. However, there was a set of significantly different amino acids concentrations between BA and NHS.
Gln significantly increased while Glu significantly decreased in the serum of NHS patients compared with BA patients in our study. Coincidentally, these amino acids participate in the urea cycle and these changes suggest impairment in the urea cycle. Urea cycle describes the conversion reactions of ammonia into urea, which help the toxic ammonia removed from the body. A lot of reports on proteomic analysis of urea cycle revealed the regulation of urea cycle enzymes had dysfunction in the liver injury such as glutamine synthetase, glutamate dehydrogenase [31], [32]. As a result, the decreased levels of Glu and increased levels of Gln were detected in serum of NHS patients. Besides, higher levels of Thr were also found in the NHS infants, presumably because of the activity of threonine-degrading enzymes decreased. Above all, these reactions mainly occur in the liver especially in mitochondria of liver cell. So this change reflects the different degree of liver cell injury. According to the previous study, the patients with hepatitis have oxidative stress, which plays a crucial role in the induction and progression of liver injury [33]. Oxidative stress produced by inflammation in the NHS infants appears to be more active than those in the BA infants and results in a more serious liver damage. Subsequently there was a different amino acids profile between the two groups.
The aim of this study was to discover new noninvasive biomarkers for BA. We quantitatively analyzed the metabolites from patients. Our results indicate that there is a metabolic profile shift of amino acid and acylcarnitine in BA from NHS. The alterations of these metabolites suggest different disturbances in metabolism between BA and NHS. They can be potentially developed into a clinically useful diagnostic tool for contributing towards an improved understanding of a disease mechanism.
Author Contributions
Conceived and designed the experiments: YZ. Performed the experiments: DZ LH. Analyzed the data: DZ YZ ZH JZ. Contributed reagents/materials/analysis tools: DZ LH ZH JZ YZ. Wrote the paper: DZ YZ JZ.
References
- 1. Hartley JL, Davenport M, Kelly DA (2009) Biliary atresia. Lancet 374: 1704–1713.
- 2. Chardot C (2006) Biliary atresia. Orphanet J Rare Dis 1: 28.
- 3. Feldman AG, Mack CL (2012) Biliary atresia: cellular dynamics and immune dysregulation. Semin Pediatr Surg 21: 192–200.
- 4. Moore SW, Zabiegaj-Zwick C, Nel E (2012) Problems related to CMV infection and biliary atresia. S Afr Med J 102: 890–892.
- 5. Szavay PO, Leonhardt J, Czech-Schmidt G, Petersen C (2002) The role of reovirus type 3 infection in an established murine model for biliary atresia. Eur J Pediatr Surg 12: 248–250.
- 6. Mack CL, Feldman AG, Sokol RJ (2012) Clues to the etiology of bile duct injury in biliary atresia. Semin Liver Dis 32: 307–316.
- 7. Baumann U, Ure B (2012) Biliary atresia. Clin Res Hepatol Gastroenterol 36: 257–259.
- 8. Kianifar HR, Tehranian S, Shojaei P, Adinehpoor Z, Sadeghi R, et al. (2013) Accuracy of hepatobiliary scintigraphy for differentiation of neonatal hepatitis from biliary atresia: systematic review and meta-analysis of the literature. Pediatr Radiol (43) 905–919.
- 9. Yang JG, Ma DQ, Peng Y, Song L, Li CL (2009) Comparison of different diagnostic methods for differentiating biliary atresia from idiopathic neonatal hepatitis. Clin Imaging 33: 439–446.
- 10. Dutta M, Joshi M, Srivastava S, Lodh I, Chakravarty B, et al. (2012) A metabonomics approach as a means for identification of potential biomarkers for early diagnosis of endometriosis. Mol Biosyst 8: 3281–3287.
- 11. Wang X, Zhang A, Sun H (2013) Power of metabolomics in diagnosis and biomarker discovery of hepatocellular carcinoma. Hepatology 57: 2072–2077.
- 12. Han X, Rozen S, Boyle SH, Hellegers C, Cheng H, et al. (2011) Metabolomics in early Alzheimer's disease: identification of altered plasma sphingolipidome using shotgun lipidomics. PLoS One 6: e21643.
- 13. Noto A, Cibecchini F, Fanos V, Mussap M (2013) NGAL and metabolomics: the single biomarker to reveal the metabolome alterations in kidney injury. Biomed Res Int 2013: 612032.
- 14. Moyer V, Freese DK, Whitington PF, Olson AD, Brewer F, et al. (2004) Guideline for the evaluation of cholestatic jaundice in infants: recommendations of the North American Society for Pediatric Gastroenterology, Hepatology and Nutrition. J Pediatr Gastroenterol Nutr 39: 115–128.
- 15. Esmaili J, Izadyar S, Karegar I, Gholamrezanezhad A (2007) Biliary atresia in infants with prolonged cholestatic jaundice: diagnostic accuracy of hepatobiliary scintigraphy. Abdom Imaging 32: 243–247.
- 16. Zytkovicz TH, Fitzgerald EF, Marsden D, Larson CA, Shih VE, et al. (2001) Tandem mass spectrometric analysis for amino, organic, and fatty acid disorders in newborn dried blood spots: a two-year summary from the New England Newborn Screening Program. Clin Chem 47: 1945–1955.
- 17. Jimenez-Rivera C, Jolin-Dahel KS, Fortinsky KJ, Gozdyra P, Benchimol EI (2013) International incidence and outcomes of biliary atresia. J Pediatr Gastroenterol Nutr 56: 344–354.
- 18. Shneider BL, Brown MB, Haber B, Whitington PF, Schwarz K, et al. (2006) A multicenter study of the outcome of biliary atresia in the United States, 1997 to 2000. J Pediatr 148: 467–474.
- 19. Zhang A, Sun H, Wang X (2012) Serum metabolomics as a novel diagnostic approach for disease: a systematic review. Anal Bioanal Chem 404: 1239–1245.
- 20. Zhou K, Lin N, Xiao Y, Wang Y, Wen J, et al. (2012) Elevated bile acids in newborns with Biliary Atresia (BA). PLoS One 7: e49270.
- 21. Mushtaq I, Logan S, Morris M, Johnson AW, Wade AM, et al. (1999) Screening of newborn infants for cholestatic hepatobiliary disease with tandem mass spectrometry. BMJ 319: 471–477.
- 22. Xu F, Zou L, Liu Y, Zhang Z, Ong CN (2011) Enhancement of the capabilities of liquid chromatography-mass spectrometry with derivatization: general principles and applications. Mass Spectrom Rev 30: 1143–1172.
- 23. Selimoglu MA, Aydogdu S, Yagci RV, Huseyinov A (2001) Plasma and liver carnitine status of children with chronic liver disease and cirrhosis. Pediatr Int 43: 391–395.
- 24. Kalhan SC, Guo L, Edmison J, Dasarathy S, McCullough AJ, et al. (2011) Plasma metabolomic profile in nonalcoholic fatty liver disease. Metabolism 60: 404–413.
- 25. Zhang A, Sun H, Han Y, Yan G, Wang X (2013) Urinary metabolic biomarker and pathway study of hepatitis B virus infected patients based on UPLC-MS system. PLoS One 8: e64381.
- 26. Ressom HW, Xiao JF, Tuli L, Varghese RS, Zhou B, et al. (2012) Utilization of metabolomics to identify serum biomarkers for hepatocellular carcinoma in patients with liver cirrhosis. Anal Chim Acta 743: 90–100.
- 27. Byrd DJ, Wiltfang A, Rodeck B, Latta A, Burdelski M, et al. (1993) The plasma amino acid profile and its relationships to standard quantities of liver function in infants and children with extrahepatic biliary atresia and preterminal liver cirrhosis. Eur J Clin Chem Clin Biochem 31: 197–204.
- 28. Krahenbuhl S, Brass EP, Hoppel CL (2000) Decreased carnitine biosynthesis in rats with secondary biliary cirrhosis. Hepatology 31: 1217–1223.
- 29. Ramsay RR, Gandour RD, van der Leij FR (2001) Molecular enzymology of carnitine transfer and transport. Biochim Biophys Acta 1546: 21–43.
- 30. Couce ML, Lopez-Suarez O, Boveda MD, Castineiras DE, Cocho JA, et al. (2013) Glutaric aciduria type I: outcome of patients with early- versus late-diagnosis. Eur J Paediatr Neurol 17: 383–389.
- 31. Fleming KE, Wanless IR (2013) Glutamine synthetase expression in activated hepatocyte progenitor cells and loss of hepatocellular expression in congestion and cirrhosis. Liver Int 33: 525–534.
- 32. Antoine DJ, Dear JW, Starkey-Lewis P, Platt V, Coyle J, et al. (2013) Mechanistic biomarkers provide early and sensitive detection of acetaminophen-induced acute liver injury at first presentation to hospital. Hepatology. 58: 777–787.
- 33. Soga T, Sugimoto M, Honma M, Mori M, Igarashi K, et al. (2011) Serum metabolomics reveals gamma-glutamyl dipeptides as biomarkers for discrimination among different forms of liver disease. J Hepatol 55: 896–905.