Metabolome and metaboproteome remodeling in nuclear reprogramming

Abstract

Nuclear reprogramming resets differentiated tissue to generate induced pluripotent stem (iPS) cells. While genomic attributes underlying reacquisition of the embryonic-like state have been delineated, less is known regarding the metabolic dynamics underscoring induction of pluripotency. Metabolomic profiling of fibroblasts vs. iPS cells demonstrated nuclear reprogramming-associated induction of glycolysis, realized through augmented utilization of glucose and accumulation of lactate. Real-time assessment unmasked downregulated mitochondrial reserve capacity and ATP turnover correlating with pluripotent induction. Reduction in oxygen consumption and acceleration of extracellular acidification rates represent high-throughput markers of the transition from oxidative to glycolytic metabolism, characterizing stemness acquisition. The bioenergetic transition was supported by proteome remodeling, whereby 441 proteins were altered between fibroblasts and derived iPS cells. Systems analysis revealed overrepresented canonical pathways and interactome-associated biological processes predicting differential metabolic behavior in response to reprogramming stimuli, including upregulation of glycolysis, purine, arginine, proline, ribonucleoside and ribonucleotide metabolism, and biopolymer and macromolecular catabolism, with concomitant downregulation of oxidative phosphorylation, phosphate metabolism regulation, and precursor biosynthesis processes, prioritizing the impact of energy metabolism within the hierarchy of nuclear reprogramming. Thus, metabolome and metaboproteome remodeling is integral for induction of pluripotency, expanding on the genetic and epigenetic requirements for cell fate manipulation.

Keywords: :

• energy metabolism
• iPS cells
• glycolysis
• mitochondria
• network biology
• oxidative phosphorylation
• proteomics
• regenerative medicine
• systems biology

Recent advances in the science of nuclear reprogramming have enabled next generation regenerative sources derived from somatic tissues. Transduction with stemness transcription factors is sufficient to direct somatic cells back to a primordial embryonic-like state.¹ This process of bioengineered reprogramming yields induced pluripotent stem (iPS) cells with the capacity to derive all three germ layers without the requirement of an embryo source. The ability to bioengineer genuine pluripotent stem cells holds significant promise to uncover mechanisms of disease, to diagnostically unravel individual variation in disease susceptibility, and to test novel therapeutic strategies.^2-5 As proof of principle, iPS cells have been successfully used as a diagnostic or therapeutic source in a number of life-threatening conditions.^6-9

Pluripotent Stringency Criteria for Bioengineered Stem Cells

To ensure acquisition of pluripotent traits in reprogrammed iPS cells, comprehensive criteria have been deployed, including morphology, gene expression, and functional characteristics.^10-12 Structurally, bona fide iPS cells resemble embryonic stem (ES) cells and share their unlimited ability to self-renew.¹⁰ These bioengineered stem cells have similar, although not identical, gene expression and epigenetic profiles (demethylation of pluripotent gene promoters, presence of bivalent domains in developmental genes, and X chromosome reactivation in cells derived from females) to their embryonic counterparts, and express key pluripotency factors (such as Nanog and Oct4) and specific cell surface markers (mouse SSEA-1).¹⁰ Validated iPS cells also meet functional pluripotency tests to ensure they can differentiate into all three embryonic germ layers. In vitro stringency tests include the use of ES cell protocols to differentiate iPS cells into defined lineages; however, due to lack of standardization, this approach in isolation does not necessarily prove pluripotency.¹³ The highest stringency test for human iPS cells is the in vivo teratoma assay, where a subcutaneous injection of cells into an immunodeficient host allows for multilineage tumor formation, and histological analysis is used to confirm induction of all three germ layers within nascent teratomas. Animal-derived iPS cells must also meet in utero criteria by recapitulating the developmental capacity of ES cells. These tests include the generation of chimeric animals with germline transmission, as well as the ultimate test of tetraploid aggregation for production of live offspring genetically identical to the transplanted iPS cells.¹³ In situ stringency tests may also be utilized to test for regenerative capacity in various disease models. A proposed minimal set of criteria that iPS cells must fulfill encompass: (1) morphological characteristics, including unlimited self-renewal; (2) expression of key pluripotent genes and downregulation of parental cell lineage specific genes; (3) transgene independence; and (4) functional differentiation validated by the highest stringency test available.^10-12 Indeed, it is not practical to assess teratoma formation in a high-throughput manner in every derived cell line. Moreover, for regenerative purposes, the isolation or selection of lines that do not form teratomas may in fact be the preferred source. Therefore identifying high-throughput surrogate markers of functional pluripotency would complement existing criteria of proficient nuclear reprogramming.¹⁴

