Time series gene expression profiling and temporal regulatory pathway analysis of BMP6 induced osteoblast differentiation and mineralization
© Luo et al; licensee BioMed Central Ltd. 2011
Received: 1 February 2011
Accepted: 23 May 2011
Published: 23 May 2011
BMP6 mediated osteoblast differentiation plays a key role in skeletal development and bone disease. Unfortunately, the signaling pathways regulated by BMP6 are largely uncharacterized due to both a lack of data and the complexity of the response.
To better characterize the signaling pathways responsive to BMP6, we conducted a time series microarray study to track BMP6 induced osteoblast differentiation and mineralization. These temporal data were analyzed using a customized gene set analysis approach to identify temporally coherent sets of genes that act downstream of BMP6. Our analysis identified BMP6 regulation of previously reported pathways, such as the TGF-beta pathway. We also identified previously unknown connections between BMP6 and pathways such as Notch signaling and the MYB and BAF57 regulatory modules. In addition, we identify a super-network of pathways that are sequentially activated following BMP6 induction.
In this work, we carried out a microarray-based temporal regulatory pathway analysis of BMP6 induced osteoblast differentiation and mineralization using GAGE method. This novel temporal analysis is more informative and powerful than the classical static pathway analysis in that: (1) it captures the interconnections between signaling pathways or functional modules and demonstrates the even higher level organization of molecular biological systems; (2) it describes the temporal perturbation patterns of each pathway or module and their dynamic roles in osteoblast differentiation. The same set of experimental and computational strategies employed in our work could be useful for studying other complex biological processes.
Osteoblasts are responsible for bone matrix production and mineralization . In concert with osteoclasts, osteoblasts coordinate bone remodeling, a physiologic process by which bone mass is maintained . Osteoblasts arise from mesenchymal stem cells (MSC) . However, the mechanism controlling this differentiation process is not well understood.
MSC differentiation to osteoblasts can be induced both in vivo and in vitro by soluble factors in the bone morphogenetic protein family (BMP) [2, 3]. Among the BMPs, BMP2, 4, 6 and 7 are the best characterized osetogenic factors . Our previous work  has shown that: (1) human MSC produce BMP6 in defined, serum-free conditions, (2) BMP6 is up-regulated under mild osteogenic stimulus (dexamethasone), (3) exogenous BMP6 potently induces osteoblast differentiation, but responses to BMP2, 4, or 7 are inconsistent and require higher doses, (4) exogenous BMP6 induces the expression or up-regulation of a set of osteoblast-related genes in human MSC, and (5) 24 hour treatment with BMP6 induce high levels of osteoblast gene expression and cause mineralization. These results established BMP6 as an endogenous regulator of human osteoblast differentiation .
Unfortunately, BMP6 signaling largely remains uncharacterized. Significant efforts, particularly a series of high throughput microarray studies [6–12] have been undertaken to uncover BMP (including BMP6) responsive genes, transcriptional programs, and their roles in osteoblast development. However, an integrated understanding of the regulatory mechanisms for osteoblast differentiation and mineralization has not been achieved. In particular, two remaining questions include: (1) what pathways and gene groups are responsible for MSC differentiation to osteoblasts in response to BMP6 stimulation? (2) what is the sequence of pathway activation events in response to BMP6 during the osteogenic induction?
To answer these two questions, we conducted a time series microarray study on BMP6 osteogenic induction and a comprehensive pathway analysis on the temporal data. We met special challenges at both experiment and data analysis levels.
At data analysis level, we developed a novel temporal pathway analysis procedure. Gene set analysis (GSA) is a well established strategy to identify pathways or gene sets associated with a particular phenotype or condition [13–16]. Temporary analysis was done on biological processes previously . However, traditional GSA methods do not directly apply without externally specified selection criteria for temporal changes. Methods like GSEA  accept continuous phenotype labels hence can correlate gene set changes to time or pre-defined temporal patterns. To our knowledge, these methods do not infer temporal pathway level changes without reference patterns. For such temporal pathway analysis to be practical, a method needs to be: (1) applicable to datasets with small or even changing sample size at each time point or condition. (2) both sensitive and selective to capture subtle yet real regulatory signals over time. To meet these challenges, we designed a special analysis procedure (Figure 1) based on our newly developed GAGE (Generally Applicable Gene-set Enrichment) method .
Using this joint experimental and computational approach, we identified both the pathways and their temporal responses to BMP6 signaling during osteoblast differentiation and mineralization.
Following our previous study , we explored BMP6 induced human MSC osteoblast gene expression and function. Our preliminary experiments showed that 8 hours BMP6 exposure was sufficient to induce expression of osteoblast differentiation marker genes in human MSCs. A minimum of 24 hours exposure was required to form mineralized matrix at 14 days after the initiation of BMP treatment. A maximal mineralization response was observed upon 96 hours of BMP6 treatment.
We designed a microarray study to explore the regulatory mechanisms underneath these phenotypic changes at different stages of human MSC osteoblast differentiation (Figure 1). Our recently developed GAGE method  was applied to infer the significantly perturbed KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways, GO (Gene Ontology) term groups, and experimentally derived gene sets (experimental sets for short) by BMP6 treatment at different times along the induction process (Figure 1). We examined these significant gene sets in details below.
Significantly perturbed KEGG pathways during BMP6 osteogenic induction
Interpretation and validation of the significant KEGG pathways inferred by GAGE.
TGF-beta signaling pathway
BMP signal targets IDs, SMAD6-7, DLX5 extremely up
Cytokine-cytokine receptor interaction
IFN target gene sets down (Table 4)
Wnt signaling pathway
Jak-STAT signaling pathway
STAT1 target genes down
MAPK signaling pathway
p53 signaling pathway
These significant KEGG pathways show different temporal perturbation patterns (Figure 2). The TGF-beta signaling pathway, Cytokine-cytokine receptor interaction, and Wnt signaling pathway are most perturbed at 8h. The Jak-STAT signaling pathway, MAPK signaling pathway and p53 signaling pathway are most perturbed at 24h. Focal adhesion and ECM-receptor interaction are most perturbed at 10d.
The overlaps in perturbed member genes between the significant KEGG pathways inferred by GAGE
Significantly perturbed GO term gene sets during BMP6 osteogenic induction
Interpretation and validation information of the significant GO term gene sets inferred by GAGE
Notch signaling pathway
insulin-like growth factor receptor binding
insulin-like growth factor binding
homophilic cell adhesion
8h-24h, 10d down
The top GO term gene sets suggest novel yet plausible regulatory processes involved in BMP6 induced osteoblast differentiation and mineralization. Some GO terms are self-explaining, including skeletal development and ossification (Figure 5b). Other less evident top GO terms are supported by the literature (Table 3). For example, the notch signaling pathway (Figure 5a) is directly involved in BMP2 induced osteoblast differentiation [11, 25], suggesting a possible connection to BMP6 too. Both insulin-like growth factor receptor binding and insulin-like growth factor binding (Figure 5a) suggest that the IGF signal is important for osteoblast differentiation and mineralization. Indeed IGF1  and IGF1R  promote mineralization and bone formation, whereas IGFBP 3-6 [28–31] sequesters IGF1 and inhibits osteoblast differentiation and mineralization. Direct connections between immune response (Figure 5b) and bone formation and metabolism have been well appreciated [32, 33], which merge into a new research area, osteoimmunology . Genes from homophilic cell adhesion  and oxidoreductase activity  (Figure 5b) sets play a role in osteoblast differentiation. Selenium (Figure 5a) deficiency is associated with osteoporosis , a disease in bone formation and metabolism.
Significantly perturbed experimentally derived gene sets during BMP6 osteogenic induction
Interpretation and validation information on the significant experimental sets inferred by GAGE
MYB transcription factor is identified as a novel regulator for osteoblast differentiation. MYB target gene set was down regulated at 8h, 24h and 10d with no change in MYB expression level. MYB transcriptional activity at protein level can be inhibited through two potential mechanisms following BMP6 treatment: the activation of Wnt signaling pathways (Table 1 and Figure 2) phosphorylates and degrades of MYB protein , and BMP/TGF-beta and Wnt responsive OVOL1 antagonizes transcriptional activation of MYB by competing for target promoter binding .
Another novel transcriptional regulator for osteoblast differentiation we predicted is BAF57, which is the regulatory subunit SWI/SNF chromatin remodeling complex . Multiple BAF57 target genes are directly related to osteoblast differentiation and function (Additional file 1: Supplementary Table 1). Indeed, SWI/SNF regulates osteoblast-specific transcription through chromatin structure modification . BMP6 treatment may target SWI/SNF to nucleus though SMAD1 signal  or p38 MAPK pathway targets SWI-SNF chromatin-remodeling complex . Interestingly, BAF57 positive target genes are up-regulated and negative targets down-regulated during BMP6 induction, which further confirms the involvement of BAF57 activity. The different timing of the positive and negative regulation likely suggests different dynamics of these actions.
Other interesting regulatory mechanisms are inferred based on the top ranking experimental sets (Table 4). Interferon beta (interferon alpha and gamma too, but not shown) positive target gene sets are down-regulated, which is consistent with Jak-STAT pathway from KEGG (Table 1). VEGF positive targets are down-regulated, likely because VEGF gene expression is down-regulated due to MYB inhibition . IRS negative targets are down-regulated is consistent with activation of IGF signal, particularly up-regulation of IGF receptor binding proteins (Table 4 and Figure 5).
Discussion and Conclusion
This is the first high throughput microarray study on BMP6 induced transcriptional program in human MSC. It covers the whole process from early to late stage osteoblast differentiation and mineralization. We conducted a comprehensive gene set analysis to identify relevant regulatory mechanisms and functional groups. We inferred a series of significant KEGG pathways, GO terms and experimental sets at different stages of BMP6 induction process. We not only showed which pathways or gene sets are significant, but also when and how they are involved in the osteoblast differentiation and mineralization. Different from common pathway analyses [13, 14, 16], our work further captures the interconnections among individual pathways or functional groups and integrate them into a whole system. Taken together, this work provides clearer mechanistic picture of osteoblast differentiation and function.
We inferred novel and coherent sets of regulatory mechanisms downstream of BMP6 signaling during osteoblast differentiation and mineralization. First, the same set of KEGG pathways are constantly involved in BMP6 induction. Their roles in osteogenic induction are clarified based on their perturbation patterns and connected to relevant discoveries in literature. These significant KEGG pathways are not separated but rather they work as a unified super regulatory system, and the pathway perturbation patterns we derived reflect a dynamic transmission process of the regulatory signal at transcriptional level along the super system. Second, a varying set of GO processes and functional groups are involved at different stage of BMP6 induced osteoblast differentiation and mineralization. These suggest novel yet plausible regulatory mechanisms, which are connected to but have not directly and explicitly introduced in literature works. Third, the most significant experimental sets suggest novel transcriptional regulators including MYB and BAF57, and regulatory pathways consistent with predictions based on KEGG and GO gene sets above.
Connections between KEGG pathways are evident as shown in the super regulatory network of pathways (Figure 4 and Table 2). Perturbations propagate along the super network at two levels: at protein level, the phosphorylation, binding, activation/inhibition events relay along pathways and transmit into interconnected pathways, as stated by KEGG graphs (Figure 3); at transcriptional level, gene expression perturbation propagates through auto-regulatory (feedback and feed forward) loops within pathways, and bridges into its neighbor pathways through the multiple component genes in common, as suggested by our pathway analysis results (Figure 4 and Table 2). Protein level transmission is faster, but transcriptional level transmission lasts longer hence ensure the long term biological effects. These longer lasting transcriptional effects are clearly evident in our data as the gene expression levels were perturbed long after the withdrawal of BMP6 treatment. The hard wired KEGG pathways and interconnections between them define how BMP6 signal triggers downstream programs (Figure 4).
There are also connections between BMP6 signal and the processes/groups represented by GO terms and experimental sets. For example, Notch signal and IGF signal are involved in the whole induction process (Figure 5a) like all significant KEGG pathways (Figure 2). It follows that these two signals should be also part of the super regulatory system and interconnected with multiple significant KEGG pathways of the network (Figure 4). This hypothesis is well supported by our data and literature: (1) Notch signal directly interacts with BMP signal. SMAD1 and NIC synergize to induce expression of HEY1 and other Notch targets [44–46]. Indeed up-regulation of HEY1 (and HEY2) requires continued BMP6 treatment by 96 hours (Additional file 1: Supplementary Figure 4c). Besides this binding-synergy at protein level, Notch ligand JAG1 expression is also up-regulated directly by BMP6 (Additional file 1: Supplementary Figure 4c). Notch also interacts with Wnt signal . (2) Growth hormone (GH) signal  as part of cytokine-cytokine receptor interaction (KEGG pathway 04060) and Jak-STAT pathways (KEGG pathway 04630) are activated and by BMP signal (Figure 4) with GHR up-regulated (Additional file 1: Supplementary Figure 4a-b), GH signal up-regulates (IGF1 and IRS1, Additional file 1: Supplementary Figure 4d) and activates IGF signal in turn . Connections between BMP signal and the two predicted transcriptional regulators, MYB and BAF57, are described above in the Results.
When examined together, we find a consistent picture emerging from the lists of significant KEGG pathways, GO terms, and experimental sets. For example, KEGG focal adhesion and ECM receptor interaction pathways (Table 1), GO homophilic cell adhesion (Table 3) and extracellular matrix structural constituent (not shown) groups consistently show the relevance of cell adhesion and extracellular matrix in osteoblast differentiation and mineralization. GO immune response groups (Table 3) echoes KEGG cytokine-cytokine receptor interaction pathway (Table 1). Similarly, significant experimental target genes sets (Table 4) closely reflected changes in the regulatory KEGG pathways (Table 1) or GO processes/groups (Table 3).
Interestingly, there are discrepancies among the clustering and the significant KEGG pathways, GO terms, and experimental sets. For example, Notch signaling is defined as both a KEGG pathway and a GO process. This KEGG pathway is not significant (not shown) but this GO process is (Figure 5). This discrepancy arises from two sources: (1) different definitions, i.e. KEGG pathways contain partially different set of genes from corresponding GO processes. While KEGG pathways tend to cover the whole homeostatic signal transmission systems even across multiple transcriptional cycles, GO usually covers one or multiple discrete steps or functional groups for a process. KEGG and GO definitions can be considered complementary and both provide valuable gene sets for our analysis. (2) GAGE  treats KEGG pathways and GO term gene sets differently: genes under a GO term are taken as a group coregulated towards a single direction, either all up or all down regulated, whereas genes in a KEGG pathway are frequently not coregulated and expression changes in both directions are counted. Timing discrepancies exist between experimental sets and corresponding KEGG pathways. For example, IFN positive target sets (only IFN beta shown) are not significant at 24-96 hours and MYB target set not at 96 hours (Figure 6) while Jak-STAT pathway and Wnt signaling pathway are significantly perturbed all the time (Figure 2). This can be explained by the fact that two-directional perturbation treatment for KEGG pathway does not account for direction or net effect of the perturbation, whether inhibited, activated or no overall effect. In the other hand, GO term analysis has no such issue, and IRS negative set and corresponding GO IGF receptor binding group are both significant all the time.
In this work, we took a systems approach in studying MSC differentiation. We combined experimental and computational work to reconstruct a unified picture of BMP6 signaling. The same set of experimental design and computational approach could be used to study other physiological and pathological processes, such as the differentiation of other cell types or tumorogenesis.
Cell culture and BMP6 osteogenic induction
Passage 5 human MSC (5 × 105) were plated in 24-well dishes and cultured for 3 days. The cells were subsequently placed in serum free media supplemented with ITS for 24 hours. BMP6 was then added for the pre-defined time periods and then removed as shown in Figure 1. To remove BMP6, cells were rinsed 2 times before fresh media without BMP6 was added. Ascorbate and b-glycerolphosphate were added 4 days after the initiation of BMP treatment. Cells were harvested at the indicated time points. To quantifying mineralization, the plates were stained with Alizarin Red S  14 days after the initiation of BMP treatment.
Microarray experiment and analysis
Human MSC cells underwent osteogenic induction with BMP6 treatment for 0 hours, 8 hours, 24 hours, and 96 hours, which correspond to four phenotypic groups, i.e. control, preosteoblast (no mineralization), (sub-maximal) mineralization, and maximal mineralization at 14 days after the initiation of BMP treatment (18 days in total). Cells were harvested at 8 hours, 24 hours, 96 hours and 10 days for microarray profiling using Affymetrix U133 plus GeneChip® platform. Assays were run in duplicate for a total of 20 arrays.
The raw data were processed by using FRAMS  with up-to-date probe set definition . GAGE  was applied to infer the most differentially expressed pathways or gene sets between the BMP6 added or withdrawn samples under comparison at different time points (Figure 1). As is shown in Figure 1, the comparisons at 8, 24 and 96 hours were between two sample conditions, whereas comparison at 10 days was between the 24 and 96 hours BMP6 groups versus 0 and 8 hours BMP6 groups.
Pathway analysis using GAGE
We use GAGE , Generally Applicable Gene set Enrichment, a novel method we developed for gene set/pathway analysis. The GAGE method is implemented in R in the "gage" package, available through Bioconductor at http://bioconductor.org/packages/release/bioc/html/gage.html. As the GAGE procedure has been described in detail in the original paper , here is a brief summary of the method.
Step 1: Gene sets separation. Gene sets are derived or collected from KEGG pathways , GO  and MSigDB  databases, as KEGG pathways, GO terms and experimental sets respectively. GAGE treats KEGG pathways differently from GO terms and experimental sets: member genes for a GO term or experimental set are taken as a group coregulated towards a single direction, either all up or all down regulated, whereas genes in a KEGG pathway are frequently not coregulated and expression changes in both directions are counted. We test for one-directional changes in GO terms because common function is frequently connected to coregulation and coexpression [52, 53], which has also been extensively used in classical GO analysis tools [13, 54, 55]. Such treatment will miss some pathway-like GO terms which are significantly perturbed in two directions, but we intended to have GO analysis as a complementary work to the KEGG pathway analysis. Nonetheless, GAGE provides two-directional test options for all types of gene sets. For GO analysis, all terms are included without differentiating categories, i.e. Biological Process, Cellular Component, Molecular Function, and without considering the hierarchical organization of the ontology tree.
Step 2: One-on-one comparisons. Instead of comparing BMP6 treated samples vs controls as two groups, GAGE does one-on-one comparison between samples from the two groups at a time. For each one-on-one comparison, log based fold changes are calculated for all genes. GAGE conducts two-sample t-test on the average fold change in specific gene sets against that for the background of the whole set. This one-on-one comparison procedure is repeated for all potential experiment-control pairs.
Step 3: Summarization. For each gene set, GAGE derives a global P-value based on a meta test on the negative log sum of multiple P-values for this set from all one-on-one comparisons between experiments and controls.
Overlap between significant pathways
Perturbation pattern visualization
For significant KEGG pathways, we generated graphs to visualize the dynamic expression perturbation at two levels: individual genes (Figure 3 and Additional file 2,3,4,5) and whole pathways (Figure 4). Gene expression level fold changes are standardized over the standard deviation of fold changes for all genes. The standardized fold changes for individual genes in KEGG pathways are visualized by using KEGGanim web tool  in Figure 3 and Additional file 2,3,4,5. We present an integrated network to show the connections between pathways and average expression perturbations for them in Figure 4. Significant KEGG pathways or closely related other pathways connected to them are represented by nodes. Connections between these pathways collected from KEGG database and graphs  are represented by arrows plus edges. The mean t-statistics from two-sample t-test from multiple one-on-one comparisons are plotted in pseudo heat colors as average perturbation magnitude for significant KEGG pathways.
This work was supported by NIH grant R01 DE017471. In addition, PJW is supported by NIH grant U54-DA-021519 and KDH is supported by NIH grants R01 AR054714 and R01 AR049682.
- Harada S, Rodan GA: Control of osteoblast function and regulation of bone mass. Nature. 2003, 423: 349-55. 10.1038/nature01660View ArticlePubMedGoogle Scholar
- Chen D, Zhao M, Mundy GR: Bone morphogenetic proteins. Growth Factors. 2004, 22: 233-41. 10.1080/08977190412331279890View ArticlePubMedGoogle Scholar
- Li X, Cao X: BMP signaling and skeletogenesis. Ann N Y Acad Sci. 2006, 1068: 26-40. 10.1196/annals.1346.006View ArticlePubMedGoogle Scholar
- Lavery KS, Swain PM, Falb D, Alaoui-Ismaili MH: BMP-2/4 and BMP-6/7 differentially utilize cell surface receptors to induce osteoblastic differentiation of human bone marrow derived mesenchymal stem cells. J Biol Chem. 2008Google Scholar
- Friedman MS, Long MW, Hankenson KD: Osteogenic differentiation of human mesenchymal stem cells is regulated by bone morphogenetic protein-6. J Cell Biochem. 2006, 98: 538-54. 10.1002/jcb.20719View ArticlePubMedGoogle Scholar
- de Jong DS, Vaes BL, Dechering KJ, Feijen A, Hendriks JM, Wehrens R, Mummery CL, van Zoelen EJ, Olijve W, Steegenga WT: Identification of novel regulators associated with early-phase osteoblast differentiation. J Bone Miner Res. 2004, 19: 947-58. 10.1359/JBMR.040216View ArticlePubMedGoogle Scholar
- Peng Y, Kang Q, Cheng H, Li X, Sun MH, Jiang W, Luu HH, Park JY, Haydon RC, He TC: Transcriptional characterization of bone morphogenetic proteins (BMPs)-mediated osteogenic signaling. J Cell Biochem. 2003, 90: 1149-65. 10.1002/jcb.10744View ArticlePubMedGoogle Scholar
- Harris SE, Guo D, Harris MA, Krishnaswamy A, Lichtler A: Transcriptional regulation of BMP-2 activated genes in osteoblasts using gene expression microarray analysis: role of Dlx2 and Dlx5 transcription factors. Front Biosci. 2003, 8: s1249-65. 10.2741/1170View ArticlePubMedGoogle Scholar
- Korchynskyi O, Dechering KJ, Sijbers AM, Olijve W, ten Dijke P: Gene array analysis of bone morphogenetic protein type I receptor-induced osteoblast differentiation. J Bone Miner Res. 2003, 18: 1177-85. 10.1359/jbmr.2003.18.7.1177View ArticlePubMedGoogle Scholar
- Kalajzic I, Staal A, Yang WP, Wu Y, Johnson SE, Feyen JH, Krueger W, Maye P, Yu F, Zhao Y, Kuo L, Gupta RR, Achenie LE, Wang HW, Shin DG, Rowe DW: Expression profile of osteoblast lineage at defined stages of differentiation. J Biol Chem. 2005, 280: 24618-26. 10.1074/jbc.M413834200View ArticlePubMedGoogle Scholar
- de Jong DS, Steegenga WT, Hendriks JM, van Zoelen EJ, Olijve W, Dechering KJ: Regulation of Notch signaling genes during BMP2-induced differentiation of osteoblast precursor cells. Biochem Biophys Res Commun. 2004, 320: 100-7. 10.1016/j.bbrc.2004.05.150View ArticlePubMedGoogle Scholar
- Balint E, Lapointe D, Drissi H, van der Meijden C, Young DW, van Wijnen AJ, Stein JL, Stein GS, Lian JB: Phenotype discovery by gene expression profiling: mapping of biological processes linked to BMP-2-mediated osteoblast differentiation. J Cell Biochem. 2003, 89: 401-26. 10.1002/jcb.10515View ArticlePubMedGoogle Scholar
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102: 15545-50. 10.1073/pnas.0506580102PubMed CentralView ArticlePubMedGoogle Scholar
- Kelder T, Conklin BR, Evelo CT, Pico AR: Finding the right questions: exploratory pathway analysis to enhance biological discovery in large datasets. PLoS Biol. 2010, 8: e1000472- 10.1371/journal.pbio.1000472PubMed CentralView ArticlePubMedGoogle Scholar
- Davies MN, Meaburn EL, Schalkwyk LC: Gene set enrichment; a problem of pathways. Brief Funct Genomics. 9: 385-90.
- Nam D, Kim SY: Gene-set approach for expression pattern analysis. Brief Bioinform. 2008, 9: 189-97. 10.1093/bib/bbn001View ArticlePubMedGoogle Scholar
- Kleemann R, van Erk M, Verschuren L, van den Hoek AM, Koek M, Wielinga PY, Jie A, Pellis L, Bobeldijk-Pastorova I, Kelder T, Toet K, Wopereis S, Cnubben N, Evelo C, van Ommen B, Kooistra T: Time-resolved and tissue-specific systems analysis of the pathogenesis of insulin resistance. PLoS One. 5: e8817-
- Luo W, Friedman MS, Shedden K, Hankenson KD, Woolf PJ: GAGE: generally applicable gene set enrichment for pathway analysis. BMC Bioinformatics. 2009, 10: 161- 10.1186/1471-2105-10-161PubMed CentralView ArticlePubMedGoogle Scholar
- Gosselet FP, Magnaldo T, Culerrier RM, Sarasin A, Ehrhart JC: BMP2 and BMP6 control p57(Kip2) expression and cell growth arrest/terminal differentiation in normal primary human epidermal keratinocytes. Cell Signal. 2007, 19: 731-9. 10.1016/j.cellsig.2006.09.006View ArticlePubMedGoogle Scholar
- Kim JB, Leucht P, Luppen CA, Park YJ, Beggs HE, Damsky CH, Helms JA: Reconciling the roles of FAK in osteoblast differentiation, osteoclast remodeling, and bone regeneration. Bone. 2007, 41: 39-51. 10.1016/j.bone.2007.01.024PubMed CentralView ArticlePubMedGoogle Scholar
- Salasznyk RM, Klees RF, Williams WA, Boskey A, Plopper GE: Focal adhesion kinase signaling pathways regulate the osteogenic differentiation of human mesenchymal stem cells. Exp Cell Res. 2007, 313: 22-37. 10.1016/j.yexcr.2006.09.013PubMed CentralView ArticlePubMedGoogle Scholar
- Ulsamer A, Ortuno MJ, Ruiz S, Susperregui AR, Osses N, Rosa JL, Ventura F: BMP-2 induces Osterix expression through up-regulation of Dlx5 and its phosphorylation by p38. J Biol Chem. 2008, 283: 3816-26.View ArticlePubMedGoogle Scholar
- Celil AB, Campbell PG: BMP-2 and insulin-like growth factor-I mediate Osterix (Osx) expression in human mesenchymal stem cells via the MAPK and protein kinase D signaling pathways. J Biol Chem. 2005, 280: 31353-9. 10.1074/jbc.M503845200View ArticlePubMedGoogle Scholar
- KEGG: Kyoto Encyclopedia of Genes and Genomes. http://www.genome.jp/kegg/
- Nobta M, Tsukazaki T, Shibata Y, Xin C, Moriishi T, Sakano S, Shindo H, Yamaguchi A: Critical regulation of bone morphogenetic protein-induced osteoblastic differentiation by Delta1/Jagged1-activated Notch1 signaling. J Biol Chem. 2005, 280: 1584-8.View ArticleGoogle Scholar
- Zhao G, Monier-Faugere MC, Langub MC, Geng Z, Nakayama T, Pike JW, Chernausek SD, Rosen CJ, Donahue LR, Malluche HH, Fagin JA, Clemens TL: Targeted overexpression of insulin-like growth factor I to osteoblasts of transgenic mice: increased trabecular bone volume without increased osteoblast proliferation. Endocrinology. 2000, 141: 2674-82. 10.1210/en.141.7.2674PubMedGoogle Scholar
- Zhang M, Xuan S, Bouxsein ML, von Stechow D, Akeno N, Faugere MC, Malluche H, Zhao G, Rosen CJ, Efstratiadis A, Clemens TL: Osteoblast-specific knockout of the insulin-like growth factor (IGF) receptor gene reveals an essential role of IGF signaling in bone matrix mineralization. J Biol Chem. 2002, 277: 44005-12. 10.1074/jbc.M208265200View ArticlePubMedGoogle Scholar
- Strohbach C, Kleinman S, Linkhart T, Amaar Y, Chen ST, Mohan S, Strong D: Potential involvement of the interaction between insulin-like growth factor binding protein (IGFBP)-6 and LIM mineralization protein (LMP)-1 in regulating osteoblast differentiation. J Cell Biochem. 2008Google Scholar
- Mukherjee A, Rotwein P: Insulin-like growth factor-binding protein-5 inhibits osteoblast differentiation and skeletal growth by blocking insulin-like growth factor actions. Mol Endocrinol. 2008, 22: 1238-50. 10.1210/me.2008-0001PubMed CentralView ArticlePubMedGoogle Scholar
- Silha JV, Mishra S, Rosen CJ, Beamer WG, Turner RT, Powell DR, Murphy LJ: Perturbations in bone formation and resorption in insulin-like growth factor binding protein-3 transgenic mice. J Bone Miner Res. 2003, 18: 1834-41. 10.1359/jbmr.2003.18.10.1834View ArticlePubMedGoogle Scholar
- Zhang M, Faugere MC, Malluche H, Rosen CJ, Chernausek SD, Clemens TL: Paracrine overexpression of IGFBP-4 in osteoblasts of transgenic mice decreases bone turnover and causes global growth retardation. J Bone Miner Res. 2003, 18: 836-43. 10.1359/jbmr.2003.18.5.836View ArticlePubMedGoogle Scholar
- Arron JR, Choi Y: Bone versus immune system. Nature. 2000, 408: 535-6. 10.1038/35046196View ArticlePubMedGoogle Scholar
- Walsh MC, Kim N, Kadono Y, Rho J, Lee SY, Lorenzo J, Choi Y: Osteoimmunology: interplay between the immune system and bone metabolism. Annu Rev Immunol. 2006, 24: 33-63. 10.1146/annurev.immunol.24.021605.090646View ArticlePubMedGoogle Scholar
- Arai F, Ohneda O, Miyamoto T, Zhang XQ, Suda T: Mesenchymal stem cells in perichondrium express activated leukocyte cell adhesion molecule and participate in bone marrow formation. J Exp Med. 2002, 195: 1549-63. 10.1084/jem.20011700PubMed CentralView ArticlePubMedGoogle Scholar
- Bennett KP, Bergeron C, Acar E, Klees RF, Vandenberg SL, Yener B, Plopper GE: Proteomics reveals multiple routes to the osteogenic phenotype in mesenchymal stem cells. BMC Genomics. 2007, 8: 380- 10.1186/1471-2164-8-380PubMed CentralView ArticlePubMedGoogle Scholar
- Ebert R, Jakob F: Selenium deficiency as a putative risk factor for osteoporosis. International Congress Series: Nutritional Aspects of Osteoporosis 2006. Proceedings of the 6th International Symposium on Nutritional Aspects of Osteoporosis, 4-6 May 2006, Lausanne, Switzerland. 2007, 1297: 158-164.View ArticleGoogle Scholar
- Kanei-Ishii C, Ninomiya-Tsuji J, Tanikawa J, Nomura T, Ishitani T, Kishida S, Kokura K, Kurahashi T, Ichikawa-Iwata E, Kim Y, Matsumoto K, Ishii S: Wnt-1 signal induces phosphorylation and degradation of c-Myb protein via TAK1, HIPK2, and NLK. Genes Dev. 2004, 18: 816-29. 10.1101/gad.1170604PubMed CentralView ArticlePubMedGoogle Scholar
- Nair M, Bilanchone V, Ortt K, Sinha S, Dai X: Ovol1 represses its own transcription by competing with transcription activator c-Myb and by recruiting histone deacetylase activity. Nucleic Acids Res. 2007, 35: 1687-97. 10.1093/nar/gkl1141PubMed CentralView ArticlePubMedGoogle Scholar
- Garcia-Pedrero JM, Kiskinis E, Parker MG, Belandia B: The SWI/SNF chromatin remodeling subunit BAF57 is a critical regulator of estrogen receptor function in breast cancer cells. J Biol Chem. 2006, 281: 22656-64. 10.1074/jbc.M602561200View ArticlePubMedGoogle Scholar
- Villagra A, Cruzat F, Carvallo L, Paredes R, Olate J, van Wijnen AJ, Stein GS, Lian JB, Stein JL, Imbalzano AN, Montecino M: Chromatin remodeling and transcriptional activity of the bone-specific osteocalcin gene require CCAAT/enhancer-binding protein beta-dependent recruitment of SWI/SNF activity. J Biol Chem. 2006, 281: 22695-706. 10.1074/jbc.M511640200View ArticlePubMedGoogle Scholar
- Colland F, Jacq X, Trouplin V, Mougin C, Groizeleau C, Hamburger A, Meil A, Wojcik J, Legrain P, Gauthier JM: Functional proteomics mapping of a human signaling pathway. Genome Res. 2004, 14: 1324-32. 10.1101/gr.2334104PubMed CentralView ArticlePubMedGoogle Scholar
- Simone C, Forcales SV, Hill DA, Imbalzano AN, Latella L, Puri PL: p38 pathway targets SWI-SNF chromatin-remodeling complex to muscle-specific loci. Nat Genet. 2004, 36: 738-43. 10.1038/ng1378View ArticlePubMedGoogle Scholar
- Lutwyche JK, Keough RA, Hunter J, Coles LS, Gonda TJ: DNA binding-independent transcriptional activation of the vascular endothelial growth factor gene (VEGF) by the Myb oncoprotein. Biochem Biophys Res Commun. 2006, 344: 1300-7. 10.1016/j.bbrc.2006.04.045View ArticlePubMedGoogle Scholar
- Itoh F, Itoh S, Goumans MJ, Valdimarsdottir G, Iso T, Dotto GP, Hamamori Y, Kedes L, Kato M, ten Dijke Pt P: Synergy and antagonism between Notch and BMP receptor signaling pathways in endothelial cells. Embo J. 2004, 23: 541-51. 10.1038/sj.emboj.7600065PubMed CentralView ArticlePubMedGoogle Scholar
- Dahlqvist C, Blokzijl A, Chapman G, Falk A, Dannaeus K, Ibanez CF, Lendahl U: Functional Notch signaling is required for BMP4-induced inhibition of myogenic differentiation. Development. 2003, 130: 6089-99. 10.1242/dev.00834View ArticlePubMedGoogle Scholar
- Takizawa T, Ochiai W, Nakashima K, Taga T: Enhanced gene activation by Notch and BMP signaling cross-talk. Nucleic Acids Res. 2003, 31: 5723-31. 10.1093/nar/gkg778PubMed CentralView ArticlePubMedGoogle Scholar
- Deregowski V, Gazzerro E, Priest L, Rydziel S, Canalis E: Notch 1 overexpression inhibits osteoblastogenesis by suppressing Wnt/beta-catenin but not bone morphogenetic protein signaling. J Biol Chem. 2006, 281: 6203-10. 10.1074/jbc.M508370200View ArticlePubMedGoogle Scholar
- Leroith D, Nissley P: Knock your SOCS off!. J Clin Invest. 2005, 115: 233-6.PubMed CentralView ArticlePubMedGoogle Scholar
- Hochreiter S, Clevert DA, Obermayer K: A new summarization method for Affymetrix probe level data. Bioinformatics. 2006, 22: 943-9. 10.1093/bioinformatics/btl033View ArticlePubMedGoogle Scholar
- Gene Ontology. http://www.geneontology.org/
- The Molecular Signature Database. http://www.broad.mit.edu/gsea/msigdb/
- Horan K, Jang C, Bailey-Serres J, Mittler R, Shelton C, Harper JF, Zhu JK, Cushman JC, Gollery M, Girke T: Annotating genes of known and unknown function by large-scale coexpression analysis. Plant Physiol. 2008, 147: 41-57. 10.1104/pp.108.117366PubMed CentralView ArticlePubMedGoogle Scholar
- Allocco DJ, Kohane IS, Butte AJ: Quantifying the relationship between co-expression, co-regulation and gene function. BMC Bioinformatics. 2004, 5: 18- 10.1186/1471-2105-5-18PubMed CentralView ArticlePubMedGoogle Scholar
- Al-Shahrour F, Diaz-Uriarte R, Dopazo J: FatiGO: a web tool for finding significant associations of Gene Ontology terms with groups of genes. Bioinformatics. 2004, 20: 578-80. 10.1093/bioinformatics/btg455View ArticlePubMedGoogle Scholar
- Zhang B, Schmoyer D, Kirov S, Snoddy J: GOTree Machine (GOTM): a web-based platform for interpreting sets of interesting genes using Gene Ontology hierarchies. BMC Bioinformatics. 2004, 5: 16- 10.1186/1471-2105-5-16PubMed CentralView ArticlePubMedGoogle Scholar
- Adler P, Reimand J, Janes J, Kolde R, Peterson H, Vilo J: KEGGanim: pathway animations for high-throughput data. Bioinformatics. 2008, 24: 588-90. 10.1093/bioinformatics/btm581View ArticlePubMedGoogle Scholar
- Hughes FJ, Collyer J, Stanfield M, Goodman SA: The effects of bone morphogenetic protein-2, -4, and -6 on differentiation of rat osteoblast cells in vitro. Endocrinology. 1995, 136: 2671-7. 10.1210/en.136.6.2671PubMedGoogle Scholar
- Takayanagi H, Sato K, Takaoka A, Taniguchi T: Interplay between interferon and other cytokine systems in bone metabolism. Immunol Rev. 2005, 208: 181-93. 10.1111/j.0105-2896.2005.00337.xView ArticlePubMedGoogle Scholar
- Fujita K, Janz S: Attenuation of WNT signaling by DKK-1 and -2 regulates BMP2-induced osteoblast differentiation and expression of OPG, RANKL and M-CSF. Mol Cancer. 2007, 6: 71- 10.1186/1476-4598-6-71PubMed CentralView ArticlePubMedGoogle Scholar
- Rawadi G, Vayssiere B, Dunn F, Baron R, Roman-Roman S: BMP-2 controls alkaline phosphatase expression and osteoblast mineralization by a Wnt autocrine loop. J Bone Miner Res. 2003, 18: 1842-53. 10.1359/jbmr.2003.18.10.1842View ArticlePubMedGoogle Scholar
- Rajan P, Panchision DM, Newell LF, McKay RD: BMPs signal alternately through a SMAD or FRAP-STAT pathway to regulate fate choice in CNS stem cells. J Cell Biol. 2003, 161: 911-21. 10.1083/jcb.200211021PubMed CentralView ArticlePubMedGoogle Scholar
- Wang X, Kua HY, Hu Y, Guo K, Zeng Q, Wu Q, Ng HH, Karsenty G, de Crombrugghe B, Yeh J, Li B: p53 functions as a negative regulator of osteoblastogenesis, osteoblast-dependent osteoclastogenesis, and bone remodeling. J Cell Biol. 2006, 172: 115-25. 10.1083/jcb.200507106PubMed CentralView ArticlePubMedGoogle Scholar
- Lengner CJ, Steinman HA, Gagnon J, Smith TW, Henderson JE, Kream BE, Stein GS, Lian JB, Jones SN: Osteoblast differentiation and skeletal development are regulated by Mdm2-p53 signaling. J Cell Biol. 2006, 172: 909-21. 10.1083/jcb.200508130PubMed CentralView ArticlePubMedGoogle Scholar
- Xiao G, Gopalakrishnan R, Jiang D, Reith E, Benson MD, Franceschi RT: Bone morphogenetic proteins, extracellular matrix, and mitogen-activated protein kinase signaling pathways are required for osteoblast-specific gene expression and differentiation in MC3T3-E1 cells. J Bone Miner Res. 2002, 17: 101-10. 10.1359/jbmr.2002.17.1.101View ArticlePubMedGoogle Scholar
- Xiao G, Wang D, Benson MD, Karsenty G, Franceschi RT: Role of the alpha2-integrin in osteoblast-specific gene expression and activation of the Osf2 transcription factor. J Biol Chem. 1998, 273: 32988-94. 10.1074/jbc.273.49.32988View ArticlePubMedGoogle Scholar