Skip to main content

Proteome characteristics of liver tissue from patients with parenteral nutrition-associated liver disease



Parenteral nutrition (PN)-associated liver disease (PNALD) is a common and life-threatening complication in patients receiving PN. However, its definitive etiology is not yet clear. Therefore, performed proteomic analyses of human liver tissue to explore the same.


Liver tissue was derived and compared across selected patients with (n = 3) /without (n = 4) PNALD via isobaric Tag for Relative and Absolute Quantitation (iTRAQ)-based quantitative proteomics. Bioinformatics analysis was performed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases to explore the mechanisms of PNALD based on differentially expressed proteins (DEPs). The essential proteins that were differentially expressed between the two groups were explored and verified by western blotting.


A total of 112 proteins were found to be differentially expressed, of which 73 were downregulated, and 39 were upregulated in the PNALD group. Bioinformatics analysis showed DEPs to be associated with mitochondrial oxidative phosphorylation (mainly involved in mitochondrial respiratory chain complex I assembly), hepatic glycolipid metabolism (involved primarily in glycogen formation and gluconeogenesis), and oxidative stress (mainly involved in antioxidant change).


Overall, our results indicated that mitochondrial energy metabolism impairment, hepatic glycolipid metabolism disorder, and excessive oxidative stress injury might explain the comprehensive mechanism underlying PNALD. Moreover, we have provided multiple potential targets for further exploring the PNALD mechanism.


Total parenteral nutrition (TPN) is a vital therapeutic measure for patients with impaired gut function, including short bowel syndrome, severe inflammatory bowel disease, or chronic idiopathic intestinal pseudo-obstruction [1,2,3]. Parenteral nutrition-associated liver disease (PNALD) is one of the most common complications of long-term parenteral nutrition (PN), which severely impairs the physical health of patients, and is the primary factor limiting long-term PN therapy.

The incidence of PNALD in infants receiving long-term parenteral nutrition is between 25 and 60%, and that in adults is between 15 and 40% [4]. In adults, a history of PNALD is characterized by elevated liver enzymes and associated steatosis, with ensuing complications such as steatohepatitis, cholestatic hepatitis, as well as fibrosis and cirrhosis [5]. Most biochemical abnormalities of the liver can be reversed by weaning at an early stage of PNALD, whereas most of the advanced pathophysiological changes are irreversible, resulting in cirrhosis, decompensated liver disease, liver failure, and liver carcinoma. In these patients, PNALD progresses to end-stage liver disease, requiring combined intestinal and liver transplant [6,7,8].

Research on the mechanistic pathways and ameliorative modalities in PNALD is a significant focus in the field of gastroenterology and hepatology. Several theories have been proposed, for instance,the components of PN directly harm the liver or that the absence of enteral nutrition disrupts the enterohepatic axis leading to liver injur y[9]. However, the pathophysiology and etiology of PNALD remain unclear, and there have been very few studies exploring the overall protein expression in liver tissue of patients with PNALD.

Proteomics is a large-scale comprehensive study of proteins, including information on protein abundance and modification, along with their interacting networks [10]. Studies regarding PNALD have been limited to the transcriptomic level in mouse models [11], and, to the best of our knowledge, neither genomic nor proteomics studies on PNALD have been conducted till date. More importantly, no prior studies have explored the use of human specimens to understand this disease. Studies are, thus, required that use liver samples from PNALD patients for proteomics research, which would be more conducive to solving clinical problems. Based on this consideration, this study aimed to explore the possible mechanisms underlying the pathogenesis of PNALD, as well as identify potential therapeutic targets by performing hepatic proteomics in patients with or without PNALD in the present study.

Materials and methods


Patients with diagnostic percutaneous or intraoperative liver biopsy were recruited from the Department of General Surgery, Jinling Hospital, Medical School of Nanjing University. Liver tissue was obtained from seven patients, with (n = 3) and without (n = 4) PNALD. Disease information and clinical characteristics were collected from the electronic database and detailed in Table 1. Subjects with PNALD had liver enzymes elevated 1.5 times the upper limit of normal, in the absence of any other cause, such as viral hepatitis or drug-induced changes [12]. All subjects gave their informed consent for inclusion before they participated in the study. The study was approved by the Ethics Committee of the Jinling hospital, Medical School of Nanjing University (2017NZGKJ-071).

Table 1 Demographic and Clinical Characteristics of Participants

Sample preparation

The tissues were ground in liquid nitrogen. One milliliter of lysis buffer (7 M urea, 4% SDS, 1x Protease Inhibitor Cocktail (Roche Ltd. Basel, Switzerland)) was added to samples, followed by sonication on ice and centrifugation at 13000 rpm for 10 min at 4 °C. The supernatant was transferred to a fresh tube.

Protein digestion and iTRAQ labeling

Determine the protein concentration of the supernatant using the BCA protein assay, and then transfer 100 μg protein per condition into a new tube and adjust to a final volume of 100 μL with 100 mM TEAB (triethylammonium bicarbonate). Add 5 μL of the 200 mM DTT and incubate sample at 55 °C for 1 h, then add 5 μL of the 375 mM iodoacetamide to the sample and incubate for 30 min protected from light at room temperature. For each sample, proteins were precipitated with ice-cold acetone, and then were redissolved in 20 μL TEAB. Then proteins were tryptically digested with sequence-grade modified trypsin (Promega, Madison, WI), and the resultant peptide mixture was labeled using chemicals from the iTRAQ reagent kit. The labeled samples were combined, desalted using C18 SPE column (Sep-Pak C18, Waters, Milford, MA) and dried in vacuo.

High-pH reverse-phase separation