Metabolic Regulation of Cell Fate

In addition to the extensively characterized genetic and epigenetic regulators of nuclear reprogramming, recent studies indicate that modulation of mitochondrial dynamics and energy metabolism may contribute to the induction and maintenance of pluripotency.^15-25 In this regard, the metabolome of iPS cells resembles that of ES cells and significantly differs from their somatic source, based upon elevated utilization of glucose and production of glycolytic end products, and downregulation of tricarboxylic acid cycle and cellular respiration metabolites.^15,16 In fact, multimodal metabolomic profiling in tandem with functional metabolic approaches revealed a bioenergetic transition from somatic oxidative metabolism to glycolysis during iPS cell derivation.^15,16 Moreover, changes in stoichiometry of the mitochondrial electron transport chain occur in the initial stage of nuclear reprogramming, and elevated glycolytic gene expression preceded induction of pluripotent genes, suggesting that a switch in energy metabolism is required to fuel nuclear reprogramming.^15,26 This metabolic transition from oxidative metabolism to glycolysis is consistent across species, cell lines, and reprogramming formulations, and somatic sources that display a greater glycolytic and lower oxidative capacity demonstrate higher efficiency of nuclear reprogramming.^15-17,27 The energetic conversion can be targeted to control the reprogramming process directly by stimulating glycolysis, either pharmacologically or by increasing glycolytic intermediates, or indirectly through hypoxic stimulation or inhibition of the p53 signaling pathway, increasing pluripotent induction.^15,16,28-39 Nuclear reprogramming can be impaired by inhibiting glycolysis or stimulating oxidative metabolism.^15,16,39 With growing evidence for remodeling of energy metabolism in cell fate decisions, acquisition of a glycolytic metabotype of stem cells may offer early criteria for establishing functional pluripotency.

Metabolic Blueprint of Pluripotency

Metabolic flux in iPS cells

Multimodal bioenergetic characterization of nuclear reprogramming has established remodeling occurring at the level of the metabolome, metabolic infrastructure, and mitochondria.^{15,16,21,27,40} Pluripotent stem cells benefit from acquiring a glycolytic metabotype associated with a faster rate of ATP generation compared with oxidative metabolism and the production of biosynthetic precursors in conjunction with the pentose phosphate pathway, thus enabling iPS cells to fuel anabolic requirements.^40-42 Agents that promote glycolysis expedite the reprogramming process, while stimulation of oxidative metabolism impairs induction of pluripotency.^15,39,43,44 Utilization of nuclear magnetic resonance spectroscopy (NMR) and mass spectrometry (MS) has enabled unprecedented coverage of small molecular weight compounds or metabolites in biological systems.^45-47 Applying such high-throughput technologies to nuclear reprogramming demonstrated that the metabolome of iPS cells segregates away from parental sources and clusters with that of ES cell counterparts, exhibiting a metabolomic extracellular footprint and intracellular fingerprint characteristic of bona fide pluripotent cells.^15,16 Monitoring extracellular metabolites in cell culture media by serial NMR measurements established a greater utilization of glucose and greater production of the glycolytic end product lactate in iPS cells compared with their parental source (Fig. 1A). Profiling of the intracellular metabolome supported an iPS cell signature distinct from the parental source and characterized by a reduction in metabolites associated with oxidative metabolism, defining the transition from cellular respiration to glycolysis.^15,16 Beyond metabolomic profiling, validated extracellular flux technology established functional monitoring of cellular bioenergetics (Fig. 1B–F). Measurement of the oxygen consumption rate (OCR) and extracellular acidification rate (ECAR) enabled real-time assessment of the balance between oxidative metabolism and glycolysis to document nuclear reprogramming-induced metabolic remodeling as indicated by a reduction in the iPS cell OCR/ECAR ratio compared with parental fibroblasts (Fig. 1C). Sequential supplementation of the ATP synthase inhibitor oligomycin, an electron transport chain uncoupler, carbonyl cyanide 4-[trifluoromethoxy]phenylhydrazone (FCCP), and the complex 1 inhibitor rotenone, dissected key parameters of mitochondrial physiology, including maximal respiration rate, reserve capacity, and ATP production (Fig. 1B). Consistent with bioenergetic remodeling and reduction in energy turnover and total ATP levels,^15-18,41,48 iPS cells were found to continuously respire near their maximal rates as indicated by lower maximal respiration (Fig. 1D) and reserve capacity (Fig. 1E), resulting in reduced ATP production (Fig. 1F) compared with the somatic source. Real-time monitoring of cellular bioenergetics thus offers a novel high-throughput platform for assessing functional pluripotency.

Figure 1. Metabolic reprogramming from oxidative metabolism in fibroblasts to glycolysis in iPS cells. Representative serial nuclear magnetic resonance-based metabolomic footprinting of cell culture media indicates elevated glucose utilization and accumulation of lactate in iPS cells (A). High-throughput simultaneous measurement of oxygen consumption rate (OCR) and extracellular acidification rate (ECAR) (B) indicates reduced oxygen utilization and greater reliance on glycolysis in iPS cells (C). Extracellular flux analysis in response to oligomycin, FCCP, and rotenone demonstrates an iPS cell specific mitochondrial profile characterized by reduced maximal respiration (D), reserve capacity (E) and ATP turnover (F). Values represent mean ± standard error, n = 8 per group. Significance was determined by Student t test. *P < 0.05 vs. MEF.

Nuclear reprogramming transforms the mitochondrial infrastructure

The cell’s central energy generators, the mitochondria, also reliably undergo multilevel reorganization during nuclear reprogramming of fibroblasts, indicating mitochondrial metamorphosis as a conserved feature of pluripotency induction.^{15,17,18,48,49} Nuclear reprogramming induces a reduction in mitochondrial density compared with the parental source, with remaining mitochondria displaying an immature ultrastructure consisting of spherical structures with undeveloped cristae, similar to those observed in ES cells (Fig. 2A).^{15,17,18,20,48,49} This is in contrast to the branched tubular shaped mitochondria with electron-dense and well-developed cristae of parental fibroblasts (Fig. 2A). Mitochondrial localization also transitions during nuclear reprogramming from widespread cytoplasmic networks to a perinuclear localization (Fig. 2A), which may reflect reorganization in bioenergetic demands from cytosolic organelles to the nucleus in support of ongoing genetic and epigenetic remodeling.^15,17,18,49 As perinuclear localization is a consistent feature of a number of diverse stem cell populations, it may provide an additional criterion for stemness.^48,50-52 Like their embryonic counterparts, iPS cells accumulate mitochondrial membrane potential-dependent probes to a greater extent than their parental fibroblasts do (Fig. 2B), which may represent a valuable marker of pluripotent induction.^15,20,53 Recent evidence indicates that inhibition of mitochondrial division significantly impairs the ability of fibroblasts to undergo nuclear reprogramming,⁵⁴ supporting an essential role for mitochondrial plasticity in nuclear reprogramming and maintenance of the pluripotent state.

Figure 2. Nuclear reprogramming induced restructuring of the mitochondrial infrastructure. The mitochondrial anatomy and prevalence is remodeled during nuclear reprogramming as resolved by electron microscopy (arrows demarcate mitochondria) (A). iPS cells demonstrate mitochondrial membrane hyperpolarization compared with parental fibroblasts based upon mitochondrial accumulation of JC-1 (B).

Nuclear reprogramming induces remodeling of the metabolic infrastructure

