Skip to main content

Changes in the intestinal expression of drug metabolism-related genes in a piglet model of parenteral nutrition



Parenteral nutrition (PN) may serve as a nutritional supportive therapy accompanied by oral medication, but the effect of PN on intestinal expression of drug metabolism-related genes remains unknown.


Twelve Bama piglets receiving PN for 14 days were used as in vivo model. Changes in intestinal drug metabolism-related genes were examined by proteomic analysis. Serum levels of fibroblast growth factor 19 (FGF19) were determined by ELISA, and the effect of FGF19 on the expression of drug metabolism-related genes was examined using murine ileum organoids.


A total of 1063 differentially expressed proteins were identified in PN group. Of note, two drug transporters (Abcb1 and Abcc2) were significantly decreased in PN group, along with two glutathione-related drug-metabolizing enzymes, glutathione peroxidase (Gpx2) and glutathione S-transferase (Gsta1). Serum FGF19 levels were dramatically reduced in PN group. Treatment with recombinant FGF19 in vitro dose-dependently up-regulated the expression of Abcb1, Abcc2, Gpx2 and Gsta1 in organoids.


Our data indicated that intestinal drug metabolism-related genes were significantly dysregulated under PN, and some of the changed genes were attributed to gut-derived FGF19.


Parenteral nutrition (PN) is a vital therapy for patients with impaired gut function who cannot tolerate enteral nutrition (EN), including severe inflammatory bowel diseases, short bowel syndrome and chronic intestinal pseudo-obstruction [1,2,3]. As a nutritional supportive therapy that may be accompanied by oral medication, we assumed that micro-environment of drug metabolism in those patients receiving PN may be different from that in those patients receiving EN. For instance, a number of studies using both PN animal models and human samples reported that genes involved in drug metabolism changed significantly in liver, suggesting that the process of drug metabolism might be disturbed by liver injury [4, 5]. Given that gut also plays a pivotal role in drug metabolism, especially the absorption of oral drugs, and that intestinal homeostasis is dramatically disrupted during PN, we assumed that intestinal expression of drug metabolism-related genes might be dysregulated, leading to impaired metabolic pathways in gut beyond liver.

Proteomic profiling is a large-scale comprehensive investigation of proteins, including information on protein abundance and modification, and their interacting networks. To the best of our knowledge, neither transcriptomic nor proteomic investigation in gut tissues have been conducted under PN background till date. Therefore, it is necessary to clarify the effect of PN on the proteome characteristics of small intestine, which may affect the pharmacokinetic profile of oral drugs. ATP-binding cassette (ABC) transporters have been reported to play a major role in protecting against dietary toxicants as well as limiting the oral bioavailability of drugs of therapeutic use [6], most notably ATP binding cassette subfamily B member 1 (Abcb1, also known as multidrug resistance protein 1) and ATP binding cassette subfamily C member 2 (Abcc2, also known as multidrug resistance-associated protein 2). Additionally, glutathione-related drug-metabolizing enzymes have been reported to play a pivotal role in biotransformation to protect against dietary oxidants and endogenously produced reactive oxygen species (ROS) in the gastrointestinal tract [7], including glutathione-s-transferases (GSTs) and glutathione peroxidases (GPXs). On the other hand, previous studies have demonstrated that gut-derived peptide hormone fibroblast growth factor 19 (FGF19), a postprandial hormone induced by Farnesoid X receptor (FXR) upon activation by bile acids, plays a fundamental role in the modulation of various metabolic processes [8]; and more importantly, reduced levels of circulating FGF19 have been implicated in the severe complications of PN use, notably PN-associated liver disease (PNALD) [9]. We hypothesized that intestinal drug metabolism-related genes might be significantly dysregulated under PN, and some of the changed genes might be attributed to gut-derived FGF19. This present study provided new insights into changes in intestinal drug metabolism-related genes under PN, and preliminary evidence showing the association of FGF19 with drug metabolism.

Material and methods