The peptide mixture was re-dissolved in buffer A (buffer A: 10 mM ammonium formate in water, pH 10.0, adjusted with ammonium hydroxide), and fractionated by high-pH separation using an Acquity UPLC system (Waters Corporation, Milford, MA) connected to a reverse-phase column (BEH C18 column, 2.1 mm × 150 mm, 1.7 μm, 300 Å, Waters Corporation, Milford, MA). The high-pH separation was performed using a linear gradient, starting from 0% B to 45% B in 45 min (B: 10 mM ammonium formate in 90% ACN, pH 10.0, adjusted with ammonium hydroxide). The column flow rate was maintained at 250 μL/min, and the column temperature was maintained at 45 °C. Twelve fractions were collected, and each was dried in a vacuum concentrator for the next step.

Low-pH nano-HPLC-MS/MS analysis

Fractions were re-suspended with 32 μl solvent C (C: water with 0.1% formic acid; D: ACN with 0.1% formic acid), separated by nanoLC, and analyzed by on-line electrospray tandem mass spectrometry. Experiments were performed on a nanoACQUITY UPLC system (Waters Corporation, Milford, MA) connected to a quadrupole-Orbitrap mass spectrometer (Q-Exactive) (Thermo Fisher Scientific, Bremen, Germany) equipped with an online nano-electrospray ion source. Eight-microliter peptide sample was loaded onto the trap column (Thermo Fisher Scientific Acclaim PepMap C18, 100 μm × 2 cm), with a flow rate of 10 μl/min for 3 min, and subsequently separated on the analytical column (Acclaim PepMap C18, 75 μm × 25 cm) with a linear gradient, from 5% D to 30% D in 105 min. The column was cleaned and re-equilibrated to initial conditions for 5 min. The column flow rate was maintained at 300 nL/min, and the column temperature was maintained at 45 °C. The electrospray voltage of 1.8 kV versus the inlet of the mass spectrometer was used.

The Q-Exactive mass spectrometer was operated in the data-dependent mode to switch automatically between MS and MS/MS acquisition. Survey full-scan MS spectra (m/z 350–1600) were acquired with a mass resolution of 70 K, followed by fifteen sequential high-energy collisional dissociations (HCD)-MS/MS scans with a resolution of 17.5 K. In all cases, one Microscan was recorded using a dynamic exclusion of 30 s. MS/MS-fixed first mass was set at 100.

Database searching and data analysis