To accommodate the metabolic transition from oxidative metabolism to glycolysis underlying induction of pluripotency, the bioenergetic infrastructure of the cell undergoes extensive remodeling at the level of the transcriptome, proteome, and mitochondria. Glycolysis and oxidative metabolism gene expression changes are consistently observed when examining iPS cells compared with the parental source, revealing a conserved feature underlying nuclear reprogramming dynamics.^16,17,19,55 In particular, genes involved in glucose uptake (SLC2A3 encoding GLUT3), glucose phosphorylation (HK3 and GCK), the final steps of glycolysis (PGAM2, ENO, PKLR, and LDH), and the non-oxidative branch of the pentose phosphate pathway (TKT and RPIA) are upregulated, while the intervening steps of glycolysis (GPI, PFK, and ALDO) are downregulated.^17,19 This suggests that glycolytic intermediates may be shunted into the pentose phosphate pathway to produce anabolic precursors for biosynthetic pathways and to maintain cellular antioxidant defenses.^17,19 Temporal assessment of glycolytic genes during reprogramming revealed the promotion of glycolytic gene expression prior to induction of pluripotent genes, indicating that metabolic remodeling may not be a consequence of reprogramming, but rather required to fuel pluripotency induction.¹⁵ Upstream regulators of gene transcription also support metabolic remodeling during reprogramming. The epigenetic methylation status of genes involved in glycolysis and oxidative metabolism of iPS cells is similar to ES cells and differs from their parental source, indicating an additional biomarker of pluripotent metabolic competence.¹⁶ Recent evidence has implicated a role for miRNAs in controlling cellular fate,^56-58 with miR-302–367 specifically implicated in pluripotent induction and maintenance, potentially by overcoming the metabolic reprogramming barrier and regulating genes involved in oxidative phosphorylation and mitochondrial biogenesis.⁵⁹

Concomitant with transcriptional remodeling of metabolic genes is the restructuring of the metaboproteome, with targeted analysis revealing that the metabolic infrastructure of iPS cells closely resembles ES cells and differs from their somatic source.¹⁵ Specifically, pluripotent cells exhibit predominant upregulation of glycolytic enzymes and changes in stoichiometry of specific mitochondrial electron transport chain complexes.^15,26,60 To examine the priority of these changes with respect to global proteomic remodeling associated with nuclear reprogramming, differential proteomic analysis was performed by two-dimensional electrophoresis together with nanoelectrospray tandem MS to detect and assign identities to significantly altered proteins (Fig. 3; Fig. S1 and Tables S1 and S2). A consistent number of protein species were observed in both fibroblast and iPS cell extracts, and within the 97 spots that differed significantly between groups, there were 441 unique proteins identified, 210 exclusively found in upregulated spots, 166 exclusive to downregulated spots, whereas 65 proteins were identified in both upregulated and downregulated spots. Ontological analysis categorized altered proteins across a spectrum of subcellular locations, including the mitochondria (Fig. S1). Transcriptional and proteomic remodeling of the metabolic infrastructure thus underlies the nuclear reprogramming-induced switch from oxidative metabolism to glycolysis; however the prioritization of these changes within the hierarchy of reprogramming remains unknown.

Figure 3. Nuclear reprogramming induced remodeling of the bioenergetic infrastructure. iPS cell-specific proteome is revealed in representative silver-stained, broad pH range (3–10) large format two-dimensional gels (A). Two-dimensional gel electrophoresis consistently resolved on the order of 550 protein species in fibroblast and iPS cell extracts, with densitometric quantification identifying a subset of 97 significantly altered protein spots (21% of the total resolved) following nuclear reprogramming, with significance determined by Student t test (P < 0.05) (B).

Systems biology prioritizes metabolic protein remodeling

