- Open Access
Transcriptional landscape of bone marrow-derived very small embryonic-like stem cells during hypoxia
Respiratory Research volume 12, Article number: 63 (2011)
Hypoxia is a ubiquitous feature of many lung diseases and elicits cell-specific responses. While the effects of hypoxia on stem cells have been examined under in vitro conditions, the consequences of in vivo oxygen deprivation have not been studied.
We investigated the effects of in vivo hypoxia on a recently characterized population of pluripotent stem cells known as very small embryonic-like stem cells (VSELs) by whole-genome expression profiling and measuring peripheral blood stem cell chemokine levels.
We found that exposure to hypoxia in mice mobilized VSELs from the bone marrow to peripheral blood, and induced a distinct genome-wide transcriptional signature. Applying a computationally-intensive methodology, we identified a hypoxia-induced gene interaction network that was functionally enriched in a diverse array of programs including organ-specific development, stress response, and wound repair. Topographic analysis of the network highlighted a number of densely connected hubs that may represent key controllers of stem cell response during hypoxia and, therefore, serve as putative targets for altering the pathophysiologic consequences of hypoxic burden.
A brief exposure to hypoxia recruits pluripotent stem cells to the peripheral circulation and actives diverse transcriptional programs that are orchestrated by a selective number of key genes.
Hypoxia is a pathophysiologic condition commonly seen in lung diseases and in many other disorders such as tissue ischemia and cancer . In the developing embryo exposure to global reduction in tissue oxygen is a critical factor in recruitment and selective differentiation of embryonic stem cells that orchestrate organogenesis . The presence of hypoxic microenvironments in adult tissues, as seen for example in solid tumors or during ischemic injury, can also lead to activation and recruitment of stem cells from specific tissue niches, including the bone marrow (BM) . Under such local hypoxic circumstances, resident BM stem cells have been shown to become mobilized and facilitate the structural and functional repair of injured organs [4–6]. The effects of global hypoxia on stem cell activation and recruitment remain poorly understood despite being a common feature of many pulmonary disorders, and increasing evidence that BM-derived circulating progenitor cells may play a causative role in acute lung injury, pulmonary hypertension, pulmonary fibrosis and COPD [7–10].
Recently, Ratajczak et al isolated a homogenous population of rare, small stem cells from adult murine BM mononuclear cells that expressed embryonic lineage markers, and were thus named very small embryonic-like stem cells (VSELs) . Pluripotency of these Sca-1+ Lin- CD45- stem cells was established by demonstrating their ability to differentiate into cells representing all three germ layers . Bone marrow-resident VSELs can be mobilized to the peripheral blood (PB) in response to specific chemoattractant gradients, including stromal derived factor-1 (SDF-1), hepatocyte growth factor (HGF), and leukemia inhibitory factor (LIF) [12, 13]. Since SDF-1 is upregulated in response to hypoxia , it seemed plausible that chemokine-dependent recruitment of VSELs can occur during tissue injury characterized by hypoxia. Indeed, recent reports have demonstrated that these pluripotent stem cells are mobilized from BM to PB following organ damage associated with significant hypoxic burden, including ischemic stroke  and myocardial infarction .
Exposure to reduced oxygen tension causes widespread perturbations in cellular transcription  that is controlled by a number of critical regulators, most prominently hypoxia inducible factors (HIFs) . While gene expression profiling of stem cells under various differentiating conditions has been undertaken , the transcriptional consequences of hypoxia in these cells is less studied. Recently, Westfall et al reported on the effects of 4% vs. 20% oxygen on gene expression of cultured human embryonic stem cells, and concluded that lower oxygen tension (4%) in vitro more accurately captures the true "physiologic" oxygen exposure in vivo by preserving the pluripotent property of these cells . Indeed, a significant difficulty with all in vitro studies of hypoxic exposure is establishing the physiologically relevant oxygen tension during exposure [3, 21] and retaining the in vivo undifferentiated state of stem cells .
To overcome these limitations, we utilized an in vivo murine model to systematically identify the transcriptional response of BM-derived VSELs to hypoxia. We hypothesized that exposing the animals to hypoxia will mobilize VSELs from their BM to peripheral circulation in a stem-cell specific chemokine gradient, and activate a spectrum of transcriptional programs that captures the diverse, pluripotent potential of these cells.
All animal experiments were performed according to protocols approved by the Institutional Animal Care and Use Committee (IACUC) of the University of Louisville. Eight-week-old, adult, male, C57BL/J6 mice were purchased from Jackson Laboratory (Bar Harbor, ME) and housed in a specific pathogen free environment.
Mice were placed in identical commercially designed chambers (Oxycycler model A44XO; Biospherix, Redfield, NY, USA) that were operated under a 12:12-h light-dark cycle (6:00 a.m to 6:00 p.m.). Gas was circulated around each of the chambers, attached tubing, and other units at 10 l/min to maintain low ambient CO2 levels. An O2 analyzer measured the O2 concentration continuously and deviations from the desired concentration were corrected by addition of N2 or O2 through solenoid valves such that the moment-to-moment desired oxygen concentration of the chamber was programmed and adjusted automatically. Ambient CO2 in the chamber was monitored periodically and maintained at < 0.01% by adjusting overall chamber basal ventilation. Humidity was measured and maintained at 40-50% by circulating the gas through a freezer and silica gel. Ambient temperature was kept at 22-24°C. Hypoxic mice were exposed to FiO2 of 8%, a level that is associated with reproducible oxyhemoglobin saturations in the 75-80% range. Control animals were exposed to circulating room air (FiO2: 21%). Animals were exposed to hypoxia or normoxia for a period of 24 h and sacrificed immediately afterwards.
Isolation of VSELs from bone marrow and peripheral blood
VSELs were isolated from murine bone marrow and peripheral blood nucleated cells using multicolor fluorescence-activated cell sorting (FACS) as previously described [4, 12, 23]. Briefly, the population of total nucleated cells (TNCs) from murine PB and BM was obtained after lysis of red blood cells using 1× BD Pharm Lyse Buffer (BD Pharmingen, San Jose, CA, USA). TNCs were subsequently stained for hematopoietic lineage markers (Lineage - lin), for Sca-1 and CD45 antigens, for 30 minutes in medium containing 2% fetal bovine serum. The following anti-mouse antibodies (BD Pharmingen, San Jose, CA, USA) were used for staining: anti-CD45R/B220 (phycoerythrin - PE; clone RA3-6B2), anti-Gr-1 (PE; clone RB6-8C5), anti-T-cell receptor-αβ (PE; clone H57-597), anti-T-cell receptor-γδ (PE; clone GL3), anti-CD11b (PE; clone M1/70), anti-Ter119 (PE; clone TER-119), and anti-Ly-6A/E (Sca-1) (biotin; clone E13-161.7, with streptavidin-conjugated to PE-Cy5) and rat anti-CD45 (allophycocyanin-Cy7; clone 30-F11). All mAbs were added at saturating concentrations. Cells were then washed, re-suspended in RPMI 1640 medium with 10% fetal bovine serum, and sorted by MoFlo cell sorter (Dako, Carpinteria, CA, USA). The Sca-1+Lin-CD45- cells (VSELs) and control Sca-1+Lin-CD45+ cells (HSCs) were isolated according to the gating and sorting strategy described before and outlined in Figure 1. For each sample run, at least 106 events were acquired, with VSELs accounting for 0.02-0.03% of the total cells (BM or PB). Given the scarcity of VSELs in peripheral blood and bone marrow, we performed six independent measurements per exposure condition (hypoxia, normoxia), each based on pooled samples from 10 mice, for a total of 120 animals.
Measurement of plasma chemokines
Mouse blood was collected from the vena cava in vacutainer tubes containing EDTA (Becton Dickinson, Franklin Lakes, NJ, USA). The collected fresh blood was centrifuged at 2000 × g for 20 min at 4°C, subsequently plasma was centrifuged for 5 minutes at 13,000 rpm to remove remaining cells and platelets, and immediately frozen at -80°C until further analysis. Levels of stem cell chemokines stromal cell derived factor-1 (SDF-1), hepatocyte growth factor (HGF), and leukemia inhibitory factor (LIF) were measured with ELISA kits (R&D Systems, Minneapolis, MN) . For the hypoxic exposure, 20 independent measurements and for the normoxic exposure, 30 independent measurements were performed. For each of these measurements, plasma was pooled from 3-4 mice.
Chemokine levels and VSEL counts between hypoxia and normoxia were performed using unpaired, two-tailed Student's t-test with P-values adjusted for unequal variances when appropriate. Box and whisker representation of the data was used to depict the upper and lower quartiles (box) and the maximum and minimum values (whiskers).
VSELs were isolated from the bone marrows of 30 additional mice exposed to hypoxia and 30 mice exposed to room air (pooled into six groups of n = 10). Total RNA was isolated using PicoPure RNA Isolation Kit (Arcturus Bioscience Inc., Mountain View, CA) and amplified using Low-RNA Input Fluorescent Linear Amplification kit (Agilent Technologies, Santa Clara, CA) with modifications. For each of the pooled groups (n = 10 mice per group), labeled cRNA was hybridized to an Agilent murine whole-genome 60-mer oligo microarray (Agilent Technologies, Santa Clara, CA) comprised of over 45,000 probes as previously described . Six independent hybridizations-3 pooled samples from mice exposed to hypoxia and 3 pooled samples from normoxic mice, were performed. After image processing, background-subtracted intensities were normalized using the quantile method .
Microarray data analysis
Two-dimensional hierarchal clustering using Pearson's correlation was performed on the entire normalized gene expression data . Differentially expressed genes in VSELs exposed to hypoxia versus normoxia were identified using a Bayesian implementation of the parametric t-test  and corrected for multiple comparisons with the Q-value method . A Q-value of < 0.05 was selected for identifying significant differential gene expression.
Gene ontology analysis
Differentially expressed genes in VSELs (Q-value < 0.05) underwent functional analysis with a web-based program, Database for Annotation, Visualization and Integrated Discovery (DAVID) and categorized using the Gene Ontology database . A permutation-based FDR analysis was employed (FDR cutoff < 5%).
By including all differentially expressed genes during hypoxia as initial seeds, we developed a gene product relational network based on published interaction databases [29–31]. To maximize the accuracy of this "interactome", we limited the connection between any two nodes to known direct relationships (e.g., protein-protein interaction). The topological architecture of the network was derived from its connectivity matrix. We tested the "scale-free" property of the network by determining if it followed a power law distribution . We identified the key network hubs by selecting differentially expressed nodes with at least 10 connections.
Effects of in vivoexposure to hypoxia on VSEL populations and stem cell chemokines
We utilized a sequential gating strategy using FACS and based on size, granularity and stem cell-specific antibody staining to purify VSELs from PB and BM of mice exposed to hypoxia and normoxic controls (Figure 1). In peripheral blood, the white blood cell counts were 9,980 ± 467/ml in normoxia and 10,503 ± 420/ml in hypoxia (p > 0.05); a minimum of 100,000 cells were counted from merging the samples from 10 mice to determine the frequency of VSEL in peripheral blood for each independent sample run. For bone marrow, cellularity was estimated from the mixing of 10 mice to include approximately 1 million cells for each independent sample run. Given the scarcity of VSELs in peripheral blood and bone marrow, we performed six independent measurements per exposure condition (hypoxia, normoxia), each based on pooled samples from 10 mice, for a total of 120 animals.
Exposure to hypoxia resulted in a significant increase in the peripheral count of VSELs relative to their BM reserve pool (Figure 2A-C). Furthermore, in vivo exposure to hypoxia was associated with a significant elevation in plasma levels of the stem cell chemoattractants SDF-1 and LIF, and borderline increase in HGF (Figure 2D-F), providing further evidence that hypoxia promotes a favorable gradient for mobilization of VSELs from BM to PB.
Transcriptional profiling of VSELs during hypoxia
Global patterns in gene expression of BM-derived VSELs exposed to hypoxia and normoxic controls were dissected using hierarchical cluster analysis and displayed as a heatmap (Figure 3). This analysis revealed that exposure to hypoxia causes a perturbation in the entire VSEL transcriptome that is distinct from the normoxic condition. Within this genome-wide expression pattern, we identified 1,388 differentially expressed genes using a Q-value cutoff of < 0.05 (Additional file 1, Table S1). The subsequent pathway and network analyses were limited to this subset of statistically significant differentially expressed genes.
Functional and network analyses of hypoxia-induced transcriptional programs in VSELs
Functional enrichment analysis of differentially expressed genes was performed using Gene Ontology analysis  and graphically displayed based the hierarchical structure of this annotation database (Figure 4). In vivo exposure to hypoxia activated a diverse array of transcriptional programs in pluripotent VSELs ranging from a number of organ-specific developmental processes (e.g., lung, kidney, skeletal, vascular, and nervous system development) to pathways involved in defense and stress response, cytokine and growth factor binding, and catalytic activity. These results provide a blueprint mapping the profound and widespread functional activation of BM-derived VSELs in response to a brief in vivo hypoxic exposure.
Given the vast transcriptional effects of hypoxia on VSELs, it is highly unlikely that a single gene will solely orchestrate this response. To better elucidate the relationships among differentially expressed genes mapped to the enriched functional groups and identify the key controllers of hypoxia-induced transcriptional programs in VSELs, we developed a gene product interaction network, or interactome, based on previously published direct interactions (Additional file 2, Figure S1). This network consisted of 424 genes (nodes) and 604 connections (edges). Topologically, the interactome possessed scale-free properties  because its degree distribution (number of nodes with a given connection, N k ) and nodal connectivity (k) were linked via a power law relationship-i.e., N k ~ k-γ (γ = 1.64, R2 = 0.98) (Additional file 2, Figure S1). We and others have demonstrated that the functional stability of complex biological networks is critically dependent on highly connected nodes or "hubs" [33, 34]. We identified the VSEL network hubs by locating 18 differentially expressed nodes with at least 10 interactions (Table 1). As shown in Figure 5A, the connectivity pattern of these hubs captures a significant proportion of the overall network structure (313 edges out of a total of 604). Furthermore, the key nodes themselves are highly interconnected implying coordinated, direct functional relationships between these hubs (Figure 5B).
The main finding of this work is that exposure to hypoxia-a common characteristic of many lung disorders, recruits pluripotent VSELs from BM to the peripheral blood and induces a distinct gene expression signature in these stem cells. Using a computationally-intensive approach, we mapped the functional and network properties of activated transcriptional programs and identified hubs that may serve as key regulators of the VSEL response to hypoxia.
To our knowledge, this is the first report on the global transcriptional effects of hypoxia on stem cells using in vivo exposure at physiologically relevant oxygen partial pressures. Utilizing a similar strategy, we recently used a murine model to investigate the consequences of cyclical hypoxia during sleep on BM-derived VSELs . The intermittent hypoxic exposure was intended to model obstructive sleep apnea and that project was therefore more narrowly defined in its focus and clinical relevance. Nevertheless, several of our current findings were also observed in the intermittent hypoxia system, including mobilization of VSELs from BM to PB, elevation of stem cell chemoattractant chemokines, and enrichment of transcriptional programs involved in development and organogenesis. Network analysis of the differentially expressed genes in VSELs also revealed that some of the key hubs were shared between the hypoxia models (intermittent, sustained), including hypoxia inducible factor 2 alpha (Hif2α), proliferator-activated receptor gamma (Pparγ), integrin beta 3 (Itgb3), and insulin growth factor binding protein 5 (Igfbp5). This observation implies that global exposure to reduced oxygen, whether continuous or intermittent, activates a number of common pathways in pluripotent stem cells. However, there were also significant differences in the transcriptional programs enriched in VSELs under sustained vs. intermittent hypoxic exposure. Firstly, a larger number of genes were differentially expressed during continuous hypoxia suggesting a more profound perturbation. Secondly, a much broader spectrum of pathways was enriched in sustained compared to intermittent hypoxia that extended beyond developmental programs and included processes involved in defense and stress response, wound healing, chemokine activity, cell adhesion and structure (Figure 4). Thirdly, many of the network hubs identified during exposure to continuous hypoxia were not differentially expressed during cyclical hypoxia (Figure 5, Table 1), including the most densely connected node, Fos (FBJ osteosarcoma oncogene), and several other key transcriptional regulators such as early growth response 1 and 2 (Egr1, Egr2), and TATA box binding protein (Tbp).
Fos and Egr1 are prototypic master regulators known as immediate early-genes (IEG) that orchestrate widespread activation of cellular gene expression in response to specific perturbations, including hypoxia, ischemia, atherosclerosis, angiogenesis, and neuronal survival [35–37]. Egr2 has been shown to govern the development and segmentation of the mouse hindbrain , and mutations in this gene have been associated with congenital hypomyelinating neuropathies in humans . Intriguingly, a recent study demonstrated that a population of Egr2-expressing precursor cells develop into the first embryonic respiratory rhythm-generating neuronal circuit , and Egr2-/- transgenic mice die from breathing irregularities resulting in severe hypoxia .
Two members of the heat shock protein 70 group, Hspa8 and Hspa5, were also densely connected network nodes whose expression increased significantly during hypoxia. Products of these genes play crucial roles in protecting cells from environmental stressors including high temperature, oxidative stress, ischemia and hypoxia . More specifically, targeting Hspa8 and Hspa5 has been demonstrated to protect cardiomyocytes  and mesenchymal stem cells from hypoxia-induced apoptosis , whereas delivering these heat shock proteins into transplanted mesenchymal stem cells rescues heart function after myocardial infarction in rats .
Another differentially upregulated hub within the interactome, Hif2α, is a member of the HIF family of basic helix-loop-helix transcription factors that control the global cellular response to oxygen deprivation. Hif2α overexpression is seen many solid tumors characterized with significant hypoxic burden including renal cell carcinoma, non small cell lung cancer and meduloblastoma, and is associated with poor outcome in these cancers . Recent reports suggest that Hif2α is a critical activator of a subpopulation of cancer cells with stem cell-like properties [46, 47], and pharmacologic and genetic targeting of this gene affects differentiation of progenitor stem cells [48, 49]. Interestingly, a single nucleotide polymorphism in Hif2α was recently demonstrated to be highly associated with adaptability to high altitude hypoxia in ethnic Tibetans .
The mobilization of VSELs from BM to peripheral circulation in response to oxygen deprivation implies the activation of organ-specific transcriptional programs under regulatory controls. Our network analysis provided an overview of the transcriptional landscape of these pluripotent stem cells during hypoxia. Since the functional integrity of such scale-free biologic networks is dependent on densely connected nodes , we postulate that these hubs comprise key regulators of hypoxia-induced gene expression in VSELs. As depicted in Figure 5, the hubs not only form the structural foundation of the interactome, but are themselves highly interconnected via direct interactions, implying a high degree of functional coordination in orchestrating the hypoxia-induced transcriptional response of stem cells.
To explore whether this hypoxia-induced transcriptional signature is unique to VSELs or common among leukocyte subpopulations, we compared our results to expression profiles of mononuclear cells (MNCs) and mesenchymal stem cells (MSCs) under hypoxia as reported by Onishi et al . While there are substantial differences between the two studies in terms of experimental protocols, hypoxic exposures, and microarray platforms, less than 1% of the differentially expressed genes in VSELs were significantly altered in MSCs or MNCs, suggesting that many of the hypoxia-induced genes in VSELs are exclusive to this subpopulation of cells. However, from a functional standpoint, differentially expressed genes in VSELs and MSCs were both enriched in several developmental processes, whereas MNCs were highly enriched in defense/immune responses and inflammation. Therefore, VSELs appear to possess transcriptional programs in common with both MSCs and MNCs-a finding that is consistent with their pluripotent properties.
Several limitations in the current study deserve mention. We have only explored the effects of hypoxia on gene expression and not assessed the post-transcriptional modifications of gene products activated during hypoxic exposure. The duration of hypoxia was short-the long term consequences of chronic oxygen deprivation on VSELs were not investigated in this study and represent a future research goal. Although bone marrow is the primary depot of VSELs, our results do not preclude the possibility that some of these stem cells were also mobilized from other tissue niches during exposure to hypoxia. Nevertheless, our expression profiling experiments were focused on BM-derived VSELs. We have not identified the tissue source of hypoxia-induced elevation in stem cells chemokines, but given the extremely low number of circulating VSELs, we believe it is unlikely that VSELs are the primary source. Although the bone marrow is a known source of these chemoattractants, hypoxic injury to various tissues can lead to the release of the chemokines into peripheral blood [5, 14, 16]. Finally, we did not follow the fate of mobilized VSELs to identify their target end organs, although given the global effects of hypoxia and the diversity of enriched developmental programs, we expect that these stem cells will be recruited to multiple sites and undergo tissue-specific differentiation.
Collectively, our results demonstrate that in vivo exposure to hypoxia in mice elicits chemoattractant gradients that promote the mobilization of pluripotent very small embryonic-like stem cells from the bone marrow to peripheral blood. A computationally-driven analysis of the transcriptome of these cells highlighted the activation of diverse transcriptional programs that appeared to be coordinated by a select number of highly connected network nodes. These hubs may represent critical targets for modulating the functional consequences of hypoxic burden on stem cells. Our approach provides a general framework for the systematic study of stem cell biology under clinically relevant pathophysiologic perturbations.
Semenza GL: Hypoxia and human disease-and the Journal of Molecular Medicine. J Mol Med. 2007, 85 (12): 1293-1294. 10.1007/s00109-007-0285-z.
Simon MC, Keith B: The role of oxygen availability in embryonic development and stem cell function. Nat Rev Mol Cell Biol. 2008, 9 (4): 285-296. 10.1038/nrm2354.
Mohyeldin A, Garzon-Muvdi T, Quinones-Hinojosa A: Oxygen in stem cell biology: a critical component of the stem cell niche. Cell Stem Cell. 2010, 7 (2): 150-161. 10.1016/j.stem.2010.07.007.
Kucia MJ, Wysoczynski M, Wu W, Zuba-Surma EK, Ratajczak J, Ratajczak MZ: Evidence that very small embryonic-like stem cells are mobilized into peripheral blood. Stem Cells. 2008, 26 (8): 2083-2092. 10.1634/stemcells.2007-0922.
Kucia M, Zhang YP, Reca R, Wysoczynski M, Machalinski B, Majka M, Ildstad ST, Ratajczak J, Shields CB, Ratajczak MZ: Cells enriched in markers of neural tissue-committed stem cells reside in the bone marrow and are mobilized into the peripheral blood following stroke. Leukemia. 2006, 20 (1): 18-28. 10.1038/sj.leu.2404011.
Schachinger V, Erbs S, Elsasser A, Haberbosch W, Hambrecht R, Holschermann H, Yu J, Corti R, Mathey DG, Hamm CW, et al: Intracoronary bone marrow-derived progenitor cells in acute myocardial infarction. N Engl J Med. 2006, 355 (12): 1210-1221. 10.1056/NEJMoa060186.
Burnham EL, Taylor WR, Quyyumi AA, Rojas M, Brigham KL, Moss M: Increased circulating endothelial progenitor cells are associated with survival in acute lung injury. Am J Respir Crit Care Med. 2005, 172 (7): 854-860. 10.1164/rccm.200410-1325OC.
Fadini GP, Schiavon M, Cantini M, Baesso I, Facco M, Miorin M, Tassinato M, de Kreutzenberg SV, Avogaro A, Agostini C: Circulating progenitor cells are reduced in patients with severe lung disease. Stem Cells. 2006, 24 (7): 1806-1813. 10.1634/stemcells.2005-0440.
Fadini GP, Schiavon M, Rea F, Avogaro A, Agostini C: Depletion of endothelial progenitor cells may link pulmonary fibrosis and pulmonary hypertension. Am J Respir Crit Care Med. 2007, 176 (7): 724-725. author reply 725
Palange P, Testa U, Huertas A, Calabro L, Antonucci R, Petrucci E, Pelosi E, Pasquini L, Satta A, Morici G, et al: Circulating haemopoietic and endothelial progenitor cells are decreased in COPD. Eur Respir J. 2006, 27 (3): 529-541. 10.1183/09031936.06.00120604.
Kucia M, Wysoczynski M, Ratajczak J, Ratajczak MZ: Identification of very small embryonic like (VSEL) stem cells in bone marrow. Cell Tissue Res. 2008, 331 (1): 125-134. 10.1007/s00441-007-0485-4.
Kucia M, Reca R, Campbell FR, Zuba-Surma E, Majka M, Ratajczak J, Ratajczak MZ: A population of very small embryonic-like (VSEL) CXCR4(+)SSEA-1(+)Oct-4+ stem cells identified in adult bone marrow. Leukemia. 2006, 20 (5): 857-869. 10.1038/sj.leu.2404171.
Kucia M, Wojakowski W, Reca R, Machalinski B, Gozdzik J, Majka M, Baran J, Ratajczak J, Ratajczak MZ: The migration of bone marrow-derived non-hematopoietic tissue-committed stem cells is regulated in an SDF-1-, HGF-, and LIF-dependent manner. Arch Immunol Ther Exp (Warsz). 2006, 54 (2): 121-135. 10.1007/s00005-006-0015-1.
Ceradini DJ, Kulkarni AR, Callaghan MJ, Tepper OM, Bastidas N, Kleinman ME, Capla JM, Galiano RD, Levine JP, Gurtner GC: Progenitor cell trafficking is regulated by hypoxic gradients through HIF-1 induction of SDF-1. Nat Med. 2004, 10 (8): 858-864. 10.1038/nm1075.
Paczkowska E, Kucia M, Koziarska D, Halasa M, Safranow K, Masiuk M, Karbicka A, Nowik M, Nowacki P, Ratajczak MZ, et al: Clinical evidence that very small embryonic-like stem cells are mobilized into peripheral blood in patients after stroke. Stroke. 2009, 40 (4): 1237-1244. 10.1161/STROKEAHA.108.535062.
Kucia M, Dawn B, Hunt G, Guo Y, Wysoczynski M, Majka M, Ratajczak J, Rezzoug F, Ildstad ST, Bolli R, et al: Cells expressing early cardiac markers reside in the bone marrow and are mobilized into the peripheral blood after myocardial infarction. Circ Res. 2004, 95 (12): 1191-1199. 10.1161/01.RES.0000150856.47324.5b.
Chi JT, Wang Z, Nuyten DS, Rodriguez EH, Schaner ME, Salim A, Wang Y, Kristensen GB, Helland A, Borresen-Dale AL, et al: Gene expression programs in response to hypoxia: cell type specificity and prognostic significance in human cancers. PLoS Med. 2006, 3 (3): e47-10.1371/journal.pmed.0030047.
Wang GL, Jiang BH, Rue EA, Semenza GL: Hypoxia-inducible factor 1 is a basic-helix-loop-helix-PAS heterodimer regulated by cellular O2 tension. Proc Natl Acad Sci USA. 1995, 92 (12): 5510-5514. 10.1073/pnas.92.12.5510.
Menicanin D, Bartold PM, Zannettino AC, Gronthos S: Genomic profiling of mesenchymal stem cells. Stem Cell Rev. 2009, 5 (1): 36-50. 10.1007/s12015-009-9056-2.
Westfall SD, Sachdev S, Das P, Hearne LB, Hannink M, Roberts RM, Ezashi T: Identification of oxygen-sensitive transcriptional programs in human embryonic stem cells. Stem Cells Dev. 2008, 17 (5): 869-881. 10.1089/scd.2007.0240.
Ivanovic Z: Hypoxia or in situ normoxia: The stem cell paradigm. J Cell Physiol. 2009, 219 (2): 271-275. 10.1002/jcp.21690.
Millman JR, Tan JH, Colton CK: The effects of low oxygen on self-renewal and differentiation of embryonic stem cells. Curr Opin Organ Transplant. 2009, 14 (6): 694-700. 10.1097/MOT.0b013e3283329d53.
Gharib SA, Dayyat E, Khalyfa A, Kim J, Clair HB, Kucia M, Gozal D: Intermittent Hypoxia Mobilizes Bone Marrow-Derived Very Small Embryonic-Like Stem Cells and Activates Developmental Transcriptional Programs in Mice. Sleep. 2010, 33 (11): 1439-1446.
Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19 (2): 185-193. 10.1093/bioinformatics/19.2.185.
Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier T, Thiagarajan M, et al: TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003, 34 (2): 374-378.
Baldi P, Long AD: A Bayesian framework for the analysis of microarray expression data: regularized t-test and statistical inferences of gene changes. Bioinformatics. 2001, 17 (6): 509-519. 10.1093/bioinformatics/17.6.509.
Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci USA. 2003, 100 (16): 9440-9445. 10.1073/pnas.1530509100.
Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4 (5): P3-10.1186/gb-2003-4-5-p3.
Aranda B, Achuthan P, Alam-Faruque Y, Armean I, Bridge A, Derow C, Feuermann M, Ghanbarian AT, Kerrien S, Khadake J, et al: The IntAct molecular interaction database in 2010. Nucleic Acids Res. 2009, D525-531. 38 Database
Ceol A, Chatr Aryamontri A, Licata L, Peluso D, Briganti L, Perfetto L, Castagnoli L, Cesareni G: MINT, the molecular interaction database: 2009 update. Nucleic Acids Res. 2009, D532-539. 38 Database
Calvano SE, Xiao W, Richards DR, Felciano RM, Baker HV, Cho RJ, Chen RO, Brownstein BH, Cobb JP, Tschoeke SK, et al: A network-based analysis of systemic inflammation in humans. Nature. 2005, 437 (7061): 1032-1037. 10.1038/nature03985.
Barabasi AL, Albert R: Emergence of scaling in random networks. Science. 1999, 286 (5439): 509-512. 10.1126/science.286.5439.509.
Luscombe NM, Babu MM, Yu H, Snyder M, Teichmann SA, Gerstein M: Genomic analysis of regulatory network dynamics reveals large topological changes. Nature. 2004, 431 (7006): 308-312. 10.1038/nature02782.
Khalyfa A, Gharib SA, Kim J, Dayyat E, Snow AB, Bhattacharjee R, Kheirandish-Gozal L, Goldman JL, Gozal D: Transcriptomic analysis identifies phosphatases as novel targets for adenotonsillar hypertrophy of pediatric obstructive sleep apnea. Am J Respir Crit Care Med. 2010, 181 (10): 1114-1120. 10.1164/rccm.200909-1398OC.
Rybnikova E, Glushchenko T, Tyulkova E, Baranova K, Samoilov M: Mild hypobaric hypoxia preconditioning up-regulates expression of transcription factors c-Fos and NGFI-A in rat neocortex and hippocampus. Neurosci Res. 2009, 65 (4): 360-366. 10.1016/j.neures.2009.08.013.
Khachigian LM: Early growth response-1 in cardiovascular pathobiology. Circ Res. 2006, 98 (2): 186-191. 10.1161/01.RES.0000200177.53882.c3.
Zhang J, Zhang D, McQuade JS, Behbehani M, Tsien JZ, Xu M: c-fos regulates neuronal excitability and survival. Nat Genet. 2002, 30 (4): 416-420. 10.1038/ng859.
Schneider-Maunoury S, Topilko P, Seitandou T, Levi G, Cohen-Tannoudji M, Pournin S, Babinet C, Charnay P: Disruption of Krox-20 results in alteration of rhombomeres 3 and 5 in the developing hindbrain. Cell. 1993, 75 (6): 1199-1214. 10.1016/0092-8674(93)90329-O.
Warner LE, Mancias P, Butler IJ, McDonald CM, Keppen L, Koob KG, Lupski JR: Mutations in the early growth response 2 (EGR2) gene are associated with hereditary myelinopathies. Nat Genet. 1998, 18 (4): 382-384. 10.1038/ng0498-382.
Thoby-Brisson M, Karlen M, Wu N, Charnay P, Champagnat J, Fortin G: Genetic identification of an embryonic parafacial oscillator coupling to the preBotzinger complex. Nat Neurosci. 2009, 12 (8): 1028-1035. 10.1038/nn.2354.
Jacquin TD, Borday V, Schneider-Maunoury S, Topilko P, Ghilini G, Kato F, Charnay P, Champagnat J: Reorganization of pontine rhythmogenic neuronal networks in Krox-20 knockout mice. Neuron. 1996, 17 (4): 747-758. 10.1016/S0896-6273(00)80206-8.
Daugaard M, Rohde M, Jaattela M: The heat shock protein 70 family: Highly homologous proteins with overlapping and distinct functions. FEBS Lett. 2007, 581 (19): 3702-3710. 10.1016/j.febslet.2007.05.039.
Hardy B, Raiter A: Peptide-binding heat shock protein GRP78 protects cardiomyocytes from hypoxia-induced apoptosis. J Mol Med. 2010
Chang W, Song BW, Lim S, Song H, Shim CY, Cha MJ, Ahn DH, Jung YG, Lee DH, Chung JH, et al: Mesenchymal stem cells pretreated with delivered Hph-1-Hsp70 protein are protected from hypoxia-mediated cell death and rescue heart functions from myocardial injury. Stem Cells. 2009, 27 (9): 2283-2292. 10.1002/stem.153.
Qing G, Simon MC: Hypoxia inducible factor-2alpha: a critical mediator of aggressive tumor phenotypes. Curr Opin Genet Dev. 2009, 19 (1): 60-66. 10.1016/j.gde.2008.12.001.
Li Z, Bao S, Wu Q, Wang H, Eyler C, Sathornsumetee S, Shi Q, Cao Y, Lathia J, McLendon RE, et al: Hypoxia-inducible factors regulate tumorigenic capacity of glioma stem cells. Cancer Cell. 2009, 15 (6): 501-513. 10.1016/j.ccr.2009.03.018.
Seidel S, Garvalov BK, Wirta V, von Stechow L, Schanzer A, Meletis K, Wolter M, Sommerlad D, Henze AT, Nister M, et al: A hypoxic niche regulates glioblastoma stem cells through hypoxia inducible factor 2 alpha. Brain. 2010, 133 (Pt 4): 983-995.
Moreno-Manzano V, Rodriguez-Jimenez FJ, Acena-Bonilla JL, Fustero-Lardies S, Erceg S, Dopazo J, Montaner D, Stojkovic M, Sanchez-Puelles JM: FM19G11, a new hypoxia-inducible factor (HIF) modulator, affects stem cell differentiation status. J Biol Chem. 2010, 285 (2): 1333-1342. 10.1074/jbc.M109.008326.
Covello KL, Kehler J, Yu H, Gordan JD, Arsham AM, Hu CJ, Labosky PA, Simon MC, Keith B: HIF-2alpha regulates Oct-4: effects of hypoxia on stem cell function, embryonic development, and tumor growth. Genes Dev. 2006, 20 (5): 557-570. 10.1101/gad.1399906.
Yi X, Liang Y, Huerta-Sanchez E, Jin X, Cuo ZX, Pool JE, Xu X, Jiang H, Vinckenbosch N, Korneliussen TS, et al: Sequencing of 50 human exomes reveals adaptation to high altitude. Science. 2010, 329 (5987): 75-78. 10.1126/science.1190371.
Ohnishi S, Yasuda T, Kitamura S, Nagaya N: Effect of hypoxia on gene expression of bone marrow-derived mesenchymal stem cells and mononuclear cells. Stem Cells. 2007, 25 (5): 1166-1177. 10.1634/stemcells.2006-0347.
This work was supported in part by the National Institutes of Health HL065270 and HL086662 (DG), and American Sleep Medicine Foundation Junior Faculty Research Award (SAG).
The authors declare that they have no competing interests.
SAG analyzed the data and wrote the manuscript. DG conceived the study, participated in its design and draft. AK participated in designing the study and collected the stem cell data and performed the transcriptional profiling experiments. MK assisted with FACS-based stem cell isolation and collection of stem cell data. EAD, JK, and HBC assisted with the stem cell isolation, animal experiments, and stem cell data measurements. All authors read and approved the final manuscript.
Sina A Gharib, Abdelnaby Khalyfa contributed equally to this work.
Electronic supplementary material
Additional file 1: Tabular list of differentially expressed genes in VSELs after hypoxic exposure. This file contains the list of differentially expressed genes (based on a Q-value cutoff less than 0.05) in VSELs exposed to hypoxia. The list includes the official gene symbols, Entrez Gene IDs, Q-values, and Log2[expression in hypoxia/expression in normoxia]. (PDF 291 KB)
Additional file 2: Gene product interaction network of differentially expressed genes in VSELs. This file contains two figures depicting a gene product interaction network of differentially expressed genes in VSELs after in vivo exposure to hypoxia. Panel (A) highlights up and downregulated genes (red and green respectively). This network has 424 genes (nodes) and 604 connections (edges). Note that each connection represents a known direct interaction between two gene products. Panel (B) demonstrates that the topology of this network is scale-free because it follows a power law distribution. N k , degree distribution; k, nodal connectivity. (PDF 2 MB)
About this article
Cite this article
Gharib, S.A., Khalyfa, A., Kucia, M.J. et al. Transcriptional landscape of bone marrow-derived very small embryonic-like stem cells during hypoxia. Respir Res 12, 63 (2011) doi:10.1186/1465-9921-12-63
- Obstructive Sleep Apnea
- Hepatocyte Growth Factor
- Pluripotent Stem Cell
- Leukemia Inhibitory Factor
- Intermittent Hypoxia