Tandem mass spectra were extracted by Proteome Discoverer software (Thermo Fisher Scientific, version Charge state deconvolution and deisotoping were not performed. All MS/MS samples were analyzed using Mascot (Matrix Science, London, UK; version 2.3). Mascot was set up to search the Uniprot-SwissProt database (Taxonomy: Homo sapiens, 20,245 entries) assuming the digestion enzyme trypsin. Mascot was searched with a fragment ion mass tolerance of 0.050 Da and a parent ion tolerance of 10.0 PPM. Carbamidomethyl of cysteine and iTRAQ 8plex of lysine and the n-terminus were specified in Mascot as fixed modifications. Oxidation of methionine and iTRAQ 8plex of tyrosine were specified in Mascot as a variable modification. Use the percolator algorithm to control peptide level false discovery rates (FDR) lower than 1%. Only unique peptides were used for protein quantification, and the method of normalization on protein median was used to correct experimental bias, the minimum number of proteins that must be observed to allow was set to 1000. Bioinformatics analysis of the identified proteins was performed, and DEPs were defined in the iTRAQ experiment according to the following criteria: unique peptides ≥1, P-value < 0.05, fold change > 1.2 or < 0.8333 [13]. DEPs were entered into the DAVID (Database for Annotation, Visualization, and Integrated Discovery) database ( for functional classification and GO enrichment analysis [14], and to determine the significant pathways according to the KEGG pathway analysis ( Protein-protein interaction (PPI) networks were generated through the STRING database (v10, Pathway analysis was conducted using Ingenuity Pathway Analysis (IPA) ( Shanghai Biotree biotech Co., Ltd) on the differentially expressed proteins.

Verification of proteins by western blot analysis

Protein samples were separated by electrophoresis on 12% SDS-PAGE gel, transferred to polyvinylidene fluoride membranes, blocked with 5% non-fat milk for 1 h, and incubated with primary antibodies against CYP2B6, DDAH1 (Thermo Scientific, Rockford, IL), FABP5 (ProteinTech, Chicago, IL, USA), CAPG, and NDUFA1 (Affinity Bioscience, USA), overnight at 4 °C. Anti-rabbit horseradish peroxidase (HRP)-conjugated antibody was used as a secondary antibody, followed by ECL substrate (Thermo Scientific, Rockford, IL) incubation, and image visualization using Tanon 5200 imaging system (Tanon, China). Gray-scale analysis of the bands was performed using ImageJ software. Data are presented as ratios of the target protein to the internal control.

Statistical analyses - Proteins with p-value < 0.05 and Fold Change < 0.83 or > 1.2 were considered as differentially expressed proteins between the two groups. Student’s t-test was used for comparison of the difference between groups.


Comprehensive identification of proteome in human liver tissue

A total of 14,307 peptides and 3337 proteins were identified by iTRAQ analysis. Mass spectrometry results showed 112 differentially expressed proteins (DEPs), of which 73 were down regulated, and 39 were upregulated in the PNALD group (Table S1).

To identify the functional classification of DEPs, this study performed gene ontology analysis according to their molecular functions (MF), biological processes (BP), and cellular components (CC) with the assistance of DAVID Bioinformatics Resources. The top 15 annotations represented in each of the three GO categories are shown in Fig. 1. Majority of enriched categories included mitochondrial and mitochondria-related proteins (Fig. 2a, Table S2; total 27, 16 downregulated and 11 upregulated). Impaired oxidative phosphorylation was the predominant process in PNALD group. Mitochondria are known to be mainly responsible for the oxidative decomposition of dextrose, fat, and protein, for providing energy. Accordingly, the glycolipid metabolism (Fig. 2b, Table S3; total 32, 28 downregulated and 4 upregulated) and amino acid metabolism (Fig. 2c, Table S4; total 22, 15 downregulated and 7 upregulated) were significantly altered between the groups.

Fig. 1
figure 1

Significant Gene Ontology analysis of differentially expressed proteins. (BP, biological process; CC, Cellular component; MF, molecular function)

Fig. 2
figure 2

Heat map of differentially expressed proteins (DEPs) in PNALD and control patient. a: Mitochondrial and mitochondrial-related proteins; b: Glycolipid metabolisms associated proteins; c: Amino acid metabolisms associated proteins. Colored boxes represent upregulation (red) and downregulation (blue) in the PNALD. The color scale shown at the upper right indicates the fold changes in protein expression of all the samples. 1, 2, 3: PNALD group; A, B, C, D: Control group

Mitochondrial oxidative phosphorylation was impaired in PNALD patients

Mitochondrial oxidative phosphorylation system is the final biochemical pathway to produce ATP and the maintenance of cell function. Mitochondrial respiratory chain NADH dehydrogenase (complex I) is the most abundant enzyme in the electron transport chain [15], and is essential for oxidative phosphorylation in mitochondria [16]. Ten subunits of NADH dehydrogenase were found to be downregulated in the PNALD group, including NDUFB11, NDUFB7, NDUFV1, NDUFA7, NDUFV2, NDUFS8, NDUFA1, NDUFS2, NDUFA12, and NDUFS1 (Table S2). Consistent with our hypothesis that mitochondrial oxidative phosphorylation in the liver of patients with PNALD might be impaired, bioinformatics analysis indicated the DEPs to be enriched in mitochondria-associated biological processes (Fig. 3), including mitochondrial electron transport (NADH to ubiquinone), response to oxidative stress, oxidation-reduction process, and ATP synthesis-coupled electron transport (Table 2). And the IPA showed that the oxidative phosphorylation was significantly inhibited (Fig. 4). These all results indicated the mitochondrial oxidative phosphorylationwas impaired in PNALD patients.

Fig. 3
figure 3

Protein-protein interaction of differentially-expressed protein for group PNALD vs. control

Table. 2 Gene Ontology analysis of mitochondria associated differentially expressed proteins
Fig. 4
figure 4

Canonical Pathway diagram from Ingenuity Pathway Analysis: Oxidative phosphorylation was significantly inhibited. Colored boxes represent upregulation (red) and downregulation (green) in the PNALD

Hepatic glycolipid metabolism disorder in patients with PNALD

Liver plays a significant role in the control of glucose homeostasis by regulating various pathways in glucose metabolism, including glycogenesis, glycogenolysis, glycolysis, and gluconeogenesis. In this study, we observed a strong enrichment of DEPs linked to metabolic enzymes. For example, phosphorylase kinases (PHKA2, PHKB, and PHKG2), which stimulate glycogen degradation [17], were downregulated. Fructose 1,6-bisphosphate 1-phosphatases (FBP1 and FBP2), the rate-limiting enzyme in gluconeogenesis [18], were downregulated. Pyruvate carboxylase (PC), which catalyzes the ATP-dependent carboxylation of pyruvate to oxaloacetate, and involved in gluconeogenesis [19, 20], was also found to be downregulated. Consistent with the function of DEPs, bioinformatics analysis indicated that glycogen metabolic process (PHKA2, PHKB, PHKG2, and GNMT) and gluconeogenesis (GOT1, FBP1, FBP2, and PC) might be damaged in PNALD. Besides, we also found some glycolipid metabolism-associated DEPs, including nicotinamide N-methyltransferase (NNMT, upregulated) [21], D-3-phosphoglycerate dehydrogenase (PHGDH, downregulated) [22], and serine/threonine-protein phosphatase (CPPED1, downregulated) [23].

Besides its role in glucose metabolism, the liver plays a pivotal role in lipid metabolism and is the hub of fatty acid metabolism and lipid circulation [24]. Estrogen sulfotransferase (SULT1E1), involved in the process of adipogenesis [25, 26], was found to be downregulated and acetyl-CoA carboxylase 2 (ACACB), which inhibits fatty acid oxidation [27], was also downregulated. Further, lipid metabolism associated DEPs, including glycine N-methyltransferase (GNMT, downregulated) [28,29,30,31,32], fatty acid-binding protein (FABP5, upregulated) [33], cathepsin S (CTSS, upregulated) [34], aldose reductase (AKR1B1, upregulated) [35], and N(G), N(G)-dimethylarginine dimethylaminohydrolase 1 (DDAH1, downregulated) were also identified [36]. These results together indicated PNALD to be related to glycolipid metabolism disorder.

Oxidative stress caused by downregulation of the antioxidant factors may generate PNALD

Previous reports had shown oxidative stress injury to be one of the significant causes of PNALD. Production of ROS is increased when mitochondrial oxidative phosphorylation is impaired and under circumstances of antioxidant defense deficiency [37]. First, the impaired mitochondrial respiratory chain complex I assembly proteins damage the normal process of oxidative phosphorylation in patients with PNALD. Second, various antioxidant factors are downregulated in PNALD. For example, peroxiredoxin-6 (PRDX6, downregulated) is a peroxiredoxin that primarily functions as an antioxidant to scavenge peroxides in biological systems [38], regucalcin (RGN, downregulated) is an antioxidant [39], delta-aminolevulinic acid dehydratase (ALAD, downregulated) is an important antioxidant enzyme, whose inhibition may result in the accumulation of its substrate d-ALA, which in turn is associated with the overproduction of ROS [40], inhibition of PHGDH (downregulated) impairs the synthesis of heme, resulting in the impairment of oxidative phosphorylation and escape of electrons to molecular oxygen generating more ROS [41], inhibition of aspartate aminotransferase (GOT1, downregulated) inhibits the synthesis of the NADPH antioxidan t[42], and DDAH1 (downregulated) deficiency significantly enhances cellular oxidative stress [43]. , Third, we found histidine metabolism-associated proteins (AMDHD1, HNMT, FTCD, HAL, and UROC1) to be downregulated in PNALD. Previous studies had reported the histidine supplementation to inhibit oxidative stress and preserve mitochondrial membrane potential as well as dehydrogenase activity [44]. Moreover, IPA results showed that the “NRF2-mediated oxidative stress response”, “EIF2 signaling”, “PFKFB4 signaling pathway”, and “Glutathione mediated detoxification” were suppressed (Fig. S2). Previous studies had reported these pathways all participated in decreasing oxidative stres s[45,46,47,48]. Combined with the above results, we can hypothesize that oxidative stress caused by downregulation of the antioxidant factors may also participate in the development of PNALD.

Validation of differential expression proteins by western blotting

In order to verify the credibility of proteomics, we randomly selected several up- and downregulated DEPs for semi-quantitative verification by western blotting, including CYP2B6, DDAH1, NDUFA1, FABP5, and CAPG (Fig. 5). As shown, expression of CYP2B6, DDAH1, and NDUFA1 was significantly downregulated, and that of FABP5 and CAPG was upregulated in the PNALD group, compared to the control group. Western blotting results were in agreement with iTRAQ proteomics results.

Fig. 5
figure 5

Confirmation of differentially expressed proteins (DEPs) by western blotting; a, c: Immunoblotting analysis of DEPs (CYP2B6, DDAH1, NDUFA1, CAPG, and FABP5) in PNALD(H-PN) and control (H-EN). b, d: Gray-Level scores. The β-actin protein was used as a control. *, Compared with the control group, P < 0.05; ***, Compared with the control group, P < 0.001


Total parenteral nutrition is a life-saving therapy for patients with chronic gastrointestinal failure [49]. However, related metabolic complications, especially PNALD, arise with prolonged PN [50]. This study represents the first-ever investigation of liver tissue proteomic profiles of patients with PNALD to identify potential biomarkers and elucidate the molecular alterations that occur following PNALD. We successfully identified 112 DEPs between PNALD and control groups. The number and interaction of DEPs indicated the mitochondrial oxidative phosphorylation damage, glycolipid metabolism disorder, and oxidative stress injury as primary biological changes occurring in PNALD. In addition, we performed western blot analysis and found the results matching perfectly with the iTRAQ proteomics data.

In our study, the cellular localization of DEPs (predicted using GO analysis) revealed most of the identified putative proteins to be localized in mitochondria, which may explain the disorders in glucose and lipid metabolism. Mitochondria are best known for harboring pathways involved in ATP synthesis through the tricarboxylic acid cycle and oxidative phosphorylation [51]. More than 80% of ATP is produced by mitochondrial oxidative phosphorylation. An earlier study had reported that TPN deteriorates hepatic mitochondrial function [52]. Consistent with this, our current study indicated that multiple subunits or components of the respiratory chain complex I assembly proteins, which are essential for oxidative phosphorylation in mitochondria [16], are downregulated in PNALD. This may affect the biosynthesis of ATP, which is necessary for the maintenance of liver function. Mean ATP levels in PN mouse model were found to be reduced as compared to those in an EN mouse, although there was no statistically significant difference between these levels [52]. Therefore, we speculated that PN might cause PNALD, at least in part, by damaging the respiratory chain complex I assembly.

Insulin resistance is one of the metabolic phenotypes found in nonalcoholic fatty liver disease but there is limited clinical evidence for insulin resistance associated with PNALD. One report in piglets had demonstrated that chronic TPN induces significant insulin resistance [29]. In agreement with the previous study, our results showed upregulated NNMT, which correlates positively with insulin resistance [53,54,55], and downregulated GNMT, deficiency of which is known to impair glucose tolerance and insulin sensitivity. Furthermore, KEGG analysis indicated the DEPs to be associated with the insulin signaling pathway and PI3K-Akt signaling pathway (Table 3). These results together suggested that insulin sensitivity was impaired in patients with PNALD. Insulin resistance was also implicated by mitochondrial dysfunction or mitochondrial oxidative stress [56]. Mitochondrial oxidative phosphorylation was hampered, proteins with antioxidant function, including PRDX6, RGN, PHGDH, ALAD, GOT1, and DDAH1, were downregulated in patients with PNALD, and production of reactive oxygen species was increased, eventually leading to oxidative stress. Consistent with these results, several studies have demonstrated the impaired capacity of antioxidants during TPN [57, 58], and indicated oxidative stress play an important role in the development of PNALD [59, 60].

Table 3 KEGG pathways of differentially expressed proteins

Our results indicated that LGALS3 (Table S1), which contributes to inflammatory injury and fibrogenesis in cholestatic liver injury [61], was upregulated, and that DEPs were associated with NOD-like receptor signaling pathway. IPA results showed that the “LPS/IL-1 mediated inhibition of RXR function” pathway was activated (Fig. S2). Inflammation triggered by macrophages induced obesity related insulin resistance [62]. Consistent with these results, previous studies had indicated that LPS-mediated macrophages and IL-1β production or LPS-activated Kupffer cell through TLR4 might be early events in the pathogenesis of PNALD [63, 64].

A limitation of this study was an insufficient number of subjects, due to the extreme rarity of PNALD and infrequent source of liver tissue. Although previous studies had reported the abnormal glycolipid metabolism and oxidative stress in TPN, this is the first study that revealed damaged mitochondrial oxidative phosphorylation in patients with PNALD and provided a new target for exploring the mechanism and treatment of PNALD.

In conclusion, we identified several biological abnormalities that occur in patients with PNALD, such as in glycolipid metabolism, mitochondrial oxidative phosphorylation, and oxidative stress. We also detected the involvement of insulin and inflammatory signaling in PNALD pathogenesis. Our results provided new insights into the changes that occur in PNALD, and further identified candidate proteins as future PNALD biomarkers or therapeutic targets.

Availability of data and materials

All data generated or analyzed during this study are included in this published article (and its supplementary information files).



Parenteral nutrition


Parenteral nutrition-associated liver disease


Isobaric Tag for Relative and Absolute Quantitation


Gene Ontology


Kyoto Encyclopedia of Genes and Genomes


Differentially expressed proteins


Total parenteral nutrition


Molecular functions


Biological processes


Cellular components


Acetyl-CoA carboxylase 2


Aldose reductase


Delta-aminolevulinic acid dehydratase


Probable imidazolonepropionase


Macrophage-capping protein


Serine/threonine-protein phosphatase


Cathepsin S


Cytochrome P450 2B6




Dimethylaminohydrolase 1


Eukaryotic translation initiation factor 2 subunit 1


Fatty acid-binding protein


Fructose-1,6-bisphosphatase 1


Fructose-1,6-bisphosphatase isozyme 2




Glycine N-methyltransferase


Aspartate aminotransferase


Aspartate aminotransferase cytoplasmic


Histidine ammonia-lyase


Histamine N-methyltransferase




Relative and Absolute Quantitation


Ingenuity Pathway Analysis






Molecular function


Nuclear factor erythroid 2-related factor 2 9.


NADH dehydrogenase [ubiquinone] 1 alpha subcomplex subunit 1


NADH dehydrogenase [ubiquinone] 1 alpha subcomplex subunit 12


NADH dehydrogenase [ubiquinone] 1 alpha subcomplex subunit 7


NADH dehydrogenase [ubiquinone] 1 beta subcomplex subunit 11


NADH dehydrogenase [ubiquinone] 1 beta subcomplex subunit 7


NADH-ubiquinone oxidoreductase 75 kDa subunit mitochondrial


NADH dehydrogenase [ubiquinone] iron-sulfur protein 2


NADH dehydrogenase [ubiquinone] iron-sulfur protein 8


NADH dehydrogenase [ubiquinone] flavoprotein 1


NADH dehydrogenase [ubiquinone] flavoprotein 2 mitochondrial


Nicotinamide N-methyltransferase


6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase 4


Pyruvate carboxylase


D-3-phosphoglycerate dehydrogenase


Phosphorylase b kinase regulatory subunit alpha liver isoform


Phosphorylase b kinase regulatory subunit beta


Phosphorylase b kinase gamma catalytic chain liver/testis isoform






Reactive oxygen species


Estrogen sulfotransferase


Toll-like receptor 4


Total parenteral nutrition


Urocanate hydratase


  1. Stanghellini V. CRDG: natural history of intestinal failure induced by chronic idiopathic intestinal pseudo-obstruction. Transplant Proc. 2010;42:15–8.

    Article  CAS  PubMed  Google Scholar 

  2. Triantafillidis JK, Papalois AE. The role of total parenteral nutrition in inflammatory bowel disease: current aspects. Scand J Gastroenterol. 2014;49:3–14.

    Article  CAS  PubMed  Google Scholar 

  3. Ekema G, Milianti S, Boroni G. Total parenteral nutrition in patients with short bowel syndrome. Minerva Pediatr. 2009;61:283–91.

    CAS  PubMed  Google Scholar 

  4. Beath SV, Kelly DA. Total parenteral nutrition-induced cholestasis: prevention and management. Clinics in Liver Disease. 2016;20:159–76.

    Article  PubMed  Google Scholar 

  5. Orso G, Mandato C, Veropalumbo C, Cecchi N, Garzi A, Vajro P. Pediatric parenteral nutrition-associated liver disease and cholestasis: Novel advances in pathomechanisms-based prevention and treatment. Dig Liver Dis. 2016;48:215–22.

    Article  CAS  PubMed  Google Scholar 

  6. Sigalet D, Boctor D, Robertson M, Lam V, Brindle M, Sarkhosh K, Driedger L, Sajedi M. Improved outcomes in paediatric intestinal failure with aggressive prevention of liver disease. Eur J Pediatr Surg. 2009;19:348–53.

    Article  CAS  PubMed  Google Scholar 

  7. Koseesirikul P, Chotinaruemol S, Ukarapol N. Incidence and risk factors of parenteral nutrition-associated liver disease in newborn infants. Pediatr Int. 2012;54:434–6.

    Article  PubMed  Google Scholar 

  8. Anez-Bustillos L, Dao DT, Baker MA, Fell GL, Puder M, Gura KM. Intravenous Fat Emulsion Formulations for the Adult and Pediatric Patient: Understanding the Differences. Nutr Clin Pract. 2016;31:596–609.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Denton C, Price A, Friend J, Manithody C, Blomenkamp K, Westrich M, Kakarla V, Phillips W, Krebs J, Munoz SA, Osei H, Jain AK. Role of the GutLiver Axis in Driving Parenteral Nutrition-Associated Injury. Children (Basel, Switzerland). 2018;5:136.

    Google Scholar 

  10. Anderson NL, Anderson NG. Proteome and proteomics: new technologies, new concepts, and new words. Electrophoresis. 1998;19:1853–61.

    Article  CAS  PubMed  Google Scholar 

  11. Le Zhan IY, Kong B, Shen J, Gorczyca L, Memon N, Buckley BT, Guo GL. Dysregulation of bile acid homeostasis in parenteral nutrition mouse model. Am J Physiol Gastrointest Liver Physiol. 2016;310:G93–G102.

    Article  PubMed  Google Scholar 

  12. Beath SV, Kelly DA. Total parenteral nutrition–induced cholestasis. Clin Liver Dis. 2016;20:159–76.

    Article  PubMed  Google Scholar 

  13. UniProt: a hub for protein information. Nucleic Acids Res. 2015;43:D204–12.

  14. Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003;4:P3.

    Article  PubMed  Google Scholar 

  15. Chaban Y, Boekema EJ, Dudkina NV. Structures of mitochondrial oxidative phosphorylation supercomplexes and mechanisms for their stabilisation. Biochimica et Biophysica Acta (BBA) - Bioenergetics. 2014;1837:418–26.

    Article  CAS  Google Scholar 

  16. Formosa LE, Dibley MG, Stroud DA, Ryan MT. Building a complex complex: assembly of mitochondrial respiratory chain complex I. Semin Cell Dev Biol. 2018;76:154–62.

  17. Skamnaki VT, Owen DJ, Noble MEM, Lowe ED, Lowe G, Oikonomakos NG, Johnson LN. Catalytic mechanism of Phosphorylase kinase probed by mutational studies†,‡. Biochemistry. 1999;38:14718–30.

    Article  CAS  PubMed  Google Scholar 

  18. Kaur R, Dahiya L, Kumar M. Fructose-1,6-bisphosphatase inhibitors: a new valid approach for management of type 2 diabetes mellitus. Eur J Med Chem. 2017;141:473–505.

    Article  CAS  PubMed  Google Scholar 

  19. Jitrapakdee S, St Maurice M, Rayment I, Cleland WW, Wallace JC, Attwood PV. Structure, mechanism and regulation of pyruvate carboxylase. Biochem J. 2008;413:369–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Kumashiro N, Beddow SA, Vatner DF, et al. Targeting Pyruvate Carboxylase Reduces Gluconeogenesis and Adiposity and Improves Insulin Resistance. Diabetes. 2013;62(7):2183–94.

  21. Pissios P. Nicotinamide N -methyltransferase: more than a vitamin B3 clearance enzyme. Trend Endocrinol Metab. 2017;28:340–53.

    Article  CAS  Google Scholar 

  22. Okabe K, Usui I, Yaku K, Hirabayashi Y, Tobe K, Nakagawa T. Deletion of PHGDH in adipocytes improves glucose intolerance in diet-induced obese mice. Biochem Biophys Res Commun. 2018;504:309–14.

    Article  CAS  PubMed  Google Scholar 

  23. Vaittinen M, Kaminska D, Kakela P, Eskelinen M, Kolehmainen M, Pihlajamaki J, Uusitupa M, Pulkkinen L. Downregulation of CPPED1 expression improves glucose metabolism in vitro in adipocytes. Diabetes. 2013;62:3747–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Nguyen P, Leray V, Diez M, Serisier S, Bloc HJL, Siliart B, Dumon H. Liver lipid metabolism. J Anim Physiol Anim Nutr. 2008;92:272–83.

    Article  CAS  Google Scholar 

  25. Wada T, Ihunnah CA, Gao J, Chai X, Zeng S, Philips BJ, Rubin JP, Marra KG, Xie W. Estrogen Sulfotransferase inhibits adipocyte differentiation. Mol Endocrinol. 2011;25:1612–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Ihunnah CA, Wada T, Philips BJ, Ravuri SK, Gibbs RB, Kirisci L, Rubin JP, Marra KG, Xie W. Estrogen Sulfotransferase/SULT1E1 Promotes Human Adipogenesis. Mol Cell Biol. 2014;34:1682–94.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Murea M, Freedman BI, Parks JS, Antinozzi PA, Elbein SC, Ma L. Lipotoxicity in diabetic nephropathy: the potential role of fatty acid oxidation. Clin J Am Soc Nephrol. 2010;5:2373–9.

    Article  CAS  PubMed  Google Scholar 

  28. Chen YM, Shiu JY, Tzeng SJ, Shih LS, Chen YJ, Lui WY, Chen PH. Characterization of glycine-N-methyltransferase-gene expression in human hepatocellular carcinoma. Int J Cancer. 1998;75:787–93.

    Article  CAS  PubMed  Google Scholar 

  29. Stoll B, Horst DA, Cui L, Chang X, Ellis KJ, Hadsell DL, Suryawan A, Kurundkar A, Maheshwari A, Davis TA, Burrin DG. Chronic parenteral nutrition induces hepatic inflammation, Steatosis, and insulin resistance in neonatal pigs. J Nutr. 2010;140:2193–200.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Liao YJ, Chen TL, Lee TS, Wang HA, Wang CK, Liao LY, Liu RS, Huang SF, Chen YM. Glycine N-methyltransferase deficiency affects Niemann-pick type C2 protein stability and regulates hepatic cholesterol homeostasis. Mol Med. 2012;18:412–22.

    Article  CAS  PubMed  Google Scholar 

  31. Liu S, Li Y, Chen Y, Chiang E, Li AF, Lee Y, Tsai T, Hsiao M, Hwang S, Chen YA. Glycine N-methyltransferase−/− mice develop chronic hepatitis and glycogen storage disease in the liver. Hepatology. 2007;46:1413–25.

    Article  CAS  PubMed  Google Scholar 

  32. Martínez-Chantar ML, Vázquez-Chantada M, Ariz U, Martínez N, Varela M, Luka Z, Capdevila A, Rodríguez J, Aransay AM, Matthiesen R, Yang H, Calvisi DF, Esteller M, Fraga M, Lu SC, Wagner C, Mato JM. Loss of the glycine N-methyltransferase gene leads to steatosis and hepatocellular carcinoma in mice. Hepatology. 2008;47:1191–9.

    Article  PubMed  CAS  Google Scholar 

  33. Hotamisligil GS, Bernlohr DA. Metabolic functions of FABPs--mechanisms and therapeutic implications. Nat Rev Endocrinol. 2015;11:592–605.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Reißhauer A, Elmer N, Liebl ME. Multiple Statistical Methods for Assessing Differential Gene Expression in Microarray Data of Diabetic Model Rats to Predict the Molecular Mechanism of Atorvastatin on Anti-Atherogenesis. Exp Clin Endocrinol Diabetes. 2013;121(05):272–9.

  35. Brown KE, Broadhurst KA, Mathahs MM, Kladney RD, Fimmel CJ, Srivastava SK, Brunt EM. Immunodetection of aldose reductase in normal and diseased human liver. Histol Histopathol. 2005;20:429–36.

    CAS  PubMed  Google Scholar 

  36. Li T, Feng R, Zhao C, Wang Y, Wang J, Liu S, Cao J, Wang H, Wang T, Guo Y, Lu Z. Dimethylarginine Dimethylaminohydrolase 1 protects against high-fat diet-induced hepatic Steatosis and insulin resistance in mice. Antioxid Redox Signal. 2017;26:598–609.

    Article  CAS  PubMed  Google Scholar 

  37. Chow J, Rahman J, Achermann JC, Dattani MT, Rahman S. Mitochondrial disease and endocrine dysfunction. Nat Rev Endocrinol. 2017;13:92–104.

    Article  CAS  PubMed  Google Scholar 

  38. Patel P, Chatterjee S. Peroxiredoxin6 in endothelial signaling. Antioxidants. 2019;8:63.

    Article  CAS  PubMed Central  Google Scholar 

  39. Yamaguchi M, Murata T. Involvement of regucalcin in lipid metabolism and diabetes. Metabolism. 2013;62:1045–51.

    Article  CAS  PubMed  Google Scholar 

  40. Pereira B, Curi R, Kokubun E, Bechara EJ. 5-aminolevulinic acid-induced alterations of oxidative metabolism in sedentary and exercise-trained rats. J Appl Physiol(Bethesda, Md. 1985). 1992;72:226–30.

    Article  CAS  Google Scholar 

  41. Vandekeere S, Dubois C, Kalucka J, Sullivan MR, García-Caballero M, Goveia J, Chen R, Diehl FF, Bar-Lev L, Souffreau J, Pircher A, Kumar S, Vinckier S, Hirabayashi Y, Furuya S, Schoonjans L, Eelen G, Ghesquière B, Keshet E, Li X. Serine synthesis via PHGDH is essential for Heme production in endothelial cells. Cell Metab. 2018;28:573–87.

    Article  CAS  PubMed  Google Scholar 

  42. Hong C, Zheng J, Li X. Inhibition of GOT1 sensitizes colorectal cancer cells to 5-fluorouracil. Cancer Chemother Pharmacol. 2017;79:835–40.

    Article  CAS  PubMed  Google Scholar 

  43. Zhao C, Li T, Han B, Yue W, Shi L, Wang H, Guo Y, Lu Z. DDAH1 deficiency promotes intracellular oxidative stress and cell apoptosis via a miR-21-dependent pathway in mouse embryonic fibroblasts. Free Radic Biol Med. 2016;92:50–60.

    Article  CAS  PubMed  Google Scholar 

  44. Ommati MM, Jamshidzadeh A, Heidari R, Sun Z, Zamiri MJ, Khodaei F, Mousapour S, Ahmadi F, Javanmard N, Shirazi Yeganeh B. Carnosine and Histidine supplementation blunt Lead-induced reproductive toxicity through Antioxidative and mitochondria-dependent mechanisms. Biol Trace Elem Res. 2019;187:151–62.

    Article  CAS  PubMed  Google Scholar 

  45. Kubben N, Zhang W, Wang L, Voss TC, Yang J, Qu J, Liu G, Misteli T. Repression of the antioxidant NRF2 pathway in premature aging. Cell. 2016;165:1361–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Dang CV. Cancer cell metabolism: there is no ROS for the weary: figure 1. Cancer Discovery. 2012;2:304–7.

    Article  CAS  PubMed  Google Scholar 

  47. Chen J. Regulation of protein synthesis by the heme-regulated eIF2α kinase: relevance to anemias. Blood. 2006;109:2693–9.

    Article  CAS  Google Scholar 

  48. Ferrucci-Da Silva C, Zhan L, Shen J, Kong B, Campbell MJ, Memon N, Hegyi T, Lu L, Guo GL. Effects of total parenteral nutrition on drug metabolism gene expression in mice. Acta Pharm Sin B. 2020;10:153–8.

    Article  PubMed  Google Scholar 

  49. Dudrick SJ, Palesty JA. Historical highlights of the development of Total parenteral nutrition. Surg Clin North Am. 2011;91:693–717.

    Article  PubMed  Google Scholar 

  50. Mutanen A, Lohi J, Heikkilä P, Jalanko H, Pakarinen MP. Liver inflammation relates to decreased Canalicular bile transporter expression in pediatric onset intestinal failure. Ann Surg. 2017;268:1.

    Google Scholar 

  51. van der Bliek AM, Sedensky MM, Morgan PG. Cell biology of the mitochondrion. Genetics. 2017;207:843–71.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  52. Katayama T, Tanaka M, Tanaka K, Asonuma K, Uemoto S, Okamura R, Utsunomiya H, Fujita S, Ueda J, Tanaka A, Ozawa K. Alterations in hepatic mitochondrial function during Total parenteral nutrition in immature rats. J Parenter Enter Nutr. 1990;14:640–5.

    Article  CAS  Google Scholar 

  53. Riederer M, Erwa W, Zimmermann R, Frank S, Zechner R. Adipose tissue as a source of nicotinamide N-methyltransferase and homocysteine. Atherosclerosis. 2009;204:412–7.

    Article  CAS  PubMed  Google Scholar 

  54. Kannt A, Pfenninger A, Teichert L, Tönjes A, Dietrich A, Schön MR, Klöting N, Blüher M. Association of nicotinamide-N-methyltransferase mRNA expression in human adipose tissue and the plasma concentration of its product, 1-methylnicotinamide, with insulin resistance. Diabetologia. 2015;58:799–808.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Kraus D, Yang Q, Kong D, Banks AS, Zhang L, Rodgers JT, Pirinen E, Pulinilkunnil TC, Gong F, Wang Y, Cen Y, Sauve AA, Asara JM, Peroni OD, Monia BP, Bhanot S, Alhonen L, Puigserver P, Kahn BB. Nicotinamide N-methyltransferase knockdown protects against diet-induced obesity. Nature. 2014;508:258–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Fazakerley DJ, Minard AY, Krycer JR, Thomas KC, Stöckli J, Harney DJ, Burchfield JG, Maghzal GJ, Caldwell ST, Hartley RC, Stocker R, Murphy MP, James DE. Mitochondrial oxidative stress causes insulin resistance without disrupting oxidative phosphorylation. J Biol Chem. 2018;293:7315–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Matilla B, Ortíz J, González P, García-Díez F, Jorquera F, Culebras JM, González-Gallego J, Tuñón MJ. Effects of parenteral nutrition supplemented with glutamine or glutamine dipeptides on liver antioxidant and detoxication systems in rats. Nutrition. 2000;16:125–8.

    Article  CAS  PubMed  Google Scholar 

  58. Sokol RJ, Taylor SF, Devereaux MW, Khandwala R, Sondheimer NJ, Shikes RH, Mierau G. TWFB: hepatic oxidant injury and glutathione depletion during total parenteral nutrition in weanling rats. Am J Phys. 1996;270:G691–700.

    CAS  Google Scholar 

  59. Staun M, Pironi L, Bozzetti F, Baxter J, Forbes A, Joly F, Jeppesen P, Moreno J, Hébuterne X, Pertkiewicz M, Mühlebach S, Shenkin A, Van Gossum A. ESPEN guidelines on parenteral nutrition: home parenteral nutrition (HPN) in adult patients. Clin Nutr. 2009;28:467–79.

    Article  PubMed  Google Scholar 

  60. Cai W, Wu J, Hong L, Xu Y, Tang Q, Shi C. Oxidative injury and hepatocyte apoptosis in total parenteral nutrition–associated liver dysfunction. J Pediatr Surg. 2006;41:1663–8.

    Article  PubMed  Google Scholar 

  61. Tian J, Yang G, Chen H, Hsu DK, Tomilov A, Olson KA, Dehnad A, Fish SR, Cortopassi G, Zhao B, Liu F, Gershwin ME, Török NJ, Jiang JX. Galectin-3 regulates inflammasome activation in cholestatic liver injury. FASEB J. 2016;30:4202–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Lauterbach MAR, Wunderlich FT. Macrophage function in obesity-induced inflammation and insulin resistance. Pflugers Arch - Eur J Physiol. 2017;469:385–96.

    Article  CAS  Google Scholar 

  63. El Kasmi KC, Vue PM, Anderson AL, Devereaux MW, Ghosh S, Balasubramaniyan N, Fillon SA, Dahrenmoeller C, Allawzi A, Woods C, McKenna S, Wright CJ, Johnson L, D Alessandro A, Reisz JA, Nozik-Grayck E, Suchy FJ, Sokol RJ. Macrophage-derived IL-1β/NF-κB signaling mediates parenteral nutrition-associated cholestasis. Nat Commun. 2018;9:1393.

  64. El Kasmi KC, Anderson AL, Devereaux MW, Fillon SA, Harris JK, Lovell MA, Finegold MJ, Sokol RJ. Toll-like receptor 4–dependent Kupffer cell activation and liver injury in a novel mouse model of parenteral nutrition and intestinal injury. Hepatology. 2012;55:1518–28.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references


We express our gratitude to all of the participants who consented to participate in this study.


This work gained support from the National Natural Science Foundation of China (General Program, 81770531), the Excellent Youth Foundation Project of Jiangsu Province (BK20170009), the Key Medical Talents Foundation of Jiangsu Province (ZDRCA2016091), the National Natural Science Foundation of China (81700518, 81900524), the Natural Science Foundation of Jiangsu Province (BK20170622), Natural Science Foundation of Shandong Province (ZR2019BH010) and the Chinese National Science Foundation (31371373).

Author information

Authors and Affiliations



W.X.Y. and T.F. designed the study. T.F., S.H.F., Z.L., and G.X.J. collected the samples. M.G. performed the experiments. M.G. and S.H.F. analyzed the data and wrote the manuscript. W.X.Y. and X.B. reviewed/edited the manuscript. W.X.Y. and X.B. supervised the research. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Bin Xue or Xinying Wang.

Ethics declarations

Ethics approval and consent to participate

All subjects gave their informed consent for inclusion before they participated in the study. The study was approved by the Ethics Committee of the Jinling hospital, Medical School of Nanjing University (2017NZGKJ-071).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1. Supplementary information:

Supplementary tables and other bioinformation analysis (Annotation Enrichment Analysis, Different Proteins Analysis, Hierarchical Clustering Analysis, KEGG Analysis and, Network Analysis and Ingenuity Pathway Analysis) were listed in the file.

Rights and permissions

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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Maitiabola, G., Tian, F., Sun, H. et al. Proteome characteristics of liver tissue from patients with parenteral nutrition-associated liver disease. Nutr Metab (Lond) 17, 43 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: