Skip to main content
  • Research article
  • Open access
  • Published:

Severe hypoxia exerts parallel and cell-specific regulation of gene expression and alternative splicing in human mesenchymal stem cells

Abstract

Background

The endosteum of the bone marrow provides a specialized hypoxic niche that may serve to preserve the integrity, pluripotency, longevity and stemness of resident mesenchymal stem cells (MSCs). To explore the molecular genetic consequences of such a niche we subjected human (h) MSCs to a pO2 of 4 mmHg and analyzed global gene expression and alternative splicing (AS) by genome-exon microarray and RT-qPCR, and phenotype by western blot and immunostaining.

Results

Out of 446 genes differentially regulated by >2.5-fold, down-regulated genes outnumbered up-regulated genes by 243:203. Exon analyses revealed 60 hypoxia-regulated AS events with splice indices (SI) >1.0 from 53 genes and a correlation between high SI and degree of transcript regulation. Parallel analyses of a publicly available AS study on human umbilical vein endothelial cells (HUVECs) showed that there was a strong cell-specific component with only 11 genes commonly regulated in hMSCs and HUVECs and 17 common differentially spliced genes. Only 3 genes were differentially responsive to hypoxia at the gene (>2.0) and AS levels in both cell types. Functional assignments revealed unique profiles of gene expression with complex regulation of differentiation, extracellular matrix, intermediate filament and metabolic marker genes. Antioxidant genes, striated muscle genes and insulin/IGF-1 signaling intermediates were down-regulated. There was a coordinate induction of 9 out of 12 acidic keratins that along with other epithelial and cell adhesion markers implies a partial mesenchymal to epithelial transition.

Conclusions

We conclude that severe hypoxia confers a quiescent phenotype in hMSCs that is reflected by both the transcriptome profile and gene-specific changes of splicosome actions. The results reveal that severe hypoxia imposes markedly different patterns of gene regulation of MSCs compared with more moderate hypoxia. This is the first study to report hypoxia-regulation of AS in stem/progenitor cells and the first molecular genetic characterization of MSC in a hypoxia-induced quiescent immobile state.

Background

The stem cell niche refers to a well-defined physiological compartment that includes cellular and acellular components and serves to integrate systemic and local signals to regulate the biology of stem cells (reviewed in [1, 2]). Like other such niches, the bone marrow provides highly specialized and heterogeneous microenvironments that determine the self-renewal, multipotency, survival and migration of residing hematopoietic and progenitor cells including mesenchymal stem cells (MSCs). Recently oxygen tension (hypoxia) has been recognized as an important component of stem cell niches that exerts control over the proliferation, differentiation and pluripotency of resident cells [35]. The oxygen tension of the endosteum, a narrow compartment of the bone marrow directly adjacent to the bone, is less than 10 mmHg, while that of the sinusoidal cavity ranges between 30–60 mmHg; therefore bone marrow cells are subject to a gradient of hypoxia the severity of which depends on their location within the niche [3, 6]. Studies of embryonic stem cells (ESCs) as well as induced pluripotential stem cells (iPSC) indicate that oxygen gradients control stem cell functions. Culture of ESCs under an aerobic pO2 of 160 mmHg causes spontaneous differentiation that is suppressed by more physiological pO2 within the range of 14–36 mmHg. Further reduction of pO2 to <10 mmHg, equivalent to the bone marrow endosteum suppresses both differentiation and proliferation of ESCs while retaining pluripotency [7, 8]. Studies on bone marrow or adipose derived MSCs have shown similarly that moderate hypoxic culture equivalent to the central BM sinusoidal niche enhances proliferation and protects against senescence while more severe hypoxia may block proliferation and induce cell death [916].

Culture of MSCs under moderate hypoxia has been shown to modulate gene expression by HIF-1/2-dependent and independent mechanisms [8, 1721]. In addition to the predicted HIF-1α target genes such as those required for anaerobic metabolism (glycolytic enzymes, glucose transporters), cell cycle (p21, p53), and angiogenesis (VEGF), moderate hypoxia was shown to mediate increased expression of Oct4 and telomerase activity of human bone marrow MSCs [16, 22, 23]. When cultured under moderate hypoxia (20–40 mmHg O2), MSCs display enhanced proliferation and migratory activity that has been attributed to increased Akt phosphorylation, expression of c-MET, VEGF, chemokine receptors CXCR4 and CXCR1, and increased phosphorylation of focal adhesion kinase [10, 18, 19, 2426]. Suppression of stem cell differentiation by hypoxia has been linked to Notch pathway signaling wherein hypoxia promotes recruitment of HIF-1α to the Notch intracellular domain and subsequently to Notch-dependent promoters thereby enhancing their expression [27]. Recently, HIF-1α was shown to regulate MSC proliferation through the enhancement of TWIST expression, which down-regulates the E2A-p21 pathway, inhibits senescence and increases proliferation [28]. Oxygen tensions between 20 and 40 mmHg enhance proliferation and pluripotency of stem cells whereas tension below 10 mmHg (<1%) inhibit proliferation and may promote apoptosis [812, 1719].

Up to 95% of all human genes are alternatively spliced [29, 30]. AS results in changes in composition of an mRNA produced from a given gene, brought about by changes in splice site choice and thence the production of proteins with different properties. AS is regulated by cell specific, developmental, and extracellular signal-regulated cues and pathways (Reviewed in [31]). Aberrant alternative splicing can cause disease and may contribute to cancer and neurodegenerative disease [32, 33]. Exon arrays allow for detection and quantification of AS on a genome-wide scale. There are currently only 2 such reports of genome-wide analyses of hypoxia-related changes in pre-mRNA splicing. One identified Lama3 as a hypoxia-related splice variant in head and neck cancers [34]. Another analyzed the effects of hypoxia on AS in human umbilical vein endothelial cells (HUVECs) and identified multiple alternative splice events [35].

Here we investigated the effects of severe hypoxia on gene expression, exon splicing, and phenotype of human (h) MSCs. The results reveal for the first time unique sets of (severe) hypoxia-activated and repressed genes, many of which differ from those reported previously for more moderate hypoxia. We report for the first time a coordinate increase in expression of acidic keratins perhaps indicating a partial mesenchymal to epithelial (epidermal) transition (MET), a decrease in insulin/IGF-1 signaling with lower phosphor-Akt, and decreased expression of anti-oxidant-related genes that suggests lower metabolism and growth compared with aerobic culture. The expression of differentiation-related markers is consistent with enhancement of osteogenic and angiogenic pathways perhaps at the expense of myogenesis and adipogenesis. We also identify a novel set of hypoxia-regulated alternatively spliced transcripts in hMSCs. To our knowledge this is the first study to report on patterns of hypoxia-mediated alternative splicing in stem cells. The results provide a molecular framework for understanding the role of severe hypoxia in preserving bone marrow progenitor cell integrity and perhaps insights into the role of hypoxia in regulating cell biology in hypoxic niche environments such as the endosteum.

Results

Isolation and characterization of human MSCs

Human bone marrow MSCs were isolated as described in Methods and used at passage 8. At this time cells were visually homogeneous, fibroblast-like and positive for the expression of mesenchymal-specific markers CD29 (99.0%) and CD166 (41.8%) and negative for the expression of hematopoietic lineage marker CD34 (0.4%). These characteristics were unaltered after exposure to hypoxia for 24 h (Figure 1A and 1B).

Figure 1
figure 1

Characterization of hMSCs cultures under air or hypoxia. (A) Morphology and (B) surface antigen profiling of normoxic (N-MSC) and hypoxic hMSCs (H-MSC). Conditions and procedures are described in Methods. Representative of n = 3.

Gene expression array

Gene expression profiles comparing normoxia and hypoxia were obtained using Agilent Human 4 × 180 K Exon and 8 × 60 K-GE microarrays as described in Methods. Only the Exon arrays are described in the present analysis and the GE arrays were used for confirmation of some gene transcripts. Hierarchical clustering of the Exon arrays confirmed high reproducibility between samples (Additional file 1: Figure S1). A robust response to hypoxia was confirmed by quantifying HIF-1α-regulated transcripts. As shown in Table 1, multiple well-characterized HIF-1α-regulated genes were represented including carbonic anhydrase (>5-fold), metallothionein (>4-fold) and VEGF-A (>4-fold). Most of these genes have been reported previously in similar high throughput analyses of MSCs exposed to hypoxia in the range 1-5% [14, 1719]. Noteworthy in our analyses are the strong inductions of leptin and insulin-like growth factor binding protein 1 transcripts, confirmed in both Exon and K-GE arrays (latter data not shown) and relatively low induction of Bnip3 and glucose transporters. Glycolytic enzyme genes are widely recognized as markers of hypoxia with at least 8 of 11 glycolytic enzymes genes responsive through the HIF-1α pathway (reviewed in [36]). Surprisingly, glycolytic gene transcripts were not represented in the 2.5-fold cut-off groupings used for our array analyses, therefore we used RT-PCR to quantify phosphofructokinase (PFK) a major rate-limiting enzyme [37], phosphoglycerate kinase (PGK) a strong hypoxia responder [38], phosphoglycerate mutase (PGAM), a moderate responder [39] and glyceraldehyde-3-phosphate dehydrogenase (GAPDH), that is regulated by hypoxia in a tissue-specific manner [4042]. As indicated in the Table, transcript levels of PFK and PGK increased, there was no significant change of GAPDH (p = 0.20, n = 7) and PGAM transcripts decreased (p = 0.02, n = 7). The results suggest mixed responses of individual genes. Glycolytic flux increased significantly under hypoxia; aerobic cultures used glucose at a rate of 2.83 ± 0.2 mg/dL/h compared with 3.33 ± 0.1 mg/dL/h for hypoxia (p > 0.05; n = 3), and there was a parallel increased rate of lactic acid generation over 24 h (data not show). Glycolytic flux is largely regulated by small molecule binding and allosteric control of PFK, the main rate-limiting enzyme and most of the pathway enzymes are present in large excess. Induction of gene expression by hypoxia is likely to be a chronic adaptation that is not required for the acute response to substrates and energy level.

Table 1 HIF-1 regulated transcripts

A total of 446 genes were found to have >2.5-fold change (p < 0.05). In agreement with previous reports on the effects of hypoxia on human MSC and endothelial cells [1719, 35] we found that more genes were down- than up-regulated (243 vs. 203). In order to confirm the microarray results RT-PCR was implemented on 3 known HIF-1α target genes in addition to the glycolytic transcripts described above, and four hypoxia-induced genes not previously identified as HIF targets. As shown in Figure 2, leptin (LEP), metallothionein-3 (MET3) and lysyl oxidase (LOX) were all confirmed to increase in a manner that reflected the array data. Similarly transcripts of non-HIF-regulated genes, keratin-16 (KRT16), Serpin peptidase inhibitor-1 (PAI-1), RAS p21 protein activator (RASA1) and Immediate Early Response 3 (IER3) were also confirmed to increase in the hypoxic samples in a manner that reflected the array data. Western blot and ELISA further confirmed that the levels of secreted hypoxia-marker genes VEGF and leptin were significantly increased by hypoxia in the spent media (Figure 3A and B). Figure 4 and Tables 2 and 3 show some of the most strongly hypoxia-regulated genes separated into functional categories. GO analyses identified 3 categories including [1] antioxidant pathways, NAD(P)H quinone-1, aldo-keto reductase family members and thioredoxin reductase were markedly down-regulated [2] extracellular matrix structural constituents in particular collagens were subject to up and down regulation suggesting dynamic changes of the ECM and [6] glucose transporters that were also subject to both up and down-regulation (see Figure 4 and Tables 2 and 3).

Figure 2
figure 2

Real-time PCR (RTqPCR) validation of differentially expressed genes. RTqPCR was implemented on the same mRNA samples used for microarray as described in Methods. Results are mean ± SEM of n = 3.

Figure 3
figure 3

Quantification of leptin and VEGF protein in MSC spent media during exposure to hypoxia. (A) Numbers in left panel refer to time (h) under hypoxia. (B) VEGF incubation time was 24 h. All results are representative of 3 separate experiments; *p < 0.05, n = 3.

Figure 4
figure 4

Heatmaps of differentially expressed genes in 2 main Gene Ontology (GO) molecular functions. Differentially expressed genes in human MSCs under hypoxia were subjected to GO analysis. A significance cut-off of p < 0.05 was used.

Table 2 Up-regulated genes and functional group allocations
Table 3 Down-regulated genes and functional group allocations

Intermediate and cytoskeleton filaments

Intermediate filaments and microfilaments provide scaffolding that allows cells to rapidly remodel in response to environmental or intracellular signals, for example to increase or decrease movement and/or cell adhesion [43, 44]. Intermediate filaments and their keratin structural components are markers of epithelium where in combination with cell junction components, desmosomes, hemidesmosomes and integrins, they contribute to the adhesion of epithelial cells to the basement membrane [43]. The keratins that compromise intermediate filaments are typically selected from 12 acidic (type 1) and 8 basic (type 2) cytokeratins that usually dimerize with specific partners from the opposite group to form organized filaments. Acidic keratins are coded on chromosome 17q whereas the basic keratins are clustered on chromosome 12q. The gene expression analysis revealed a remarkable increase in the expression of 9 out of 12 acidic keratins and decreased transcript levels of 3 (K2, K4, K5) out of 8 basic keratins with no change of the other members (Tables 2 and 3). To determine whether these changes in gene expression were reflected at the protein and cell structural levels, we measured K16 protein by western blot and the organization of intermediate and microfilaments respectively by keratin and F-actin immunostaining. As shown in Figure 5, K16 levels were increased by 2 ± 0.1-fold (p < 0.05) in MSC after exposure to hypoxia and this was associated with a change in the organization of intermediate filaments and altered cellular morphology (Figure 5B). The changes in cell shape and organization of cytoplasmic filaments indicated by keratin staining were also reflected by similar cell images imaged by F-actin immunostaining (Figure 5C). Other proteins related to intermediate filaments structure and function that were also increased by hypoxia include periplakin and desmoplakin, junctional adhesion molecule 2 (JAM2), adipocyte-specific adhesion molecule (ASAM) adhesion molecule with Ig-like domain 2 (AMIGO2) and podoplanin (PDPN) (see Table 2). In addition, six integrins (α-1,-3,-5,-6,-7, β-1) were increased by >2-fold. These changes may reflect a partial mesenchymal to epithelial transition (MET) that is a reversal of epithelial to mesenchymal transition (EMT) a well-established pathway associated with development and carcinogenesis. MET is predicted to involve reduced cell mobility and proliferative quiescence [45]. Interestingly hypoxia also increased the transcript levels of the HGF receptor c-MET that is normally expressed by cells of epithelial origin (Table 2).

Figure 5
figure 5

Keratin induction and reorganization of intermediate- and micro-filaments under hypoxia. (A) Western blot of Keratin-16 expression in normoxic (N-MSC) and hypoxic (H-MSC) human MSCs. (B) Keratin immunostaining and (C) F-actin immunostaining of MSCs cultured for 24 h under normoxia (N-MSC) or hypoxia (H-MSC). Western blots and immunostaining are described in Methods. All results are representative of 3 separate experiments; *p < 0.05, n = 3.

Regulation of differentiation

MSCs are pluripotent with the capacity to differentiate into adipocytes, osteocytes, chondrocytes, endothelium, neurons and muscle as well as hematopoietic cells. It has been reported that moderate hypoxia >10 mmHg increased [19, 46] or decreased [23, 47, 48] adipogenic differentiation and increased [10, 19, 47, 48] or decreased [20, 23] osteogenic differentiation of MSCs. In our studies, combined GO and manual analyses revealed trends of hypoxia-regulated differentiation markers that suggest repression of adipogenesis and myogenesis in favor of osteogenesis, angiogenesis and hematopoiesis. Leptin, interleukin-11 (Il-11), growth differentiation factor-6 (GDF6), and MMP11, all strongly induced by hypoxia, are secreted factors that favor osteogenesis over adipogenesis [4951]. Collagens are major constituents of the bone matrix and are highly expressed in MSCs. Type 1 collagen, the main osteogenic collagen was not amongst the regulated genes, however, multiple other osteogenic markers genes including aggrecan (ACAN), sparc/osteonectin (SPOCK1), osteoblast cadherin (CDH11) and Wnt pathway members Frizzled homologue-8 (FZD8) and Secreted frizzled related protein 4 (SFRP4) were induced in parallel with decreased expression of the negative osteogenic regulator osteocrin (Tables 2 and 3) [52, 53]. Four of the most strongly induced genes including early growth response 2 (EGR2 > 6-fold) [54], interleukin 11 (Il-11 > 5-fold), growth differentiation factor 6 (GDF6 > 5-fold) and hemopoietic cell kinase (HCK ~5-fold) as well as activin-A receptor (2-fold [55]) are associated with hematopoiesis. Transcript levels of multiple genes associated with angiogenesis were increased by hypoxia [56] while there were marked decreases of multiple striated muscle markers (Tables 2 and 3). The coordinate decrease of 3 Inhibitor of DNA binding factors ID1, ID2, and ID4 (Table 3) is also consistent with a switch in differentiation patterns. Taken together the results support inductions of multiple marker genes for osteogenesis, angiogenesis, and hematopoiesis by severe hypoxia that may predispose cells to these lineages.

Proliferation and metabolism

As discussed above, oxygen tensions above 15 mmHg enhance proliferative and migratory potential of MSCs while tensions below10 mmHg decrease proliferation [812, 1719]. In agreement with this we found that proliferation of hMSCs was decreased under a pO2 of 4 mmHg and the cultures became stationary after 3 days (data not shown). Enhanced proliferation under moderate hypoxia has been attributed to down-regulation of the p21-Ras pathway and increased activity of PI3-kinase-Akt while increased migration was attributed to increased expression of c-MET, VEGF, CXCR4 and CXCR1 (reviewed in [16]). Whereas we also observed increased expression of c-MET and VEGF-A we found no change in expression of CXCR receptors, increased expression of p21-Ras and decreased activity of the PI3-kinase-Akt pathway (see Tables 2 and 3). PI3-kinase-Akt is a central regulator of cell growth and survival. We found that the expression of both IGF-1 and PI3-kinase was significantly decreased under hypoxia while multiple IGF-1BPs were increased (Tables 2 and 3). To determine whether these changes were reflected by parallel changes in pathway activity we measured the phosphorylation of Akt-Thr308 after culture under normoxia or hypoxia. As shown in Figure 6, phosphor-Akt-Thr308 levels were significantly lower after hypoxia. This contrasts with the effects of moderate hypoxia where Akt phosphorylation is increased [26]. Other down-regulated transcripts related to growth and survival included senescence marker galactosidase beta-1-like (−3.0) and death associated kinase-2 (−2.4). There were also decreases of several interleukins and the cell migration cytokine CCL2, also known as MCP-1.

Figure 6
figure 6

Western blot and quantification of phosphor-Akt expression in normoxic (N-MSC) and hypoxic (H-MSC) human MSCs (A and B). Western blot procedures are described in Methods. Akt-P-Thr308 quantification was by NIH image using total Akt as loading control; *p < 0.05, n = 3.

Glucose metabolism

During hypoxic exposure, cells switch from oxidative metabolism to anaerobic glycolysis for energy production. Glycolysis is less efficient than oxidative phosphorylation and more glucose is required to sustain the same level of cell function. Glycolytic enzyme and glucose transporter genes are regulated by HIF-1α and previous studies have reported their induction by moderate hypoxia in the range of 20–40 mmHg [13, 14, 1719]. We found increased transcript levels of at least 2 key glycolytic pathway genes, PFK and PGK, decrease levels of PGAM, no change of GAPDH and decreased transcripts of glucose-6-phosphate dehydrogenase (G6PDH) and pyruvate dehydrogenase PDH; the latter two enzymes are involved in the pentose phosphate pathway (PPP) and acetyl-CoA production respectively. In addition, out of 7 regulated glucose transporters 4 were induced and 3 repressed by hypoxia. The results indicate mixed responses of glucose metabolizing genes. Decreased flux through the PPP is predicted to decrease antioxidant capacity through lower NADPH production, an effect consistent with the down-regulation of other anti-oxidant pathways discussed above.

Hypoxia-related AS events

Exon Microarrays identified 53 genes that were subjected to hypoxia-dependent AS; the arrays also revealed a correlation between high splice index and genes that responded the most strongly at the level of gene transcripts. As shown in Figure 7, eight of the most strongly induced genes and 5 of those that were most strongly repressed also had the highest splicing indices. This suggests a possible mechanistic link between gene regulation and control of AS by hypoxia. LEP, IL-11, IGFBP1, TEK, CA9, LOX4, HCK and EGR2 were each induced by more than 5-fold by hypoxia whereas EFNA3, CORO7, FER1L5, MYH2 and ACAT1 were each repressed (preferentially expressed in air) by greater than 4-fold under hypoxia. Four of the induced genes are HIF-1 targets and EFNA3 a down-regulated transcript, is an inhibitor of angiogenesis and target of the HIF-1-regulated microRNA-210 [57]. GO analyses indicated that the predicted functions of genes with hypoxia-regulated AS included cell adhesion, migration, apoptosis, angiogenesis and oxidation-reduction (Table 4). In terms of exon use, 8 of the identified AS events are described in the human genome annotation (5 C-terminus, 2 cassette-exon, 1 alt-3′ split site), and 16 involve constitutive exons. The remainder are novel AS events.

Figure 7
figure 7

Splicing maps of highly regulated genes with high splice indices. Genes were selected based on highest positive (LEP, IL11, IGFBP1, TEK, CA9, LOX4, HCK, ERG2) and negative (EFNA3, CORO7, FER1L5, MYH2, ACTA1) fold change of gene expression hypoxia vs. normoxia. Bar graphs indicate hypoxic/normoxic differential exonic expression levels (n = 3).

Table 4 GO enrichment classification of alternatively spliced genes after hypoxia

We chose two alternatively spliced genes, ALDH3A2 and NDRG4 for further analysis and qPCR confirmation; ALDH3A2 transcripts were decreased by hypoxia whereas NDRG4 were not regulated. ALDH3A2 is a member of the ALDH superfamily of NAD(P) + −dependent enzymes that catalyze the oxidation of a wide variety of aliphatic and aromatic aldehydes. ALDH3A2 is also known as fatty-aldehyde dehydrogenase (FALDH) because of its role in protecting against lipid peroxidation [58]; loss of ALDH3A2 activity is the cause of Sjogren-Larsson syndrome [59]. As shown in Figure 8A, hypoxia/normoxia dictate the use of AS to generate alternative isoforms of the ALDH3A2 gene. Exon array analysis revealed that the AS event occurred in the last exon, while PCR validation confirmed that the unique exon 10 was not expressed in hypoxic MSCs. This indicates preferentially expression of the M2 isoform under hypoxia (Figure 8B). Previous work has shown that hypoxia can repress the expression of both constitutive and induced ALDH3 isoforms [60] and it is also known that ALDH3A2 gene expression is positively regulated by PPARα [61]. Our finding that PPARα expression is also decreased by hypoxia in parallel with ALDH3A2 (Table 3) suggests a possible mechanism for the transcriptional repression of ALDH3A2 by hypoxia. Other work has shown that AS of the ALDH3A2 gene may determine its subcellular localization with physiological implications for function [62]. Our results are the first to show that hypoxia regulates AS of the FALDH gene.

Figure 8
figure 8

Oxygen-dependent differential splicing of the ALDH3A2 and NDRG4 genes. (A) Expression of splice isoforms confirmed by RT-PCR of genes with prior evidence of AS. (B) Annotation of each alternative isoform (M2 and M1, top graphic), exon structure (middle graphic) and expression profiles (bottom graphic) for ALDH3A2. Light blue boxes indicate down-regulation for hypoxic (H) versus normoxic (N) MSCs; gray boxes indicate no significant change; white boxes indicate no probe detection above expression threshold. The green exon (exon10) indicates the unique exon of isoform M1. Exon expression values (log2) are displayed for both H-MSC (black data points) and N-MSC (red data points), ranked in order of genomic position on the x-axis. The various isoforms can be validated in the PCR result although there is no probe annotated to the isoform-specific exon. (C) Exon structure (top graphic) modified domain of AS (middle graphic) and expression profiles (bottom graphic) for NDRG4. Light red boxes indicate up-regulation for hypoxic (H) versus normoxic MSCs (N) MSCs; gray boxes indicate no significant change; white boxes indicate no probe detection above expression threshold. A Modified or disrupted conserved domain of differentially expressed exon is from the CDD database (NCBI). Exon expression values (log2) are displayed for both H-MSC (black data points) and N-MSC (red data points), ranked in order of genomic position on the x-axis.

The N-myc downstream-regulated gene (NDRG) family belongs to the alpha/beta hydrolase super family of genes that regulate survival, growth and differentiation of host cells. The NDRG4 isoform is preferentially expressed in brain and heart and its overexpression has been linked with an aggressive behavior of meningioma tumors [63]. AS variants of NDRG4 include heart (H) and two brain (B- and Bvar) isoforms [64]. As shown in Figure 8C, exon 5 of the NDRG4 gene was subject to hypoxia-dependent AS according to the exon microarray, and this was validated by qPCR. Alternatively spliced transcripts of the NDRG genes have been described although the functional significance is unknown [64, 65].

Tissue-specific regulation of gene expression and AS by hypoxia

Weigand et al. [35] recently reported results of their global analyses of hypoxia-regulated AS using human umbilical vein endothelial cells (HUVECs). Because our results as well as those of others suggest that many of the effects of hypoxia are tissue-specific we compared our gene expression and AS array data with that of Weigand et al. These results are shown in Supplement Figures 2, 3 and 4. Out of 446 and 113 genes regulated >2.5-fold by hypoxia in hMSCs and HUVECs respectively only 9 were commonly regulated including ANGPTL4, VEGFA, STC, SLC2A3, PGF, EGLN3, ADM, BHLHE40 and NQO1 (Additional file 2: Figure S2). Out of 541 and 239 gene isoforms created by hypoxia-regulated AS in hMSCs and HUVECs respectively, only 17 isoforms were common (Additional file 3: Figure S3). As shown in Additional file 4: Figure S4 comparing all of the genes in hMSCs and HUVECs regulated by hypoxia at the levels of gene expression and AS only 3 genes were common. These included Egln3, a prolyl hydroxylase involved in HIF-1 proteolytic targeting, stanniocalcin-1 (STC) a pro-angiogenic anti-apoptosis gene product [66, 67] (both up-regulated) and NAD(P) Quinone oxidoreductase 1 (NQO1) a Phase II redox detoxification enzyme [68], (down-regulated). It should be noted that Weigand et al. exposed HUVECs to 1% oxygen for 48 h compared with our 0.5% for 24 h therefore the differences may be due in part to different oxygen tension as well as cell type.

Discussion

Our studies suggest that severity of hypoxia within a range of 4–10 mmHg is a critical determinant of global gene expression that has important implications for the biology and function of human MSCs. Transcriptional profiling revealed altered patterns of gene expression particularly involving glucose metabolism, insulin/IGF-1 signaling, intermediate filaments, extracellular matrix, anti-oxidant enzymes, and differentiation-markers. The results are consistent with the quiescent, immobile state conferred by severe hypoxia along with reduced oxidative stress and a switch in favor of osteogenic, angiogenic and perhaps hematopoietic programs over those of adipogenesis and myogenesis. The gene expression profiles are distinct in many respects from those described previously for MSCs subjected to more moderate hypoxia that mediates enhanced proliferation [9, 10, 1719].

Salient features of the present study that define the role of severe hypoxia and distinguish it from moderate hypoxia include the following: [1] Mixed responses of genes involved in glucose metabolism including the pentose phosphate and glycolysis pathways, and moderately induced or decreased levels of glucose transporter gene transcripts. Although glycolytic flux increased acutely under hypoxia, a lower than a maximal induction of all genes may be physiologically beneficial in the long term for a closed environment such as the endosteum where high glycolytic activity would generate localized acidosis and cell death. An analogous condition may be represented in anoxic turtles where the PGAM gene is also repressed and may contribute to the hypometabolic state that is necessary for the turtle to survive extended anoxia [69]. Notable also is the weak induction of Bnip3, a programmed cell death protein normally strongly induced by hypoxia through HIF-1α and co-regulated by hypoxia and acidosis [70]. Also noteworthy is the strong induction of carbonic anhydrase (CA9) suggesting a pre-emptive adaptation to acidosis [2]. IGF-1 and PI3-kinase transcripts were decreased while IGFBP-1, −3 and −5 transcripts were increased, and this was associated with depressed phosphorylation of Akt-Thr308 (Figure 6). The IGF-PI3-kinase pathway regulates cell growth and survival through key targets including mTOR and p70-S6-kinase, and its depressed activity by severe hypoxia is consistent with reduced growth and metabolism. PI3-kinase signaling and phosphor-Akt were reported to be activated in MSCs by moderate hypoxia [26] therefore the reduced activity by severe hypoxia again distinguishes the two conditions [6]. Acidic keratins were markedly increased by hypoxia and this correlated with a more polar organization of intermediate filaments and cell elongation. The acidic keratin genes are clustered on chromosome 17q and the coordinately increase levels of multiple transcripts suggests a common regulatory mechanism. Regulation of keratin expression by the HIF pathway has not been reported. It is also noteworthy that whereas acidic keratin transcripts were all induced, 3 basic keratin gene transcripts decreased. Cytosolic keratins are markers of epithelium and the strong inductions are reminiscent of mesenchymal to epithelial transition (MET) that is associated with reduced rates of growth and migration [4345]. Such a partial transition is again consistent with the quiescent immobile state predicted by the severely hypoxic environment of the endosteum. Changes in expression of multiple collagen and actin genes provides additional evidence for dynamic changes in cytoskeleton and extracellular matrix that contribute to cell growth and mobility [3]. Reduced expression of antioxidant pathway genes in particular NAD(P)H quinone-1 and aldo-keto reductase family members is consistent with lower oxidative stress created by severe hypoxia and may be associated with lower electron leakage and free radical production by mitochondria. Lower ROS and oxidative stress may also contribute to reduced proliferation [4]. Enhanced expression of osteogenic and angiogenic marker genes and depressed expression of adipogenic and myogenic markers is consistent with some but not all previous studies. Mayer et al. [48] and Hung et al. [10] found that osteogenesis was induced in MSCs by an oxygen tension of 10 mmHg. Hung et al. also reported that adipogenesis was reduced by 10 mmHg whereas Fink et al. [46] and Kato et al. [16] reported the opposite effect with increased adipogenesis under the same pO2. The differences may be due to incubation conditions and whether the cell are exposed to reoxygenation during treatments; Hung et al. [12] reported that the effects of hypoxia on proliferation and differentiation were fully reversed by reoxygenation.

The induction of numerous markers of angiogenesis is consistent with many previous reports describing angiogenic activation by hypoxia. The selective decrease of multiple myogenic gene markers by severe hypoxia is a novel observation of this study.

This is the second report to describe global changes in the alternative splicing of genes mediated by hypoxia and the first to describe such an effect in stem cells. There was only a small amount of overlap of hypoxia-regulated alternatively spliced genes in HUVECs versus MSCs with only 17 commonly regulated exons. It is noteworthy that multiple genes that were highly regulated at the transcript level in hMSCs were also subject to differential exon inclusion suggesting possible common regulatory factors in both pathways. Regulators of AS include specific RNA binding proteins and transcription factors (reviewed in [32]). It is known that transcription factor binding can influence AS perhaps by influencing the concentration of direct AS regulators within the transcription complex, or by altering the rate of RNA polymerase II elongation, leading indirectly to AS. Extracellular signals can also modify AS events by changing the activity of regulators for example, hnRNP-A1 (heterogeneous nuclear ribonucleoprotein A1) an RNA binding protein involved in nuclear pre-RNA processing that inhibits the inclusion of multiple alternative exons becomes phosphorylated upon osmotic shock resulting in cytoplasmic accumulation with consequent effects on AS [71]. HIF-1 was shown to regulate differential splicing of the LDH gene in shrimps [72], as well as AS of the hTERT gene that regulates telomerase activity in mammalian cells [73]. Of the most strongly and dually (transcripts and AS) regulated genes LEP, IGFBP1, CA9 and LOXL4 are all HIF-1α regulated. Therefore HIF-1α may contribute to the AS of dependent genes providing dual and coordinated regulation of transcription and AS by hypoxia. It is difficult to assess the physiological significance of most of the AS responding genes without information on the functional consequences, however it was possible to group the genes into functional categories of cell survival, proliferation, mobility and differentiation (Table 4).

Conclusion

We report on molecular genetic and phenotypic changes conferred on hMSC by severe hypoxia. The gene expression changes are largely distinct from those reported previously for more moderate hypoxia that support enhanced proliferation, and the results are consistent with a quiescent, immobile phenotype with reduced metabolic activity and lower oxidative stress. Hypoxia-mediated AS may contribute importantly to gene regulation and protein function during adaptation to a severely hypoxic environment such as that imposed by the bone marrow endosteum.

Methods

Isolation, culture and characterization of human MSC

Normal human bone marrow aspirates were obtained with written consent from healthy donors in accordance with the Declaration of Helsinki and with the approval of the Human subjects Ethics Committee of Second Affiliated Hospital of Zhejiang University. MSCs from 3 such donors were cultured as described previously [10, 12, 19, 74]. Cell surface markers were profiled using a BD FACS CantoTM II Flow Cytometry System after 3–5 passages as described previously [75] with the following human specific monoclonal antibodies: CD29-phycoerythrin (PE) (eBioscience, San Diego, CA, USA), CD34-PE (MACS, Miltenyi Biotec, Auburn, CA, USA) and CD166-PE (BD Biosciences Pharmingen, San Diego, CA, USA), respectively.

Hypoxia

MSCs were plated at 1 × 105 cells/cm2 in complete culture medium and incubated under hypoxia (0.5% O2, 5% CO2) or normoxia (21% O2, 5% CO2) for 24 hours using a ProOX Model C21 system (BioSpherix, Redfield, NY, USA).

RNA extraction

Total RNA from normoxic and hypoxic MSCs (3 independent donor samples each) was extracted using a Kit from Biochain, (Hayward, CA, USA), according to the manufacturer’s instructions. The RNA quality was assessed by formaldehyde agarose gel electrophoresis and quantified using a spectrophotometer (Nanodrop, Wilmington, DE, USA).

RNA amplification and labeling

Gene expression microarray: RNA was amplified, reverse transcribed and labeled as described previously [36]. Briefly, 1 μg of total RNA was amplified using a Message AmpTM II RNA Amplification kit (Life Technologies, Austin, TX, USA). The RNA was reverse transcribed in the presence of cy3-dCTP or cy5-dCTP using Klenow enzyme. For exon microarray, RNA from 3 separate replicate samples of normoxic or hypoxia MSCs was amplified using a Low Input Quick Amp WT Labeling Kit from Agilent.

Array hybridization and data acquisition

Hybridization, scanning and washing were performed on Agilent’s Microarray Platform according to standard protocols. Raw data were acquired using an Agilent DNA microarray scanner and Agilent feature extraction software. All data is MIAME compliant and raw files from the 6 arrays (n = 3 per group) are deposited in the GEO database at NCBI (GEO accession # GSE55875). In addition, we downloaded from the NCBI GEO database the 6 raw data files (n = 3 per group) from the AS study by Weigand et al. [35].

Gene and AS analyses

All raw text files (from hMSC Agilent arrays) and CEL files (from HUVEC Affymetrix arrays) were imported into GeneSpring GX 11 software (Silicon Genetics, Redwood City, CA) for either global gene analysis or alternative splicing. A total of four analytical experiments were performed with 3 replicate arrays for each condition. Normalized expression values were calculated by the Robust Multi-array Average (RMA) method. The resultant signal information was analyzed using one-way analysis of variance (ANOVA) (p < 0.05), assuming normality and equal variances. Multiple Testing Correction (MTC) of p-values by Benjamini Hochberg screened out >90% of regulated genes (33/446) including many of the known HIF-1α-regulated genes such as metallothionein (MT3), carbonic anhydrase XII (CA12; 4.6-fold), (VEGFA; 4.6-fold) and placental growth factor, (PGF; 3.28-fold) (Table 1), as well as many of the gene transcripts that were subsequently confirmed by QPCR including phosphoglycerate kinase (PGK; 7.1-fold), phosphofructokinase (PFK; 3.2-fold), RAS p21 protein activator (RASA1; 2.9-fold) and Immediate Early Response 3 (IER3; 2.5-fold). MTC is known to screen out false positives as well as many true positives, therefore as in our previous micro-array studies [74, 76] it was not employed here. The GeneSpring Cross Gene Error model was applied; this model determines the likelihood of observing a specific fold change to the likelihood of observing a fold measurement by the 50.0th percentile of all measurements in the sample. The average value of expression level for each gene across the samples is set to 1.0 and the resulting normalized signal value plotted for each sample. Lists of differentially expressed genes from different experiments were compared within GeneSpring and displayed as Venn diagrams to show overlapping and non-overlapping genes. Heatmaps and graphs were also generated within GeneSpring.

Conserved domain function annotation of AS exons

To identify protein domains of genes modified by AS, a conserved domain database (CCD) from NCBI web services was used. By default, domain predictions are derived by comparing two protein isoform sequences: one that aligns to the alternative exon and another in which the exon is absent from the corresponding mRNA sequences (competitive isoform).

Validation of the differentially expressed genes

Differentially expressed genes-of-interest were selected based on function and confirmed using quantitative real-time (RT-q) PCR. 2 μg of total RNA respectively from normoxic and hypoxic MSCs was reverse transcribed with oligo (dT)18 using M-MLV reverse transcriptase (TAKARA, Japan) in a volume of 40 μL. Following reverse transcription, 1 μl of this cDNA mixture was employed for a qPCR program of 40 cycles (10 s at 95°C/32 s at 60°C) with SYBR®Premix Ex TaqTM. Data were generated using Perfect Real Time (TAKARA, Japan) and Real-Time PCR (Applied Biosystems, Foster City, CA, USA). Data were analyzed by the 2-∆∆Ct method and results are shown as fold change relative to control.

Confirmation of alternative exon expression

24 alternative exon sequences were selected for confirmation by RT-qPCR. The sequences of the AS exon and its neighboring exons were acquired from the Human Genome UCSC Genome Database. If AS of the last exon of a gene was indicated, the last two-exon sequences of this gene were accessed. Constitutive exon-specific primers were designed and after reverse transcription, qPCR products were separated in 1% agarose gels supplemented with ethidium bromide, and visualized by UV.

Enzyme-linked immunosorbent assay (ELISA)

Leptin and VEGF concentrations in culture medium of normoxic and hypoxic MSCs were measured by ELISA (R&D System, Minneapolis, MN, USA), according to the manufacturer’s instructions. Assays were conducted on a SoftMax Pro® multiplate reader (Molecular Devices, Inc., Downingtown, PA, USA).

Western blot

Our western blot procedures are described in detail elsewhere [77, 78]. Briefly, equal amounts of protein were separated by 10-15% SDS polyacrylamide gel electrophoresis and electro-transferred onto an Immobilon-P Transfer Membrane (Millipore, Billerica, MA, USA). Membranes were blocked with 5% BSA in TBS-T and incubated with VEGF (Santa Cruz Biotechnology, Santa Cruz, CA, USA), leptin, α-tubulin, or Akt (R&D System), antibodies overnight at 4°C. Horseradish peroxide–conjugated secondary antibodies were hybridized by standard procedures. β-actin was used as loading control.

F-actin and keratin intermediate filament staining

Normoxic and hypoxic MSCs were examined for actin filaments organization using Alexa Fluor 555 phalloidin (Invitrogen). Briefly, MSCs were fixed in 4% formaldehyde for 10 min, permeabilized with 0.1% Triton X-100 for 5 min, and blocked with PBS containing 1% BSA for 25 min. Cells were stained with diluted phalloidin in PBS for 20 min at room temperature and viewed under a fluorescence microscope.

Statistical analysis

One-way ANOVA was used to compare experimental groups. Data are expressed as mean ± standard deviation (SD), and a p value < 0.05 was considered as statistically significant.

Abbreviations

AS:

Alternative splicing

MSC:

Mesenchymal stem cell

MET:

Mesenchymal to epidermal transition

VEGF:

Vascular endothelial growth factor.

References

  1. Jones DL, Wagers AJ: No place like home: anatomy and function of the stem cell niche. Nat Rev Mol Cell Biol. 2008, 9: 11-21. 10.1038/nrm2319.

    Article  CAS  PubMed  Google Scholar 

  2. Scadden DT: The stem-cell niche as an entity of action. Nature. 2006, 441: 1075-1079. 10.1038/nature04957.

    Article  CAS  PubMed  Google Scholar 

  3. Eliasson P, Rehn M, Hammar P, Larsson P, Sirenko O, Flippin LA, Cammenga J, Jönsson JI: Hypoxia mediates low cell-cycle activity and increases the proportion of long-term reconstituting hematopoietic stem cells during in vitro culture. Exp Hematol. 2010, 38: 301-310. 10.1016/j.exphem.2010.01.005.

    Article  CAS  PubMed  Google Scholar 

  4. Panchision DM: The role of oxygen in regulating neural stem cells in development and disease. J Cell Physiol. 2009, 220: 562-568. 10.1002/jcp.21812.

    Article  CAS  PubMed  Google Scholar 

  5. Silván U, Díez-Torre A, Arluzea J, Andrade R, Silió M, Aréchaga J: Hypoxia and pluripotency in embryonic and embryonal carcinoma stem cell biology. Differentiation. 2009, 78: 159-168. 10.1016/j.diff.2009.06.002.

    Article  PubMed  Google Scholar 

  6. Eliasson P, Jonsson JI: The hematopoietic stem cell niche: low in oxygen but a nice place to be. J Cell Physiol. 2010, 222: 17-22. 10.1002/jcp.21908.

    Article  CAS  PubMed  Google Scholar 

  7. Ezashi T, Das P, Roberts RM: Low O2 tensions and the prevention of differentiation of hES cells. Proc Natl Acad Sci. 2005, 102: 4783-4788. 10.1073/pnas.0501283102.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  8. Jeong CH, Lee HJ, Cha JH, Kim JH, Kim KR, Kim JH, Yoon DK, Kim KW: Hypoxia-inducible factor-1 alpha inhibits self-renewal of mouse embryonic stem cells in Vitro via negative regulation of the leukemia inhibitory factor-STAT3 pathway. J Biol Chem. 2007, 282: 13672-13679. 10.1074/jbc.M700534200.

    Article  CAS  PubMed  Google Scholar 

  9. Basciano L, Nemos C, Foliguet B, de Isla N, de Carvalho M, Tran N, Dalloul A: Long term culture of mesenchymal stem cells in hypoxia promotes a genetic program maintaining their undifferentiated and multipotent status. BMC Cell Biol. 2011, 12: 12-10.1186/1471-2121-12-12.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  10. Hung SP, Ho JH, Shih YR, Lo T, Lee OK: Hypoxia promotes proliferation and osteogenic differentiation potentials of human mesenchymal stem cells. J Orthop Res. 2012, 30: 260-266. 10.1002/jor.21517.

    Article  PubMed  Google Scholar 

  11. Follmar KE, Decroos FC, Prichard HL, Wang HT, Erdmann D, Olbrich KC: Effects of glutamine, glucose, and oxygen concentration on the metabolism and proliferation of rabbit adipose-derived stem cells. Tissue Eng. 2006, 12: 3525-3533. 10.1089/ten.2006.12.3525.

    Article  CAS  PubMed  Google Scholar 

  12. Hung SC, Pochampally RR, Hsu SC, Sanchez C, Chen SC, Spees J, Prockop DJ: Short-term exposure of multipotent stromal cells to low oxygen increases their expression of CX3CR1 and CXCR4 and their engraftment in vivo. PLoS One. 2007, 2 (5): e416-10.1371/journal.pone.0000416.

    Article  PubMed Central  PubMed  Google Scholar 

  13. Tsai CC, Yew TL, Yang DC, Huang WH, Hung SC: Benefits of hypoxic culture on bone marrow multipotent stromal cells. Am J Blood Res. 2012, 2: 148-159.

    CAS  PubMed Central  PubMed  Google Scholar 

  14. Martin-Rendon E, Hale SJ, Ryan D, Baban D, Forde SP, Roubelakis M, Sweeney D, Moukayed M, Harris AL, Davies K, Watt SM: Transcriptional profiling of human cord blood CD133+ and cultured bone marrow mesenchymal stem cells in response to hypoxia (1.5%). Stem Cells. 2007, 25: 1003-1012. 10.1634/stemcells.2006-0398.

    Article  CAS  PubMed  Google Scholar 

  15. Nazmul H, Mohammad Tariqur R, Noor Hayaty Abu K, Aied Mohammed A: Hypoxic culture conditions as a solution for mesenchymal stem cell based regenerative therapy. Sci World J. 2013, 2013: 632972-

    Google Scholar 

  16. Jin Y, Kato T, Furu M, Nasu A, Kajita Y, Mitsui H, Ueda M, Aoyama T, Nakayama T, Nakamura T, Toguchida J: Mesenchymal stem cells cultured under hypoxia escape from senescence via down-regulation of p16 and extracellular signal regulated kinase. Biochem Biophys Res Commun. 2010, 39: 1471-1476.

    Article  Google Scholar 

  17. 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: 1166-1177. 10.1634/stemcells.2006-0347.

    Article  CAS  PubMed  Google Scholar 

  18. Nekanti U, Dastidar S, Venugopal P, Totey S, Ta M: Increased proliferation and analysis of differential gene expression in human Wharton’s jelly-derived mesenchymal stromal cells under hypoxia. Int J Biol Sci. 2010, 6: 499-512.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  19. Grayson WL, Zhao F, Bunnell B, Ma T: Hypoxia enhances proliferation and tissue formation of human mesenchymal stem cells. Biochem Biophys Res Commun. 2007, 358: 948-953. 10.1016/j.bbrc.2007.05.054.

    Article  CAS  PubMed  Google Scholar 

  20. Salim A, Nacamuli RP, Morgan EF, Giaccia AJ, Longaker MT: Transient changes in oxygen tension inhibit osteogenic differentiation and Runx2 expression in osteoblasts. J Biol Chem. 2004, 279: 40007-40016. 10.1074/jbc.M403715200.

    Article  CAS  PubMed  Google Scholar 

  21. Robins JC, Akeno N, Mukherjee A, Dalal RR, Aronow BJ, Koopman P, Clemens TL: Hypoxia induces chondrocyte-specific gene expression in mesenchymal cells in association with transcriptional activation of Sox9. Bone. 2005, 37: 313-322. 10.1016/j.bone.2005.04.040.

    Article  CAS  PubMed  Google Scholar 

  22. Tsai CC, Su PF, Huang YF, Yew TL, Hung SC: Oct4 and Nanog directly regulate Dnmt1 to maintain self-renewal and undifferentiated state in mesenchymal stem cells. Mol Cell. 2012, 47: 169-182. 10.1016/j.molcel.2012.06.020.

    Article  CAS  PubMed  Google Scholar 

  23. Fehrer C, Brunauer R, Laschober G, Unterluggauer H, Reitinger S, Kloss F, Gülly C, Gassner R, Lepperdinger G: Reduced oxygen tension attenuates differentiation capacity of human mesenchymal stem cells and prolongs their lifespan. Aging Cell. 2007, 6: 745-757. 10.1111/j.1474-9726.2007.00336.x.

    Article  CAS  PubMed  Google Scholar 

  24. Hu X, Wei L, Taylor TM, Wei J, Zhou X, Wang JA, Yu SP: Hypoxic preconditioning enhances bone marrow mesenchymal stem cell migration via Kv2.1 channel and FAK activation. Am J Physiol Cell Physiol. 2011, 301: C362-C372. 10.1152/ajpcell.00013.2010.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  25. Lee SH, Lee YJ, Song CH, Ahn YK, Han HJ: Role of FAK phosphorylation in hypoxia-induced hMSCS migration: involvement of VEGF as well as MAPKS and eNOS pathways. Am J Physiol Cell Physiol. 2010, 298: C847-C856. 10.1152/ajpcell.00418.2009.

    Article  CAS  PubMed  Google Scholar 

  26. Rosová I, Dao M, Capoccia B, Link D, Nolta JA: Hypoxic preconditioning results in increased motility and improved therapeutic potential of human mesenchymal stem cells. Stem Cells. 2008, 26: 2173-2182. 10.1634/stemcells.2007-1104.

    Article  PubMed Central  PubMed  Google Scholar 

  27. Gustafsson MV, Zheng X, Pereira T, Gradin K, Jin S, Lundkvist J, Ruas JL, Poellinger L, Lendahl U, Bondesson M: Hypoxia requires Notch signaling to maintain the undifferentiated cell state. Dev Cell. 2005, 9: 617-628. 10.1016/j.devcel.2005.09.010.

    Article  CAS  PubMed  Google Scholar 

  28. Tsai CC, Chen YJ, Yew TL, Chen LL, Wang JY, Chiu CH, Hung SC: Hypoxia inhibits senescence and maintains mesenchymal stem cell properties through down-regulation of E2A-p21 by HIF-TWIST. Blood. 2011, 117: 459-469. 10.1182/blood-2010-05-287508.

    Article  CAS  PubMed  Google Scholar 

  29. Pan Q, Shai O, Lee LJ, Frey BJ, Blencowe BJ: Deep surveying of alternative splicing complexity in the human transcriptome by high-throughput sequencing. Nat Genet. 2008, 40: 1413-1415. 10.1038/ng.259.

    Article  CAS  PubMed  Google Scholar 

  30. Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, Kingsmore SF, Schroth GP, Burge CB: Alternative isoform regulation in human tissue transcriptomes. Nature. 2008, 456: 470-476. 10.1038/nature07509.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  31. David CJ, Manley JL: The search for alternative splicing regulators: new approaches offer a path to a splicing code. Genes Dev. 2008, 22: 279-285. 10.1101/gad.1643108.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  32. David CJ, Manley JL: Alternative pre-mRNA splicing regulation in cancer: pathways and programs unhinged. Genes Dev. 2010, 24: 2343-2364. 10.1101/gad.1973010.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  33. Mills JD, Janitz M: Alternative splicing of mRNA in the molecular pathology of neurodegenerative diseases. Neurobiol Aging. 2011, 33: 1012 e1011-1024.

    Google Scholar 

  34. Moller-Levet CS, Betts GN, Harris AL, Homer JJ, West CM, Miller CJ: Exon array analysis of head and neck cancers identifies a hypoxia related splice variant of LAMA3 associated with a poor prognosis. PLoS Comput Biol. 2009, 5: e1000571-10.1371/journal.pcbi.1000571.

    Article  PubMed Central  PubMed  Google Scholar 

  35. Weigand JE, Boeckel JN, Gellert P, Dimmeler S: Hypoxia-induced alternative splicing in endothelial cells. PLoS One. 2012, 7: e42697-10.1371/journal.pone.0042697.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  36. Semenza GL: HIF-1 mediates metabolic responses to intratumoral hypoxia and oncogenic mutations. J Clin Invest. 2013, 123: 3664-3671. 10.1172/JCI67230.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  37. Minchenko O, Opentanova I, Minchenko D, Ogura T, Esumi H: Hypoxia induces transcription of 6-phosphofructo-2-kinase/fructose-2,6-biphosphatase-4 gene via hypoxia-inducible factor-1alpha activation. FEBS Lett. 2004, 576: 14-20. 10.1016/j.febslet.2004.08.053.

    Article  CAS  PubMed  Google Scholar 

  38. Li H, Ko HP, Whitlock JP: Induction of phosphoglycerate kinase 1 gene expression by hypoxia. Roles of Arnt and HIF1alpha. J Biol Chem. 1996, 271: 21262-21267. 10.1074/jbc.271.35.21262.

    Article  CAS  PubMed  Google Scholar 

  39. Takahashi Y, Takahashi S, Yoshimi T, Miura T: Hypoxia-induced expression of phosphoglycerate mutase B in fibroblasts. Eur J Biochem. 1998, 254: 497-504. 10.1046/j.1432-1327.1998.2540497.x.

    Article  CAS  PubMed  Google Scholar 

  40. Lu S, Gu X, Hoestje S, Epner DE: Identification of an additional hypoxia responsive element in the glyceraldehyde-3-phosphate dehydrogenase gene promoter. Biochim Biophys Acta. 2002, 1574: 152-156. 10.1016/S0167-4781(01)00359-1.

    Article  CAS  PubMed  Google Scholar 

  41. Said HM, Hagemann C, Stojic J, Schoemig B, Vince GH, Flentje M, Roosen K, Vordermark D: GAPDH is not regulated in human glioblastoma under hypoxic conditions. BMC Mol Biol. 2007, 8: 55-10.1186/1471-2199-8-55.

    Article  PubMed Central  PubMed  Google Scholar 

  42. Higashimura Y, Nakajima Y, Yamaji R, Harada N, Shibasaki F, Nakano Y, Inui H: Up-regulation of glyceraldehyde-3-phosphate dehydrogenase gene expression by HIF-1 activity depending on Sp1 in hypoxic breast cancer cells. Arch Biochem Biophys. 2011, 509: 1-8. 10.1016/j.abb.2011.02.011.

    Article  CAS  PubMed  Google Scholar 

  43. Chung BM, Rotty JD, Coulombe PA: Networking galore: intermediate filaments and cell migration. Curr Opin Cell Biol. 2013, 25: 600-612. 10.1016/j.ceb.2013.06.008.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  44. Windoffer R, Beil M, Magin TM, Leube RE: Cytoskeleton in motion: the dynamics of keratin intermediate filaments in epithelia. J Cell Biol. 2011, 194: 669-678. 10.1083/jcb.201008095.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  45. Lamouille S, Subramanyam D, Blelloch R, Derynck R: Regulation of epithelial- mesenchymal and mesenchymal-epithelial transitions by microRNAs. Curr Opin Cell Biol. 2013, 25: 200-207. 10.1016/j.ceb.2013.01.008.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  46. Fink T, Abildtrup L, Fogd K, Abdallah BM, Kassem M, Ebbesen P, Zachar V: Induction of adipocyte-like phenotype in human mesenchymal stem cells by hypoxia. Stem Cells. 2004, 22: 1346-1355. 10.1634/stemcells.2004-0038.

    Article  PubMed  Google Scholar 

  47. Wagegg M, Gaber T, Lohanatha FL, Hahne M, Strehl C, Fangradt M, Tran CL, Schonbeck K, Hoff P, Ode A, Perka C, Duda GN, Buttgereit F: Hypoxia promotes osteogenesis but suppresses adipogenesis of human mesenchymal stromal cells in a hypoxia-inducible factor-1 dependent manner. PLoS One. 2012, 7 (9): e46483-10.1371/journal.pone.0046483. doi:10.1371/journal.pone.0046483

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  48. Mayer H, Bertram H, Lindenmaier W, Korff T, Weber H, Weich H: Vascular endothelial growth factor (VEGF-A) expression in human mesenchymal stem cells: autocrine and paracrine role on osteoblastic and endothelial differentiation. J Cell Biochem. 2005, 95: 827-839. 10.1002/jcb.20462.

    Article  CAS  PubMed  Google Scholar 

  49. Avecilla ST, Hattori K, Heissig B, Tejada R, Liao F, Shido K, Jin DK, Dias S, Zhang F, Hartman TE, Hackett NR, Crystal RG, Witte L, Hicklin DJ, Bohlen P, Eaton D, Lyden D, De Sauvage F, Rafii S: Chemokine-mediated interaction of hematopoietic progenitors with the bone marrow vascular niche is required for thrombopoiesis. Nat Med. 2004, 10: 64-71. 10.1038/nm973.

    Article  CAS  PubMed  Google Scholar 

  50. Turner KJ, Neben S, Weich N, Schaub RG, Goldman SJ: The role of recombinant interleukin 11 in megakaryocytopoiesis. Stem Cells. 1996, 14 (Suppl 1): 53-61.

    Article  PubMed  Google Scholar 

  51. Motrescu ER, Rio MC: Cancer cells, adipocytes and matrix metalloproteinase 11: a vicious tumor progression cycle. Biol Chem. 2008, 389: 1037-1041.

    Article  CAS  PubMed  Google Scholar 

  52. Eijken M, Meijer IM, Westbroek I, Koedam M, Chiba H, Uitterlinden AG, Pols HA, van Leeuwen JP: Wnt signaling acts and is regulated in a human osteoblast differentiation dependent manner. J Cell Biochem. 2008, 104 (2): 568-579. 10.1002/jcb.21651.

    Article  CAS  PubMed  Google Scholar 

  53. Yerges LM, Zhang Y, Cauley JA, Kammerer CM, Nestlerode CS, Wheeler VW, Patrick AL, Bunker CH, Moffett SP, Ferrell RE, Zmuda JM: Functional characterization of genetic variation in the Frizzled 1 (FZD1) promoter and association with bone phenotypes: more to the LRP5 story?. J Bone Miner Res. 2009, 24: 87-96. 10.1359/jbmr.080816.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  54. Li S, Miao T, Sebastian M, Bhullar P, Ghaffari E, Liu M, Symonds AL, Wang P: The transcription factors Egr2 and Egr3 are essential for the control of inflammation and antigen-induced proliferation of B and T cells. Immunity. 2012, 37: 685-696. 10.1016/j.immuni.2012.08.001.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  55. Shav-Tal Y, Zipori D: The role of activin a in regulation of hemopoiesis. Activin A expression is induced during hemopoietic differentiation. Stem Cells. 2002, 20: 493-500. 10.1634/stemcells.20-6-493.

    Article  CAS  PubMed  Google Scholar 

  56. Sadanandam A, Rosenbaugh EG, Singh S, Varney M, Singh RK: Semaphorin 5A promotes angiogenesis by increasing endothelial cell proliferation, migration, and decreasing apoptosis. Microvasc Res. 2010, 79: 1-9. 10.1016/j.mvr.2009.10.005.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  57. Hu S, Huang M, Li Z, Jia F, Ghosh Z, Lijkwan MA, Fasanaro P, Sun N, Wang X, Martelli F, Robbins RC, Wu JC: MicroRNA-210 as a novel therapy for treatment of ischemic heart disease. Circulation. 2010, 122 (11 Suppl): S124-S131.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  58. Demozay D, Mas JC, Rocchi S, Van Obberghen E: FALDH reverses the deleterious action of oxidative stress induced by lipid peroxidation product 4-hydroxynonenal on insulin signaling in 3T3-L1 adipocytes. Diabetes. 2008, 57: 1216-1226. 10.2337/db07-0389.

    Article  CAS  PubMed  Google Scholar 

  59. De Laurenzi V, Rogers GR, Hamrock DJ, Marekov LN, Steinert PM, Compton JG, Markova N, Rizzo WB: Sjögren-Larsson syndrome is caused by mutations in the fatty aldehyde dehydrogenase gene. Nat Genet. 1996, 12: 52-57. 10.1038/ng0196-52.

    Article  CAS  PubMed  Google Scholar 

  60. Reisdorph R, Lindahl R: Hypoxia exerts cell-type-specific effects on expression of the class 3 aldehyde dehydrogenase gene. Biochem Biophys Res Commun. 1998, 249: 709-712. 10.1006/bbrc.1998.9216.

    Article  CAS  PubMed  Google Scholar 

  61. Motojima K, Hirai T: Peroxisome proliferator-activated receptor alpha plays a vital role in inducing a detoxification system against plant compounds with crosstalk with other xenobiotic nuclear receptors. FEBS J. 2006, 273: 292-300. 10.1111/j.1742-4658.2005.05060.x.

    Article  CAS  PubMed  Google Scholar 

  62. Ashibe B, Hirai T, Higashi K, Sekimizu K, Motojima K: Dual subcellular localization in the endoplasmic reticulum and peroxisomes and a vital role in protecting against oxidative stress of fatty aldehyde dehydrogenase are achieved by alternative splicing. J Biol Chem. 2007, 282: 20763-20773. 10.1074/jbc.M611853200.

    Article  CAS  PubMed  Google Scholar 

  63. Kotipatruni RP, Ferraro DJ, Ren X, Vanderwaal RP, Thotala DK, Hallahan DE, Jaboin JJ: NDRG4, the N-Myc downstream regulated gene, is important for cell survival, tumor invasion and angiogenesis in meningiomas. Integr Biol (Camb). 2012, 4: 1185-1197. 10.1039/c2ib20168b.

    Article  CAS  Google Scholar 

  64. Zhou RH, Kokame K, Tsukamoto Y, Yutani C, Kato H, Miyata T: Characterization of the human NDRG gene family: a newly identified member, NDRG4, is specifically expressed in brain and heart. Genomics. 2001, 73: 86-97. 10.1006/geno.2000.6496.

    Article  CAS  PubMed  Google Scholar 

  65. Svensson E, Vidovic K, Olofsson T, Vallon-Christersson J, Borg A, Gullberg U: The Wilms’ tumor gene 1 (WT1) induces expression of the N-myc downstream regulated gene 2 (NDRG2). DNA Cell Biol. 2007, 26: 589-597. 10.1089/dna.2007.0586.

    Article  CAS  PubMed  Google Scholar 

  66. Kim SJ, Ko JH, Yun JH, Kim JA, Kim TE, Lee HJ, Kim SH, Park KH, Oh JY: Stanniocalcin-1 protects retinal ganglion cells by inhibiting apoptosis and oxidative damage. PLoS One. 2013, 8 (5): e63749-10.1371/journal.pone.0063749. doi:10.1371/journal.pone.0063749

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  67. Law AY, Wong CK: Stanniocalcin-1 and −2 promote angiogenic sprouting in HUVECs via VEGF/VEGFR2 and angiopoietin signaling pathways. Mol Cell Endocrinol. 2013, 374: 73-81. 10.1016/j.mce.2013.04.024.

    Article  CAS  PubMed  Google Scholar 

  68. Harada S, Fujii C, Hayashi A, Ohkoshi N: An association between idiopathic Parkinson’s disease and polymorphisms of phase II detoxification enzymes: glutathione S-transferase M1 and quinone oxidoreductase 1 and 2. Biochem Biophys Res Commun. 2001, 288: 887-892. 10.1006/bbrc.2001.5868.

    Article  CAS  PubMed  Google Scholar 

  69. Zhang J, Biggar KK, Storey KB: Regulation of p53 by reversible post-transcriptional and post-translational mechanisms in liver and skeletal muscle of an anoxia tolerant turtle. Trachemys Scr Elegans Gene. 2013, 513: 147-155.

    CAS  Google Scholar 

  70. Kubasiak LA, Hernandez OM, Bishopric NH, Webster KA: Hypoxia and acidosis activate cardiac myocyte death through the Bcl-2 family protein BNIP3. Proc Natl Acad Sci. 2002, 99: 12825-12830. 10.1073/pnas.202474099.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  71. Allemand E, Guil S, Myers M, Moscat J, Caceres JF, Krainer AR: Regulation of heterogenous nuclear ribonucleoprotein A1 transport by phosphorylation in cells stressed by osmotic shock. Proc Natl Acad Sci. 2005, 102: 3605-3610. 10.1073/pnas.0409889102.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  72. Soñanez-Organis JG, Rodriguez-Armenta M, Leal-Rubio B, Peregrino-Uriarte AB, Gómez-Jiménez S, Yepiz-Plascencia G: Alternative splicing generates two lactate dehydrogenase subunits differentially expressed during hypoxia via HIF-1 in the shrimp Litopenaeus vannamei. Biochimie. 2012, 94: 1250-1260. 10.1016/j.biochi.2012.02.015.

    Article  PubMed  Google Scholar 

  73. Anderson CJ, Hoare SF, Ashcroft M, Bilsland AE, Keith WN: Hypoxic regulation of telomerase gene expression by transcriptional and post-transcriptional mechanisms. Oncogene. 2006, 25: 61-69.

    CAS  PubMed  Google Scholar 

  74. Wilson A, Shehadeh L, Yu H, Webster KA: Age-related changes in the molecular genetic profiles of murine bone marrow mesenchymal stem cells. BMC Genomics. 2010, 11: 229-10.1186/1471-2164-11-229.

    Article  PubMed Central  PubMed  Google Scholar 

  75. Wu R, Gu B, Zhao X, Tan Z, Chen L, Zhu J, Zhang M: Derivation of multipotent nestin(+)/CD271 (−)/STRO-1 (−) mesenchymal-like precursors from human embryonic stem cells in chemically defined conditions. Hum Cell. 2013, 26: 19-27. 10.1007/s13577-011-0022-3.

    Article  CAS  PubMed  Google Scholar 

  76. Shehadeh LA, Webster KA, Hare JM, Vazquez-Padron RI: Dynamic regulation of vascular myosin light chain (MYL9) with injury and aging. PLoS One. 2011, 6 (10): e25855-10.1371/journal.pone.0025855.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  77. Wei J, Wang W, Chopra I, Li HF, Dougherty CJ, Adi J, Adi N, Wang H, Webster KA: c-Jun N-terminal kinase (JNK-1) confers protection against brief but not extended ischemia during acute myocardial infarction. J Biol Chem. 2011, 286: 13995-14006. 10.1074/jbc.M110.211334.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  78. Shao H, Xu Q, Wu Q, Ma Q, Salgueiro L, Wang J, Eton D, Webster KA, Yu H: Defective CXCR4 expression in aged bone marrow cells impairs vascular regeneration. J Cell Mol Med. 2011, 15: 2046-2056. 10.1111/j.1582-4934.2010.01231.x.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by grants from the National High Technology Research and Development Program (“863” Program) of China (2011AA020102), National Natural Science Foundation of China (31000650; 31101052; 31171418), Zhejiang Provincial Natural Science Foundation of China (Y2090308; Y2110254), Florida State (3KN05 to LAS), and the National Institute of Health (K01AG040468 to LAS and RO1 HL072924-01 to KAW).

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Keith A Webster or Jian’an Wang.

Additional information

Competing interest

The authors declare that they have no competing interests financial or non-financial related to the content of this article.

Authors’ contributions

XH: conception, design of experiments and manuscript writing; RW: implement experiments, data analysis, manuscript writing; LS: data analysis, bioinformatics, statistics, manuscript writing; LZ, XH, FG, XL: implement experiments; HY: design of experiments and manuscript writing; KAW: conception, design of experiments, data analysis, manuscript writing, corresponding author; JW, conception, design of experiments, data analysis, manuscript writing, corresponding author. All authors read and approved the final manuscript.

Xinyang Hu, Rongrong Wu, Lina A Shehadeh contributed equally to this work.

Electronic supplementary material

12864_2013_6014_MOESM1_ESM.tiff

Additional file 1: Figure S1: Differential MSC gene expression. Heatmap of the top differentially expressed genes by at least 2.0-fold (p < .01) in hMSCs under hypoxia relative to normoxia. (TIFF 6 MB)

12864_2013_6014_MOESM2_ESM.tiff

Additional file 2: Figure S2: Venn diagram and heatmap of differentially expressed genes that overlap in hMSCs and HUVECs under hypoxia. (A). Differentially expressed genes determined by our analyses of HUVECs and human MSCs under hypoxia were compared to find common differentially expressed transcripts. A 2.0 fold change and p < 0.01 significance cut-offs were used. (B). Selected genes, including VEGF-A, from the 9 transcripts overlapping in HUVEC and hMSCs, are shown by heatmap displaying expression levels in normoxic and hypoxic MSCs. Color bar shown is Log2. (TIFF 15 MB)

12864_2013_6014_MOESM3_ESM.jpeg

Additional file 3: Figure S3: Venn diagram and graphs of differentially expressed isoforms overlapping in hMSCs and HUVECs under hypoxia. (A). Differentially expressed isoforms determined by our analysis of HUVECs and human MSCs under hypoxia were compared to find common differentially expressed exons. A 1.0 splicing index and p < 0.01 significance cutoffs were used. (B). Exonic expression of 17 isoforms overlapping in HUVECs and hMSCs are shown in the graphs displaying the hypoxic vs. normoxic exonic expression levels. (JPEG 802 KB)

12864_2013_6014_MOESM4_ESM.tiff

Additional file 4: Figure S4: Venn Diagrams and heatmap of differentially expressed genes overlapping genes in hMSCs versus HUVECs. Differentially expressed genes and isoforms determined by our analysis of HUVECs and human MSCs under hypoxia were compared to find common or unique differentially expressed genes and isoforms. A 1.0 splicing index, 2.0 fold change, and p < 0.01 significance cut-offs were used. (TIFF 12 MB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( https://creativecommons.org/publicdomain/zero/1.0/ ) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hu, X., Wu, R., Shehadeh, L.A. et al. Severe hypoxia exerts parallel and cell-specific regulation of gene expression and alternative splicing in human mesenchymal stem cells. BMC Genomics 15, 303 (2014). https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2164-15-303

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2164-15-303

Keywords