Bioenergetic requirements of hematopoietic stem cells and pluripotent stem cells (PSCs) vary with lineage fate, and cellular adaptations rely largely on substrate (glucose/glutamine) availability and mitochondrial function to balance tricarboxylic acid (TCA)-derived anabolic and redox-regulated antioxidant functions. Heme synthesis and degradation converge in a linear pathway that utilizes TCA cycle-derived carbon in cataplerotic reactions of tetrapyrrole biosynthesis, terminated by NAD(P)H-dependent biliverdin reductases (IXα, BLVRA and IXβ, BLVRB) that lead to bilirubin generation and cellular antioxidant functions. We now demonstrate that PSCs with targeted deletion of BLVRB display physiologically defective antioxidant activity and cellular viability, associated with a glutamine-restricted defect in TCA entry that was computationally predicted using gene/metabolite topological network analysis and subsequently validated by bioenergetic and isotopomeric studies. Defective BLVRB-regulated glutamine utilization was accompanied by exaggerated glycolytic accumulation of the rate-limiting hexokinase reaction product glucose-6-phosphate. BLVRB-deficient embryoid body formation (a critical size parameter of early lineage fate potential) demonstrated enhanced sensitivity to the pentose phosphate pathway (PPP) inhibitor 6-aminonicotinamide with no differences in the glycolytic pathway inhibitor 2-deoxyglucose. These collective data place heme catabolism in a crucial pathway of glutamine-regulated bioenergetic metabolism and suggest that early stages of lineage fate potential require glutamine anaplerotic functions and an intact PPP, which are, in part, regulated by BLVRB activity. In principle, BLVRB inhibition represents an alternative strategy for modulating cellular glutamine utilization with consequences for cancer and hematopoietic metabolism.
Bioenergetic requirements of hematopoietic stem cells (HSCs) and pluripotent stem cells (PSC) vary with lineage fate, and cellular adaptations rely largely on substrate (glucose/glutamine) availability and mitochondrial function to balance TCA (tricarboxylic acid)-derived anabolic requirements, redox homeostasis, and reactive oxygen species (ROS) . These adaptations are especially relevant to PSCs which maintain a balance between metabolically quiescent, non-cycling (ROSlow) cells in the hypoxic bone marrow niche and a metabolically active (ROShigh) phenotype associated with the switch to aerobic (oxidative phosphorylation) metabolism and proliferation. ROS accumulation promotes lineage-restricted hematopoietic development , although ROS generation exceeding the capacity of cellular antioxidant mechanisms promotes senescence and apoptosis. Mitochondria and NADPH oxidase are the primary sources of ROS generation , and oxidation of glucose and glutamine via the mitochondrial TCA cycle provides the critical source of biosynthetic precursors required for cellular proliferation. Quiescent stem cells largely depend on glycolysis for ATP synthesis, and glutamine oxidation is indispensable for survival of PSCs , serves as a crucial mitochondrial substrate for cancer cells , and regulates human HSC lineage specification . Similarly, since the NAD+/NADH ratio is regulated by glycolytic and mitochondrial activities that change dramatically during differentiation or reprogramming , the NAD+/NADH redox state may have a role in driving PSC fate.
Heme (protoporphyrin IX complexed with iron) serves as an indispensable cofactor for all aerobic cells, providing the prosthetic group for phylogenetically distinct hemoproteins functioning in the electron transport chain (ETC), gas exchange, or as light-sensing plant or bacterial phytochromes . Heme biosynthesis is initiated by the condensation of TCA-derived succinyl CoA and glycine, whose generation from serine provides 1-carbon units coupled to the folate cycle for purine and thymidine biosynthesis . The cytotoxicity of free heme is minimized by a two-step catabolic reaction initiated by heme oxygenases (HMOX1 and HMOX2), with sequential generation of biliverdin (BV) and bilirubin (BR) tetrapyrroles in a major pathway functioning in cellular defenses against oxidative and nitrosative stress [10,11]. BV to BR derivatization is regulated by two non-redundant BV reductases (BLVRA and BLVRB), differing by their isomeric specificity in NAD(P)H-dependent oxidation/reduction (redox) coupling. BLVRA retains specificity for the predominant BV found in adults (BVIXα) with less efficient utilization of BV IXβ, IXγ, or IXδ [12,13], while BLVRB is promiscuous, catalyzing the NAD(P)H-dependent reduction of non-IXα BVs , several flavins , and pyrroloquinoline quinones . Thus, heme synthesis and degradation converge in a linear pathway that utilizes TCA cycle-derived carbon in cataplerotic reactions ultimately linked to BR generation and cellular antioxidant functions .
Preferential reliance on glycolytic pathways is a common bioenergetic feature of PSC and cancer metabolism despite low efficient ATP generation [1,5], and both PSCs and cancer cells reside in hypoxic niches and rely on glutamine metabolism for survival and growth [4,5]. The overlapping glucose/glutamine requirements for both heme biosynthesis and PSC substrate utilization suggested that dysregulated heme metabolism would affect bioenergetic utilization. We focused on BLVRB because of its acknowledged activity in early fetal development , coupled with recent evidence for a redox-regulated function governing hematopoietic lineage fate .
Materials and methods
All regents were of the highest purity commercially available. BV was purchased from Frontier Scientific, Inc. (Logan, UT); BR, NADPH, and FMN (flavin mononucleotide) were purchased from Sigma–Aldrich. High-performance liquid chromatography (HPLC) grade acetonitrile, methanol, and water were purchased from Burdick & Jackson (Morristown, NJ).
Bioethical studies using induced PSCs
Studies using induced PSCs were approved by the Stony Brook IRB (Institutional Review Board) in accordance with the Declaration of Helsinki .
Genome editing in PSCs
Induced PSCs derived from CD34+ human umbilical cords (NCRM1) were obtained from the NIH Center for Regenerative Medicine and propagated on matrigel-coated surfaces under normoxic conditions in feeder- and serum-free cell culture medium (mTeSR, Stem Cell Technologies, Vancouver, BC). Single-guide RNA (sgRNA) oligonucleotides were designed using the clustered regularly interspaced short palindromic repeats (CRISPR) Design Tool (http://crispr.mit.edu/) to minimize off-target effects and selected to precede a 5′-NGG protospacer-adjacent motif (PAM). sgRNAs used in the present study were synthesized to target BLVRB exon 3 nucleotides 430–449 (RefSeq ID NM_000713.2) contiguous to Serine 111, sequence as follows: sgRNA1: 5′-TCATGGTGTGGACAAGGTCG-3′; control sgRNA: 5′-GTAGCGAACGTGTCCGGCGT-3′. Oligonucleotides were annealed and cloned into the BbsI site of the Streptococcus pyogenes Cas9 (SpCas9) nuclease (plasmid PX458; Addgene, Cambridge, MA), thereby allowing for transfection monitoring using green fluorescent protein (GFP) detection. NCRM1 transfection was completed using Lipofectamine 3000 (Life Technologies), and GFP+ cells were sorted and collected by flow cytometry for expansion and further characterization. Putative off-target effects were excluded by RNASeq and detailed sequence alignments of the top 50 sites predicted for the guide sequence (http://crispr.mit.edu).
PSC growth curves and embryoid body formation
PSCs plated at a density of 5 × 103 cells/well were propagated in mTesR media and harvested at specific time points using Accutase (Innovative Cell Technologies) for quantification. Phosphatidylserine-expressing (apoptotic) cells were detected using 5 μl Annexin-V-FITC for 15 min at 25°C in the dark. Samples (both counts and Annexin V-positivity) were immediately analyzed using a BD LSRFortessa flow cytometer (BD Biosciences) using the 488 nm laser to detect Annexin-V-FITC (490/525 nm), or absolute cell quantification using a fixed flow rate of 100 μl/min. The effect of 2-deoxyglucose (2-DG, Sigma–Aldrich) or 6-AN (6-aminonicotinamide, Cayman Chemical) on cell growth was evaluated in an identical fashion, modified for supplementation at 2-day intervals. Embryoid bodies (EBs) were generated and quantified as previously described . Briefly, cells were initially dissociated using Accutase (Stem Cell Technologies), followed by pelleting and resuspension in mTeSR medium containing 10 μM ROCK inhibitor (R&D Systems) and various concentrations (or not) of 6-AN or 2-DG. Viable cells were counted using trypan blue exclusion and equivalent numbers (1 × 105 cells) were plated by gentle centrifugation (300×g for 3 min) into AggreWell plates (Stem Cell Technologies), followed by 37°C incubation to allow for EB formation. EB formation was confirmed by immunofluorescence staining using the pluripotency marker Oct3/4 (R&D Systems); quantification and microscopic visualization were completed at Day 3, and individual EB size was quantified using the Image J software.
Statistical comparisons were completed using ANOVA or t-test, modified for bootstrap resampling to compare relative EB size change in the presence of 6-AN; for all comparisons, a value of P < 0.05 was applied for establishing statistical significance.
Details of other experimental procedures are described in Supplementary Materials and Methods.
BLVRB-deficient ( BLVRB−/−) PSCs
Targeted BLVRB disruption in PSCs was completed by genome editing using PSCs derived from human umbilical cord CD34+ hematopoietic cells (Figure 1A and Supplementary Figure S1). The guide RNA (gRNA) strand (or control gRNA) designed to target BLVRB exon 3 contiguous to the active-site serine111  was introduced into a plasmid encoding the single-strand nickase Cas9n and GFP, followed by flow cytometric enrichment of GFP+ cells. Sequence analysis of six expanded subclones confirmed the presence of distinct exon 3 mutations, and two independent subclones (BLVRB−/−CL19 and BLVRB−/−CL20) encompassing a 6-bp BLVRB deletion had no evidence for genetic aberrations (in/dels) as established both by RNASeq and focused off-target prediction algorithms . Unsupervised hierarchical clustering using the limited subset of differentially expressed genes (N = 42, P-value < 0.001 to minimize false discovery rate) confirmed genetic co-segregation of both BLVRB−/− PSCs from either naive or control PSCs (Figure 1B and Supplementary Figure S1). Only two genes (HIST1H4F and MAGEA4) displayed >2-fold difference between BLVRB−/−CL19 and BLVRB−/−CL20 (Supplementary Table S1); importantly, both BLVRB−/−CL19 and BLVRB−/−CL20 demonstrated concordant down-regulation of both HIST1H4F (∼67-fold) and MAGEA4 (∼40-fold) compared with control and naive iPSCs, establishing the relative genetic homogeneity compared with control or naive iPSCs.
Generation and characterization of
Immunoblot analysis confirmed the absence of BLVRB protein in both clones, with no evidence that BLVRB-deficient PSCs had reciprocal induction of BLVRA protein or enzymatic activity (Figure 1C,D). HMOXs have potent effects on regulation of cellular integrity and function, and we confirmed that BLVRB-deficient PSCs had no statistically significant alterations in inducible HMOX1 or constitutive HMOX2 transcripts (Figure 1E). Similarly, flow cytometric analyses using pluripotency markers Tra-1-60 and SSEA4 were nearly identical between the genotypes, with modest loss of SSEA3 and CD117 (c-Kit proto-oncogene) in BLVRB-deficient PSCs (more evident in BLVRB−/−CL20); nonetheless, there was no evidence for exaggerated differentiation as established using KDR (vascular endothelial growth factor receptor 2) as the hemato-endothelial marker (Supplementary Figure S1).
BLVRB maintains a requisite function in PSC oxidant neutralization
Both BVs and BRs retain signaling and metabolic functions [12,13] and a BV–BR redox cycle has been proposed as an amplifiable potent buffer against ROS and oxidant stress [10,11]. Virtually, all mammalian BVs are generated by regioselective cleavage at the α-meso carbon (BV IXα) , and across the genotypes, we were unable to detect secreted BV IXβ isomers in any of the samples [sensitivity 1 ng/μl, although secreted BV IXα isomers were readily detectable and similar across all PSCs (Figure 1F)]. Thus, PSCs (like adult tissues) appear to have limited capacity for non-BV IXα generation under steady-state conditions . Immunofluorescent microscopy using a cell-permeant fluorogenic probe to quantify ROS in live cells demonstrated exaggerated ROS accumulation in both BLVRB−/−CL19 and BLVRB−/−CL20 compared with control PSCs (Figure 1G). The greatest ROS accumulation was most evident along the metabolically active cellular leading edge, although a statistically significant increase was evident throughout the culture in both BLVRB−/− PSCs. Defective ROS handling was associated with time-dependent loss of cell number in both BLVRB−/−CL19 and BLVRB−/−CL20 PSCs compared with control (Figure 1H), results which were largely accounted for by apoptotic cell death as documented by cell-surface phosphatidylserine exposure using Annexin V binding (Figure 1I). These collective data established that BLVRB maintained a physiologically significant and requisite role in metabolic regulation of PSC antioxidant homeostasis and cellular viability. The lack of identifiable defects affecting isomeric BV IXα/IXβ accumulation suggested that the mechanism(s) leading to these phenotypic differences extended beyond those directly associated with BLVRB-dependent redox coupling.
Perturbed metabolic networks in
Since interconnected metabolic networks are highly adaptive to single-gene perturbations , we developed a genome-scale systems biology approach to better delineate pathway perturbations causally implicated in BLVRB-deficient bioenergetic responses. Across the genotypes, we quantified the steady-state accumulation of 290 intracellular metabolites using multiple reaction monitoring and microcapillary liquid chromatography/tandem mass spectrometry (LC–MS/MS, N = 6 determinations/genotype) and applied the Pearson correlation coefficient (r) to sequentially identify BLVRB-associated metabolites across the PSCs (r ≥ 0.8), followed by a similar approach to identify the genetic subset (r ≥ 0.8) correlated with at least one of these metabolites (Figure 2A). Hierarchical clustering of this extended correlation matrix (N = 45 metabolites and 1932 genes) delineated gene/metabolite pairs co-segregating into three predominant biclusters (P < 2.2 × 10−16), whose metabolic enrichment and topological pathway analyses identified a limited number of perturbed pathways (N = 9 using P < 0.05), collectively involved in nucleotide biosynthesis (purine and pyrimidine), redox homeostasis (pentose phosphate pathway [PPP]), and bioenergetics/biosynthesis (TCA cycle) (Figure 2B and Supplementary Table S2). Delineation of biological connectivity was further elucidated by constructing a comprehensive network encompassing the gene/metabolite subsets (Figure 2C). The resulting network self-organized into distinct metabolic modules centrally anchored by BLVRB. Near-neighbor connectivity was closest for the known BLVRB substrate FMN , along with one TCA cycle intermediate fumarate and maleic acid (hydrogenated succinic acid).
Computational studies identify
BLVRB-regulated gene/metabolic networks.
The collective pathway topology and network structure placed BLVRB at the intersection of converging pathways that utilize glucose and glutamine for biosynthesis and bioenergetics, along with the maintenance of cellular reducing capacity via PPP-mediated glycolysis. Indeed, while glycolytic entry into the oxidative component of the PPP maintains redox homeostasis by generating NAD(P)H reducing equivalents, the downstream non-oxidative branch of the cycle generates ribose-5-phosphate, the essential nucleotide component used for purine/pyrimidine synthesis and riboneogenesis . Comparison of glycolytic metabolites across the genotypes confirmed nearly 9-fold accumulation of glucose-6-phosphate (G6P) in BLVRB−/− PSCs compared with control PSCs (P < 1 × 10−5), associated with statistically decreased levels of all downstream glycolytic metabolites (Figure 2D). G6P maintains the intracellular glucose pool, and its generation by the rate-limiting hexokinase reaction has been implicated in the switch to aerobic glycolysis that occurs during pluripotency reprogramming , consistent with a phenotypic switch of exaggerated glucose utilization in BLVRB−/− PSCs. Similar but less dramatic statistically different metabolic changes were evident in BLVRB−/− PPP (three of seven metabolites) and TCA (3 of 10 metabolites) pathways (Figure 2D), two of which (α-ketoglutarate and succinyl CoA) have roles in the maintenance of stem cell pluripotency . Furthermore, while the NADP+/NADPH and glutathione (reduced/oxidized) ratios were unchanged across the genotypes, BLVRB−/− PSCs demonstrated a diminished NAD+/NADH ratio (Figure 2E–G), results collectively suggest a discrete bioenergetics defect with an adaptive glycolytic phenotype.
BLVRB−/− PSCs have a bioenergetic defect involving mitochondrial function and the ETC
The TCA cycle provides the primary source of mitochondrial NADH for the ETC and cellular respiration, prompting more detailed evaluation of mitochondrial function. MitoTracker green fluorescence demonstrated no obvious differences across the PSCs (Figure 3A), results which were extended using cell-permeant TMRE (tetramethylrhodamine ethyl ester) fluorescence as a functional marker of the mitochondrial transmembrane potential (ΔΨm). Across the PSCs, we saw no evidence that substrate (i.e. glucose [Glc] or glutamine [Gln]) restriction affected TMRE fluorescence; furthermore, mitochondrial depolarization using the protonophore carbonyl cyanide m-chlorophenyl hydrazine (CCCP) was identical across all the PSCs, results collectively excluding a discrete mitochondrial defect exacerbated by substrate availability. We then compared basal and maximal oxygen consumption rates (OCRs) across the PSCs, modified for the presence or absence of Glc or Gln as primary TCA anaplerotic substrate(s). In the absence of either substrate, basal (oligomycin-sensitive) OCRs were similar, confirming that the rate of respiration used to drive ATP synthesis was stable across the PSCs (Figure 3B). These results contrasted to those using Gln as the primary fuel, wherein control PSCs displayed an increased OCR that was blunted or diminished in both BLVRB−/−CL19 and BLVRB−/−CL20 PSCs. Across the genotypes, the availability of both Gln + Glc rescued the Gln response, confirming preferential and sufficient utilization of Glc for basal bioenergetics . Since basal OCRs do not adequately reflect PSC response to situations of increased energy demand, we determined protonophore (FCCP, carbonyl cyanide-p-trifluoromethoxyphenylhydrazone)-stimulated respiration as a measure of mitochondrial substrate oxidation and maximal respiratory capacity . In the absence of either substrate (–Gln, –Glc), both BLVRB−/−CL19 and BLVRB−/−CL20 PSCs displayed defective maximal OCRs, with a persistent defect in Gln-mediated OCRs compared with control PSCs (most pronounced in BLVRB−/−CL20 PSCs); across the genotypes, Glc supplementation rescued the Gln defect (Figure 3C). Defective Gln-mediated basal and maximal OCRs were consistent with bioenergetic uncoupling of ATP synthesis and oxygen consumption , results which were underscored by striking loss of Gln-mediated proton leak-linked respiration in BLVRB−/−CL19 and BLVRB−/−CL20 PSCs (Figure 3D). These collective analyses were most consistent with a discrete defect involving Gln entry into the TCA cycle as the primary bioenergetics defect, largely (if not fully) reversible in the presence of Glc supplementation under either basal or maximal conditions.
Bioenergetic defects in
BLVRB-deficient stem cells.
Defective glutamine TCA entry in
Glutamine retains key metabolic functions by donating its nitrogen and carbon into various growth-promoting pathways , while its oxidation to α-ketoglutarate (α-KG) provides for TCA entry with effects on PSC maintenance  and survival  (Figure 4A). Metabolic flux quantification using [U-13C5] glutamine demonstrated nearly identical fractional uptake between control and BLVRB−/− PSCs, with comparable fractional labeling (and conversion) to [U-13C5] glutamate (Figure 4B). These results sharply contrasted to glutamate-derived α-KG where both BLVRB−/−CL19 and BLVRB−/−CL20 demonstrated statistically significant decreased accumulation of M+5 (i.e. labeled) isotopomers with decreased M+0 (endogenous) isotopomer in BLVRB−/−CL19. Expectedly, diminished α-KG availability resulted in decreased accumulation of downstream metabolites fumarate and malate (M+0 and M+4) isotopomers with diminished succinate M+4 in BLVRB−/−CL20 (Figure 4D–G), results confirming the presence of a discrete Gln defect in TCA entry and rewired glutamate metabolism. In contrast, basal metabolic flux using [U6-13C] glucose revealed high-level and comparable (>97%) fractional labeling of 13C lactate across the genotypes (Figure 4G), confirming an intact and saturated glycolytic metabolism that was comparable between control and BLVRB−/− PSCs and unaffected by defective TCA glutamine entry.
Glutamine-restricted TCA defect in
BLVRB-deficient stem cells.
BLVRB maintains PPP integrity during EB formation
Defective glutamate TCA entry suggested that the metabolic requirements of BLVRB−/− PSCs were adapted for differential glucose utilization, either for bioenergetics through the glycolytic pathway or for NAD(P)H reductive capacity and ribosylated nucleotides through PPP. In the presence of the glycolytic pathway inhibitor 2-DG or the PPP pathway inhibitor 6-AN (a 6-phosphogluconate dehydrogenase inhibitor ), both control and BLVRB−/− PSCs demonstrated comparable dose-dependent growth inhibition (Figure 5A,B), results highlighting a requisite role for both pathways in cell growth independent of BLVRB expression and glutamate anaplerotic functions. We then compared these results with EB three-dimensional aggregate formation, a definable endpoint of early ES cell differentiation that recapitulates embryonic gastrulation and functions as a critical size parameter of lineage fate potential [30,31]. In the absence of any inhibitors, BLVRB−/−CL19 and BLVRB−/−CL20 EBs displayed smaller size distributions compared with control EBs, establishing a requisite role for BLVRB-dependent glutamine anaplerotic functions in molecular regulation of early lineage fate potential (Figure 5C). Furthermore, 2-DG supplementation at a concentration (0.5 mM) maintaining >75% fractional growth survival in suspension cultures across the genotypes had a profound effect on EB formation (Figure 5D), inhibiting EB formation across the genotypes and establishing a requisite role for the glycolytic pathway in EB lineage fate potential unrelated to glutamine TCA entry. In contrast, 50 μM 6-AN which results in complete inhibition of cell growth across the genotypes in suspension culture (vide supra) did not inhibit EB formation; rather, 6-AN supplementation modulated EB size distribution across the genotypes, with an exaggerated effect on both BLVRB−/− PSCs (Figure 5C,D). Indeed, whereas the relative percent size changes for BLVRB−/−CL19 and BLVRB−/−CL20 in the presence of 6-AN supplementation are not significantly different, the percent change in each is significantly greater than control (P-value < 0.000001), confirming exaggerated 6-AN sensitivity across both BLVRB−/− PSCs. Both 6-AN and 2-DG pathway inhibitor effects were readily evident by the microscopic appearance of BLVRB−/− EBs, with dissolution of EB formation in the presence of 2-DG and morphologically smaller and less compact EB structure with 6-AN supplementation.
BLVRB-regulated glycolytic and PPP pathway functions.
Our data establish that PSCs require BLVRB reductase activity for physiologically relevant antioxidant activity and cellular viability, accompanied by a bioenergetics defect with diminished glutamate TCA entry. Dysregulated TCA glutamine utilization in BLVRB-deficient PSCs links heme generation and degradation with anaplerotic entry into the TCA cycle. TCA-derived succinate functions as one component (with glycine) of the rate-limiting ALAS (aminolevulinic acid synthase) reaction catalyzing the first step of tetrapyrrole biosynthesis and presumably explains the shunted glutamine response as an adaptive means of limiting heme generation. Secondary effects would lead to decreased ETC proteins (many of which functionalize with heme), acquired ETC defects with ROS accumulation , and dysfunctional TCA cycle NAD+/NADH generation. Shunted glutamine was confirmed by isotopomeric tracings at the initial TCA entry point (α-ketoglutarate) and resulted in glutamine-restricted defects in oxygen consumption under both basal and maximal respiration.
These data are consistent with those previously demonstrating the critical role of heme generation/degradation pathways in rescuing the phenotype of the TCA cycle enzyme fumarate hydratase (FH) . FH deficiency leads to truncation of the TCA cycle, fumarate accumulation, and hereditary leiomyomatosis and renal cell cancer ; the extended (glutamine to BR) heme generation/degradation pathway enables FH-deficient cells to use accumulated TCA cycle metabolites as the carbon source and for partial mitochondrial NADH production . Indeed, genetic adaptations of FH-deficient cells specifically demonstrate dysregulated BLVRB expression to the exclusion of any effects on BLVRA, an observation that underscores the crucial role of BLVRB reductase activity in bioenergetic metabolism linked to glutamine oxidation and TCA activity. Approaches directed at limiting glutamine utilization have been proposed as strategies for cancer therapeutics and for preparation of purified sources of terminally committed cells free of undifferentiated PSCs ; in theory, the shunted glutamine metabolism in BLVRB-deficient PSCs suggests an alternative approach of modulating cellular glutamine utilization by inhibiting BLVRB activity [33,34].
Loss of antioxidant function in BLVRB−/− PSCs is consistent with previously reported functions for BLVRA in various cell lines [11,35]. Unlike the BLVRA homolog whose pleiotropic cellular effects are mediated by multidomain features , BLVRB structure is generally featureless with a function (and therefore an antioxidant effect) presumably limited to its substrate/NAD(P)H redox activity . Nonetheless, BLVRB−/− PSCs retain normal BLVRA activity, display no differential accumulation of BV IXα isomers, and undetectable BV Ixβ isomers, supporting the working model that BLVRB's antioxidant function occurs (in part) independently of its redox activity, and partially (if not fully) by a requisite function in TCA glutamate entry. Indeed, despite the enrichment of non-BV (BR) Ixα isomers (i.e. IXβ, IXγ, and IXδ) in early fetal development , identification of these isomers is rare in mammals and largely restricted to distinct HMOX variants found in bacteria such as Pseudomonas aeruginosa . When integrated with our data, generation of these isomers appears to occur during a later (committed) stage of embryogenic development beyond that of pluripotency. These data support prior speculation that BLVRB may utilize alternative substrate(s) for its redox coupling and antioxidant effect(s) ; indeed, our computational network model places FMN most proximate to BLVRB (Figure 2), although additional yet-to-be substrates cannot be excluded.
Energy metabolism switches from glycolysis to oxidative phosphorylation with cellular differentiation, although preferential reliance on glycolytic pathways is a common bioenergetic feature of PSCs despite low efficient ATP generation [1,5]. The glutamine bioenergetic defect was largely reversible using glucose as a substrate, and indeed, there was striking nearly 10-fold increased accumulation of the rate-limiting hexokinase metabolite G6P in BLVRB−/− PSCs, consistent with enhanced glycolysis for bioenergetics or for redox homeostasis through PPP. Expectedly, both PSC proliferation and EB formation were inhibited by the glycolytic pathway inhibitor 2-DG irrespective of BLVRB expression, highlighting functional importance of the adaptive glycolytic phenotype known to accompany pluripotency reprogramming . Interestingly, dichotomous results were seen in the presence of PPP inhibition, the bifurcating arm of the glycolytic pathway largely responsible for generation of NADPH reducing equivalents and ribose-5-phosphate that is required during riboneogenesis . Indeed, studies in Saccharomyces cerevisiae suggest that a thermodynamically driven route of ribose production may be uncoupled from NADPH generation contingent on cellular requirements, an observation that may explain the differential sensitivity to 6-AN during proliferative vs. early-stage EB formation, whose size defines the critical first developmental stage regulating lineage fate [30,31]. Furthermore, the observation that BLVRB−/− EB formation displayed enhanced 6-AN sensitivity suggests that BLVRB plays a requisite role in providing PPP support during early stages of EB formation-associated lineage fate potential. BLVRB preferentially uses NADPH for its redox coupling , although its ability to directly regulate cellular NADPH levels is unestablished. Alternatively, BLVRB-regulated PPP support could be mediated by glutamine's function as a nitrogen donor for purine/pyrimidine bases during the biosynthetic transition to multicellular growth and differentiation. While the precise mechanism remains under investigation, our data suggest that glycolytic and PPP pathways maintain requisite functions in proliferative capacity of PSCs independently of glutamine anaplerotic functions, although early stages of lineage fate potential require glutamine anaplerotic functions and an intact PPP pathway which is, in part, regulated by BLVRB activity.
Intracellular levels of α-KG maintain the pluripotency of embryonic stem cells, and diminished levels are sufficient to regulate chromatin modifications including histone methylation and PSC fate . The diminished α-KG levels in BLVRB-deficient PSCs establish that BLVRB expression and function is responsible for the maintenance of the intracellular α-KG pool. Indeed, BLVRB-deficient PSCs demonstrate striking down-regulation of the histone cluster family member HIST1H4F, consistent with the initiation of an epigenetic switch in the transition towards lineage commitment. Previous data have demonstrated that glutamine-supported nucleotide biosynthesis is required for hematopoietic (erythroid) lineage specification, with a shift to myelomonocytic differentiation in the setting of limiting glutamine . Interestingly, BLVRB transcript demonstrates ∼40-fold induction during erythropoiesis and a restricted expression during early stages of megakaryocytopoiesis . These collective observations provide a metabolic mechanism whereby BLVRB-regulated glutamine metabolism during critical stages of HSC development could modulate lineage fate unrelated to its redox-regulated antioxidant function.
aminolevulinic acid synthase
carbonyl cyanide m-chlorophenyl hydrazine
clustered regularly interspaced short palindromic repeats
electron transport chain
green fluorescent protein
hematopoietic stem cells
oxygen consumption rates
pluripotent stem cells
pentose phosphate pathway
reactive oxygen species
Streptococcus pyogenes Cas9
tetramethylrhodamine ethyl ester
Z.L. and W.F.B. conceptualized this study. Z.L., S.W., D.W., W.Z., G.D.G., and W.F.B. performed methodology of this work. L.E.M., D.V.G., and N.M.N. investigated this study. Z.L. and W.F.B. involved in writing the original draft and review and editing. N.M.N., D.V.G., G.D.G., and W.F.B. acquired funding of this work. Z.L., G.D.G., and W.F.B. collected resources for this work. Z.L., G.D.G., and W.F.B. supervised this study.
This work was supported by grants from the National Institutes of Health NIH/NCRR 1 S10 RR023680-1 (to the Stony Brook Proteomics Facility); NIH/HL091939 (W.F.B.); NIH/HL129545 (N.M.N.), NIH/AG059277 (D.V.G.), and the New York State Stem Cell Foundation (C026716, to W.F.B.). CRISPR plasmid pSpCas9(BB)-2A-GFP (PX458) was provided by F.Z.
We thank Wan-yi Yen for assistance with image capture and analysis, Irina Zaitseva for protoeomic studies, Emily Montal for assistance with 13C tracing analyses, and Christian Ruiz for assistance with oxygen consumption and 13C tracing analyses.
The Authors declare that there are no competing interests associated with the manuscript.