Cell culture reagents for ileum organoids were purchased from STEMCELL (#06005, USA). Recombinant human FGF19 protein was purchased from Abcam (#ab200246, USA). GSH/GSSG kit was purchased from Beyotime (#S0053, China). All other chemicals were purchased from Sigma (China).

PN piglet model

Twelve 2-day-old Guangxi Bama minipigs (0.7–1 kg) were housed in a heat- and light-controlled room with 12 h light/dark cycle, randomized to EN and PN group (n = 6/group). Surgery procedure and PN formula was identical to our previous work [10]. Briefly, BTPU-040 catheters (Instech, USA) were used for cannulation via external jugular vein. For PN group, the piglets were parenterally fed with PN solution containing 240 ml/kg/day of protein (13 g), lipid (5 g) and carbohydrate (25 g); Detailed formula shown in Table 1. PN piglets received 50% of the total requirement of the PN solution during the first 24 h. EN pigs were fed with milk powder containing soy protein and glucose, providing digestion energy of 4,050 kcal/kg and 48 g/kg in 240 mL of warm water. All animals were sacrificed at 14 days after surgery. This study was approved by Xin Hua Hospital Animal Use Committee (XHEC-F-2020–008).

Table 1 Components of parenteral nutrition solution

Sample collection

Ileal tissues were isolated about 10 cm proximal to the ileocecal valve. The segments were divided into several parts (1–2 cm each) for intestinal morphology assessment, proteome analysis, RNA isolation and glutathione measurement. All the samples were stored at − 80 °C until further use.

Intestinal morphology assessment

Formalin-fixed and paraffin-embedded tissue Sects. (5 μm) were stained using standard immunohistochemistry procedures. Villus height and crypt depth were measured in 20+ well oriented, full-length crypt/villus units per specimen and averaged. Data were analyzed with LAS AF LITE image processing software (Leica, Germany).

Proteome analysis

Ileum protein extracts were generated by homogenizing 20 mg tissue in 200 μL lysis buffer (4% SDS, 100 mM DTT, 150 mM Tris–HCl pH 8.0), followed by tandem mass tag (TMT)-labelling and proteome analysis using LC–MS/MS. TMT-labelled peptides were loaded onto a the C18-reversed phase column (12 cm long, 75 μm ID, 3 μm) in buffer A (2% acetonitrile and 0.1% Formic acid) and separated with a linear gradient of buffer B (90% acetonitrile and 0.1% Formic acid) at a flow rate of 300 nL/min over 60 min as follows: 0–2 min, linear gradient from 2 to 5% buffer B; 2–42 min, linear gradient from 5 to 20% buffer B; 42–50 min, linear gradient from 20 to 35% buffer B; 50–52 min, linear gradient from 35 to 90% buffer B; 52–60 min, buffer B maintained at 90%. MS data was acquired using a data-dependent top 15 method dynamically choosing the most abundant precursor ions from the survey scan (300–1800 m/z) for high-energy collisional dissociations (HCD) fragmentation.Dynamic exclusion duration was 30 s. Survey scans were acquired at a resolution of 70,000 at m/z 200 and resolution for HCD spectra was set to 17,500 at m/z 200. Normalized collision energy was 30 eV. Normalized collision energy was 30.

Database searches and bioinformatics analysis

The resulting LC–MS/MS raw files were imported into MaxQuant software (version for data interpretation and protein identification against the database Uniprot_Hordeum-vulgare_201747-20180125 (downloaded on 25/01/2018, and including 201747 protein sequences). Differentially expressed proteins (DEPs) were defined in the TMT experiment according to the following criteria: unique peptides ≥ 1, p-value < 0.05, fold change > 1.2 or < 0.83 [4]. Bioinformatics analysis were then performed, including Gene ontology (GO)-term classification (, and Protein–protein interaction (PPI) networks (v10,

Quantitative real-time PCR (qPCR)

Total RNA was extracted using Trizol (life technology, USA), and cDNA was synthesized with RNA isolation plus kit (Takara, Japan) according to the manufacturer’s protocol. Quantitative real-time PCR was performed using PIKO96 (Thermo, Germany). Primers used in this study are listed in Table 2.

Table 2 Primers for quantitative RT-PCR analysis

Glutathione measurement

Oxidized glutathione (GSSG) and reduced glutathione (GSH) were measured using a glutathione kit (Beyotime, China). Briefly, 10 mg ileum tissue obtained from 10 cm proximal to the ileocecal valve was homogenated in 0.5 ml cold GSH assay buffer, followed by centrifugation at 100,000 g for 10 min. GSH and GSSG levels were determined according to manufacturer’s instruction, and were normalized to tissue weight.

FGF19 levels

Blood was taken via portal vein, and serum FGF19 levels in portal vein were determined using an ELISA kit according to manufacturer’s protocol (#ELP-FGF19, Raybiotech, USA).

Culture and treatment of ileum organoids

Isolation and culture of primary organoids were conducted following the method described previously [11]. Briefly, primary crypts were released from the ileum mucosa of a C57BL wild-type mouse (8w, male). Isolated crypts were then mixed with Matrigel (Corning Inc, Corning, NY) and cultured in murine intestinal organoid growth medium (#06005, STEMCELL). For FGF19 treatment, organoids were treated with recombinant human FGF19 (20 and 50 ng/ml) at 1-day after seeding for 6 days. Morphological changes were examined under phasecontrast microscopy.

Transport assay in organoids

For Abcb1 function assay, Rhodamine 123 (Rho 123) assay was performed as described previously [12]. Briefly, organoids were incubated with growth medium containing Rho 123 (10 μM) in the dark at 37 °C for 1 h. Afterwards, these organoids were washed five times with PBS and then incubated in culture medium to allow Rho 123 efflux in the dark at 37 °C for 1 h, with or without Verapamil (100 μM, a specific inhibitor for Abcb1). A 100 μl aliquot of sample was taken for the determination of Rho 123 content using a fluorescence spectrophotometer (SpectraMax, CA, USA) under ex/em = 485 nm/520 nm. Alternatively, organoids were also imaged using Leica DMI6000B fluorescence microscopy with LAS AF LITE image processing software (Leica, Germany).

Statistical analysis

Data are expressed as means ± SD, statistically analyzed and plotted using GraphPad Prism 8.0 (GraphPad Software Inc., San Diego, CA, USA). Variables were analyzed either by Student’s t-test for two groups or ANOVA analysis for multiple groups, with statistical significance when p < 0.05.


Proteome characteristics of ileum tissue from piglet model of PN

In ileum tissues, a total of 40,682 peptides and 6534 proteins were identified by TMT analysis, and 1063 proteins were significantly changed in PN group (Fig. 1A), including 540 down-regulated proteins and 523 up-regulated proteins (Fig. 1B). In GO-term analysis, a number of DEPs involved in drug binding (MF #7/10) were identified by annotation clustering (Fig. 1C). PPI analysis revealed that interaction between DEPs could be roughly clustered to three fields (Fig. 1D), including RNA RECOGNITION & IMMUNE RESPONSE (red), DRUG BINDING & METABOLISM (blue) and DNA BINDING & CHROMATIN MODULATION (green).

Fig. 1
figure 1

Proteome characteristics of ileum tissue from piglet model of PN. A Volcano plot depicting the protein changes induced by PN treatment. Plots are up-regulated (red) or down-regulated (green) proteins with fold change > 1.2 or < 0.83 (p < 0.05). B Heat map of differentially expressed proteins (DEPs). C Gene Ontology analysis of DEPs. BP biological process; CC Cellular component; MF molecular function D Protein–protein interaction of DEPs. Three key interactions were identified, including RNA recognition & immune response (red), drug binding & metabolism (blue) and DNA binding & chromatin modulation (n = 3/group)

Effect of PN on ABC transporters and glutathione-related genes

For ABC transporters, Abcb1 (0.77-fold vs EN) and Abcc2 (0.62-fold vs EN) significantly decreased in PN group. For glutathione peroxidases, Gpx2 significantly decreased in PN group (0.67-fold vs EN), while Gpx8 increased (1.3-fold vs EN). For glutathione transferases, Gsta1 (0.53-fold vs EN) and Gstm3 (0.81-fold vs EN) significantly decreased in PN group, while Gstt1 (1.26-fold vs EN) increased (Fig. 2A). Consistently, qPCR data indicated that mRNA changes basically mirrored protein changes except Gpx8 and Gstm3 (Fig. 2B).

Fig. 2
figure 2

Changes in the intestinal expression of ABC transporters and Glutathione-related genes in vivo. A Relative protein levels determined by LC–MS. Ileum tissues were isolated from 14d-PN piglets. Data are presented as mean ± SD. **p < 0.01, ***p < 0.001 (n = 3/group). B Relative mRNA levels determined by qPCR using ∆∆Ct. Data are presented as mean ± SD. **p < 0.01, ***p < 0.001 (n = 6/group)

Effect of PN on glutathione levels in ileum

In addition to glutathione-related genes, glutathione levels also changed significantly in PN group. As shown, GSH levels were 1.42 ± 0.20 μmol/g tissue and 0.72 ± 0.16 μmol/g tissue in EN and PN group (p < 0.05), respectively (Fig. 3A). However, GSSG levels were slightly reduced without statistical differences (Fig. 3B), and consequently there were no significant differences in GSH/GSSG ratio between groups (Fig. 3C).

Fig. 3
figure 3

Total levels of GSH, GSSG and GSH/GSSG ratio in ileum tissue. Ileum tissues were isolated from 14d-PN piglets. Data are presented as mean ± SD. *p < 0.05 (n = 6/group)

Changes in morphological features of ileum and FGF19 levels

After 14d treatment, PN group exhibited significant villus atrophy in ileum tissues (Fig. 4A), characterized by decreased villus height and crypt depth comparing to EN group. Villus height in EN and PN group was 286 ± 65 μm and 180 ± 48 μm (p < 0.05), respectively. Crypt depth in EN and PN group was 143 ± 16 μm and 73 ± 24 μm (p < 0.05), respectively (Fig. 4B). Additionally, serum FGF19 levels in EN and PN group was 11.0 ± 2.1 ng/ml and 5.6 ± 1.5 ng/ml (p < 0.05), respectively (Fig. 4C).

Fig. 4
figure 4

Effect of PN on ileum epithelium and serum FGF19 levels. A Representative images showing villus atrophy in the ileum epithelium. Scale bar = 100 μm. B Quantification of villus height and crypt depth in (A). Data are presented as mean ± SD. *p < 0.05 (n = 6/group). C Serum levels of FGF19 in portal vein. Data are presented as mean ± SD. *p < 0.05 (n = 6/group)

Effect of FGF19 on the expression of drug metabolism-related genes

In order to address whether those changes in drug metabolism-related genes were attributed to reduced levels of FGF19, we then introduced murine ileum organoids as in vitro model. Firstly, we confirmed that FGF19 on the dosage of 50 ng/ml induced no significant impacts on organoid expansion and differentiation. As shown, no significant changes were observed in morphological features (average organoid number/field and average crypt number/organoid) (Fig. 5A, B). Of note, most of those genes changed in vivo were also dramatically regulated by FGF19 in vitro. For ABC transporters, both Abcb1 (Abcb1a and Abcb1b) and Abcc2 were up-regulated by FGF19 treatment in a dose-dependent manner. For glutathione-related genes, Gpx2 and Gsta1 were markedly up-regulated by FGF19 treatment, while Gpx8, Gstm3 and Gstt1 were unchanged. (Fig. 5C). Taken together with above-mentioned in vivo data, these findings suggest changes in Abcb1, Abcc2, Gpx2 and Gsta1 under PN were probably attributed to reduced levels of FGF19.

Fig. 5
figure 5

Effect of FGF19 on ABC transporters Glutathione-related genes in ileum organoids. A Representative images of organoids treated with FGF19 (50 ng/ml). B Expansion and differentiation of organoids described in (A). Data are presented as mean ± SD C Relative mRNA levels of ABC transporters and Glutathione-related genes. Data are presented as mean ± SD from three independent experiments. *p < 0.05, **p < 0.01

Effect of FGF19 on the efflux function of Abcb1

As shown, fluorescence signal of Rho 123 was enriched within the lumen domain of organoids, which was significantly enhanced in FGF19 group, and was suppressed in Verapamil group (Fig. 6A). Quantitatively, average signal intensity in FGF19 group was increased by about 30% compared to UT (untreated) group (p < 0.05, Fig. 6B). Consistently, released Rho 123 levels in medium were approximately 2.4 nM/60 min and 3.2 nM/60 min in UT and FGF19 group, respectively (p < 0.05, Fig. 6C). Collectively, these results suggested that FGF19 treatment can facilitate the efflux function of Abcb1 in vitro.

Fig. 6
figure 6

Effect of FGF19 on Abcb1 efflux function. A Efflux transport of Rho 123 through Abcb1. Organoids were treated with FGF19 (50 ng/ml) for 6 days. Verapamil (100 uM) served as positive control. Scale bar = 200 μm B Fluorescence intensity in organoids described in (A). At least 20 organoids were analyzed in each group. Data are presented as mean ± SD. *p < 0.05, **p < 0.01. C Rho 123 in medium. Released Rho123 in culture medium was determined as described in Methods section. Data are presented as mean ± SD. **p < 0.01. Three independent experiments were performed that presented similar results


According to the United States Healthcare Cost and Utilization Project, it is estimated that over 375,000 patients rely on PN each year in the United States, a number that has almost doubled in the past 20 years [13]. Given that PN mostly serves as a supportive role in nutrient feeding for those patients who are simultaneously receiving oral medication, whether administration of PN can affect drug metabolism raised our great interest. In this study, we introduced piglet model instead of rodent models, because it has more anatomical, physiological, immunological and metabolic similarities with human than any other rodent model [14]. For instance, the pig is more appropriate for modeling liver function and metabolism, as it has hepatic features similar to that of humans and, unlike the rat lack of a gallbladder [15]. The piglet has also been proven to be a useful representation of the human neonate when studying lipid nutrition, including the effect of long-chain n-3 polyunsaturated fatty acids on protein metabolism in the neonate during growth [16]. Overall, the prominence of the piglet in studies concerning human neonatal nutrition reflects the thought that the pig is most similar to humans compared to other animals and therefore the preferred model.

By using ileal tissues, this study represented the first-ever investigation of proteomic profiles in gut to elucidate the whole picture of DEPs by PN treatment. Notably, GO-term analysis revealed a number of DEPs involved in drug binding, namely TOP2A, DHFR, PDE4D, ALB, ATP1B1, DCK (Fig. 1C). Some of these genes are implicated in drug resistance that may significantly affect clinic outcomes. For instance, DHFR (dihydrofolate reductase, 1.45-fold increase in PN) that catalyzes tetrahydrofolate regeneration by reduction of dihydrofolate, is a key factor involved in the resistance of methotrexate [17]; TOP2A (DNA topoisomerase 2a, 1.66-fold increase in PN) that controls and alters the topologic states of DNA during transcription, has been implicated in the resistance of etoposide [18]. Additionally, PPI analysis suggested strong connections among drug binding, RNA recognition and DNA modulation (Fig. 1D), indicating that the delivery of nucleic acid drugs (e.g. anti-viral agents or anti-tumor agents) might be influenced by PN. For instance, more attention should be paid to the use of antibiotics for the management of intestinal infection for pediatric patients receiving PN.

For ABC transporters, we found both Abcb1 and Abcc2 were significantly decreased after PN administration (Fig. 2). The protein product of Abcb1 (P-glycoprotein) is an ATP-dependent drug efflux pump for xenobiotic compounds with broad substrate specificity. Interestingly, anti-viral agents, most of which are nucleic acid drugs that might be influenced by PN as mentioned above, frequently behave as Abcb1 substrates. On the other hand, the protein product of Abcc2 (Mrp2) serves as an efflux pump for some anionic drug conjugates including glucuronides, sulfates and glutathiones, contributing to elimination of a wide range of drugs, including, antihypertensives (e.g., olmesartan and temocaprilate), antineoplastics (e.g., methotrexate, cisplatin, doxorubicin and irinotecan), antiretrovirals (e.g., adefovir, lopinavir and saquinavir), and antibiotics (e.g., ampicillin and azithromycin) [19]. Therefore, we speculate that if the expression of these above-mentioned efflux transporters is decreased in the patients receiving PN, there would be an increase in bioavailability of certain drugs, which would increase drug concentrations in the body leading to potential risks of drug-drug interaction.

On the other hand, excessive endogenous or exogenous toxins may persist in the gut lumen due to gastrointestinal hypomotility under PN, contributing to the development of gut pathologies such as inflammation. Therefore, biotransformation and drug-metabolizing enzymes in gut can be as important as those in liver to metabolic detoxication. In this study, we found that a number of glutathione-related genes like Gpx2 and Gsta1 were decreased in PN group (Fig. 2), as well as the absolute content of GSH (Fig. 3). Gsta1 belongs a member of a family of enzymes that add GSH to target electrophilic compounds, including carcinogens, therapeutic drugs, environmental toxins, and products of oxidative stress [20]; Gpx2 is a member of glutathione peroxidase family which catalyzes the reduction of organic hydro peroxides and hydrogen peroxide by GSH, and thereby protect cells against oxidative damage [21]. GSH is a ubiquitous tripeptide (-glu-cys-gly). Reduced GSH and its oxidized form, GSSG, is the major thiol redox system in cells, of crucial importance for cellular function. In gut, the maintenance of intestinal mucosal GSH status is crucial for detoxication reactions, as sufficient GSH levels can enhance intestinal peroxide metabolism and thus drive peroxide uptake from gut lumen and reduce peroxide output into lymph [22]. Collectively, we speculate that reduction in GSH and GSH-related enzymes by PN may not only impair elimination of cellular ROS, but also facilitate the absorption and distribution of oxidant compounds, toxins or lipophilic drug derivatives via lymph circulation. Nevertheless, despite the absolute content of GSH was obviously reduced, GSH/GSSG ratio was unchanged, suggesting that oxidative injury in gut was less severe than liver, where GSH/GSSG ratio was significantly decreased as described elsewhere [5].

One of the key findings in this study is the role of FGF19 in the dysregulation of these drug metabolism-related genes, a well-studied postprandial hormone involved in the regulation of bile acid homeostasis and energy metabolism. First of all, FGF19 signals can activate FGFR4/β-klotho receptor complex to repress hepatic cholesterol 7α-hydroxylase (CYP7A1), thereby limiting the synthesis of bile acids [23]. In addition, FGF19 can prevent lipid accumulation by inducing fatty acid oxidation and suppressing acetyl coenzyme A carboxylase 2 (ACC2), a key gene involved in fatty acid synthesis [24]. Furthermore, FGF19 can regulate glucose metabolism through stimulation of glycogen synthesis and inhibition of gluconeogenesis via inactivation of cAMP regulatory element-binding protein (CREB) and proliferator-activated receptor g coactivator-1α (PGC-1α) [25]. Therefore, the beneficial impact of FGF19 on glucose, lipid and bile acid homeostasis raise the possibility to pursue FGF19 as a therapeutic target for treating various metabolic disorders including diabetes, non-alcoholic fatty liver disease and cholestatic liver diseases [26, 27]. Given that PN is a condition commonly featured by long-term deprivation of oral feeding along with disrupted bile acid homeostasis, we hypothesized that changes in those above-mentioned genes may be, at least partially attributed to FGF19. For in vivo experiments, we examined serum FGF19 levels in portal vein and morphological changes in ileum where FGF19 is predominantly produced. As expected, PN induced significant villus atrophy leading to reduced production of FGF19 (Fig. 4). For in vitro experiments, we examined the effect of FGF19 on the expression of drug metabolism-related genes by using ileum organoids. One aspect needs comment here. The dosage of FGF19 here was refered to Kim et al.’s work in which the effect of 50 ng/ml of recombinant human FGF19 was clarified on mouse ileum organoids [28]. As FGF19 signaling is also involved in cell survival and proliferation in addition to a wide range of metabolic processes including cholesterol, lipid and glucose, it’s important to discriminate whether these changes in drug metabolism-related genes were resulted from cell proliferation or differentiation. In this study, we found a number of drug-related genes were significantly regulated by FGF19 at 50 ng/ml, including Abcb1a/b, Abcc2, Gpx2 and Gsta1. Given no morphological changes observed (Fig. 5), we assumed that these changes in the expression of drug metabolism-related genes were attributed to the specific effect of FGF19 on gene regulation, other than non-specific effect on cell survival or proliferation. Finally, we confirmed that the efflux function of Abcb1 was substantially facilitated by FGF19 (Fig. 6), indicating that the regulation of Abcb1 by FGF19 was evident not only on gene levels, but also on protein levels.

Finally, two limitations of this study need comment, 1) the small number of replicates in the proteomic analysis, and 2) no data from other intestinal segments except ileum. Although previous studies have reported impaired intestinal homeostasis under PN, this present study is the first one into abnormal expression pattern of intestinal drug metabolism-related genes. In conclusion, this study indicated that PN can induce significant changes in the expression of drug metabolism-related genes in ileum, which may affect the pharmacokinetc process of certain drugs in clinical practice. Some of the changed genes may be attributed to gut-derived FGF19, including glutathione-related genes (Gsta1 and Gpx2) and ABC transporters (Abcb1 and Abcc2).

Availability of data and materials

All data generated or analyzed during this study are included in this published article.



Parenteral nutrition


Enteral nutrition


Differentially expressed proteins


Tandem mass tag


High-energy collisional dissociations


Quantitative real-time PCR


Gene ontology


Protein–protein interaction


Oxidized glutathione




Fibroblast growth factor 19

ABC transporters:

ATP-binding cassette transporters


Farnesoid x receptor


Reactive oxygen species


ATP binding cassette subfamily B member 1


ATP binding cassette subfamily C member 2




Glutathione peroxidase


Rhodamine 123


Dihydrofolate reductase


DNA topoisomerase 2a


Cholesterol 7α-hydroxylase


Acetyl coenzyme A carboxylase 2


CAMP regulatory element-binding protein


Proliferator-activated receptor g coactivator-1α


  1. Ko D, Yang HB, Youn J, Kim HY. Clinical outcomes of pediatric chronic intestinal pseudo-obstruction. J Clin Med. 2021;10:2376.

    Article  CAS  Google Scholar 

  2. Forbes A, Escher J, Hebuterne X, Klek S, Krznaric Z, Schneider S, Shamir R, Stardelova K, Wierdsma N, Wiskin AE, Bischoff SC. ESPEN guideline: clinical nutrition in inflammatory bowel disease. Clin Nutr. 2017;36:321–47.

    Article  Google Scholar 

  3. Baxter KJ, Srinivasan JK, Ziegler TR, Dhere T, Ricketts RR, Durham MM. Pediatric short bowel syndrome and subsequent development of inflammatory bowel disease: an illustrative case and literature review. Pediatr Surg Int. 2017;33:731–6.

    Article  Google Scholar 

  4. Maitiabola G, Tian F, Sun H, Zhang L, Gao X, Xue B, Wang X. Proteome characteristics of liver tissue from patients with parenteral nutrition-associated liver disease. Nutr Metab (Lond). 2020;17:43.

    Article  CAS  Google Scholar 

  5. Ferrucci-Da SC, 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  Google Scholar 

  6. Zhang J, Wang H, Fan Y, Yu Z, You G. Regulation of organic anion transporters: role in physiology, pathophysiology, and drug elimination. Pharmacol Ther. 2021;217:107647.

    Article  CAS  Google Scholar 

  7. Walsh ME, Shi Y, Van Remmen H. The effects of dietary restriction on oxidative stress in rodents. Free Radical Biol Med. 2014;66:88–99.

    Article  CAS  Google Scholar 

  8. Madnawat H, Welu AL, Gilbert EJ, Taylor DB, Jain S, Manithody C, Blomenkamp K, Jain AK. Mechanisms of parenteral nutrition-associated liver and gut injury. Nutr Clin Pract. 2020;35:63–71.

    Article  Google Scholar 

  9. van Erpecum KJ, Schaap FG. Intestinal failure to produce FGF19: A culprit in intestinal failure-associated liver disease? J Hepatol. 2015;62:1231–3.

    Article  Google Scholar 

  10. Liu Y, Xiao Y, Chen S, Tian X, Wang W, Wang Y, Cai W. The farnesoid X receptor agonist tropifexor prevents liver damage in parenteral nutrition-fed neonatal piglets. J Pediatr Gastroenterol Nutr. 2021;73:e11–9.

    Article  CAS  Google Scholar 

  11. Liu Y, Gong Z, Zhou J, Yan J, Cai W. Lin 28A/Occludin axis: An aberrantly activated pathway in intestinal epithelial cells leading to impaired barrier function under total parenteral nutrition. FASEB J. 2021;35:e21189.

    CAS  PubMed  Google Scholar 

  12. Zhao J, Zeng Z, Sun J, Zhang Y, Li D, Zhang X, Liu M, Wang X. A novel model of P-glycoprotein inhibitor screening using human small intestinal organoids. Basic Clin Pharmacol Toxicol. 2017;120:250–5.

    Article  CAS  Google Scholar 

  13. Pimiento JM, Evans DC, Tyler R, Barrocas A, Hernandez B, Araujo-Torres K, Guenter P. Value of nutrition support therapy in patients with gastrointestinal malignancies: a narrative review and health economic analysis of impact on clinical outcomes in the United States. J Gastrointest Oncol. 2021;12:864–73.

    Article  Google Scholar 

  14. Puiman P, Stoll B. Animal models to study neonatal nutrition in humans. Curr Opin Clin Nutr Metab Care. 2008;11:601–6.

    Article  CAS  Google Scholar 

  15. Shulman RJ. The piglet can be used to study the effects of parenteral and enteral nutrition on body composition. J Nutr. 1993;123:395–8.

    Article  CAS  Google Scholar 

  16. Odle J. New insights into the utilization of medium-chain triglycerides by the neonate: observations from a piglet model. J Nutr. 1997;127:1061–7.

    Article  CAS  Google Scholar 

  17. Bertino JR. Cancer research: from folate antagonism to molecular targets. Best Pract Res Clin Haematol. 2009;22:577–82.

    Article  Google Scholar 

  18. Wong N, Yeo W, Wong WL, Wong NL, Chan KY, Mo FK, Koh J, Chan SL, Chan AT, Lai PB, Ching AK, Tong JH, Ng HK, Johnson PJ, To KF. TOP2A overexpression in hepatocellular carcinoma correlates with early age onset, shorter patients survival and chemoresistance. Int J Cancer. 2009;124:644–52.

    Article  CAS  Google Scholar 

  19. Chen Z, Shi T, Zhang L, Zhu P, Deng M, Huang C, Hu T, Jiang L, Li J. Mammalian drug efflux transporters of the ATP binding cassette (ABC) family in multidrug resistance: a review of the past decade. Cancer Lett. 2016;370:153–64.

    Article  CAS  Google Scholar 

  20. Bogaards JJ, Venekamp JC, van Bladeren PJ. Stereoselective conjugation of prostaglandin A2 and prostaglandin J2 with glutathione, catalyzed by the human glutathione S-transferases A1–1, A2–2, M1a–1a, and P1–1. Chem Res Toxicol. 1997;10:310–7.

    Article  CAS  Google Scholar 

  21. Brigelius-Flohe R, Maiorino M. Glutathione peroxidases. Biochim Biophys Acta. 2013;1830:3289–303.

    Article  CAS  Google Scholar 

  22. Aw TY. Determinants of intestinal detoxication of lipid hydroperoxides. Free Radical Res. 1998;28:637–46.

    Article  CAS  Google Scholar 

  23. Holt JA, Luo G, Billin AN, Bisi J, McNeill YY, Kozarsky KF, Donahee M, Wang DY, Mansfield TA, Kliewer SA, Goodwin B, Jones SA. Definition of a novel growth factor-dependent signal cascade for the suppression of bile acid biosynthesis. Genes Dev. 2003;17:1581–91.

    Article  CAS  Google Scholar 

  24. Fu L, John LM, Adams SH, Yu XX, Tomlinson E, Renz M, Williams PM, Soriano R, Corpuz R, Moffat B, Vandlen R, Simmons L, Foster J, Stephan JP, Tsai SP, Stewart TA. Fibroblast growth factor 19 increases metabolic rate and reverses dietary and leptin-deficient diabetes. Endocrinology. 2004;145:2594–603.

    Article  CAS  Google Scholar 

  25. Potthoff MJ, Boney-Montoya J, Choi M, He T, Sunny NE, Satapati S, Suino-Powell K, Xu HE, Gerard RD, Finck BN, Burgess SC, Mangelsdorf DJ, Kliewer SA. FGF15/19 regulates hepatic glucose metabolism by inhibiting the CREB-PGC-1alpha pathway. Cell Metab. 2011;13:729–38.

    Article  CAS  Google Scholar 

  26. Jahn D, Rau M, Hermanns HM, Geier A. Mechanisms of enterohepatic fibroblast growth factor 15/19 signaling in health and disease. Cytokine Growth Factor Rev. 2015;26:625–35.

    Article  CAS  Google Scholar 

  27. Mazzetti M, Marconi G, Mancinelli M, Benedetti A, Marzioni M, Maroni L. The management of cholestatic liver diseases: current therapies and emerging new possibilities. J Clin Med. 2021;10:1763.

    Article  CAS  Google Scholar 

  28. Kim YC, Byun S, Seok S, Guo G, Xu HE, Kemper B, Kemper JK. Small heterodimer partner and fibroblast growth factor 19 inhibit expression of NPC1L1 in mouse intestine and cholesterol absorption. Gastroenterology. 2019;156:1052–65.

    Article  CAS  Google Scholar 

Download references


We thank Dr. Yang Liu for the technical support.


This work was supported by National Natural Science Foundation of China (81900465 and 82170909); Medicine and Engineering Cross Research Foundation of Shanghai Jiao Tong University (YG2019QNA42).

Author information

Authors and Affiliations



Li-Na Dai and Yu-Ling Zhao performed in vitro and in vivo experiments; Lu Jiang contributed to data analysis. Jun-Kai Yan designed research and drafted manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jun-Kai Yan.

Ethics declarations

Ethics approval and consent to participate

This study was approved by Xin Hua Hospital Animal Use Committee.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

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

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

Dai, LN., Zhao, YL., Jiang, L. et al. Changes in the intestinal expression of drug metabolism-related genes in a piglet model of parenteral nutrition. Nutr Metab (Lond) 19, 18 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: