Pantothenic acid is a precursor of two coenzymes, coenzyme A (CoA) and acyl-carrier-protein. The coenzymes of this vitamin participate in various metabolic reactions, such as glucose, fatty acids and amino acids entering into energy-yielding tricarboxylic acid (TCA) cycle, fatty acid oxidation and synthesis, cholesterol synthesis, acetylcholine synthesis, and heme synthesis etc. [1, 2]. Its importance is highlighted by the adverse effects of pantothenic acid deficiency (PAD) in mammals such as rats, cats, and pigs, including growth depression, skin lesions, diarrhea, loss of hair [2–6]. Also in poultry studies, PAD results in growth retardation, poor feathering, dermatosis, and high mortality in chicks, turkeys, geese, and ducks [7–13]. It has been demonstrated extensively that pantothenic acid can keep the structure of intestine integrity and maintain the intestinal function of animals [3, 14–16]. Intestinal hypofunction was shown to be a major consequence of PAD in rats, dogs, cats, chicks, and fish, such as intestinal ulceration, diarrhea, and colitis [4, 5, 17–21]. Rats deficient in pantothenic acid exhibit duodenitis and duodenal ulcers [18], as well as duodenal changes including eventual atrophy of crypts, diminution in size of villi [16]. In addition, previous studies in fish showed that PAD decreased intestinal digestive and absorptive capacities by reducing the activities of both intestinal brush border enzymes and digestive enzymes [14, 21].
In addition, pantothenic acid can protect cell membrane against peroxidative damage by increasing glutathione content [22–24]. Previous studies have shown that dietary PAD could lead to oxidative stress in geese [9] and fish [21]. It is proposed that oxidative stress induced by PAD leads to intestinal injury and hypofunction. So far, the detailed mechanisms of growth depression and intestinal hypofunction due to PAD are still unclear. In order to understand the underlying mechanisms, we established a PAD duck model, identified an overview of underlying processes and the extend of alterations of intestinal mucosa using a proteomic approach.
The mortality of PAD ducks was greater than the control (CON) birds (P < 0.001; Table 1). In comparison with the CON group, average daily weight gain (ADG) and average daily feed intake (ADFI) were declined in the PAD group, while feed conversion ratio (FCR) was increased (P < 0.001; Table 1).
Variable | PAD | CON | SEM | P-value |
---|---|---|---|---|
Mortality (%) | 68.7a | 0b | 3.16 | < 0.001 |
ADG (g/d) | 10.6b | 31.3a | 1.13 | < 0.001 |
ADFI (g/d) | 16.8b | 40.1b | 1.30 | < 0.001 |
FCR (g:g) | 1.66a | 1.28b | 0.03 | < 0.001 |
ADG, average daily weight gain; ADFI, average daily feed intake; FCR, feed conversion ratio; and SEM, standard error of the mean
a, b Mean values with unlike superscript letters were significantly different (P < 0.05). Data were analyzed by the Student’s t test. Data are expressed as mean and SEM (n = 8)
Plasma pantothenic acid concentration decreased by 85%, while plasma alkaline phosphatase (ALP) activity decreased by approximately 70% as a result of PAD when compared to the CON group (P < 0.001; Table 2). Compared to the CON group, plasma glucose and insulin contents were declined in the PAD group (P < 0.01; Table 2). Plasma glucagon did not differ between the PAD and CON group (P > 0.05; Table 2). Compared to the CON group, plasma malondialdehyde (MDA) content was increased, while plasma total superoxide dismutase (T-SOD) activity was decreased in the PAD group (P < 0.001; Table 2).
Variable | PAD | CON | SEM | P-value |
---|---|---|---|---|
Pantothenic acid (nmol/L) | 121b | 800a | 54.5 | < 0.001 |
Glucose (mmol/L) | 7.33b | 9.91a | 0.14 | < 0.001 |
ALP (U/L) | 236b | 740a | 27.7 | < 0.001 |
Insulin (μIU/mL) | 6.11b | 8.84a | 0.38 | 0.008 |
Glucagon (pg/mL) | 147 | 147 | 5.72 | 0.971 |
MDA | 6.47a | 4.18b | 0.21 | < 0.001 |
T-SOD | 80.2b | 103a | 2.24 | 0.020 |
ALP, alkaline phosphatase; MDA, malondialdehyde; T-SOD, total superoxide dismutase; SEM, standard error of the mean
a, b Mean values with unlike superscript letters were significantly different (P < 0.05). Data were analyzed by the Student’s t test. Data are expressed as mean and SEM (n = 8)
Dietary PAD resulted in morphological alterations of the small intestine of Pekin ducks. Compared to the CON group, PAD reduced villus height and villus surface area of duodenum, jejunum, and ileum (P < 0.001; Table 3). Villus width and crypt depth of both duodenum and jejunum were not affected by PAD in ducks (P > 0.05; Table 3). Compared to the CON group, PAD decreased crypt depth of ileum in ducks (P < 0.001; Table 3) but not villus width of ileum (P > 0.05; Table 3).
Items | PAD | CON | SEM | P-value |
---|---|---|---|---|
Duodenum | ||||
Villus height (μm) | 666b | 893a | 27.7 | < 0.001 |
Villus width (μm) | 108 | 130 | 5.60 | 0.093 |
Crypt depth (μm) | 175 | 214 | 9.90 | 0.153 |
Villus surface area (mm2) | 0.23b | 0.36a | 0.018 | 0.001 |
Jejunum | ||||
Villus height (μm) | 467b | 678a | 42.4 | 0.033 |
Villus width (μm) | 121 | 134 | 7.20 | 0.249 |
Crypt depth (μm) | 123 | 158 | 8.30 | 0.076 |
Villus surface area (mm2) | 0.18b | 0.28a | 0.019 | 0.010 |
Ileum | ||||
Villus height (μm) | 305b | 534a | 32.4 | < 0.001 |
Villus width (μm) | 118 | 122 | 3.40 | 0.618 |
Crypt depth (μm) | 126b | 194a | 10.5 | 0.002 |
Villus surface area (mm2) | 0.11b | 0.20a | 0.014 | < 0.001 |
SEM, standard error of the mean
a, b Mean values with unlike superscript letters were significantly different (P < 0.05). Data were analyzed by the Student’s t test. Data are expressed as mean and SEM (n = 8)
A total of 22,973 peptide belonging to 3345 proteins were identified in the duodenum mucosa of two groups. Comparisons of the relative abundance of proteins from mucosa of PAD ducks with those of CON ducks showed that a total of 421 proteins showed a fold change (FC) > 1.5, of which 198 proteins were up-regulated and 223 proteins were down-regulated. The complete list of proteins altered by PAD is presented in Additional file 1.
All the differentially expressed proteins in mucosa in response to PAD were used to conduct GO categories of biological process, cellular component, and molecular function, and pathway analysis. As shown in Fig. 1, the top 15 enriched terms under biological process included oxidation-reduction process, small molecule metabolic process, generation of precursor metabolites and energy, purine nucleoside triphosphate metabolic process, monocarboxylic acid metabolic process, nucleoside triphosphate metabolic process, carboxylic acid metabolic process, ATP metabolic process, oxoacid metabolic process, ribonucleoside monophosphate metabolic process, organic acid metabolic process, nucleoside monophosphate metabolic process, ribonucleoside triphosphate metabolic process, coenzyme metabolic process, and cellular respiration. The top 15 enriched terms under cellular component included extracellular exosome, extracellular vesicle, extracellular organelle, vesicle, mitochondrion, cytoplasmic part, myelin sheath, cytoplasm, brush border, adherens junction, organelle inner membrane, focal adhesion, cell-substrate adherens junction, mitochondrial inner membrane, and cluster of actin-based cell projections. The top 15 enriched terms under molecular function included cell adhesion molecule binding, oxidoreductase activity, coenzyme binding, identical protein binding, cadherin binding, actin filament binding, actin binding, aminopeptidase activity, RNA binding, fatty-acyl-CoA binding, exopeptidase activity, cytoskeletal protein binding, hydro-lyase activity, glutathione transferase activity, and S100 protein binding.
Based on the pathway analysis by Kyoto Encyclopedia of Genes and Genomes (KEGG), the significantly affected pathways were glycolysis and gluconeogenesis, amino acid metabolism, cori cycle, fatty acid beta oxidation, striated muscle contraction, TCA cycle, trans-sulfuration pathway, PPAR signaling pathway, cytoplasmic ribosomal proteins, vitamin A and carotenoid metabolism, oxidative stress, regulation of actin cytoskeleton, glutathione metabolism, and oxidative phosphorylation (Fig. 2). Among them, energy metabolism (glycolysis and gluconeogenesis, fatty acid beta oxidation, TCA cycle, and oxidative phosphorylation), regulation of actin cytoskeleton, and oxidative stress were selected in the present study and the proteins involved in these processes are listed in Table 4.
UniProtKB ID | Protein name | Short name | Fold change* | P-Value |
---|---|---|---|---|
Glycolysis and gluconeogenesis | ||||
U3IHG8 | Fructose-bisphosphate aldolase | ALDOB | −5.47 | 5.90E-09 |
U3IVG9 | Hexokinase domain containing 1 | HKDC1 | −4.34 | 1.13E-13 |
U3IR52 | Alpha-enolase | ENO1 | −3.75 | 6.19E-08 |
U3J2H8 | Fructose-bisphosphatase 1 | FBP1 | −2.42 | 1.35E-10 |
U3J1L1 | Glyceraldehyde-3-phosphate dehydrogenase | GAPDH | −3.02 | 9.49E-06 |
U3I0F9 | Pyruvate kinase | PKM | −3.28 | 1.94E-07 |
U3IZA5 | ATP-dependent 6-phosphofructokinase | PFKP | −2.33 | 3.93E-05 |
U3I8D8 | Triosephosphate isomerase | TPI1 | −2.93 | 9.45E-05 |
U3ILF5 | Phosphoglycerate kinase | PGK1 | −2.01 | 1.71E-06 |
U3IR48 | Dihydrolipoyl dehydrogenase | DLD | −1.66 | 2.09E-02 |
U3IEW2 | Pyruvate dehydrogenase E1 beta subunit | PDHB | − 1.58 | 1.79E-02 |
U3I939 | Fructose-bisphosphate aldolase | ALDOA | 1.89 | 1.34E-03 |
U3IR68 | Hexokinase 1 | HK1 | 3.31 | 1.12E-02 |
Fatty acid beta oxidation | ||||
U3J4Z9 | Acyl-CoA synthetase long chain family member 5 | ACSL5 | −3.06 | 6.74E-09 |
U3I9A1 | Acyl-CoA dehydrogenase family member 11 | ACAD11 | −1.96 | 1.34E-02 |
U3ITA9 | Medium-chain specific acyl-CoA dehydrogenase | ACADM | −1.85 | 1.79E-05 |
U3IHS8 | Carnitine O-acetyltransferase | CRAT | −1.78 | 6.54E-03 |
U3IR48 | Dihydrolipoyl dehydrogenase | DLD | −1.66 | 2.09E-02 |
U3I6S1 | Hydroxyacyl-CoA dehydrogenase trifunctional multienzyme complex subunit beta | HADHB | −1.65 | 1.60E-03 |
U3IEF4 | Enoyl-CoA delta isomerase 2 | ECI2 | −1.65 | 4.61E-02 |
U3J928 | Acyl-coenzyme A oxidase | ACOX1 | −1.51 | 1.68E-05 |
U3IDQ1 | Acyl-coenzyme A oxidase | ACOX2 | −1.74 | 3.50E-05 |
TCA cycle | ||||
U3IR48 | Dihydrolipoyl dehydrogenase | DLD | −1.66 | 2.09E-02 |
R0J775 | Aconitase 1 (Fragment) | ACO1 | −1.64 | 2.75E-10 |
U3IC15 | Aconitate hydratase, mitochondrial | ACO2 | −2.13 | 9.22E-11 |
R0JXM5 | Malate dehydrogenase (Fragment) | MDH1 | −1.78 | 2.25E-05 |
U3IA60 | Malate dehydrogenase | MDH2 | −2.89 | 4.45E-08 |
U3IEW2 | Pyruvate dehydrogenase E1 beta subunit | PDHB | −1.58 | 1.79E-02 |
U3J597 | Isocitrate dehydrogenase [NADP] | IDH1 | 1.58 | 5.05E-06 |
Oxidative phosphorylation | ||||
U3J532 | NADH:ubiquinone oxidoreductase subunit A5 | NDUFA5 | −1.58 | 1.94E-02 |
R0LLX6 | NADH dehydrogenase [ubiquinone] 1 alpha subcomplex subunit 6 (Fragment) | NDUFA6 | −1.56 | 4.97E-02 |
U3J175 | ATP synthase, H+ transporting, mitochondrial Fo complex subunit B1 | ATP5F1 | −1.80 | 2.12E-03 |
R0LYJ7 | ATP synthase subunit d, mitochondrial (Fragment) | ATP5H | −1.79 | 2.69E-04 |
R0LIL9 | ATP synthase subunit O, mitochondrial (Fragment) | ATP5O | −1.60 | 2.16E-03 |
Regulation of actin cytoskeleton | ||||
U3IFN5 | Villin 1 | VIL1 | −3.12 | 1.13E-08 |
U3IY96 | Ezrin | EZR | −2.00 | 1.61E-05 |
U3IRY0 | Vimentin | VIM | 2.66 | 5.54E-04 |
U3J6G2 | Tropomyosin 3 | TPM3 | 4.51 | 1.36E-02 |
U3I7J8 | Tropomyosin alpha-1 chain | TPM1 | 9.58 | 6.76E-04 |
U3IA79 | Myosin light chain 1 | MYL1 | 11.47 | 2.50E-06 |
R0LM85 | Myosin light chain kinase (Fragment) | MYLK | 2.88 | 1.83E-02 |
U3I4I9 | Actinin alpha 1 | ACTN1 | 2.13 | 5.89E-07 |
U3IZ83 | Fibronectin 1 | FN1 | 2.28 | 4.31E-02 |
U3IR26 | Vinculin | VCL | 1.87 | 4.54E-11 |
U3I935 | Moesin | MSN | 2.07 | 5.48E-03 |
Oxidative stress | ||||
U3J0T0 | Amine oxidase | MAOA | −2.56 | 5.04E-08 |
A0A172QNN4 | Catalase | CAT | −1.84 | 1.49E-05 |
R0JZP2 | Microsomal glutathione S-transferase 1 (Fragment) | MGST1 | −1.57 | 2.66E-02 |
U3I5T1 | Glutamate-cysteine ligase catalytic subunit | GCLC | 1.57 | 9.62E-03 |
Intestinal absorption | ||||
U3IRP7 | Solute carrier family 2 member 2 | SLC2A2 | −5.11 | 3.17E-02 |
U3IFN5 | Villin 1 | VIL1 | −3.12 | 1.13E-08 |
U3IY96 | Ezrin | EZR | −2.00 | 1.61E-05 |
U3IUS6 | Monoacylglycerol O-acyltransferase 2 | MOGAT2 | −2.28 | 5.35E-03 |
TCA, tricarboxylic acid
* Fold change is expressed as the ratio of the pantothenic acid-deficient to the control group. For the down-regulated proteins, the fold change was transformed to the corresponding negative value
Of the proteins enriched in glycolysis and gluconeogenesis, eleven proteins were down-regulated (fructose-bisphosphate aldolase (ALDOB), hexokinase domain containing 1 (HKDC1), alpha-enolase (ENO1), fructose-bisphosphatase 1 (FBP1), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), pyruvate kinase (PKM), ATP-dependent 6-phosphofructokinase (PFKP), triosephosphate isomerase (TPI1), phosphoglycerate kinase (PGK1), dihydrolipoyl dehydrogenase (DLD), and pyruvate dehydrogenase E1 beta subunit (PDHB)), while two proteins were up-regulated (fructose-bisphosphate aldolase (ALDOA) and hexokinase 1 (HK1)).
Nine proteins involved in fatty acid beta oxidation were all down-regulated in the PAD group, including acyl-CoA synthetase long chain family member 5 (ACSL5), acyl-CoA dehydrogenase family member 11 (ACAD11), medium-chain specific acyl-CoA dehydrogenase (ACADM), carnitine O-acetyltransferase (CRAT), DLD, hydroxyacyl-CoA dehydrogenase trifunctional multienzyme complex subunit beta (HADHB), enoyl-CoA delta isomerase 2 (ECI2), acyl-coenzyme A oxidase (ACOX1), and acyl-coenzyme A oxidase (ACOX2).
Of the proteins enriched in TCA cycle, six proteins were down-regulated (DLD, aconitase 1 (Fragment) (ACO1), aconitate hydratase, mitochondrial (ACO2), malate dehydrogenase (fragment) (MDH1), malate dehydrogenase (MDH2), and pyruvate dehydrogenase E1 beta subunit (PDHB)), while one protein was up-regulated (isocitrate dehydrogenase [NADP] (IDH1)).
Five proteins were involved in the oxidative phosphorylation, which were all down-regulated in the PAD group (NADH:ubiquinone oxidoreductase subunit A5 (NDUFA5), NADH dehydrogenase [ubiquinone] 1 alpha subcomplex subunit 6 (Fragment) (NDUFA6), ATP synthase, H+ transporting, mitochondrial Fo complex subunit B1 (ATP5F1), ATP synthase subunit d, mitochondrial (Fragment) (ATP5H), and ATP synthase subunit O, mitochondrial (Fragment) (ATP5O)).
Of the proteins involved in regulation of actin cytoskeleton, two proteins were down-regulated (Villin 1 (VIL1) and Ezrin (EZR)), while nine proteins were up-regulated (vimentin (VIM), tropomyosin 3 (TPM3), tropomyosin alpha-1 chain (TPM1), myosin light chain 1 (MYL1), myosin light chain kinase (Fragment) (MYLK), actinin alpha 1 (ACTN1), fibronectin 1 (FN1), vinculin (VCL), and moesin (MSN)).
Of the proteins involved in oxidative stress, three proteins were all down-regulated (amine oxidase (MAOA), catalase (CAT), and microsomal glutathione S-transferase 1 (Fragment) (MGST1)), while one protein was up-regulated (glutamate-cysteine ligase catalytic subunit (GCLC)).
Besides, we found that four proteins were involved in intestinal absorption, which were all down-regulated by PAD, including solute carrier family 2 member 2 (SLC2A2), VIL1, EZR, and monoacylglycerol O-acyltransferase 2 (MOGAT2).
To validate the isobaric tags for relative and absolute quantification (iTRAQ) results, the abundance of ACADM and GAPDH were analyzed by Western blot. Compared to the CON group, the protein expressions of ACADM and GAPDH were decreased in the PAD group (Fig. 3), which were consistent with the iTRAQ results.
Previous studies reported that dietary PAD causes deficiency symptoms such as growth retardation, dermatosis, diarrhea, and even death in both mammals [2–6] and poultry [7–13, 19, 25]. In ducks, growth retardation and excessive exudate from the eyes were exhibited in response to PAD [13]. Consistent with previous studies, the PAD ducks showed growth depression, exudate on eyelids, poor feathering, and a high mortality rate in the present study. Furthermore, poor pantothenic acid status was observed in PAD ducks indicated by reduced plasma pantothenic acid concentration, which is a sensitive biomarker for pantothenic acid status [21, 26, 27]. These results indicate that a severe PAD animal model was successfully established.
It has been shown that pantothenic acid can maintain intestinal function by keeping the structure of intestine integrity of animals [3, 14–16]. Insufficient of pantothenic acid would cause intestinal hypofunction in rats, dogs, cats, chicks and fish [4, 5, 17–21]. And a feature common to experimental PAD in different animals is intestinal ulceration [4, 5, 17–20]. Duodenitis and duodenal ulcers [18], as well as duodenal changes including eventual atrophy of crypts, diminution in size of villi, were observed in PAD rats [16]. In agreement with the previous studies, we observed small intestine morphological changes indicated by the reduced villus height and villus surface area in PAD ducks in the present study, indicating damage to the intestinal epithelium.
Furthermore, PAD leads to abnormalities in carbohydrate metabolism. Low fasting blood glucose levels and increased sensitivity to insulin were found in PAD rats and dogs [28–32]. In agreement with previous studies, PAD caused fasting hypoglycemia and decreased plasma insulin level in ducks in the present study, indicating abnormal glucose metabolism.
Together, dietary PAD of ducks resulted in growth retardation, alterations of intestinal morphology and function, and abnormal glucose metabolism. However, limited data are available currently concerning the molecular mechanisms behind. Therefore, we used a proteomic approach, for the first time, to investigate the metabolic disorder of small intestine induced by PAD to explain intestinal hypofunction and growth depression. Proteomic analysis revealed 421 differentially expressed proteins in the mucosa of PAD ducks compared to CON birds, indicating an important impact of pantothenic acid on intestinal function. The differentially expressed proteins are mainly enriched in glycolysis and gluconeogenesis, fatty acid beta oxidation, oxidative phosphorylation, TCA cycle, intestinal absorption, regulation of actin cytoskeleton, and oxidative stress. These processes probably underlie the intestinal hypofunction and growth depression induced by PAD. Notably, the CoA-binding proteins were reduced in the intestinal mucosa of PAD ducks, such as ACSL5, ACAD11, ACADM, HADHB, ECI2, ACOX1, and ACOX2. This finding supports the hypothesis that CoA-binding proteins may be depressed as a result of dietary PAD.
A total of thirteen proteins participating in the glycolysis and gluconeogenesis pathway were altered by PAD. Of these, twelve proteins were involved in glycolysis, two enhanced (ALDOA and HK1) and ten diminished (ALDOB, HKDC1, ENO1, GAPDH, PKM, PFKP, TPI1, PGK1, DLD, and PDHB). HKDC1 and HK1 are two isozymes of hexokinases, which mediate the initial step of glycolysis by catalyzing phosphorylation of D-glucose to D-glucose 6-phosphate [33]. PFKP catalyzes the phosphorylation of fructose 6-phosphate to fructose 1,6-bisphosphate by ATP [34]. ALDOA and ALDOB are two isoforms of aldolase family which are located in skeletal muscle and liver tissue respectively, cleaves fructose-1,6-bisphosphate to triose phosphates [35]. TPI1 catalyzes the interconversion between dihydroxyacetone phosphate and D-glyceraldehyde-3-phosphate [36]. GAPDH catalyzes the conversion of glyceraldehyde 3-phosphate into 1,3-diphosphoglycerate [37]. PGK1 catalyzes one of the two ATP producing reactions in the glycolytic pathway via the reversible conversion of 1,3-diphosphoglycerate to 3-phosphoglycerate [38]. ENO1 catalyzes the formation of phosphoenolpyruvate from 2-phosphoglycerate [39]. PKM catalyzes the second ATP generation reaction in the glycolytic pathway via conversion of phosphoenolpyruvate to pyruvate [40]. PDHB is a subunit of pyruvate dehydrogenase (E1). As the E3 component of pyruvate dehydrogenase complex, DLD oxidizes dihydrolipoic acid to lipoic acid [41]. Pyruvate dehydrogenase complex converts pyruvate to acetyl-CoA [42, 43]. Therefore, ten out of twelve proteins (ALDOB, HKDC1, ENO1, GAPDH, PKM, PFKP, TPI1, PGK1, DLD, and PDHB) were down-regulated in the PAD group suggests that the rate of glycolysis in the intestinal mucosa may be impaired. Besides, one protein (FBP1) involved in gluconeogenesis process was down-regulated by PAD, indicating a depressed gluconeogenesis. FBP1, a rate-limiting enzyme in gluconeogenesis, catalyzes the hydrolysis of fructose 1,6-bisphosphate to fructose 6-phosphate.
In the present study, PAD down-regulated 9 proteins involved in fatty acid beta oxidation, including ACSL5, ACAD11, ACADM, CRAT, DLD, HADHB, ECI2, ACOX1, and ACOX2. Notably, among these, ACSL5, ACAD11, ACADM, HADHB, ECI2, ACOX1, and ACOX2 are belong to CoA-binding proteins, which were all reduced due to PAD. ACSL5 belongs to the acyl-CoA synthetase family, catalyzing free fatty acids into fatty acyl-CoA esters, which play key roles in lipid biosynthesis and fatty acid degradation [44]. ACAD11 and ACADM are two members of fatty acyl-CoA dehydrogenases that catalyze the first step in each cycle of fatty acid beta oxidation [45]. CRAT catalyzes the reversible transfer of an acetyl group from acyl-CoA to carnitine [46]. ECI2 is involved in the beta oxidation of unsaturated fatty acids, catalyzing the formation of 2-trans-enoyl-CoA from 3-cis or trans-enoyl-CoA [47]. HADHB catalyzes the final step of beta-oxidation, in which 3-ketoacyl CoA is cleaved by the thiol group of another molecule of Coenzyme A [48]. ACOX1 and ACOX2 catalyze the first and rate-limiting step of peroxisomal fatty acid beta oxidation, the desaturation of acyl-CoAs to 2-trans-enoyl-CoAs. ACOX1 catalyzes medium to very long straight-chain fatty acids [49], while ACOX2 catalyze the CoA-esters of very long-chain fatty acids, branched-chain fatty acids and the C27-bile acid intermediates [50]. The down-regulation of all these proteins involved in the fatty acid beta oxidation process in PAD ducks may indicate that this process is impaired. This implication is supported by previous finding in rats that PAD reduced CoA and short-chain acyl-CoA contents [51], as well as hepatic peroxisomal fatty acid beta oxidation [6].
PAD down-regulated six proteins involved in the TCA cycle, including DLD, ACO1, ACO2, MDH1, MDH2, and PDHB, and up-regulated one protein, IDH1. PDHB is a subunit of pyruvate dehydrogenase (E1). As the E3 component of pyruvate dehydrogenase complex and α-ketoglutarate dehydrogenase complex, DLD oxidizes dihydrolipoic acid to lipoic acid [41]. Pyruvate dehydrogenase complex converts pyruvate to acetyl-CoA [42, 43]. α-ketoglutarate dehydrogenase catalyzes the oxidation of α-ketoglutarate to form succinyl-CoA, which is a rate-limiting enzyme of the TCA cycle [52, 53]. ACO1 and ACO2 are two isozymes of aconitate hydratase located in cytoplasmic and mitochondrial respectively, which catalyzes the isomerization of citrate to isocitrate via cis-aconitate. MDH1 and MDH2 are two isozymes of malate dehydrogenase located in cytoplasmic and mitochondrial respectively, which catalyzes the reversible interconversion of malate and oxaloacetate [54]. IDH1 is a subunit of isocitrate dehydrogenase, catalyzing the oxidative decarboxylation of isocitrate into α-ketoglutarate [55]. Six out of seven proteins (DLD, ACO1, ACO2, MDH1, MDH2, and PDHB) were down-regulated due to PAD, indicating a decreased TCA cycle.
PAD down-regulated 5 proteins involved in oxidative phosphorylation, including NDUFA5, NDUFA6, ATP5F1, ATP5H, and ATP5PO. NDUFA5 and NDUFA6, the subunits of complex I, are involved in complex I assembly [56, 57]. ATP5F1, ATP5H, and ATP5PO, the subunits of complex V, are involved in complex V assembly [58–60]. The down-regulated expression of proteins involved in oxidative phosphorylation, probably indicates that this process is impaired by PAD.
Together, our duodenum mucosal proteomic analysis revealed glycolysis and gluconeogenesis, fatty acid beta oxidation, TCA cycle, and oxidative phosphorylation are probably impaired in response to PAD, consequently resulting in insufficient energy generation in the small intestine and subsequent growth retardation.
Eleven proteins altered by PAD have actin-binding domains and play direct roles in the organization of structure of the cytoskeleton, including two down-regulated proteins (VIL1 and EZR) and nine up-regulated proteins (ACTN1, FN1, MSN, MYL1, MYLK, VCL, TPM3, TPM1, and VIM). VIL1 and EZR (also known as villin 2; VIL2) are microvillar proteins in intestinal epithelial cells [61]. VIL1 is an epithelial cell-specific actin-binding protein that regulates cell migration, cell death, cell morphology, and epithelial-to-mesenchy-mal transition [62, 63]. EZR is critical for the de novo lumen formation and expansion during villus morphogenesis, and EZR absence resulted in abnormal villus morphogenesis [64, 65].
Actinin is a component of stress fibers and links the cytoskeleton to adherens-type junctions. As one isoform of actinin, ACTN1 plays a major role in cell migration and adhesion [66]. Besides, ACTN1 can directly binds to VCL, and the two proteins cooperate to organize the cytoskeleton at adhesion junctions [67]. FN1, a glycoprotein component of the extracellular matrix, has a role in cell adhesion and migration [68]. VIM induces changes in cell shape, motility, and adhesion during the epithelial to mesenchymal transition [69]. MYLK induces contraction of the perijunctional actomyosin ring through myosin II regulatory light chain phosphorylation, and thereby increases intestinal epithelial permeability when activated [70]. TPM3 and TPM1 are two isoforms of tropomyosins which play important roles in the regulation of assembly, stability, and motility of the intestinal epithelial cells [71]. Collectively, the over expression of nine proteins (ACTN1, FN1, MSN, MYL1, MYLK, VCL, TPM3, TPM1, and VIM) and down-regulation of two proteins (VIL1 and EZR) in response to PAD probably indicates orchestrated regulation of actin cytoskeleton dynamics and the negative impact on intestinal integrity. Remodeling of the cytoskeleton is fundamental in proliferation, apoptosis, cell invasion and metastasis [72]. Therefore, these altered proteins involved in regulation of actin cytoskeleton in response to PAD probably resulted in morphological changes of small intestine, such as atrophy of their intestinal villus. Furthermore, it is reported that ATP depletion uncouples the gate and fence functions of the tight junction and induces actin network dissolution of epithelial cells [73]. The alterations of intestinal morphology and actin cytoskeleton due to PAD in the present study may be attributed to ATP depletion indicated by impaired glycolysis, fatty acid beta oxidation, TCA cycle, and oxidative phosphorylation processes.
PAD down-regulated three proteins involved in the oxidative stress, including MAOA, CAT, and MGST1, and up-regulated one protein, GCLC. MAOA metabolizes dopamine to dihydroxyphenylacetic acid and H2O2, a potential source of reactive oxygen species [74]. Both CAT and MGST1 take part in the oxidative stress defense as its scavenging of H2O2 [75, 76]. GCLC is a subunit of glutamate cysteine ligase which catalyzes the rate-limiting step in reduced glutathione (GSH) synthesis. And glutamate cysteine ligase is often activated to increase cellular GSH content in response to oxidative stress [77]. Collectively, the reduction of MAOA, CAT, and MGST1 in PAD ducks, as well as the enhanced GCLC, indicates small intestinal oxidative stress was induced. This is supported by the results of increased plasma MDA content and decreased T-SOD activity in the present study and is consistent with previous studies in geese [9] and fish [21]. It has been shown that pantothenic acid protects cells against oxidative stress by promoting glutathione levels and cellular repair mechanisms [22–24]. Furthermore, PAD induced oxidative stress in ducks, which may be associated with intestinal injury and morphological alterations.
A novel and important finding of this study is that four proteins (SLC2A2, VIL1, EZR, and MOGAT2) involved in intestinal absorption were down-regulated in duodenum mucosa as a result of PAD. Intestinal SLC2A2 (GLUT2) is known as a means to transfer glucose and fructose from the lumen to the bloodstream and, thereby, to provide sugar to tissue. And intestinal SLC2A2 deletion in mice induced glucose malabsorption visualized by the delay in the distribution of oral sugar in tissues, as well as decreased microvillus length and body weight gain [78]. It is reported pantothenic acid appears to be part of a glucose carrier system [79], therefore, PAD may direct reduce SLC2A2 and resulted in abnormal glucose absorption and hypoglycemia in ducks. Villin (VIL1 and EZR) is also involved in the absorptive and secretory function of epithelial cells by modulating F-actin polymerization/depolymerization. Specifically, villin-depleted mice showed a reduction in intestinal glucose absorption [80]. Decreased protein expression of SLC2A2, VIL1, and EZR due to PAD probably impair glucose absorption system, which may provide a possible explanation for fasting hypoglycemia that is seen. MOGAT2 plays a central role in absorption of dietary fat in the small intestine by catalyzing the resynthesis of triacylglycerol in enterocytes [81]. MGAT2 deficient specifically in the small intestine showed a delay in fat absorption in mice [82]. Therefore, a reduction of SLC2A2, VIL1, EZR, and MOGAT2 suggests that glucose and fat malabsorption in the small intestine may be induced by PAD, which probably leads to growth depression observed. This is in line with the previous findings in fish that PAD decreased the digestive and absorptive capacities indicated by the reduced the activities of both intestinal brush border enzymes and digestive enzymes [14, 21].
PAD caused growth retardation, fasting hypoglycemia, morphological alterations of the small intestine, and oxidative stress in ducks. We have performed a proteomic profiling further to investigate the effects of PAD on duodenum mucosal proteome of ducks. The results indicated that PAD may suppress energy generation processes such as glycolysis and gluconeogenesis, fatty acid beta oxidation, TCA cycle, and oxidative phosphorylation, leading to impaired ATP generation. Furthermore, PAD may induce glucose and fat malabsorption in the small intestine indicated by four diminished proteins involved. Besides, PAD probably leads to negative impacts on intestinal integrity and morphology because most of proteins involved in regulation of actin cytoskeleton were up-regulated. To be sum, PAD causes intestinal hypofunction and growth depression probably by impairing glycolysis and gluconeogenesis, fatty acid beta oxidation, TCA cycle, oxidative phosphorylation, actin cytoskeleton, and intestinal absorption processes. These findings add to our understanding of the mechanisms of intestinal hypofunction due to PAD.
All experimental procedures with ducks were performed according to the Guidelines for Experimental Animals established by the Ministry of Science and Technology (Beijing, China). Ethical approval on animal survival was given by the animal ethics committee of the Institute of Animal Sciences (IAS), Chinese Academy of Agricultural Sciences (CAAS, Beijing, China) with the following reference number: IASCAAS-2019–19.
A total of 128 one-day-old male white Pekin ducks (Anas platyrhynchos) were assigned randomly to one of two dietary treatments of 8 replicate pens with 8 birds per pen. The ducks were either fed with a PAD diet or a pantothenic acid adequate (control, CON) diet. From hatch to 16 days of age, all ducks were reared in pens in a temperature-controlled room with feed and water ad libitum. And the light was continuously on.
The pantothenic acid-deficient basal diet was formulated in line with NRC (1994) [83], containing 4.65 mg pantothenic acid /kg of diet (Table 5). To produce the PAD diet and control diet, the basal diets were supplemented with 0 and 8 mg/kg of crystalline calcium pantothenate respectively (purity, 99%; Xinfu Technology Co. Ltd., Hangzhou, China). The pantothenic acid concentration for the control diet (12.65 mg/kg) met the recommendation provided by NRC (1994) [83].
Item | Value |
---|---|
Ingredient, % | |
Corn | 79.7 |
Soy isolate protein | 16.0 |
Limestone | 1.0 |
Dicalcium phosphate | 1.6 |
Vitamin and trace mineral premixa | 1.0 |
Sodium chloride | 0.3 |
DL-Methionine | 0.3 |
L-Lysine·HCl | 0.1 |
Total | 100.0 |
Calculated composition | |
Metabolizable energyb, MJ/kg | 13.35 |
Crude protein | 20.39 |
Calcium | 0.93 |
Nonphytate phosphorus | 0.43 |
Lysine | 1.17 |
Methionine | 0.57 |
Methionine + cysteine | 0.80 |
Threonine | 0.77 |
Tryptophan | 0.19 |
Arginine | 1.38 |
Pantothenic acidc, mg/kg | 4.65 |
a Supplied per kilogram of total diet: Cu (CuSO4•5H2O), 10 mg; Fe (FeSO4•7H2O), 60 mg; Zn (ZnO), 60 mg; Mn (MnSO4•H2O), 80 mg; Se (NaSeO3), 0.3 mg; I (KI), 0.2 mg; choline chloride, 1000 mg; vitamin A (retinyl acetate), 10,000 IU; vitamin D3 (Cholcalciferol), 3000 IU; vitamin E (DL-α-tocopheryl acetate), 20 IU; vitamin K3 (menadione sodium bisulfate), 2 mg; thiamin (thiamin mononitrate), 2 mg; riboflavin, 10 mg; pyridoxine hydrochloride, 4 mg; cobalamin, 0.02 mg; nicotinic acid, 50 mg; folic acid, 1 mg; biotin, 0.2 mg
b The value is calculated according to the AME of ducks (Ministry of Agriculture of China, 2012)
c The value was based on high performance liquid chromatography coupled with triple quadrupole mass spectrometry
At 16 days of age, after overnight fasting, duck weight and feed intake from each pen were recorded to calculate ADG, ADFI, and FCR. ADFI and FCR were all corrected for mortality. Two ducks from each pen were randomly selected for sampling as follows. Blood was taken via jugular vein to harvest plasma, stored at − 20 °C until further analysis. Thereafter, these selected ducks were sacrificed by CO2 inhalation and intestinal tissues were obtained immediately. The intestinal sections were divided into the duodenum, jejunum, and ileum. And they were rinsed with physiological saline and then cut into 1 cm length segments, fixed in 10% neutral formalin used for histological analysis. Duodenum mucosa from the remaining segment was obtained as described previously [84], snap frozen in liquid nitrogen, and stored at − 80 °C until further analysis.
Feed and plasma pantothenic acid concentrations were measured by HPLC (Agilent 1290) coupled with triple quadrupole mass spectrometry (Agilent 6470) according to the methods described previously [85]. Before LC/MS analysis, feed and plasma samples were pretreated according to the methods described previously [86, 87]. The peak was identified by the pure authentic standards (Sigma-Aldrich, St. Louis, MO, USA).
Plasma glucose and ALP were determined by commercial kits following the manufacturer’s instructions (BioSino Bio-technology and Science Inc., Beijing, China). Plasma MDA, T-SOD, insulin, and glucagon were measured by commercial kits following the manufacturer’s protocols (Nanjing Jiancheng Institute of Bioengineering, Nanjing, Jiangsu, China).
The duodenum, jejunum, and ileum sections were embedded in paraffin and transversely sectioned in (4 μm thick) and stained with hematoxylin and eosin following deparaffinization and dehydration. Intestinal tissues and structures were observed using a BH2 Olympus microscope (Olympus, Tokyo, Japan) and analyzed using an image analysis system (Olympus 6.0). Villus height, villus width, crypt depth, and villus surface area were assessed following the method as described previously [84].
Three individual duodenum mucosa samples were randomly chosen from each group to conduct the iTRAQ assays. Proteins were extracted and digested as described previously [88]. The digested samples were labelled following the manufacturer’s protocols with an 8-plex iTRAQ kit (AB SCIEX, Foster City, USA). The PAD samples were labelled with iTRAQ tags 113, 114, and 115, while the CON samples were labelled with tags 116, 117, and 118. The labelled samples were combined and fractionated into 20 fractions by HPLC (DINOEX Ultimate 3000 BioRS, Thermo Fisher, Waltham, MA, USA). We performed LC-electrospray ionization-MS/MS analysis on a Triple TOF 5600 plus system (AB SCIEX, Framingham, USA). The collected raw MS/MS data were searched according to the Uniprot-Swissport Anas platyrhynchos Database UniProt_Mallard_8839 using ProteinPilot Software (version 4.5, AB SCIEX). To minimize false-positive results, we counted only peptides at the 95% CI and a false discovery rate < 0.01, containing at least one unique peptide. For protein quantitation, a differentially expressed protein containing at least two unique spectra was used with a FC value > 1.5 or < − 1.5 with P < 0.05 between the PAD and CON groups. For functional annotation, the differentially expressed proteins induced by PAD were performed KEGG pathway and Gene Ontology (GO) enrichment analysis using ClueGo software as described previously [88].
Western blot analysis of two proteins, ACADM and GAPDH, were performed following the method as described [89]. Primary antibodies (1 μg/ml) against ACADM (ab92461; Abcam) and GAPDH (HX1828; Huaxingbio) were used. Histone H3 (A2348; ABclonal) served as a loading control.
For results of growth performance, plasma parameters, and intestinal mucosal histomorphology, data were analyzed using the Student’s t test procedures of SAS software (SAS Institute Inc., 2011). The replicate pen of 8 ducks for growth performance or two ducks for other indices served as the experimental unit. The variability in the data was expressed as the standard error of the means (SEM). Differences between means were considered statistically significant at P < 0.05.
ADG | average daily weight gain |
ADFI | average daily feed intake |
ALP | alkaline phosphatase |
CoA | coenzyme A |
CON | control |
FC | fold change |
FCR | feed conversion ratio |
GO | Gene Ontology |
HPLC | high performance liquid chromatography |
iTRAQ | isobaric tags for relative and absolute quantification |
KEGG | Kyoto Encyclopedia of Genes and Genomes |
MDA | malondialdehyde |
PAD | pantothenic acid deficiency |
SEM | standard error of the mean |
TCA | tricarboxylic acid |
T-SOD | total superoxide dismutase |
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
We thank the members of our laboratory for their helpful comments on the manuscript. We thank Wuhan GeneCreate Biological Engineering Co., Ltd. for performing the mass spectrometric analysis.
JT performed the study, was involved in all aspects of analysis and drafted the manuscript. YF, MX, ZZ and SH were involved in experimental design. BZ, YW, ZG performed data analysis and sample analysis. JT, YF, SL and SH were involved in the data interpretation and writing. All authors read and approved the final manuscript.
This work was supported by the earmarked fund for the China Agriculture Research System (CARS-42), the science and technology innovation project of the Chinese Academy of Agricultural Sciences (CXGC-IAS-09). The funding bodys had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
The data sets supporting the conclusions of this article are included within the article and its additional file. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE (https://www.ebi.ac.uk/pride/archive/) partner repository with the data set identifier PXD026710.
All experimental procedures with ducks were performed according to the Guidelines for Experimental Animals established by the Ministry of Science and Technology (Beijing, China). Ethical approval on animal survival was given by the animal ethics committee of the Institute of Animal Sciences (IAS), Chinese Academy of Agricultural Sciences (CAAS, Beijing, China) with the following reference number: IASCAAS-2019–19.
Not applicable.
All authors declare that they have no conflict of interests.
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
15.
16.
17.
18.
19.
20.
21.
22.
23.
24.
25.
26.
27.
28.
29.
30.
31.
32.
33.
34.
35.
36.
37.
38.
39.
40.
41.
42.
43.
44.
45.
46.
47.
48.
49.
50.
51.
52.
53.
54.
55.
56.
57.
58.
59.
60.
61.
62.
63.
64.
65.
66.
67.
68.
69.
70.
71.
72.
73.
74.
75.
76.
77.
78.
79.
80.
81.
82.
83.
84.
85.
86.
87.
88.
89.