Network analysis was performed to obtain the distinguishing characteristics of upregulated and downregulated subproteomes.⁶¹ Ingenuity Pathway Analysis interaction mapping clustered upregulated proteins within a single fully connected 457 node, 3069 edge network (Fig. 4; Fig. S2), while downregulated proteins were encompassed by a separate 392 node, 3237 edge network (Fig. 5; Fig. S2), both exhibiting scale-free topology. Bioinformatic interrogation of the upregulated subproteome extracted 8 overrepresented canonical pathways (Fig. 6A), half of which were metabolic processes: purine metabolism, glycolysis, arginine/proline metabolism, and pyruvate metabolism. Eight canonical pathways were also prioritized from the downregulated subproteome, including pyruvate metabolism and oxidative phosphorylation (Fig. 6B). Network prioritization of gene ontology (GO) biological processes reflected similarities with prioritized canonical pathways and revealed that metabolic functions related mainly to purine metabolism and glycolysis were uniquely represented in the upregulated network (Fig. 7A), while oxidative phosphorylation, phosphate metabolism regulation, precursor biosynthesis processes, and RNA-related metabolic functions were unique to the downregulated network (Fig. 7B). Appearance of glycolysis only in the upregulated network and oxidative phosphorylation only in the downregulated counterpart (demarcated by arrows in Fig. 7A and B) is consistent with observed functional bioenergetic differences between somatic and iPS cells. Proteomic deconvolution combined with in silico network generation and interrogation of proteome-associated canonical pathways and network-associated biological processes thus provides a systems-based assessment of the metabolic reprogramming associated with pluripotent induction.

Figure 4. Pathway analysis interaction mapping clustered the upregulated subproteome into a fully connected network of 457 nodes and 3069 edges. Network degree distribution exhibits scale-free topology as defined by the relationship between node degree (k) vs. node degree distribution (P[k]), the proportion of all nodes at each specified degree (see Fig. S2A).^68,69

Figure 5. Pathway analysis interaction mapping clustered the downregulated subproteome into a fully connected network of 392 nodes and 3237 edges. Network degree distribution exhibits scale-free topology as defined by the relationship between node degree (k) vs. node degree distribution (P[k]), the proportion of all nodes at each specified degree (see Fig. S2B).^68,69

Figure 6. Pathway analysis interaction mapping prioritizes canonical metabolic processes within the upregulated and downregulated iPS subproteomes. Ingenuity Pathway Analysis identified a number of overrepresented (P < 0.001) metabolic canonical pathways within the upregulated (A) and downregulated proteins arising from somatic cell reprogramming (B).

Figure 7. Gene ontology biological process prioritization within the upregulated and downregulated iPS cell subproteome networks. Unique to the upregulated network were 36 significantly overrepresented (P < 0.001) GO processes, including catabolic/metabolic functions mainly related to purine metabolism and glycolysis, cell cycle functions, and processes related to cell organization/biogenesis, kinase signaling, and an assortment of cell regulation functions (A). Unique to the downregulated network were 34 significantly overrepresented processes, associated with oxidative phosphorylation and several phosphate metabolism regulation and RNA related metabolic functions, as well as cytoskeletal structure/morphology and organization processes, cell localization, and transport (B).

In summary, the metabolic transition from oxidative metabolism to glycolysis is a consistent feature of the nuclear reprogramming process resulting in an ES cell-like metabolome that characterizes generated iPS cells. The present study delineated metabolome and metaboproteome remodeling underlying the nuclear reprogramming of somatic tissue into iPS cells. Indeed, acquiring a glycolytic metabotype may surpass a bioenergetic barrier that cells must achieve to undergo successful reprogramming and establish metabolic competence for regenerative applications. As such, metabolic profiling represents a novel stringency criterion of bioengineered pluripotency.

Materials and Methods

Cell derivation, protein extraction, and quantification

Mouse embryonic fibroblasts were obtained from embryos at 14.5 d post coitum. Resulting fibroblasts were plated for 2–3 passages to full confluency. Transduced fibroblasts were cultured in ES cell maintenance medium, consisting of Dulbecco’s modified Eagle medium (Millipore) supplemented with pyruvate (Lonza), L-glutamine (Invitrogen), non-essential amino acids (Mediatech), 2-mercaptoethanol (Sigma-Aldrich), 15% fetal bovine serum (Invitrogen), and LIF (Millipore).¹⁵ Transfer vectors were generated with pSIN-SEW based vector, pSIN-CSGWdlNotI, with full-length human OCT3/4, SOX2, KLF4, and c-MYC cDNAs (Open Biosystems).⁶ Clonogenic expansion produced reprogrammed cell lines that were maintained in ES cell maintenance medium.

Live cell metabolic analysis

