- Open Access
Maternal diet deficient in riboflavin induces embryonic death associated with alterations in the hepatic proteome of duck embryos
Nutrition & Metabolismvolume 16, Article number: 19 (2019)
Maternal riboflavin deficiency (RD) induces embryonic death in poultry. The underlying mechanisms, however, remain to be established and an overview of molecular alterations at the protein level is still lacking. We investigated embryonic hepatic proteome changes induced by maternal RD to explain embryonic death.
A total of 80 45-week-old breeding female ducks were divided into two groups of 40 birds each, and all birds were raised individually for 8 weeks. All the female ducks received either a RD or a riboflavin adequate (control, CON) diet, which supplemented the basal diet with 0 or 10 mg riboflavin /kg of diet respectively.
The riboflavin concentrations of maternal plasma and egg yolk, as well as egg hatchability declined markedly in the RD group compared to those in the CON group after 2 weeks, and declined further over time. The hepatic proteome of E13 viable embryos from 8-week fertile eggs showed that 223 proteins were upregulated and 366 proteins were downregulated (> 1.5-fold change) in the RD group compared to those in the CON group. Pathway analysis showed that differentially expressed proteins were mainly enriched in the fatty acid beta-oxidation, electron transport chain (ETC), and tricarboxylic acid (TCA) cycle. Specifically, all the proteins involved in the fatty acid beta-oxidation and ETC, as well as six out of seven proteins involved in the TCA cycle, were diminished in the RD group, indicating that these processes could be impaired by RD.
Maternal RD leads to embryonic death of offspring and is associated with impaired energy generation processes, indicated by a number of downregulated proteins involved in the fatty acid beta-oxidation, ETC, and TCA cycle in the hepatic of duck embryos. These findings contribute to our understanding of the mechanisms of liver metabolic disorders due to maternal RD.
Riboflavin (vitamin B2) is an essential precursor of flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD). A number of flavin-dependent proteins that utilize FMN or/and FAD, so-called flavoproteins, participate in a range of redox reactions in the tricarboxylic acid (TCA) cycle, fatty acid beta-oxidation, amino acid degradation, and electron transport chain (ETC), among others [1, 2]. Due to its involvement in primary metabolic pathways, maternal riboflavin deficiency (RD) causes abnormal development of offspring. In mammals, insufficient riboflavin supplementation in rat dams during lactation markedly decreased riboflavin concentrations in the liver, carcass, and milk , and reduced the body weight of offspring by 20% . In poultry, maternal RD in laying hens decreased the riboflavin concentrations in egg yolk and resulted in embryonic death [5, 6]. Inherited RD chicken embryos have delayed and aberrant feather development (clubbed down), fatty livers, and extensive cutaneous hemorrhaging, which suddenly die at mid-incubation [7,8,9,10]. An explanation of the sudden death in RD embryos is due to the inhibition of the various critical biological pathways, leading to energy depletion, lipid accumulation in the liver, and severe hypoglycemia [8, 9]. Of these, the beta-oxidation of lipids was severely impaired in RD embryos, as indicated by the reduction in the activity of acyl-CoA dehydrogenase, a flavin-dependent enzyme, and the accumulation of lipids and fatty acid oxidation intermediates in the embryonic liver [7, 9]. Cogburn et al. (2018) recently confirmed that RD in mid-stage embryos leads to a reduced expression of numerous genes involved in the beta-oxidation of lipids and energy depletion at the transcriptional level, including medium chain acyl-CoA dehydrogenase transcripts . However, the underlying mechanisms of embryonic death caused by RD are unknown at the protein level-the functional relevant integration level. Here, we used a proteomic approach to investigate the effects of maternal RD on the hepatic protein levels of offspring.
Materials and methods
Animals ethics statement
This study, which complied with institutional and national guidelines for the care and use of animals, was approved by the Committee of Animal Experiments of the Institute of Animal Sciences, Chinese Academy of Agricultural Sciences. All efforts were made to minimize animal suffering.
Animals and housing
This study was designed to investigate the effects of maternal RD on the hepatic protein levels of offspring using a proteomic approach. A total of 80 45-week-old breeding female white Pekin ducks (Anas platyrhynchos) were obtained from the Pekin duck breeding center (Chinese Academy of Agricultural Sciences) and randomly divided into two treatment groups with 10 replicates (4 birds per replicate) for each group. All female ducks were initially adapted for 2 weeks using a conventional corn-soybean meal diet supplemented with 10 mg riboflavin /kg of diet. After that, the birds received two experimental diets for 8 weeks that supplemented the basal diet with 0 or 10 mg riboflavin /kg of diet. The female ducks were fertilized by male Pekin ducks of the same age fed the experimental diet supplemented the basal diet with 10 mg riboflavin /kg of diet during the experimental period. The birds in the two treatments were housed in individual plastic cages in an environmentally controlled room and received ad libitum access to water and feed. Eighteen hours of light were provided daily from 04:00 to 22:00 h.
All the female ducks were raised with common diets from hatch to 45 weeks of age, and all nutrients met the recommendations for ducks as established in the Nutrient Requirements of Meat-type Ducks of China . The basal diet during the experimental period was riboflavin-deficient containing 1.48 mg free riboflavin/kg of diet (Table 1). The RD diet and CON diet were produced from the basal diet, which was supplemented with 0 or 10 mg crystalline riboflavin/kg diet (purity, 99%; Sigma Aldrich, MO, USA). Except for the riboflavin content of the basal experimental diet, all nutrients met the recommendations for laying ducks .
All female ducks were weighed at the start and end of the experiment. One duck from each replicate was randomly selected to collect whole blood from a wing vein twice weekly. Blood samples were collected into heparin sodium-anticoagulant tubes and centrifuged at 1520×g for 10 min to obtain plasma. Plasma was stored at − 20 °C until assayed for free riboflavin concentration. Meanwhile, eggs were collected and weighed every day during the experimental period, and egg production was recorded. Twelve eggs per group were selected randomly from different ducks on the last day of every two week period and were broken to collect yolk samples, which were then stored at − 20 °C until riboflavin analysis using HPLC. Sixteen eggs from each replicate (four eggs from each duck) were selected to measure fertility and hatchability in the incubator at weekly intervals from 1 to 8 weeks of the experiment.
During the incubation of 7-week eggs, all eggs from the RD group were examined for livability of the embryo at seven days and followed by every three days by means of transmitted light. During the incubation of 8-week eggs, one viable embryo at three embryonic ages (E13, E20, and E27) from each replicate was selected randomly. The embryos were sacrificed by cervical dislocation, then liver samples were collected and frozen in liquid nitrogen immediately, and subsequently stored at − 80 °C for further analyses.
The riboflavin concentration in feed, plasma, and egg yolk were determined by reversed-phase high performance liquid chromatography (HPLC) according to the methods described previously [12, 13]. Before HPLC analysis, feed and plasma samples were prepared according to the method described previously [14, 15], while egg yolk samples were prepared according to the method for animal tissue described previously . The peak was identified and quantified by analysis of authentic standard (Sigma Aldrich).
Total lipids were extracted by homogenizing minced liver tissue samples in chloroform-methanol (2:1) as described previously . The extracts were evaporated under a stream of nitrogen, weighed, and resuspended in chloroform-methanol (2:1) containing 0.01% butyrated hydroxytoluene. Measurements of total lipids in the liver were performed as described previously . Aliquots were dried and resuspended in 1-butanol for analysis of triglyceride (TG) using commercial kits according to the manufacturer’s instructions (BioSino Bio-technology and Science Inc., Beijing, China).
A total of six liver samples (three biological replicates per group) were used to conduct the isobaric tags for relative and absolute quantification (iTRAQ) assays. Each liver sample was ground in liquid nitrogen. The grinded powder was lysed in a solution containing 200 μl L3 buffer (50 mM Tris-Cl, pH 8, 8 M urea, 2 M thiourea, 2 M EDTA, 1 × protease inhibitors cocktails), 800 μl of ice-cold acetone, and 10 mM DTT. The suspensions were incubated at − 20 °C for 2 h. The precipitate pellets were obtained via centrifugation at 12,000×g for 20 min at 4 °C and subsequently resuspended in 800 μl of ice-cold acetone and 10 mM DTT. The suspensions were further centrifuged at 12,000×g for 20 min at 4 °C to collect the precipitated pellets, and then vacuum dried. The dried precipitated pellets were dissolved in 200 μl L3 buffer. Subsequently, total protein concentration was measured using the Bradford assay.
For each sample, 100 μg of protein was reduced, alkylated, and digested with trypsin according to the manufacturer’s protocol (Applied Biosystems, Framingham, MA, USA). Each digested sample was labelled with iTRAQ 8-plex reagents (AB Sciex, Foster City, USA) according to the manufacturer’s instructions. The RD samples were labelled with iTRAQ tags 113, 114, and 115, and the CON samples were labelled with tags 116, 117, and 118. Labelled samples were mixed and fractionated into 20 fractions by HPLC (DINOEX Ultimate 3000 BioRS, Thermo Fisher, Waltham, MA, USA) using a Durashell C18 column (5 μm, 100 Å, 4.6 × 250 mm). LC-electrospray ionization-MS/MS analysis was carried out with a Triple TOF 5600 plus system (AB SCIEX, Framingham, USA). The original MS/MS file data for identification and quantitation were analysed against the database UniProt_Mallard_8839 using ProteinPilot Software version 4.0 (AB SCIEX). To minimize the false discovery rate, a threshold for protein identification was applied. Only unique peptides whose confidences are more than 95% were contained in the iTRAQ labelling quantification and used for further analysis.
For analysis of the proteomic results, the relative expressions of identified proteins were based on the ratio of the reporter ions of the peptides between the two groups (RD vs CON). A protein was considered differentially expressed when the protein had both a fold change (FC) greater than 1.5 and a P-value less than 0.05.
To enrich the differentially expressed proteins with respect to specific functional terms, the protein lists were analysed using ClueGo software (http://www.ici.upmc.fr/cluego/) with the Gene Ontology (GO) database (release date: February 2018). A pathway enrichment analysis of the differentially expressed proteins  was performed using ClueGo software and applying database from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (release date: February 2018).
All the data were analysed using a T-test in SAS software (SAS Institute Inc., 2003). 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. The relative expression of identified proteins was based on the ratio of the reporter ions of the peptides in the RD group to the CON group. Fold changes (FC) were used to determine whether identified proteins were enhanced (FC > 1.5) or diminished (FC < − 1.5) by RD. A fold change (FC) of 1.5 and P < 0.05 was set as the threshold to identify differentially proteins induced by RD.
Maternal growth performance and embryonic growth
The final body weight of maternal ducks fed the RD diet was not different from the CON group (initial body weight: 3880 ± 264 g in the RD group vs 3863 ± 301 g in the CON group; final body weight: 3546 ± 265 g in the RD group vs 3604 ± 297 g in the CON group). Egg weight, egg production, and egg fertility of maternal ducks fed the RD diet were not different from those birds in the CON group from 1 to 8 weeks (P > 0.05, Additional files 1, 2, 3).
The hatchability of eggs from the maternal ducks in the RD group dropped dramatically to approximately half of that in the controls after feeding the RD diet for 2 weeks (P < 0.05, Fig. 1a). Subsequently, from 3 to 5 weeks of the experiment, the hatchability of eggs declined further as riboflavin depletion last longer in the maternal ducks. Strikingly, the hatchability dropped to approximately zero after 6 weeks of riboflavin depletion in maternal ducks.
The embryonic mortality of all fertile RD eggs was approximately zero after 6 weeks of riboflavin depletion in maternal ducks. The curve on the time distribution of 164 dead embryos from the eggs at 7 weeks is shown in Fig. 1b. This curve displays two critical periods with the peaks on the seventh and nineteenth days of incubation, while the former was the major peak of embryonic mortality. Accumulated 54.2, 73.2, and 92.7% of RD embryos were dead within 7, 13, and 19 days of incubation, respectively, indicating that the majority of RD embryos were dead by day 13 of incubation.
Riboflavin concentrations in maternal plasma and egg yolk
Figure 1c shows that, the RD diet reduced the riboflavin concentration in the maternal plasma by 40% of the value present in birds fed the control diet at 2 weeks of the experiment (P < 0.05). Afterward, maternal plasma riboflavin concentration decreased gradually to only 21.6% of the control values at 8 weeks of the experiment (P < 0.05). Figure 1d shows that, the riboflavin concentration in egg yolk was also decreased by RD in the maternal diet by 28.7% at the end of 2 weeks compared with the controls. Then, it was reduced further to only 42% of the control at the end of 4 weeks. Yolk riboflavin depletion was complete at 4 weeks in the RD group, stabilizing at an average of 3.7 μg riboflavin/g yolk.
Total lipids and TG content in maternal plasma and liver, egg yolk, and embryo
Total lipids in the maternal liver, egg yolk, and embryonic liver was not affected by the riboflavin concentration in the diet (P > 0.05, Table 2). The TG content of the maternal liver and egg yolk from the animals fed the RD diet was not significantly different from the controls (P > 0.05, Table 2). However, the TG content in the liver of E27 embryos from the RD group was greater than that in the CON group (P < 0.05, Table 2), even though the E20 embryo was not affected (P > 0.05, Table 2).
Proteomic analysis of the embryonic liver
Using iTRAQ analysis, a total of 25,617 peptide spectral matches were identified, and 3801 proteins were identified in the E13 embryonic livers of two groups. Comparisons of the relative abundance of proteins from the embryonic livers of maternal ducks fed the RD diet with those fed the control diet showed that 223 proteins were enhanced and 366 diminished. The complete list of proteins regulated by maternal RD is presented in Additional file 4.
We performed GO categories of biological process, cellular component, and molecular function, and pathway analysis on the set of 589 differentially expressed proteins in livers from the RD group compared with those in livers from the CON group. As shown in Fig. 2, an enrichment analysis of the significantly enriched under the terms carboxylic acid metabolic process, oxoacid metabolic process, organic acid metabolic process, small molecule metabolic process, oxidation-reduction process, organonitrogen compound metabolic process, small molecule catabolic process, single-organism catabolic process, organic acid catabolic process, carboxylic acid catabolic process, cellular amino acid metabolic process, monocarboxylic acid metabolic process, organic substance catabolic process, alpha-amino acid metabolic process, and small molecule biosynthetic process. The top 15 enriched terms under cellular component included extracellular exosome, extracellular vesicle, extracellular organelle, vesicle, mitochondrion, mitochondrial part, mitochondrial inner membrane, organelle inner membrane, mitochondrial matrix, mitochondrial envelope, mitochondrial membrane, cytosol, myelin sheath, organelle envelope, and envelope. The top 15 enriched terms under molecular function included coenzyme binding, anion binding, nucleotide binding, nucleoside phosphate binding, cell adhesion molecule binding, cadherin binding, identical protein binding, flavin adenine dinucleotide binding, oxidoreductase activity, RNA binding, NAD binding, protein homodimerization activity, NADP binding, intramolecular oxidoreductase activity, and vitamin binding.
The pathway analysis by KEGG on differentially expressed proteins suggested that the significantly affected pathways were amino acid metabolism, fatty acid beta-oxidation, glycolysis and gluconeogenesis, TCA cycle and deficiency of pyruvate dehydrogenase complex, mitochondrial LC-fatty acid beta-oxidation, peroxisomal beta-oxidation of tetracosanoyl-CoA, TCA cycle, PPAR signaling pathway, Cori cycle, tryptophan metabolism, cytoplasmic ribosomal proteins, electron transport chain, and synthesis and degradation of ketone bodies (Fig. 3). Since riboflavin deficiency mainly affected energy generation pathways, such as the TCA cycle, fatty acid beta-oxidation, and ETC, the proteins associated with these processes are listed in Table 3. Notably, some of these proteins are flavin-containing enzymes such as electron transfer flavoprotein subunit alpha (ETFA), flavoprotein-ubiquinone oxidoreductase (ETFDH), short-chain specific acyl-CoA dehydrogenase (ACADS), Acyl-CoA dehydrogenase long chain (ACADL), acyl-CoA dehydrogenase family member 9 (ACAD9), succinate dehydrogenase [ubiquinone] flavoprotein subunit (SDHA), dihydrolipoyl dehydrogenase (DLD), NADH dehydrogenase [ubiquinone] flavoprotein 1 (NDUFV1), NADH dehydrogenase [ubiquinone] 1 alpha subcomplex subunit 9 (NDUFA9), and NADH dehydrogenase [ubiquinone] 1 alpha subcomplex subunit 10 (NDUFA10), which were all diminished in the RD group.
Of the proteins in the enriched TCA cycle, six proteins were downregulated (SDHA, succinate dehydrogenase [ubiquinone] iron-sulfur subunit (SDHB), DLD, succinate--CoA ligase [ADP/GDP-forming] subunit alpha (SUCLG1), dihydrolipoamide S-succinyltransferase (DLST), and fumarate hydratase (FH)), while one protein was upregulated (isocitrate dehydrogenase [NAD] subunit, mitochondrial (IDH3A)). Sixteen proteins were involved in fatty acid beta-oxidation (carnitine O-palmitoyltransferase 1 (CPT1A), carnitine palmitoyltransferase 2 (CPT2), DLD, ETFDH, long-chain-fatty-acid--CoA ligase 1 (ACSL1), ACADL, long-chain-fatty-acid--CoA ligase 5 (ACSL5), solute carrier family 25 member 20 (SLC25A20), ETFA, ACADS, 2,4-dienoyl-CoA reductase (DECR1), trifunctional enzyme subunit alpha (HADHA), carnitine palmitoyltransferase 1B (CPT1B), ACAD9, enoyl-CoA delta isomerase 1 (ECI1), and carnitine O-acetyltransferase (CRAT)), which were all diminished in the RD group. Fourteen proteins were involved in the electron transport chain (SDHA, SDHB, NADH-ubiquinone oxidoreductase 75 kDa subunit (NDUFS1), ATP synthase subunit O (ATP5PO), NADH:ubiquinone oxidoreductase core subunit V2 (NDUFV2), ATP synthase peripheral stalk-membrane subunit b (ATP5PB), ATP synthase F1 subunit gamma (ATP5F1C), NDUFA9, NDUFV1, NDUFA10, NADH dehydrogenase [ubiquinone] 1 alpha subcomplex subunit 12 (NDUFA12), ACAD9, Ubiquinol-cytochrome c reductase core protein 2 (UQCRC2), and NADH:ubiquinone oxidoreductase core subunit S3 (NDUFS3)), which were all diminished in the RD group.
Previous studies described inadequate riboflavin in the diet of laying hens resulted in low egg hatchability and abnormal embryonic development [5, 6]. Similarly, inherited RD chicken embryos presented delayed development, fatty livers, and extensive cutaneous hemorrhaging, and suddenly death at mid-incubation . Recently, Cogburn et al. (2018) showed that RD in mid-stage embryos leads to the reduced expression of numerous genes involved in the beta-oxidation of lipids and energy depletion at the transcriptional level . However, limited data are available currently at the protein level-the functional relevant integration level concerning the molecular mechanisms behind induced by maternal RD. The aim of this study was to investigate the underlying mechanisms of embryonic death induced by maternal RD using a proteomic approach. Proteomic analysis revealed that maternal RD mainly affected proteins involved in the TCA cycle, fatty acid beta-oxidation, and ETC processes, suggesting that these processes are associated with embryonic death. Notably, all flavoproteins or subunits of flavin-dependent enzymes among the RD altered proteins, such as DLD, SDHA, ETFDH, ACADS, ACADL, ETFA, NDUFA9, NDUFV1, NDUFA10, and ACAD9, were greatly reduced in the liver from RD embryos. This finding is consistent with the hypothesis that flavoprotein expression may be downregulated due to a reduced supply of riboflavin in the diet and is in line with previous animal and human studies [8, 20, 21].
The riboflavin concentrations in maternal plasma and egg yolk dropped dramatically in the RD group, indicating that these ducks were riboflavin deficient. The current study demonstrated that maternal diet deficient in riboflavin caused embryonic death in ducks, and the major peak of embryonic mortality was at the early stage of incubation. This finding confirms data from laying hens fed a diet with inadequate riboflavin showed depressed hatchability [5, 6], and is in accordance with the observation that the excessive depletion of riboflavin reserves in hens caused a major peak of embryonic mortality on the fourth day . In the present study, the majority of RD embryos (73.2%) were nonviable within 13 days of incubation, confirming previous findings inherited RD chicken embryos suddenly died mid-incubation (days 13–15) . Therefore, we used a proteomic approach to investigate the metabolic disorders of E13 embryonic liver tissue induced by maternal RD to explain embryonic death. Proteomics analysis revealed 589 differentially expressed proteins in the livers of RD embryos compared to those that were adequately supplied with riboflavin, indicating an important impact of riboflavin on embryonic development. The identical proteins are mainly enriched in the TCA cycle, fatty acid beta-oxidation, and ETC processes based on KEGG analysis.
RD downregulated six proteins involved in the TCA cycle, including DLD, SDHA, SDHB, SUCLG1, DLST, and FH, and upregulated one protein, IDH3A. The results obtained in previous human and animal studies [20, 21], showing that DLD was diminished in the skeletal muscle of riboflavin-responsive multiple acyl-CoA dehydrogenase deficiency (RR-MAD) patient and in RD ducks, consistent with our results. DLD (E3) is a common component of pyruvate and α-ketoglutarate dehydrogenase complexes, converting dihydrolipoic acid and NAD+ into lipoic acid and NADH . Pyruvate dehydrogenase catalyses the oxidative decarboxylation of pyruvate to acetyl-CoA, thereby linking glycolysis to the TCA cycle and fatty acid synthesis [24, 25]. DLST (E2) is also a common component of α-ketoglutarate dehydrogenase, catalysing the conversion of α-ketoglutarate to succinyl-CoA and NADH, which is a rate-limiting enzyme of the TCA cycle [26, 27]. SDHA and SDHB are two subunits of succinate dehydrogenase complex, which catalyses the oxidation of succinate to fumarate . SUCLG1 is a subunit of the heterodimeric enzyme succinate coenzyme A ligase, which catalyses the conversion of succinyl CoA and ADP or GDP to succinate and ATP or GTP . FH catalyses the reversible hydration/dehydration of fumarate to malate . IDH3A is a subunit of isocitrate dehydrogenase, which catalyzes the decarboxylation of isocitrate into alpha-ketoglutarate . Six out of seven proteins (DLD, SDHA, SDHB, SUCLG1, DLST, and FH) were depressed in the RD group, which probably indicates a decreased liver TCA cycle.
RD downregulated 16 proteins involved in fatty acid beta-oxidation, including ETFA, ETFDH, CPT1A, CPT1B, CPT2, ACADS, ACAD9, ACADL, DLD, ACSL1, ACSL5, DECR1, SLC25A20, HADHA, ECI1, and CRAT. Previous studies showed that ACADS, ACAD9, and ETFDH were downregulated in the RD group in humans and animals, which is consistent with our results [20, 21]. CPT1A and CPT1B are two isoforms of carnitine palmitoyltransferase 1 (CPT1), an outer membrane protein that catalyses activated fatty acids into acylcarnitines, which is the first committed and regulated step in mitochondrial fatty acid oxidation [32, 33]. CPT2, an inner membrane protein, catalyses the formation of acyl-CoA from acylcarnitine and CoA . CRAT catalyses the reversible transfer of acyl-CoA from carnitine to free CoA . SLC25A20 facilitates the transfer of acylcarnitine esters in exchange for free carnitine across the mitochondrial membrane . ACADS, ACADL, and ACAD9 belong to the family of fatty acyl-CoA dehydrogenases that catalyse the initial rate-limiting step of the beta-oxidation cycle . HADHA converts medium- and long-chain 2-enoyl-CoA compounds into 3-ketoacyl-CoA. ETFA and ETFDH accept electrons from multiple acetyl-CoA dehydrogenases and subsequently transfer these to ETC . ACSL1 and ACSL5 belong to the acyl-CoA synthetase family, which plays an important role in fatty acid catabolism and de novo lipid synthesis, catalysing the initial fatty acid activation by forming a thioester with CoA . ECI1 and DECR1 have important roles in the metabolism of unsaturated fatty acids in beta oxidation. ECI1 is an auxiliary enzyme in the beta oxidation of unsaturated fatty acids that converts 3-cis or trans-enoyl-CoA to 2-trans-enoyl-CoA . DECR1 participates in the metabolism of (poly)unsaturated fatty enoyl-CoA esters with double bonds in both even- and odd-numbered positions, and catalyses the NADP-dependent reduction of 2,4-dienoyl-CoA to yield trans-3-enoyl-CoA . The decreased expression of these proteins involved in the fatty acid beta-oxidation process may imply that mitochondrial fatty acid beta-oxidation is impaired by RD. Furthermore, downregulation of these proteins involved in the fatty acid beta-oxidation process due to RD could very well explain the observed elevation of liver TG levels in E27 embryos. This explanation is supported by previous findings in humans and animals [9, 20, 21, 42,43,44].
RD downregulated 14 proteins involved in the ETC process, including SDHA, SDHB, NDUFS1, ATP5PO, NDUFV2, ATP5PB, ATP5F1C, NDUFA9, NDUFV1, NDUFA10, NDUFA12, ACAD9, UQCRC2, and NDUFS3. NDUFS1, NDUFS3, NDUFV1, NDUFV2, NDUFA9, NDUFA10, and NDUFA12 are seven subunits of complex I, which play a direct role in complex I assembly . ACAD9 not only plays a physiological role in the beta oxidation of fatty acids, but also serves as an assembly factor for mitochondrial respiratory chain complex I . The reduction of these proteins in RD embryos is in line with the role of riboflavin in complex I assembly . As two subunits of complex II, SDHA and SDHB not only play an important role in the TCA cycle but also feed electrons to the respiratory chain ubiquinone pool . A previous study showed that SDHA mutations caused a complex II deficiency . UQCRC2 is a subunit of complex III, which is required for the assembly and stabilization of the complex . ATP5PO, ATP5PB, and ATP5F1C are three subunits of complex V, which play a direct role in complex V assembly [50, 51]. The downregulated expression of proteins involved in the ETC process, including complex I, complex II, complex III, and complex V, likely indicates that mitochondrial oxidative phosphorylation is impaired by RD, which was supported by previous studies [8, 21, 44, 52].
Our proteomic analysis revealed an impairment of fatty acid beta-oxidation, TCA cycle, and ETC processes in RD embryos. Together, these impaired processes may lead to insufficient ATP production and subsequent embryonic death.
Maternal RD causes embryonic death and abnormal development in ducks. Our analysis of embryonic liver proteomics provides the first global view of the protein level responses to maternal RD and illustrates the impairment of vital biological processes. Proteomic analysis showed that RD mainly diminished the expression of proteins involved in the TCA cycle, fatty acid beta-oxidation, and ETC processes in the livers of embryos, indicating that these processes were impaired and thus energy generation was reduced, which may lead to embryonic death. This finding adds to our understanding of the mechanisms underlying embryonic death as a result of maternal RD.
Electron transport chain
Flavin adenine dinucleotide
High performance liquid chromatography
Isobaric tags for relative and absolute quantification
Kyoto Encyclopedia of Genes and Genomes database
The standard error of the means
Powers HJ. Riboflavin (vitamin B-2) and health. Am J Clin Nutr. 2003;77:1352–60.
Lienhart W-D, Gudipati V, Macheroux P. The human flavoproteome. Arch Biochem Biophys. 2013;535:150–62.
Kirchgessner M, Hirschvogl G, Eder K, RothMaier DA. The effect of dietary riboflavin supplementation during pregnancy and lactation on riboflavin concentrations of milk, liver and carcass of lactating rats. J Anim Physiol an N. 1997;77:44–52.
Roth-Maier DA, Kettler SI, Hirschvogl G, Kirchgessner M. Investigations on the effect of maternal dietary riboflavin supplementation during gravidity and lactation on riboflavin concentrations in blood, liver, and total body of rat offspring. J Anim Physiol an N. 2000;83:158–67.
Squires MW, Naber EC. Vitamin profiles of eggs as indicators of nutritional status in the laying hen: riboflavin study. Poultry Sci. 1993;72:483–94.
Naber EC, Squires MW. Research note: early detection of the absence of a vitamin premix in layer diets by egg albumen riboflavin analysis. Poultry Sci. 1993;72:1989–93.
Abrams VA, Han CC. Riboflavin-deficient chicken embryos: hypoglycemia without dicarboxylic aciduria. Comp Biochem Phys B. 1995;111:233–41.
Cogburn LA, Smarsh DN, Wang X, Trakooljul N, Carré W, White HB. Transcriptional profiling of liver in riboflavin-deficient chicken embryos explains impaired lipid utilization, energy depletion, massive hemorrhaging, and delayed feathering. BMC Genomics. 2018;19:177.
White HB III. Sudden death of chicken embryos with hereditary riboflavin deficiency. J Nutr. 1996;126:1303S–17S.
Maw AJG. Inherited riboflavin deficiency in chicken eggs. Poultry Sci. 1954;33:216–7.
Ministry of Agriculture of China. Nutrient Requirements of Meat-type Ducks of China. Beijing, China: China Agriculture Press; 2012.
Tang J, Xie M, Yang J, Wen Z, Zhu Y, Huang W, Hou S. Riboflavin requirements of white Pekin ducks from hatch to 21 d of age. Br Poult Sci. 2013;54:407–11.
Tang J, Wen Z, Guo Z, Huang W, Guo Y, Xie M, Hou S. Dietary riboflavin supplementation improve the growth performance and antioxidant status of starter white Pekin ducks fed a corn-soybean meal diets. Livest Sci. 2014;170:131–6.
Britton NL, Riter KL, Smallidge RL, Hillebrandt J. Reversed-phase liquid chromatographic determination of riboflavin in feeds. J AOAC Int. 2003;86:197–201.
Petteys BJ, Frank EL. Rapid determination of vitamin B 2 (riboflavin) in plasma by HPLC. Clin Chim Acta. 2011;412:38–43.
Batey DW, Eckhert CD. Identification of FAD, FMN, and riboflavin in the retina by microextraction and high-performance liquid chromatography. Anal Biochem. 1990;188:164–7.
Folch J, Lees M, Sloane-Stanley G. The method of lipid estimation. J Biol Chem. 1957;228:497.
Wen Z, Tang J, Hou S, Guo Y, Huang W, Xie M. Choline requirements of White Pekin ducks from hatch to 21 days of age. Poult Sci. 2014;93:1–6.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25:25–9.
Gianazza E, Vergani L, Wait R, Brizio C, Brambilla D, Begum S, Giancaspero TA, Conserva F, Eberini I, Bufano D. Coordinated and reversible reduction of enzymes involved in terminal oxidative metabolism in skeletal muscle mitochondria from a riboflavin-responsive, multiple acyl-CoA dehydrogenase deficiency patient. Electrophoresis. 2006;27:1182–98.
Tang J, Hegeman MA, Hu J, Xie M, Shi W, Jiang Y, Vd B, Guo Y, Hou S, Keijer J. Severe riboflavin deficiency induces alterations in the hepatic proteome of starter Pekin ducks. Br J Nutr. 2017;118:641–50.
Romanoff AL, Bauernfeind JC. Influence of riboflavin-deficiency in eggs on embryonic development (gallus domesticus). Anat Rec. 1942;82:11–23.
Carothers DJ, Pons G, Patel MS. Dihydrolipoamide dehydrogenase: functional similarities and divergent evolution of the pyridine nucleotide-disulfide oxidoreductases. Arch Biochem Biophys. 1989;268:409–25.
Sugden MC, Holness MJ. Recent advances in mechanisms regulating glucose oxidation at the level of the pyruvate dehydrogenase complex by PDKs. Am J Physiol Endocrinol Metab. 2003;284:E855–62.
Holness M, Sugden M. Regulation of pyruvate dehydrogenase complex activity by reversible phosphorylation. Biochem Soc Trans. 2003;31:1143–51.
Gibson GE, Blass JP, Beal MF, Bunik V. The α-ketoglutarate-dehydrogenase complex. Mol Neurobiol. 2005;31:43–63.
Zündorf G, Kahlert S, Bunik V, Reiser G. α-Ketoglutarate dehydrogenase contributes to production of reactive oxygen species in glutamate-stimulated hippocampal neurons in situ. Neuroscience. 2009;158:610–6.
Rustin P, Munnich A, Rötig A. Succinate dehydrogenase and human diseases: new insights into a well-known enzyme. Eur J Hum Genet. 2002;10:289–91.
Carrozzo R, Verrigni D, Rasmussen M, de Coo R, Amartino H, Bianchi M, Buhas D, Mesli S, Naess K, Born AP, et al. Succinate-CoA ligase deficiency due to mutations in SUCLA2 and SUCLG1: phenotype and genotype correlations in 71 patients. J Inherit Metab Dis. 2016;39:243–52.
King A, Selak MA, Gottlieb E. Succinate dehydrogenase and fumarate hydratase: linking mitochondrial dysfunction and cancer. Oncogene. 2006;25:4675–82.
Huh TL, Kim YO, Oh IU, Song BJ, Inazawa J. Assignment of the human mitochondrial NAD+ −specific isocitrate dehydrogenase α subunit (IDH3A) gene to 15q25.1→q25.2 by in situ hybridization. Genomics. 1996;32:295–6.
Lee K, Kerner J, Hoppel CL. Mitochondrial carnitine palmitoyltransferase 1a (CPT1a) is part of an outer membrane fatty acid transfer complex. J Biol Chem. 2011;286:25655–62.
van der Leij FR, Cox KB, Jackson VN, Huijkman NC, Bartelds B, Kuipers JR, Dijkhuizen T, Terpstra P, Wood PA, Zammit VA, Price NT. Structural and functional genomics of the CPT1B gene for muscle-type carnitine palmitoyltransferase I in mammals. J Biol Chem. 2002;277:26994–7005.
Bonnefont J, Djouadi F, Prip-Buus C, Gobin S, Munnich A, Bastin J. Carnitine palmitoyltransferases 1 and 2: biochemical, molecular and medical aspects. Mol Asp Med. 2004;25:495–520.
Godárová A, Litzlbauer E, Brunner S, Agu AC, Lohninger A, Hofbauer R. L -carnitine regulates mRNA expression levels of the carnitine acyltransferases – CPT1A, CPT2, and CRAT. Monatsh Chem. 2005;136:1349–63.
Tachibana K, Takeuchi K, Inada H, Yamasaki D, Ishimoto K, Tanaka T, Hamakubo T, Sakai J, Kodama T, Doi T. Regulation of the human SLC25A20 expression by peroxisome proliferator-activated receptor alpha in human hepatoblastoma cells. Biochem Bioph Res Co. 2009;389:501–5.
Ghisla S, Thorpe C. Acyl-CoA dehydrogenases. Eur J Biochem. 2004;271:494–508.
Olsen RK, Andresen BS, Christensen E, Bross P, Skovby F, Gregersen N. Clear relationship between ETF/ETFDH genotype and phenotype in patients with multiple acyl-CoA dehydrogenation deficiency. Hum Mutat. 2003;22:12–23.
Watkins PA, Maiguel D, Jia Z, Pevsner J. Evidence for 26 distinct acyl-coenzyme a synthetase genes in the human genome. J Lipid Res. 2007;48:2736–50.
Gurvitz A, Mursula AM, Firzinger A, Hamilton B, Kilpeläinen SH, Hartig A, Ruis H, Hiltunen JK, Rottensteiner H. Peroxisomal Δ3-cis-Δ2-trans-enoyl-CoA isomerase encoded by ECI1 is required for growth of the yeast Saccharomyces cerevisiae on unsaturated fatty acids. J Biol Chem. 1998;273:31366–74.
Helander HM, Koivuranta KT, Horelli-Kuitunen N, Palvimo JJ, Palotie A, Hiltunen JK. Molecular cloning and characterization of the human mitochondrial 2,4-dienoyl-CoA reductase gene (DECR). Genomics. 1997;46:112–9.
Olpin S, Bates C. Lipid metabolism in riboflavin-deficient rats: 2. Mitochondrial fatty acid oxidation and the microsomal desaturation pathway. Br J Nutr. 1982;47:589–96.
Sakurai T, Miyazawa S, Furuta S, Hashimoto T. Riboflavin deficiency and β-oxidation systems in rat liver. Lipids. 1982;17:598–604.
Brijlal S, Lakshmi A, Bamji M. Mitochondrial oxidative metabolism during respiratory infection in riboflavin deficient mice. J Nutr Biochem. 1999;10:728–32.
Lazarou M, McKenzie M, Ohtake A, Thorburn DR, Ryan MT. Analysis of the assembly profiles for mitochondrial- and nuclear-DNA-encoded subunits into complex I. Mol Cell Biol. 2007;27:4228–37.
Schiff M, Haberberger B, Xia C, Mohsen A-W, Goetzman ES, Wang Y, Uppala R, Zhang Y, Karunanidhi A, Prabhu D, et al. Complex I assembly function and fatty acid oxidation enzyme activity of ACAD9 both contribute to disease severity in ACAD9 deficiency. Hum Mol Genet. 2015;24:3238–47.
Grad LI, Lemire BD. Riboflavin enhances the assembly of mitochondrial cytochrome c oxidase in C. Elegans NADH-ubiquinone oxidoreductase mutants. Biochim Biophys Acta. 2006;1757:115–22.
Parfait B, Chretien D, Rötig A, Marsac C, Munnich A, Rustin P. Compound heterozygous mutations in the flavoprotein gene of the respiratory chain complex II in a patient with Leigh syndrome. Hum Genet. 2000;106:236–43.
Miyake N, Yano S, Sakai C, Hatakeyama H, Matsushima Y, Shiina M, Watanabe Y, Bartley J, Abdenur JE, Wang RY, et al. Mitochondrial complex III deficiency caused by a homozygous UQCRC2 mutation presenting with neonatal-onset recurrent metabolic decompensation. Hum Mutat. 2013;34:446–52.
Rönn T, Poulsen P, Tuomi T, Isomaa B, Groop L, Vaag A, Ling C. Genetic variation in ATP5O is associated with skeletal muscle ATP50 mRNA expression and glucose uptake in young twins. PLoS One. 2009;4:e4793.
Senior AE, Sashi N, Joachim W. The molecular mechanism of ATP synthesis by F1F0-ATP synthase. Biochim Biophys Acta. 2002;1553:188–211.
Xin Z, Pu L, Gao W, Wang Y, Wei J, Shi T, Yao Z, Guo C. Riboflavin deficiency induces a significant change in proteomic profiles in HepG2 cells. Sci Rep-UK. 2017;7:45861.
We want to thank the National Funds for helping support this research. Funders had no role in the study design, analysis, or writing of this article.
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), and the Basic Research Program of the Institute of Animal Sciences of the Chinese Academy of Agricultural Sciences (2018-YWF-YTS-8).
Availability of data and materials
All data generated or analysed during the current study are available from the corresponding author on reasonable request. The datasets supporting the conclusions of this article are included within the article.
Ethics approval and consent to participate
All animal experiments were reviewed and approved by the Institutional Animal Care and Use Committee of the Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, and performed in accordance with the “Guidelines for Experimental Animals” of the Ministry of Science and Technology (Beijing, China).
Consent for publication
All authors approve the submission of this manuscript and declare no conflict of interest. The manuscript has not been published previously, and not under consideration for publication elsewhere.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Egg weight of ducks in the riboflavin-deficient (RD) group and the control (CON) group. (DOCX 14 kb)
Egg production of ducks in the riboflavin-deficient (RD) group and the control (CON) group. (DOCX 15 kb)
Egg fertility rate of ducks in the riboflavin-deficient (RD) group and the control (CON) group. (DOCX 14 kb)
Differentially expressed proteins in E13 embryonic liver from maternal ducks after 8-week riboflavin depletion. (DOCX 108 kb)