Vitamin B2 enables regulation of fasting glucose availability
Abstract
Flavin adenine dinucleotide (FAD) interacts with flavoproteins to mediate oxidation-reduction reactions required for cellular energy demands. Not surprisingly, mutations that alter FAD binding to flavoproteins cause rare inborn errors of metabolism (IEMs) that disrupt liver function and render fasting intolerance, hepatic steatosis, and lipodystrophy. In our study, depleting FAD pools in mice with a vitamin B2-deficient diet (B2D) caused phenotypes associated with organic acidemias and other IEMs, including reduced body weight, hypoglycemia, and fatty liver disease. Integrated discovery approaches revealed B2D tempered fasting activation of target genes for the nuclear receptor PPARα, including those required for gluconeogenesis. We also found PPARα knockdown in the liver recapitulated B2D effects on glucose excursion and fatty liver disease in mice. Finally, treatment with the PPARα agonist fenofibrate activated the integrated stress response and refilled amino acid substrates to rescue fasting glucose availability and overcome B2D phenotypes. These findings identify metabolic responses to FAD availability and nominate strategies for the management of organic acidemias and other rare IEMs.
Editor's evaluation
This paper provides important findings on the metabolic roles of vitamin B2 in the liver that will be of broad interest. Convincing data establish the effects of vitamin B2 deficiency on body composition, energy expenditure, and glucose metabolism.
https://doi.org/10.7554/eLife.84077.sa0Introduction
Flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD) serve as essential cofactors for diverse proteins that mediate oxidation-reduction (redox) reactions, transcriptional regulation, and metabolism (Powers, 2003). In particular, FAD supports the activity of flavoproteins that enable the electron transport chain (ETC), the tricarboxylic acid (TCA) cycle, and fatty acid oxidation (FAO). Along these lines, mutations occurring in more than 50% of human flavoproteins cause inborn errors of metabolism (IEMs) with heterogeneous clinical presentations frequently characterized by organic acidemia, fasting intolerance, and fatty liver disease (Balasubramaniam et al., 2019). Compared to the more well-understood roles of nicotinamide adenine dinucleotide (NAD), the physiological relevance of FAD has remained largely ignored. This gap in knowledge slows the pursuit of therapeutic strategies to treat IEM and leaves fundamental energy balance roles for FAD unresolved.
The liver coordinates whole-body metabolism during fasting by releasing glucose and other fuels to spare the brain and survive low-nutrient conditions. Nuclear receptors sense and receive the signals of nutrient abundance in the liver to perform precise regulation of genes that ultimately maintain the energetic needs of the FAO and amino acid catabolism pathways serving gluconeogenesis (Scholtes and Giguère, 2022). Among nuclear receptors, peroxisome proliferator-activated receptor a (PPARα) modulates fasting responses and pivots hepatocytes toward conservation and recycled substrates to sustain energy production. NAD can activate nuclear receptor co-activators in the liver for gluconeogenesis (Rodgers et al., 2005). However, the role of FAD as a bridging molecule for metabolic reactions and transcriptional regulation of energy balance are almost completely unknown. Synthetic PPARα agonists promote FAO by directing the activity of pathways for balanced lipid metabolism in the liver, which consequently supported a series of FDA-approved fibrate drugs for the treatment of hypertriglyceridemia (Jackevicius et al., 2011). The lipid-lowering properties of bezafibrate and fenofibrate motivated preclinical studies (Steele et al., 2020; Waskowicz et al., 2019; Yavarow et al., 2020) that form ongoing efforts to overcome the limited armament of therapies for IEMs.
Nutrition impacts FAD availability for the organism and dietary riboflavin (vitamin B2) supplies the backbone for all FAD and FMN synthesis. Riboflavin deficiency gives rise to abnormal development and energy balance disorders. Low riboflavin also occurs in more broadly observed diseases including cancer (Kabat et al., 2008; Powers, 2005), HIV (Beach et al., 1992; Fouty et al., 1998), and cardiovascular disease (Tzoulaki et al., 2012). Thus, models that expand how FAD requirements form the regulatory environment for metabolic homeostasis provide opportunities for preclinical experiments and studies of fundamental nutrient-sensing and disease mechanisms. Here, we define the outcomes of vitamin B2 depletion resembling IEMs of flavoprotein deficiency and transcriptional pathways that preserve glucose availability in fasted mice.
Results
Riboflavin deficiency alters body composition and energy expenditure
In mammals, diet furnishes vitamin B2 to synthesize all the FAD for electron transfer in the mitochondria and redox reactions required for cellular homeostasis (Powers, 2003). Among key metabolic organs, ad libitum FAD levels were highest in the liver, heart, and kidney (Figure 1a). To study how FAD depletion influences energy balance, we exposed male mice to vitamin B2-deficient (B2D) or control diets for 4 wk and performed metabolic phenotyping (Figure 1b). We found that 99% B2D was sufficient to reduce liver FAD levels by 70% (Figure 1c). Moreover, B2D significantly blunted weight gain (Figure 1d and e), and body composition measurements showed B2D also reduced fat mass (Figure 1f). When we examined the contribution of B2 to energy expenditure, we were surprised the stunted body weight phenotype of B2D did not arise from higher oxygen consumption or canonical markers of brown adipose tissue thermogenesis (Figure 1—figure supplement 1). B2D strongly reduced oxygen consumption (Figure 1g), but substrate preference was unchanged (Figure 1h). Animals moved less during nighttime measurements (Figure 1i). In absolute terms, B2D-fed mice consumed the same amount of food as controls (Figure 1j). Female mice exposed to B2D also showed reduced body weight gain and locomotory activity with trends toward lower oxygen consumption in the night phase (Figure 1—figure supplement 2). In contrast to B2D, reducing riboflavin in the diet by 90% did not impact liver FAD levels nor influence energy balance (Figure 1—figure supplement 3). These results identify B2 requirements that make FAD available for energy expenditure requirements and body weight in male and female mice.
FAD is required for hepatic glucose production during fasting
Common clinical phenotypes of flavoprotein and nutritional riboflavin depletion disorders include fasting intolerance and hypoglycemia derived from impaired liver glucose production and metabolic flexibility (Houten et al., 2016). We found that liver FAD displayed circadian accumulation (Patel et al., 2012) coinciding with the onset of gluconeogenesis (ZT12) that occurs before the active period for mice (Figure 2a). Follow-up experiments also demonstrated re-feeding mice after 16 hr of fasting raised blood glucose levels (Figure 2b), but depleted liver FAD concentrations (Figure 2c). We next sought to understand whether the changes in FAD that occur in the mouse liver during B2D affected gluconeogenesis using multiple in vivo methods. After 4 wk of B2D, we subjected mice to a pyruvate tolerance test (PTT) after an overnight fast. Consistent with impaired gluconeogenesis from pyruvate, blood glucose concentrations were lower in B2D conditions compared with control diets at all times during the PTT (Figure 2d). Further plasma analysis established trends toward higher fasting concentrations of lactate, triglycerides (TG), free fatty acids (FFAs), and ketone bodies (Figure 2—source data 2). To further determine the effects of B2D on gluconeogenesis in vivo, we measured basal liver glucose production using an intravenous infusion of radiolabeled glucose in awake mice (Chopra et al., 2008; Saha et al., 2004). In a liver-specific way, B2D suppressed in vivo hepatic basal glucose production inferred from 3H-glucose infusions into fasted mice (Figure 2e). Steady-state basal glucose production rates were reached within 1 hr of infusion (Figure 2e’) and B2D reduced average rates by more than 50% (Figure 2e”). These data indicate that FAD depletion directly affects liver glucose metabolism.
-
Figure 2—source data 1
- https://cdn.elifesciences.org/articles/84077/elife-84077-fig2-data1-v1.zip
-
Figure 2—source data 2
- https://cdn.elifesciences.org/articles/84077/elife-84077-fig2-data2-v1.docx
-
Figure 2—source data 3
- https://cdn.elifesciences.org/articles/84077/elife-84077-fig2-data3-v1.docx
To explore molecular outcomes of riboflavin depletion, we used RNA-seq to identify biologically cohesive gene programs of B2D in the liver. Consistent with known roles of FAD in FAO, pathways related to fatty acid catabolism and lipid oxidation were strongly repressed in response to B2D (Figure 2f). Among the B2D-repressed genes annotated to the GO fatty acid catabolism pathway, we noted several encoded flavoproteins whose loss-of-function cause rare organic acidemias, including Sqor (Friederich et al., 2020) and Ivd (Vockley and Ensenauer, 2006). To further investigate the effect of B2D on flavoprotein gene expression, we curated a set of 117 genes encoding flavoproteins that require FAD or FMN for activity. Reflecting a specific impact of riboflavin deficiency, flavoprotein genes were enriched among B2D-repressed genes (Figure 2g) but not B2D-induced genes.
IEMs that arise from mutations in genes encoding mitochondria FAD transfer enzymes for FAO are identified by elevated organic acids in the blood, lipodystrophy, and fatty liver disease (Balasubramaniam et al., 2019). Lipidomics analysis in the liver identified phospholipid (phosphatidylethanolamine [PE], phosphatidylcholine [PC], and lyso-PC) and diacylglycerol (DG) as robustly attenuated (Figure 2h) in B2D-exposed mice compared to controls. Hepatic steatosis (Figure 2i) also accompanied B2D effects, presumably due to increased TG (Figure 2—source data 3). Accordingly, we hypothesized that altered expression of genes in the liver of B2D relative to normal diet controls reflected non-alcoholic fatty liver disease (NAFLD). To test this hypothesis, we performed high-confidence transcriptional target (HCT) intersection analysis (Ochsner et al., 2019) to identify signaling nodes with significant regulatory footprints among liver B2D-induced or -repressed genes. We validated the HCT footprint of B2D against mouse genes from the International Mouse Phenotyping Consortium, whose knockout increased liver TG (Figure 2j). From this approach, we identified strong enrichment of genes induced by B2D with metabolic transcription factor knockouts that cause fatty liver disease, such as Ppara (Cotter et al., 2014; Kersten et al., 1999; Montagner et al., 2016), Nr1h4 (Sinal et al., 2000), and Nr0b2 (Huang et al., 2007). Using other computational approaches to expand upon the underpinnings of macrosteatosis caused by B2D, we retrieved nodes previously shown to contain significant transcriptional footprints within genes differentially expressed in clinical NAFLD (Bissig-Choisat et al., 2021). B2D-induced genes consisted of footprints for transcription factor nodes active in NAFLD (Figure 2k), including HNF family members (Xu et al., 2021) and SREBP1 (Shimano et al., 1997). Collectively, our unbiased approach converged metabolic phenotypes and regulatory networks of B2D with those driving NAFLD and macrosteatosis observed in organic acidemias.
PPARα activity maintains liver FAD pools for fasting glucose availability
Analysis of RNA-seq and ChIP-seq data (Lee et al., 2014; Oshida et al., 2015) discovered strongly enriched PPARα binding near promoter regions for 43 out of 121 putative flavoproteins (p=1.19 × 10–11, hypergeometric test) and supported direct coupling of FAD availability and nuclear receptor activity in the liver. Using DNA motif analysis, we also found canonical PPARα target genes among the genes downregulated by B2D (Figure 3a). Likewise, we identified a set of BioGRID-curated interaction partners of PPARα that was enriched among nodes with strong footprints in B2D-repressed genes (Figure 3a). Furthermore, B2D reduced PPARα-regulated flavoprotein genes (Figure 3b). The observation that patterns of B2D-sensitive gene expression in the liver overlap with PPARα regulatory footprints indicated a convergent role for PPARα and riboflavin in the transcriptional control of gluconeogenic responses.
To explore the physiological intersections between PPARα and B2D, we deleted Ppara in the liver by administering adeno-associated virus (AAV) expressing either Cre recombinase or a GFP under a liver-specific promoter to mice (AAV-TBG-Cre or AAV-TBG-GFP, respectively) harboring a floxed Ppara allele (Pparaflox/flox). Cre-mediated excision of the loxP sites depletes Ppara in adult mice (Ppara ΔHep) and allows us to determine how PPARα mediates B2D responses in mice (Figure 3c). During 4 wk of diet intervention, Ppara ΔHep showed modest body weight gain that was slightly greater than B2D in magnitude. Combined B2D and knockdown of Ppara completely abolished body weight gain (Figure 3d and e) and reduced fat mass (Figure 3f) relative to control diet conditions. We also learned Ppara ΔHep largely recapitulated B2D effects on blood glucose excursion curves during PTTs (Figure 3g and h). Consistent with previous studies (Stec et al., 2019), Ppara ΔHep showed greater fat accumulation in the liver relative to Pparaflox/flox on control diet. B2D strongly increased liver lipid across genotypes (Figure 3i), demonstrating parallel effects of B2D and Ppara ΔHep on liver energy balance. Additionally, FAD levels were lowered to similar extents by Ppara ΔHep or B2D relative to Pparaflox/flox on control diet (Figure 3j). We also phenotyped male Ppara whole-body knockout (pKO) mice exposed to control or B2D for 1 mo (Figure 3—figure supplement 1). As expected (Cotter et al., 2014), we learned pKO largely negated effects of B2D on body weight gain. PTT demonstrated B2D and pKO decreased the conversion of glucose release relative to control diets and wild-type controls.
To extend our studies of gene regulation, we profiled Ppara and the PPARα target gene Acox1. Ppara levels were reduced in Ppara ΔHep, along with Acox1 (Figure 3k). Furthermore, the blunted expression levels of Acox1 to diet and genetic interventions (Ppara ΔHep) confirm B2D confers diminished PPARα activity in the mouse liver. In a complementary approach, we performed experiments to determine how FAD depletion affects PPARα transcriptional activity using the riboflavin analog and competitive inhibitor lumiflavin (Figure 3l). We found lumiflavin reduced PPRE-luciferase activity in the presence of the PPARA agonist WY-14643, suggesting intact FAD pools contribute to PPARα transcriptional activity. Altogether, these findings suggest that PPARα sustains FAD levels and requires riboflavin to direct glucogenic responses in the liver.
PPARα activation by fenofibrate rescues liver glucose production even when FAD cannot be generated from diet
The FDA-approved fibrate drugs act selectively on PPARα to lower blood lipids and treat hypertriglyceridemia (Bougarne et al., 2018). PPARα agonists also show promise for the treatment of some mitochondrial disorders (Steele et al., 2020). The convergence of B2D effects on PPARα regulation of gene expression and metabolic responses suggested fenofibrate treatment may restore some metabolic competency in animals on riboflavin-deficient diet. To explore this possibility, we administered fenofibrate after 7 wk of B2D exposure (Figure 4a). Daily gavage with fenofibrate (300 mg/kg) for 2 wk while maintaining mice on diet interventions significantly reduced body weight gain under B2D conditions (Figure 4b and c). At the end of the experiment, fasted mice received fenofibrate 2 hr before blood glucose measurements. Fenofibrate increased blood glucose in both groups of mice far above pre-gavage levels (Figure 4d). When liver histology was examined, we did not detect hepatic steatosis in B2D mice that received fenofibrate (Figure 4e).
Given the ability of PPARα to regulate flavoprotein gene expression, we pursued additional RNA-Seq studies to understand the mechanisms that allowed fenofibrate to rescue hypoglycemia in B2D. Consistent with restoration of the flavoprotein transcriptome in response to PPARα activation, we found the number of flavoprotein genes induced by B2D+ fenofibrate more than doubled when compared to B2D alone (Figure 4f). Given the improvement in hepatic steatosis after PPARα activation, we hypothesized the B2D+ fenofibrate treatment caused inversion of the alignments between B2D and NAFLD transcription networks. Using this approach, we found gene footprints enriched by B2D+ fenofibrate and those depleted in clinical NAFLD converged (Figure 4g), including the GABP transcriptional program inactivated in inflammatory liver diseases (Niopek et al., 2017). These unbiased approaches strengthen the notion that PPARα activation overcomes fatty liver and hypoglycemia phenotypes imposed by B2D.
Altered sphingolipid pools and respiratory chain efficiency in B2D
Fat and protein metabolism produces substrates for the synthesis of sphingolipids, such as ceramides and dihydroceramides, whose tissue accrual associates with severity of fatty liver disease (Luukkonen et al., 2016) and mitochondrial dysfunction (Hammerschmidt et al., 2019; Park et al., 2016). To determine whether B2D or B2D+ fenofibrate also changed the sphingolipid composition of the mouse liver, we measured a battery of sphingolipids and detected varying levels of multiple species altered by diet and fenofibrate treatments (Figure 5—figure supplement 1). Pool sizes for deoxysphingolipids that require alanine condensation with palmitoyl-CoA trended higher with B2D (Figure 5a). This finding may derive from the serine availability during hypoxic stress and higher NADH levels (Yang et al., 2020). Consistent with this idea, B2D conditions reduced the NAD/NADH ratio in the liver (Figure 5b). Notably, deoxysphingolipids impair mitochondrial function (Alecu et al., 2017; Muthusamy et al., 2020) and lead to an energy deficit, especially in the liver, where fasting requires elevated demand for biomass.
Pathogenic variants in genes for riboflavin transport and metabolism that deplete FAD impair electron transfer from flavoenzymes in the ETC to coenzyme Q10 and, ultimately, the energy requirements for gluconeogenesis and fasting tolerance (Rinaldo et al., 2002). During prolonged fasting, beta-oxidation and proteolysis produce high fluxes of electrons that flow through the ETC. Coenzyme Q10 collects and converges electron flow on complex III via complex II (oxidizing succinate into fumarate) or complex I. In this way, the coenzyme Q10 pool accommodates variable electron fluxes and manages the mitochondria redox environment. Analysis of FAD and FMN (Figure 5c) confirmed fenofibrate was incapable of reconciling cofactor pools and lacked meaningful impacts on relative levels of some oxidative phosphorylation proteins (Figure 5d). In contrast to the reduced NAD/NADH, we found B2D treatments accumulated oxidized coenzyme Q10 and rodent-biased coenzyme Q9 (Figure 5e), suggesting higher Q/QH2, more negative free energy for complexes I/II, and compensation for the defects in mitochondrial energy efficiency associated with B2D (Satapati et al., 2015).
Fenofibrate activates the integrated stress response in B2D
We next measured concentrations of carnitines, amino acid, and hydrophilic metabolites in the liver using mass spectroscopy across the B2D experiments. B2D altered the steady-state levels of valine, methionine, phenylalanine, and caused accumulation of short-chain C5 carnitines that reflect incomplete beta-oxidation of fatty acids (Figure 6a). Upper arms of glucose metabolism (3PG/2PG) showed lower activities coupled with higher lactate and buildup of metabolites above pyruvate oxidation (G3P, PEP). Metaboanalyst (Pang et al., 2021) revealed complete removal of dietary riboflavin caused metabolite changes in the liver that enriched for organic acidemias and inborn errors of the TCA cycle (Figure 6b). Moreover, B2D caused accumulation of oxidized TCA cycle metabolic intermediates, including fumarate, malate, and 2HG. These findings were compatible with gene sets (Figure 6c) associated with hypoxia and epithelial-mesenchymal transition (Sciacovelli et al., 2016; Ward et al., 2010).
PPARα activity favors conversion of proteins to provide amino acids as substrates for anabolic processes (Kersten et al., 2001). Our phenotyping analysis of fenofibrate suggested alternative carbon sources might supply the substrates to support glucose production during B2D (Figure 4c), including pyruvate, which recovered blood glucose in B2D to pre-gavage control levels (Figure 6d). Amino acids, such as alanine and serine, are also significant contributors to de novo synthesis of glucose. In line with this idea, gluconeogenic amino acids (serine, glutamate, histidine, alanine) showed selective and unique accumulation in the combined B2D and fenofibrate treatments (Figure 6a). Steady-state levels of other anaplerotic amino acids that replenish the TCA cycle (Figure 6a), isoleucine, and threonine, also increased at the expense of reduced ketosis and β-hydroxybutyrate (Figure 2—source data 3). Similarly, B2D+ fenofibrate elevated levels of TCA cycle metabolites citrate and aconitate and increased the fumarate to succinate ratio. Importantly, we noted moderated levels of carnitines enriched in organic acidemias (C5 and C6). These conditions suggest PPARα activation inhibits cataplerosis of TCA cycle intermediates.
The constellation of phenotypic effects resulting from riboflavin depletion and fenofibrate suggested global shifts in gene expression and metabolism. RNA-seq indicated B2D+ fenofibrate precipitated a response with elements of greater amino acid and glucose metabolism (Figure 6e). We also observed a gene signature for the integrated stress response (ISR) further supported by increased expression of the master transcription factor regulator Atf4, as well as its key target genes (Psat1, Gls, Fgf21, Gpt2, Ddit4) during combined B2D+ fenofibrate treatments. Mechanistically, ISR activation occurs through phosphorylation of eIF2α and ATF4 translation, which mediates transcription of target genes to resolve the ISR and regain amino acid homeostasis (Harding et al., 2003; Ye et al., 2012). In line with this idea, B2D+ fenofibrate achieved robust GCN2 expression, eIF2α phosphorylation, and higher ATF4 levels relative to any other treatment (Figure 6f). Together, our integrated transcription and mass spectrometry analyses favor a model in which PPARα activation precipitates concerted activation of the ISR, which in turn overcomes the fasting intolerance imposed by B2D (Figure 6g).
Discussion
Results presented here reveal that vitamin B2 provides critical substrates for glucose availability during fasting in mice. Mammals cannot synthesize vitamin B2, so diet remains the only available source of FAD and FMN (Powers, 2003). Our results show that 4 wk of 90% B2D did not influence whole-body metabolic phenotypes. Based on data reported by Fenton and Cowgill, 1947, mice require about 4 mg riboflavin/kg diet for normal growth. However, mouse chow generally contains almost twice the daily riboflavin requirement (Anonymous, 1977). Most tissues store very little of the vitamin as riboflavin. The majority of bioavailable riboflavin is converted to FAD and excess is excreted in urine as FMN (Rivlin, 1970). Riboflavin turnover rates in rodents are also relatively slow at approximately 16 d (Yang and McCormick, 1967), which underscores sufficient FAD concentrations remain in the liver after 4 wk of 90% B2D. Therefore, the historical literature suggests 10% of the required level of riboflavin may be sufficient for growth and metabolism in the mouse during our 4-week monitoring period.
Our work sheds light on how the liver copes with severe metabolic crises and the FAO disorders caused by flavoprotein disruption or FAD depletion during prolonged fasting when both beta-oxidation and gluconeogenesis are concomitantly activated. We speculate these FAO disorders rely on conservation responses when the mitochondria are starved of FAD and FMN required for ETC activity. B2D reduces fat stores and slows metabolism for the energy conservation functions that adapt the animal to survive FAD and FMN depletion. Despite lower body weight gain when compared to control diet, B2D causes reduced energy expenditure. IEMs that arise from mutations in genes encoding mitochondrial FAD transfer enzymes for FAO are often identified by ectopic fat distribution and liver disease (Balasubramaniam et al., 2019). People with long-chain fatty acid oxidation disorders (LC-FAODs), including mutations in the FAD-dependent protein ACADVL, also show similar energy expenditure phenotypes as B2D. Indeed, energy expenditure is significantly lower among subjects with LC-FAODs compared to the reference population in the National Health and Nutrition Examination Survey (DeLany et al., 2023). These observations suggest that responses to B2D are relevant to studying IEMs defined by mutations in FAD-dependent metabolic enzymes.
Interestingly, the metabolic phenotype of B2D resists conditions unfavorable for mitochondrial function, including hypoxia, suggesting that these changes select for survival. Increased reliance upon glycolysis occurs in mitochondrial disease (Robinson, 2006) and hypoxia responses activate glycolytic enzymes that allow energy production when the mitochondria are starved of oxygen as a substrate for oxidative phosphorylation. Likewise, our data support the idea that the stress of B2D gives rise to an environment for fenofibrate and PPARα activation to co-opt the ISR, adapt to TCA cycle dysfunction, and engage hypoxia enzymes to reconcile anaplerosis. Elevated TCA cycle intermediates, such as fumarate, then act to stabilize antioxidant transcription factors and protect against oxidative stress in the liver (Ashrafian et al., 2012) during the stress of IEMs modeled in our study.
IEMs modeled in our study frequently present NAFLD-like phenotypes that contribute to fasting intolerance (Rinaldo et al., 2002). In our model, B2D alters lipid profiles and gene expression patterns that converge B2D with more common NAFLD phenotypes. The liver dominates mass-specific metabolic rates (Rolfe and Brown, 1997) and, for this reason, the NAFLD caused by B2D likely reflects a combination of lower hepatic fatty acid oxidation contributing to the accumulation of liver triglycerides and other complex lipid species observed in obesity (Shi et al., 2022). Our lipidomics also revealed that B2D caused accumulation of deoxysphingolipids, which also become more abundant in NAFLD from incomplete fat oxidation and accrual of toxic intermediates (Gai et al., 2020). These findings are particularly relevant for discovering new biomarkers for fatty liver disease.
The energetic requirements of fasting dictate substrate oxidation and electron transport. B2D restricts mitochondrial function, which causes reductive pressures that likely exacerbate ROS formation in part through proton leak. PPARα activation prevents oxidative stress (Ip et al., 2004) by merging lower NAD/NADH and oxidized Q, which, in turn, lowers the free energy of electron transport through complexes I and II. These data are consistent with inhibition of complex II, resulting in a more oxidized Q pool and reveals an important adaptation that allows the liver to overcome riboflavin deficiency (Treberg et al., 2011).
Inefficient amino acid availability, or increased requirements of amino acids to maintain gluconeogenesis, activate ATF4 and the ISR. ATF4 is the principal downstream effector of the ISR, whose regulation becomes altered in human and rodent NAFLD (Puri et al., 2008; Seo et al., 2009). Integrated metabolomic and RNA-seq studies demonstrated that fenofibrate unveils the ISR to increase the abundance of anaplerotic amino acids. We are unaware of other studies that observe coincident activation of PPARα and ATF4 to drive adaptive responses to liver stress. However, induction of shared PPARα and ATF4 targets, including Fgf21, may contribute to hepatoprotection from lipotoxic lipid accumulation (Montagner et al., 2016). It will be interesting in the future to determine why B2D exposes a vulnerability to fenofibrate that engages selective genome regulation by PPARα and ATF4 for recovering glucose production in settings of fasting intolerance.
This study is unique in its comprehensive approach to understand the complex metabolic consequences of FAD depletion and riboflavin auxotrophy in mouse models. Nevertheless, we acknowledge some limitations of the study. While the effect of B2D may act through different pathways and tissues, we demonstrated some of these effects may be mediated through disturbance of nuclear receptor activity and altered glucose availability in the liver. A limitation of our study is that we have focused on how PPARα requires vitamin B2 in the mouse liver to produce glucose, and it is possible that other fasting responsive transcription factors (Bideyan et al., 2021) are sensitive to B2D. One important caveat of these experiments is that the kidney also contributes to glucose production during fasting (Joseph et al., 2000). Although our in vivo experiments do not explore whether the kidney reconciled any liver gluconeogenic deficiency, our observations reveal fundamental vitamin requirements to source the liver with the FAD and FMN pools necessary for energy balance.
There is evidence fibrates generate beneficial effects in IEMs, including improvement of antioxidant defenses and resolution of inflammation (Seminotti et al., 2023). While further studies are needed, we describe allostatic outcomes of PPARα activation that overcome bioenergetic costs of fasting. Rare diseases of flavoprotein mutations, including organic acidemias, cause substantial morbidity and have no cure. Therefore, understanding how nuclear receptor regulation of flavoprotein function and FAD pool distribution surmounts hypoglycemia and fatty liver is valuable for describing fundamental physiology and implementing future therapeutic strategies.
Materials and methods
Mice and housing conditions
Request a detailed protocolAll animal procedures were approved by the Institutional Animal Care and Use Committee of Baylor College of Medicine (AN-6411). All mice were housed in a barrier-specific pathogen-free animal facility with 12 hr dark-light cycle and free access to water and food. C57BL/6J wild-type mice (RRID:IMSR_JAX:000664) were obtained from the BCM Center for Comparative Medicine and global Ppara-/- mice (RRID:IMSR_JAX:008154) were generated previously (Lee et al., 1995). In all experiments, male mice were randomly placed on control or riboflavin deficient diet starting at 4 wk of age. Riboflavin-deficient and matched control diets (Xu et al., 2018) were provided by Research Diets: 90% control, D10012G; 90% riboflavin-deficient, D12030102; 99% control, A18041301; or 99% riboflavin-deficient, A19080901. All diets were isocaloric, and amino acids were kept constant (Supplementary files 1 and 2).
AAV TBG administration
Request a detailed protocolFor AAV-TBG injections, 1-month old Pparaflox/flox mice (Brocker et al., 2017) were injected with AAV8-TBG-iCre-WPRE (1 × 1011 vc/ml, Addgene# 192878) or AAV8-TBG-EGFP-WPRE (1 × 1011 vc/ml, Addgene #192880) into the tail vein in a total volume of 100 μl saline. These viruses were provided by the Gene Vector Core at Baylor College of Medicine. Mice were allowed to recover for 2 wk to permit for viral expression and recombination of floxed Ppara alleles before being used for B2D experiments.
Fenofibrate gavage
Request a detailed protocolWe used malnutrition (van Zutphen et al., 2016), Ppara knockout (Montagner et al., 2016), and cancer cachexia (Goncalves et al., 2018) studies to design fenofibrate interventions. 0.5% methylcellulose solution was prepared by heating 150 ml water with 2.5 g 400 cP methylcellulose (Sigma, #M0262) added with stirring. Chilled water was added (350 ml) and stirred overnight at 4°C. Fenofibrate-gavage solution was made with 112.5 mg fenofibrate (Sigma, #F6020) in 1.5 ml 0.5% methylcellulose solution. Mice were gavaged daily at 300 mg/kg.
Pyruvate tolerance test (PTT)
Request a detailed protocolTo determine pyruvate tolerance, mice were fasted for 16 hr, and Na-pyruvate was administered (1 g/kg body weight) by intraperitoneal injection. Blood glucose levels were monitored at 0, 15, 30, 60, and 120 min by a Freestyle Glucose Monitoring System (Abbott Laboratories).
Basal glucose production rate in live mice
Request a detailed protocolBasal glucose production was measured as described previously (Saha et al., 2004). In brief, mice were anesthetized, and a midline neck incision was made to expose the jugular vein. A microcannula was inserted into the jugular vein, threaded into the right atrium, and anchored at the venotomy site. Mice were allowed to recover for 4–5 d with ad libitum access to water and food. Following an overnight fast, the conscious mice received a bolus dose (10 uCi) and then a constant rate intravenous infusion (0.1 uCi/min) of chromatography-purified [3-3H]-glucose using a syringe infusion pump for 120 min through surgically implanted catheter in jugular vein. For determination of basal glucose production, blood samples were collected after 0, 60, 90, and 120 min of labeled glucose infusion to calculate the basal glucose production rates using plasma glucose-specific activity of tritiated labeled glucose (3-3H). Steady states were reached within 1 hr of infusion. Plasma was separated and deproteinized using equal volumes of barium hydroxide and zinc sulfate, dried to remove 3H2O, resuspended in water, and counted in scintillation fluid using a liquid scintillation counter (Beckman Instruments, Palo Alto, CA). We calculated basal whole-body glucose production (mg/kg/min) by dividing the [3-3H] glucose infusion rate by the plasma glucose-specific activity corrected to body weight (Wall et al., 1957).
Indirect calorimetry
Request a detailed protocolWild-type mice were maintained on experimental diets and housed at room temperature in Comprehensive Lab Animal Monitoring System Home Cages (CLAMS-HC, Columbus Instruments). Oxygen consumption, CO2 emission, energy expenditure, food and water intake, and activity were measured for 5 d (BCM Mouse Metabolic and Phenotyping Core). Mouse body composition was examined by magnetic resonance imaging (Echo Medical Systems) before indirect calorimetry.
Histology
Request a detailed protocolSections of liver tissue were frozen in Tissue-Tek OCT compound (4583; Sakura Finetek USA), and neutral lipids stained with Oil Red O. Formalin-fixed paraffin-embedded tissue sections were stained with hematoxylin and eosin (H/E). Images were captured on a Nikon Ci-L Brightfield microscope.
Serum and lipid assays
Request a detailed protocolFasted serum was used to measure serum lactate (K607; Biovision), ALT (TR71121; Thermo Scientific), AST (TR70121; Thermo Scientific), beta-hydroxybutyrate (Biovision K632), and serum-free fatty acids (#sfa-1; Zen-Bio).
Liver FAD measurements
Request a detailed protocol10–15 mg of liver tissue was deproteinated (K808; Biovision), followed by measurement of FAD through a colorimetric assay (K357; Biovision). FAD concentration was standardized to input tissue weight.
Hepatic TG and cholesterol
Request a detailed protocolBoth serum and tissue samples were analyzed for triglycerides (Triglyceride reagent TR22421; Thermo Scientific) and cholesterol (Total Cholesterol Reagent TR13421; Thermo Scientific). Hepatic TGs and cholesterol were assayed as described previously (Kim et al., 2019). Briefly, liver homogenates were mixed with a 1:2 chloroform:methanol solution followed by isolation of the lipid-rich chloroform layer (modified Folch method).
Immunoblotting
Request a detailed protocolCell and tissue lysates were prepared in Protein Extraction Reagent (Thermo Fisher) supplemented with Halt Protease and Phosphatase Inhibitor Cocktail (Thermo Fisher). Western blotting was performed with whole-cell lysates run on 4–12% Bis-Tris NuPage gels (Life Technologies) and transferred onto Immobilon-P Transfer Membranes (Millipore), followed by antibody incubation. Immunoreactive bands were visualized by chemiluminescence. Antibodies used in this study are listed in Supplementary file 3.
RNA extraction and RNA-seq analysis
Request a detailed protocolTotal liver RNA was extracted using the QIAGEN RNeasy Plus Mini kit (74034; QIAGEN). Sample quality was confirmed on an Agilent 2100 Bioanalyzer (Agilent). mRNA library preparation and RNA sequencing were performed by Novogene. mRNA libraries were prepared with NEBNext Ultra RNA Library Prep Kit for Illumina (NEB) and size selection for libraries was performed using AMPure XP system (Beckman Coulter), followed by library purity analysis. Libraries were sequenced on NovaSeq PE150 and reads mapped to the UCSC mouse reference genome mm10 using STAR. FeatureCounts was used to calculate the expression level as reads per kilobase per million (RPKM). DESeq2 calculated differentially expressed genes with p-values adjusted using Benjamini and Hochberg’s method for controlling the false discovery rate (FDR). Genes with significant differential expression were determined by p<0.05. Gene Set Enrichment Analysis was performed with the Molecular Signatures Database, and -log10(p-value) calculated for Hallmark gene sets.
qPCR
Request a detailed protocolcDNA was synthesized from total RNA using qScript (QuantBio #95048-100). Relative mRNA expression was measured with SsoAdvanced Universal Probes Supermix reactions (Bio-Rad #175284) read out with a QuantStudio 3 real-time PCR system (Applied Biosystems). TATA-box binding protein (Tbp) was the invariant control. Roche Universal Probe Gene Expression Assays were used as previously described (Koh et al., 2018). The following primer and probe sets were used to detect the following genes: Ppara: tcgagttcatgcaagtttcg (F), tccctcctggcttctctagg (R), Roche Universal Probe Library probe #1; Acox1: caccattgccattcgataca (F), tgcgtctgaaaatccaaaatc (R), Roche Universal Probe Library probe #106; Ucp1: ggcctctacgactcagtcca (F), taagccggctgagatcttgt (R), Roche Universal Probe Library probe #34; Cidea: aaaccatgaccgaagtagcc (F), aggccagttgtgatgactaagac (R), Roche Universal Probe Library probe #66; Prdm16: acaggcaggctaagaaccag (F), cgtggagaggagtgtcttcag (R), Roche Universal Probe Library probe #56; Pparg2: gaaagacaacggacaaatcacc (F), gggggtgatatgtttgaacttg (R), Roche Universal Probe Library probe #7; Cpt1b: gagtgactggtgggaagaatatg (F), gctgcttgcacatttgtgtt (R), Roche Universal Probe Library probe #92; Tbp: cggtcgcgtcattttctc (F), gggttatcttcacacaccatga (R), and Roche Universal Probe Library probe #51.
Luciferase assays
Request a detailed protocolWe used HepG2 cells purchased directly by our lab for the studies. All cells are tested monthly to ensure an absence of mycoplasma contamination. HepG2 cells (ATCC, HB-8065, RRID:CVCL_002) were maintained in the following media: Eagle’s Minimum Essential Media supplemented with 10% FBS, and 1% penicillin/streptomycin antibiotics. We used Lipofectamine 2000 to transiently express pcDNA3.1-PPARA (Addgene #169019), pCMV β-gal (ATCC) and PPRE-luciferase fusion plasmids (Addgene #1015). 20 hr later, cells were treated with lumiflavin (Cayman) or 10 μM WY14643 (Cayman) for an additional 24 hr. 0.1% DMSO served as the vehicle treatment. Reporter gene activity was detected using the Promega Luciferase Assay Kit. Relative luminescence units were normalized to β-galactosidase activity (Sigma).
Consensome and high-confidence transcriptional (HCT) target intersections
Request a detailed protocolTranscription factor footprint analysis and consensome enrichments were performed (Ochsner et al., 2019). For transcription factors, node and node family consensomes are gene lists ranked according to the strength of their regulatory relationship with upstream signaling pathway nodes derived from independent publicly archived transcriptomic or ChIP-Seq datasets. Genes in the 95th percentile of a given node consensome were designated high-confidence transcriptional targets (HCTs) for that node and used as the input for the HCT intersection analysis using the Bioconductor GeneOverlap analysis package implemented in R. For both consensome and HCT intersection analysis, p-values were adjusted for multiple testing using the method of Benjamini and Hochberg to control the FDR as implemented with the p. adjust function in R, to generate q values. Evidence for a transcriptional regulatory relationship between a node and a gene set was inferred from a larger intersection between the gene set and HCTs for a given node or node family than would be expected by chance after FDR correction (q < 0.05). The HCT intersection analysis code has been deposited in the SPP GitHub account at https://github.com/signaling-pathways-project/ominer/ (copy archived at signaling-pathways-project, 2023).
Metabolomics
Targeted measurements of hepatic carnitines, fatty acids, lipids species, CoAs, glycolysis, and TCA metabolites were carried out by the BCM Dan L Duncan Cancer Center CPRIT Cancer Proteomics and Metabolomics Core. Parallel analysis of lipids and ceramides was performed at the University of Utah. All samples were processed and analyzed as described previously (Chaurasia et al., 2019; Kettner et al., 2016).
Reagents
High-performance liquid chromatography grade and mass spectrometry grade reagents were used: acetonitrile, methanol, and water (Burdick & Jackson); formic acid, ammonium acetate, and internal standards (Sigma-Aldrich); MS grade lipid standards (Avanti Polar Lipids).
Internal standards and quality control
Request a detailed protocolTo assess overall process reproducibility, mouse pooled liver or serum samples were run along with the experimental samples. A number of internal standards, including injection standards, process standards, and alignment standards, were used to assure QA/QC targets to control for experimental variability. Aliquots (200 µl) of 10 mM solutions of isotopically labeled standards were mixed and diluted up to 8000 µl (final concentration 0.25 mM) and aliquoted into a final volume of 20 µl. The aliquots were dried and stored at –80°C until further analysis. To monitor instrument performance, 20 µl of a matrix-free mixture of the internal standards were reconstituted in 100 µl of methanol:water (50:50) and analyzed by multiple reaction monitoring (MRM). Metabolite extraction from the samples was monitored using pooled mouse serum or liver samples and spiked internal standards. The median coefficient of variation (CV) value for the internal standard compounds was 5%. To address overall process variability, metabolomic studies included a set of nine experimental sample technical replicates, which were spaced evenly among the injections for each day. The matrix-free internal standards and serum and liver samples were analyzed twice daily.
Separation of CoAs and carnitines
Request a detailed protocolTargeted profiling for CoAs and carnitines in electrospray ionization-positive mode by the RP chromatographic method employed a gradient containing water (solvent A) and acetonitrile (ACN, solvent B, with both solvents containing 0.1% formic acid). Separation of metabolites was performed on a Zorbax Eclipse XDBC18 column (50 × 4.6 mm i.d.; 1.8 μm, Agilent Technologies) maintained at 37°C. The gradient conditions were 0–6 min in 2% B, 6.5 min in 30% B, 7 min in 90% B, 12 min in 95% B; followed by re-equilibration to the initial conditions.
Separation of glycolysis and TCA metabolites
Request a detailed protocolGlycolysis and TCA metabolites were separated by normal phase chromatography using solvents containing water (solvent A), solvent A modified by the addition of 5 mM ammonium acetate (pH 9.9), and 100% acetonitrile (solvent B). The binary pump flow rate was 0.2 ml/min with a gradient spanning 80% B to 2% B over 20 min, 2% B to 80% B for 5 min, and 80% B for 13 min. The flow rate was gradually increased during the separation from 0.2 ml/min (0–20 min) to 0.3 ml/min (20–25 min), and then 0.35 ml/min (25–30 min), 0.4 ml/min (30–37.99 min), and finally to 0.2 ml/min (5 min). Metabolites were separated on a Luna Amino (NH2) column (4 µm, 100 A 2.1 × 150 mm, Phenominex) maintained in a temperature-controlled chamber (37°C). All the columns used in this study were washed and reconditioned after every 50 injections.
Separation of fatty acids
Request a detailed protocolTargeted profiling for fatty acids employed a reverse phase chromatographic method by a gradient containing water (solvent A) with 10 mM ammonium acetate (pH 8) and 100% methanol (solvent B) on a Luna Phyenyl Hexyl column (3 µm, 2 × 150 mm; Phenominex, CA) maintained at 40°C. The binary pump flow rate was 0.2 ml/min with a gradient spanning 40% B to 50% B over 8 min, 50% B to 67% B over 5 min, hold 67% B for 9 min, 67% B to 100% B over 1 min, hold 100% B for 6 min, 100% B to 40% B over 1 min, and hold 40% B for 7 min.
Liquid chromatography/mass spectrometry (LC/MS) analyses
Request a detailed protocolThe chromatographic separation of non-lipid metabolites was performed using either reverse phase separation or normal phase online with the unbiased profiling platform based on a 1290 SL Rapid resolution LC and a 6490 triple Quadrupole mass spectrometer (Agilent Technologies, Santa Clara, CA). Lipidomics required a Shimadzu CTO-20A Nexera X2 UHPLC coupled with TripleTOF 5600 equipped with a Turbo VTM ion source (AB Sciex, Concord, Canada). Using a dual electrospray ionization source, the samples were independently examined in both positive and negative ionization modes. The data acquisition during the analysis was controlled using the Mass Hunter workstation data acquisition software.
Lipidomics
Request a detailed protocolMouse liver lipids were extracted using a modified Bligh-Dyer method. Briefly, 50 mg of crushed tissue sample from mouse whole liver was used. A 2:2:2 volume ratio of water/methanol/dichloromethane was used for lipid extract at room temperature after spiking internal standards 17:0 LPC, 17:0 PC, 17:0 PE, 17:0 PG, 17:0 ceramide, 17:0 SM, 17:0 PS, 17:0 PA, 17:0 TAG, 17:0 MAG, DAG 16:0/18:1, CE 17:0. The organic layer was collected, followed by a complete drying procedure under nitrogen. Before MS analysis, the dried extract was resuspended in 100 μl of Buffer B (10:5:85 acetonitrile/water/isopropyl alcohol) containing 10 mM NH4OAc and subjected to reverse-phase chromatography and LC/MS. Internal standards prepared in chloroform/methanol/water (100 pmol/μl) were LPC 17:0/0:0, PG 17:0/17:0, PE 17:0/17:0, PC 17:0/17:0, TAG 17:0/17:0/17:0, SM 18:1/17:0, MAG 17:0, DAG 16:0/18:1, CE 17:0, ceramide 18:1/17:0, PA 17:0, PI 17:0/20:4, and PS 17:0/17:0.
For lipid separation, 5 ml of the lipid extract was injected into a 1.8 mm particle 50 × 2.1 mm Acquity HSS UPLC T3 column (Waters). The column heater temperature was set at 55°C. For chromatographic elution, a linear gradient was used over a 20 min total run time, with 60% Solvent A (acetonitrile/water [40:60, v/v] with 10 mM ammonium acetate) and 40% Solvent B (acetonitrile/water/isopropanol [10:5:85 v/v] with 10 mM ammonium acetate) gradient in the first 10 min. The gradient was ramped linearly to 100% Solvent B for 7 min. Then the system was switched back to 60% Solvent B and 40% Solvent A for 3 min. A 0.4 ml/min flow rate was used at an injection volume of 5 μl. The column was equilibrated for 3 min and run at a flow rate of 0.4 ml/min for a total run time of 20 min. TOF MS survey scans (150 ms) and 15 MS/MS scans with a total duty cycle time of 2.4 s were performed. The mass range in both modes was 50–1200 m/z. The acquisition of MS/MS spectra by data-dependent acquisition (DDA) function of the Analyst TF software (AB Sciex).
The raw data in .mgf format were converted using ProteoWizard software. The NIST MS PepSearch Program was used to search the converted files against LipidBlast libraries. The m/z width was determined via the mass accuracy of internal standards at 0.001 for positive mode and 0.005 for a negative mode at an overall mass error of less than 2 ppm. The minimum match factor at 400 was set for the downstream data processing. The MS/MS identification results from all the files were combined using an in-house software tool to create a library for quantification. The raw data files were searched against this in-house library of known lipids with mass and retention time using Multiquant 1.1.0.26 (ABsciex). The lipid species identified in the positive or negative ion modes were analyzed separately using relative abundance of peak spectra for the downstream analyses. The identified lipids were quantified by normalizing against their respective internal standard.
Ceramides and lipids
Request a detailed protocolLipid extracts are separated on an Acquity UPLC CSH C18 1.7 μm 2.1 × 50 mm column maintained at 60°C connected to an Agilent HiP 1290 Sampler, Agilent 1290 Infinity pump, equipped with an Agilent 1290 Flex Cube and Agilent 6490 triple quadrupole (QQQ) mass spectrometer. In positive ion mode, sphingolipids are detected using dynamic multiple reaction monitoring (dMRM). Source gas temperature is set to 210°C, with an N2 flow of 11 l/min and a nebulizer pressure of 30 psi. Sheath gas temperature is 400°C, sheath gas (N2) flow of 12 l/min, capillary voltage is 4000 V, nozzle voltage 500 V, high-pressure RF 190 V, and low-pressure RF is 120 V. Injection volume is 2 μl, and the samples are analyzed in a randomized order with the pooled QC sample injection eight times throughout the sample queue. Mobile phase A consists of ACN:H2O (60:40 v/v) in 10 mM ammonium formate and 0.1% formic acid, and mobile phase B consists of IPA:ACN:H2O (90:9:1 v/v) in 10 mM ammonium formate and 0.1% formic acid. The five chromatography gradient starts at 15% mobile phase B, increases to 30% B over 1 min, increases to 60% B from 1 to 2 min, increases to 80% B from 2 to 10 min, and increases to 99% B from 10 to 10.2 min where it is held until 14 min. Post-time is 5 min, and the flow rate is 0.35 ml/min throughout. Collision energies and cell accelerator voltages were optimized using sphingolipid standards with dMRM transitions as [M+H]+→[m/z = 284.3] for dihydroceramides, [M+H]+→[m/z = 287.3] for isotope-labeled dihydroceramides, [M-H2O+H]+→[m/z = 264.2] for ceramides, [MH2O+H]+→[m/z = 267.2] for isotope-labeled ceramides and [M+H]+→[M-H2O+H] for all targets. Sphingolipids and ceramides without available standards are identified based on HR-LC/MS, quasi-molecular ions, and characteristic product ions. Their retention times are either taken from HR-LC/MS data or inferred from the available standards. Results from LC-MS experiments are collected using Agilent Mass Hunter Workstation and analyzed using the software package Agilent Mass Hunter Quant B.07.00. Ceramide and lipid species are quantitated based on peak area ratios to the standards added to the extracts.
Analysis of FAD and FMN by IC-HRMS
Request a detailed protocolTo determine the relative abundance of FAD and FMN in mouse liver tissue, extracts were prepared and analyzed by high-resolution mass spectrometry (HRMS). Approximately 20–30 mg of tissue were pulverized in liquid nitrogen then homogenized with a Precellys Tissue Homogenizer. Metabolites were extracted using 80/20 (v/v) methanol/water with 0.1% ammonium hydroxide. Samples were centrifuged at 17,000 × g for 5 min at 4°C, supernatants were transferred to clean tubes, followed by evaporation under vacuum. Samples were reconstituted in deionized water, then 10 μl was injected into a Thermo Scientific Dionex ICS-5000+ capillary ion chromatography (IC) system containing a Thermo IonPac AS11 250 × 2 mm 4 μm column. IC flow rate was 360 µl/min (at 30°C), and the gradient conditions are as follows: initial 1 mM KOH, increased to 35 mM at 25 min, increased to 99 mM at 39 min, and held at 99 mM for 10 min. The total run time was 50 min. To increase desolvation for better sensitivity, methanol was delivered by an external pump and combined with the eluent via a low dead volume mixing tee. Data were acquired using a Thermo Orbitrap Fusion Tribrid Mass Spectrometer under ESI negative mode and imported to Thermo Trace Finder software for final analysis. Relative abundance was normalized by tissue weight.
Analysis of reduced and oxidized coenzymes by triple quadruple LC-MS
Request a detailed protocolTo determine the relative abundance of ubiquinone (oxidized CoQ10), ubiquinol (reduced CoQ10), ubiquinone-9 (CoQ9), and ubiquinol-9 (reduced CoQ9) in mouse liver samples, extracts were prepared and analyzed by Thermo Scientific TSQ triple quadrupole mass spectrometer coupled with a Dionex UltiMate 3000 HPLC system. Approximately 20–30 mg of tissue were pulverized in liquid nitrogen then homogenized with a Precellys Tissue Homogenizer. Coenzymes were extracted with 500 µl ice-cold 100% isopropanol. Tissue extracts were vortexed, centrifuged at 17,000 × g for 5 min at 4°C, and supernatants were transferred to clean autosampler vials. The mobile phase was methanol containing 5 mM ammonium formate. Separations of CoQ9, CoQ10, reduced-CoQ9, and reduced-CoQ10 were achieved on a Kinetex 2.6 µm C18 100 Å, 100 × 4.6 mm column. The flow rate was 400 µl/min at 35°C. The mass spectrometer was operated in the MRM positive ion electrospray mode with the following transitions. CoQ10/oxidized: m/z 863.7→197.1 CoQ10/reduced: m/z 882.7→197.1, CoQ9/oxidized: m/z 795.6 →197.1, and CoQ9/reduced: m/z 814.7 →197.1. Raw data files were imported to Thermo Trace Finder software for final analysis. Relative abundance was normalized by tissue weight.
Quantification and statistical analysis
Request a detailed protocolAll measurements were taken from distinct biological samples. Unless otherwise noted, all statistical analyses were performed using GraphPad Prism (version 9) and tests described in the figure legends. In the case of multiple groups, a one- or two-way ANOVA with post-hoc tests were used to determine statistical significance. When only two groups were compared, nonparametric Mann–Whitney tests were used to determine statistical significance. Gene expression and metabolomic heatmaps were plotted as Z-scores using R(4.0.3) and ComplexHeatmap (2.6.2). The species-by-species t-test was applied for metabolomics data to identify the top differentially regulated metabolites that passed the nominal threshold p-values. For multiple comparisons, the Sidak procedure was used for false discovery rate (FDR) correction. Statistical analysis of energy balance was performed by ANCOVA with lean body mass as a co-variate (Mina et al., 2018). No statistical method was used to predetermine sample size. Unblinded analysis of histology was performed by the investigators. All data are expressed as mean ± SEM, unless otherwise specified.
Materials availability
Request a detailed protocolThis study did not generate new unique reagents. The authors declare that reagents utilized are available upon reasonable request to the corresponding author.
Data availability
All data generated or analyzed during this study are included in the manuscript and supporting files. RNA-sequencing data have been deposited in GEO under accession code GSE206200. Numerical data used to generate the figures are attached as source data files to each figure. Figure 5 - Source Data 2 and Figure 6 - Source Data 2 contain Western blot scans.
-
NCBI Gene Expression OmnibusID GSE206200. Vitamin B2 enables PPARa regulation of fasting glucose availability.
References
-
Localization of 1-Deoxysphingolipids to mitochondria induces mitochondrial dysfunctionJournal of Lipid Research 58:42–59.https://doi.org/10.1194/jlr.M068676
-
Report of the American Institute of Nurtition ad hoc Committee on Standards for Nutritional StudiesThe Journal of Nutrition 107:1340–1348.https://doi.org/10.1093/jn/107.7.1340
-
Disorders of riboflavin metabolismJournal of Inherited Metabolic Disease 42:608–619.https://doi.org/10.1002/jimd.12058
-
Hepatic transcriptional responses to fasting and feedingGenes & Development 35:635–657.https://doi.org/10.1101/gad.348340.121
-
Molecular actions of PPARα in lipid metabolism and inflammationEndocrine Reviews 39:760–802.https://doi.org/10.1210/er.2018-00064
-
Hepatocyte-specific PPARA expression exclusively promotes agonist-induced cell proliferation without influence from Nonparenchymal cellsAmerican Journal of Physiology. Gastrointestinal and Liver Physiology 312:G283–G299.https://doi.org/10.1152/ajpgi.00205.2016
-
Impairments of hepatic Gluconeogenesis and Ketogenesis in PPARα-deficient neonatal miceAmerican Journal of Physiology. Endocrinology and Metabolism 307:E176–E185.https://doi.org/10.1152/ajpendo.00087.2014
-
Resting and total energy expenditure of patients with long-chain fatty acid oxidation disorders (LC-Faods)Molecular Genetics and Metabolism 138:107519.https://doi.org/10.1016/j.ymgme.2023.107519
-
The nutrition of the Mouse; a difference in the riboflavin requirements of two highly inbred strainsThe Journal of Nutrition 34:273–283.https://doi.org/10.1093/jn/34.3.273
-
Pathogenic variants in SQOR Encoding sulfide:Quinone Oxidoreductase are a potentially treatable cause of Leigh diseaseJournal of Inherited Metabolic Disease 43:1024–1036.https://doi.org/10.1002/jimd.12232
-
Dietary intake of selected B vitamins in relation to risk of major cancers in womenBritish Journal of Cancer 99:816–821.https://doi.org/10.1038/sj.bjc.6604540
-
Peroxisome Proliferator-activated receptor alpha mediates the adaptive response to fastingThe Journal of Clinical Investigation 103:1489–1498.https://doi.org/10.1172/JCI6223
-
Metaboanalyst 5.0: Narrowing the gap between raw spectra and functional insightsNucleic Acids Research 49:W388–W396.https://doi.org/10.1093/nar/gkab382
-
A role for Ceramides, but not Sphingomyelins, as antagonists of insulin signaling and mitochondrial metabolism in C2C12 MyotubesThe Journal of Biological Chemistry 291:23978–23988.https://doi.org/10.1074/jbc.M116.737684
-
Riboflavin (vitamin B-2) and healthThe American Journal of Clinical Nutrition 77:1352–1360.https://doi.org/10.1093/ajcn/77.6.1352
-
Interaction among folate, riboflavin, genotype, and cancer, with reference to colorectal and Cervical cancerThe Journal of Nutrition 135:2960S–2966S.https://doi.org/10.1093/jn/135.12.2960S
-
Fatty acid oxidation disordersAnnual Review of Physiology 64:477–502.https://doi.org/10.1146/annurev.physiol.64.082201.154705
-
Riboflavin metabolismThe New England Journal of Medicine 283:463–472.https://doi.org/10.1056/NEJM197008272830906
-
Lactic Acidemia and mitochondrial diseaseMolecular Genetics and Metabolism 89:3–13.https://doi.org/10.1016/j.ymgme.2006.05.015
-
Metabolic adaptations in the absence of Perilipin: Increased beta-oxidation and decreased hepatic glucose production associated with peripheral insulin resistance but normal glucose tolerance in Perilipin-null miceThe Journal of Biological Chemistry 279:35150–35158.https://doi.org/10.1074/jbc.M405499200
-
Mitochondrial metabolism mediates oxidative stress and inflammation in fatty liverThe Journal of Clinical Investigation 125:4447–4462.https://doi.org/10.1172/JCI82204
-
Transcriptional control of energy metabolism by nuclear receptorsNature Reviews. Molecular Cell Biology 23:750–770.https://doi.org/10.1038/s41580-022-00486-7
-
Isoform 1C of Sterol regulatory element binding protein is less active than Isoform 1A in livers of transgenic mice and in cultured cellsThe Journal of Clinical Investigation 99:846–854.https://doi.org/10.1172/JCI119248
-
SoftwareSPP Ominer U, version swh:1:rev:d95703fb6c298f3ed68e96bf8f35aa02f730dbb7Software Heritage.
-
Loss of hepatic PPARα promotes inflammation and serum hyperlipidemia in diet-induced obesityAmerican Journal of Physiology. Regulatory, Integrative and Comparative Physiology 317:R733–R745.https://doi.org/10.1152/ajpregu.00153.2019
-
Metabolic effects of Bezafibrate in mitochondrial diseaseEMBO Molecular Medicine 12:e11589.https://doi.org/10.15252/emmm.201911589
-
Evidence for two sites of superoxide production by mitochondrial NADH-Ubiquinone Oxidoreductase (complex I)The Journal of Biological Chemistry 286:27103–27110.https://doi.org/10.1074/jbc.M111.252502
-
A nutrient-wide Association study on blood pressureCirculation 126:2456–2464.https://doi.org/10.1161/CIRCULATIONAHA.112.114058
-
Isovaleric Acidemia: New aspects of genetic and Phenotypic heterogeneityAmerican Journal of Medical Genetics. Part C, Seminars in Medical Genetics 142C:95–103.https://doi.org/10.1002/ajmg.c.30089
-
Effect of insulin on utilization and production of circulating glucoseThe American Journal of Physiology 189:43–50.https://doi.org/10.1152/ajplegacy.1957.189.1.43
-
Bezafibrate induces Autophagy and improves hepatic lipid metabolism in Glycogen storage disease type IAHuman Molecular Genetics 28:143–154.https://doi.org/10.1093/hmg/ddy343
-
Degradation and excretion of riboflavin in the ratThe Journal of Nutrition 93:445–453.https://doi.org/10.1093/jn/93.4.445
-
Serine catabolism feeds NADH when respiration is impairedCell Metabolism 31:809–821.https://doi.org/10.1016/j.cmet.2020.02.017
Article and author information
Author details
Funding
American Diabetes Association (1-18-IBS-105)
- Sean M Hartig
National Institute of Diabetes and Digestive and Kidney Diseases (R01DK114356)
- Sean M Hartig
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was funded by the Nancy Chang, PhD Award for Research Excellence at Baylor College of Medicine, American Diabetes Association #1-18-IBS-105 and NIH R01DK114356 (to SMH). We thank the Gene Vector Core at Baylor College of Medicine for vectors and support. Other BCM core services received support from NCI P30CA125123: Human Tissue Acquisition and Pathology, Metabolomics, and Integrated Microscopy.
Ethics
This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All of the animals were handled according to approved institutional animal care and use committee (IACUC) protocols (AN-6411) of the Baylor College of Medicine. Every effort was made to minimize suffering.
Copyright
© 2023, Masschelin et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,326
- views
-
- 137
- downloads
-
- 3
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cell Biology
The β-catenin-dependent canonical Wnt signaling is pivotal in organ development, tissue homeostasis, and cancer. Here, we identified an upstream enhancer of Ctnnb1 – the coding gene for β-catenin, named ieCtnnb1 (intestinal enhancer of Ctnnb1), which is crucial for intestinal homeostasis. ieCtnnb1 is predominantly active in the base of small intestinal crypts and throughout the epithelia of large intestine. Knockout of ieCtnnb1 led to a reduction in Ctnnb1 transcription, compromising the canonical Wnt signaling in intestinal crypts. Single-cell sequencing revealed that ieCtnnb1 knockout altered epithelial compositions and potentially compromised functions of small intestinal crypts. While deletion of ieCtnnb1 hampered epithelial turnovers in physiologic conditions, it prevented occurrence and progression of Wnt/β-catenin-driven colorectal cancers. Human ieCTNNB1 drove reporter gene expression in a pattern highly similar to mouse ieCtnnb1. ieCTNNB1 contains a single-nucleotide polymorphism associated with CTNNB1 expression levels in human gastrointestinal epithelia. The enhancer activity of ieCTNNB1 in colorectal cancer tissues was stronger than that in adjacent normal tissues. HNF4α and phosphorylated CREB1 were identified as key trans-factors binding to ieCTNNB1 and regulating CTNNB1 transcription. Together, these findings unveil an enhancer-dependent mechanism controlling the dosage of Wnt signaling and homeostasis in intestinal epithelia.
-
- Cell Biology
- Stem Cells and Regenerative Medicine
Our understanding of the transitions of human embryonic stem cells between distinct stages of pluripotency relies predominantly on regulation by transcriptional and epigenetic programs with limited insight on the role of established morphological changes. We report remodeling of the actin cytoskeleton of human embryonic stem cells (hESCs) as they transition from primed to naïve pluripotency which includes assembly of a ring of contractile actin filaments encapsulating colonies of naïve hESCs. Activity of the Arp2/3 complex is required for the actin ring, to establish uniform cell mechanics within naïve colonies, promote nuclear translocation of the Hippo pathway effectors YAP and TAZ, and effective transition to naïve pluripotency. RNA-sequencing analysis confirms that Arp2/3 complex activity regulates Hippo signaling in hESCs, and impaired naïve pluripotency with inhibited Arp2/3 complex activity is rescued by expressing a constitutively active, nuclear-localized YAP-S127A. Moreover, expression of YAP-S127A partially restores the actin filament fence with Arp2/3 complex inhibition, suggesting that actin filament remodeling is both upstream and downstream of YAP activity. These new findings on the cell biology of hESCs reveal a mechanism for cytoskeletal dynamics coordinating cell mechanics to regulate gene expression and facilitate transitions between pluripotency states.