Oxygen consumption rates and extracellular acidification rates were measured using an XF24 Extracellular Flux Analyzer (Seahorse Biosciences). In brief, cells were plated into wells of an XF24 cell culture microplate and maintained overnight to ensure 80% confluence. Prior to assay, plates were equilibrated in unbuffered XF assay medium supplemented with 25 mM glucose, 2 mM glutamax, 1 mM sodium pyruvate, 1× nonessential amino acids, and 1% fetal bovine serum in the absence of CO₂ for 1 h.⁶² Mitochondrial processes were interrogated by serial addition of oligomycin (0.5 μg/ml), FCCP (1 μM), and rotenone (0.5 μM). Each plotted value is the mean of at least 10 replicates and is normalized to baseline oxygen consumption and to total protein quantified using a Bradford assay (Bio-Rad).

Two-dimensional gel electrophoresis and protein identification

Fully confluent cultured cells (MEF, iPS; three 10 cm dishes per cell line) were washed extensively with phosphate buffered saline (10 × 10 mL each). Protein was then extracted in situ from adherent cells in 500 μL lysis buffer (7 M urea, 2 M thiourea, 2% [w/v] 3-[(3-cholamidopropyl)dimethylammonio]-1-propanesulfonate), and quantified in triplicate by Bio-Rad protein assay using the microassay procedure with a bovine γ-globulin standard. As described previously,^63,64 protein extracts (100 μg) were resolved in the first dimension by immobilized pH gradients (pH 3–10) and in the second dimension by SDS-PAGE (12.5%) using a Protean^® II XL system (Bio-Rad). Two-dimensional gels were silver stained, digitized for spot image analysis and normalized by total spot intensity using Bio-Rad PDQuest v.7.4.0.^63,64 Altered protein species were isolated, destained, and prepared for nanoelectrospray linear ion trap tandem MS by reduction, alkylation, tryptic digestion, peptide extraction, and drying.⁶⁵ Peptides were separated on ProteoPep C18 PicoFrit^™ nanoflow column (New Objective) using an Eksigent nanoHPLC system (MDS Sciex) coupled to a hybrid LTQ-Orbitrap mass spectrometer (Thermo Fisher Scientific). Eluted peptide ions were continuously monitored between 375–1600 m/z, with automatic switching to MS/MS collision-induced dissociation mode for ions exceeding an intensity of 8000. Protein identities were assigned by matching multiple peptide spectra to theoretical tryptic fragments in Swiss-Prot (v.53.0) using Mascot^™ v.2.2.⁶⁶

Network derivation and bioinformatic interrogation

Upregulated and downregulated functional networks associated with differentially expressed proteins were identified using Ingenuity Pathway Analysis and Ingenuity Pathways Knowledge Base (Ingenuity^® Systems, http://www.ingenuity.com), as previously described.⁶⁴ Ingenuity was also used to identify significantly overrepresented canonical pathways associated with the altered proteome. Composite networks, constructed by merging Ingenuity functional subnetworks, were depicted using Cytoscape 2.8.3, with topological properties characterized as an undirected network using Network Analyzer.⁶⁷ Interrogation of composite networks was performed with the Cytoscape module BiNGO 2.3 (Bioinformatic Network Gene Ontology), to assess gene ontology (GO) biological process enrichment for interpretation of composite network functional interdependencies. Individual categories were considered for overrepresentation by BiNGO and Ingenuity on the basis of hypergeometric distributions, comparing proportional occurrence within the network relative to overall occurrence in the entire proteome.

Statistical analysis

Unless otherwise indicated, comparison between groups was performed using a Student t test of variables with 95% confidence intervals, with data expressed as mean ± standard error. For BiNGO, a hypergeometric distribution with Benjamini and Hochberg false discovery rate correction was utilized to determine significant overrepresentation of GO biological processes.

Acknowledgments

Supported by National Institutes of Health, American Heart Association, Canadian Institutes of Health Research, Fondation Leducq, Marriott Foundation, and Mayo Clinic Center for Regenerative Medicine.

Disclosure of Potential Conflicts of Interest

No potential conflicts of interest were disclosed.

Supplemental Materials

Supplemental materials may be found here:
http://www.landesbioscience.com/journals/cc/article/25509