TGFβ-induced switch from adipogenic to osteogenic differentiation of human mesenchymal stem cells: identification of drug targets for prevention of fat cell differentiation
© The Author(s). 2016
Received: 7 January 2016
Accepted: 25 July 2016
Published: 26 August 2016
Patients suffering from osteoporosis show an increased number of adipocytes in their bone marrow, concomitant with a reduction in the pool of human mesenchymal stem cells (hMSCs) that are able to differentiate into osteoblasts, thus leading to suppressed osteogenesis.
In order to be able to interfere with this process, we have investigated in-vitro culture conditions whereby adipogenic differentiation of hMSCs is impaired and osteogenic differentiation is promoted. By means of gene expression microarray analysis, we have investigated genes which are potential targets for prevention of fat cell differentiation.
Our data show that BMP2 promotes both adipogenic and osteogenic differentiation of hMSCs, while transforming growth factor beta (TGFβ) inhibits differentiation into both lineages. However, when cells are cultured under adipogenic differentiation conditions, which contain cAMP-enhancing agents such as IBMX of PGE2, TGFβ promotes osteogenic differentiation, while at the same time inhibiting adipogenic differentiation. Gene expression and immunoblot analysis indicated that IBMX-induced suppression of HDAC5 levels plays an important role in the inhibitory effect of TGFβ on osteogenic differentiation. By means of gene expression microarray analysis, we have investigated genes which are downregulated by TGFβ under adipogenic differentiation conditions and may therefore be potential targets for prevention of fat cell differentiation. We thus identified nine genes for which FDA-approved drugs are available. Our results show that drugs directed against the nuclear hormone receptor PPARG, the metalloproteinase ADAMTS5, and the aldo-keto reductase AKR1B10 inhibit adipogenic differentiation in a dose-dependent manner, although in contrast to TGFβ they do not appear to promote osteogenic differentiation.
The approach chosen in this study has resulted in the identification of new targets for inhibition of fat cell differentiation, which may not only be relevant for prevention of osteoporosis, but also of obesity.
Human mesenchymal stem cells (hMSCs) from bone marrow have the ability to differentiate into cells from multiple lineages, including osteoblasts and adipocytes. The commitment of hMSCs towards either the osteogenic or adipogenic lineage depends on the local availability of growth factors and hormones, which are able to activate lineage-specific transcriptional regulators . Patients suffering from osteoporosis show an increased number of adipocytes in their bone marrow, concomitant with a reduction in the pool of hMSCs that are able to differentiate into osteoblasts, thus leading to suppressed osteogenesis [2, 3]. It is still unclear to what extent this age-related increase in differentiation of hMSCs towards adipocytes results from intrinsic changes in the stem cells or from alterations in the microenvironment of the bone marrow [1, 4]. These observations have recently stirred increasing interest in anabolic therapies for osteoporosis, whereby osteogenic differentiation of hMSCs is stimulated by preventing adipogenic differentiation [5, 6].
Most information about the signaling pathways that are required for osteogenic and adipogenic differentiation of hMSCs has come from in-vitro studies, whereby cells are treated with specific combinations of growth factors and hormones . Differentiation into both lineages requires treatment of monolayer cells with dexamethasone (DEX) and is enhanced by the presence of bone morphogenetic proteins (BMPs). Osteogenic differentiation is obtained by the additional presence of ascorbate and β-glycerophosphate, whereas adipogenic differentiation requires treatment with insulin, 3-isobutyl-1-methylxanthine (IBMX), and the PPARG activator rosiglitazone. Our previous data have indicated that under these experimental conditions at least 75 % of the hMSCs differentiate into either osteoblasts or adipocytes .
Activation of the RUNX2 nuclear transcription factor appears to be essential for the osteogenic pathway of hMSCs, while the adipogenic pathway requires the transcriptional activity of the PPARG nuclear hormone receptor . Evidence has been presented that transcriptional regulators promoting differentiation into one of these lineages actively suppress differentiation into the other lineage . It has therefore been postulated that a reciprocal relationship exists between the osteogenic and adipogenic pathways, implying that impaired adipogenic differentiation of hMSCs may result in enhanced osteogenic differentiation .
Multiple regulators are known that affect the choice between the osteogenic and adipogenic lineage. Most notably, activation of the WNT/β-catenin pathway promotes osteogenic differentiation and inhibits adipogenic differentiation of hMSCs . BMP enhances the outgrowth of both osteoblasts and adipocytes, while the related cytokine transforming growth factor beta (TGFβ) inhibits both osteogenic and adipogenic differentiation, at least when tested under these in-vitro conditions [10, 11]. Moreover, Kim et al.  have shown that cAMP-activated protein kinases regulate the differentiation choice of hMSCs between osteogenesis and adipogenesis.
In order to find drugs affecting the choice between osteogenic and adipogenic differentiation of hMSCs, we have optimized the adipogenic culture conditions in such a way that, within the same culture, a fraction of the cells differentiates into osteoblasts and another fraction into adipocytes. Our results show that addition of TGFβ to these cultures fully blocks adipogenic differentiation but, surprisingly, enhances osteogenic differentiation. Analysis of the individual components in the culture medium showed that the presence of the phosphodiesterase inhibitor IBMX, which stabilizes cAMP levels in the cell, converted TGFβ from an inhibitor to an enhancer of osteogenic differentiation.
Based on these observations we have set up a gene expression microarray experiment to identity genes that under these optimized adipogenic culture conditions are downregulated by TGFβ. The genes identified in this way seem to play an important role in adipogenesis, since inhibitors of the corresponding proteins were found to be able to block adipogenic differentiation of hMSCs. Some of these inhibitors have already received FDA approval for the treatment of various diseases, and may therefore be good candidates for therapeutic drug repurposing in order to treat patients suffering from such diseases as obesity and osteoporosis.
Culture and differentiation of hMSCs
hMSCs harvested from normal human bone marrow were purchased from Lonza (Walkersville, MD, USA) at passage 2. Cells were expanded for no more than five passages in “mesenchymal stem cell growth medium” (MSCGM; Lonza) at 37 °C in a humidified atmosphere containing 7.5 % CO2. Studies were performed with hMSCs from three different donors, encoded 5F0138, 6F4085, and 7F3458.
For differentiation experiments, 4.0 × 104 cells per cm2 were seeded in high-glucose-containing Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10 % fetal bovine serum (a selected lot from Lonza), 100 U/ml penicillin, and 100 μg/ml streptomycin. This medium will be further referred to as “proliferation medium” (PM). The next day, cells were switched to either osteogenic or adipogenic differentiation medium. Cells treated with PM were used as negative controls. Osteogenic differentiation medium (ODM) was composed of PM supplemented with 10−7 M DEX, 0.2 mM ascorbate, and 10 mM β-glycerophosphate. Adipogenic differentiation medium (ADM) was composed of PM supplemented with 10−6 M DEX, 10 μg/ml insulin (R&D Systems, Minneapolis, MN, USA), 10−7 M rosiglitazone (Sigma-Aldrich, St. Louis, MO, USA), and 500 μM IBMX (Sigma-Aldrich). Unless indicated otherwise, recombinant human TGFβ1 and BMP2 (both from R&D Systems) were used at concentrations of 2 ng/ml and 125 ng/ml, respectively. Prostaglandin E2 (PGE2), Batimastat, and Zopolrestat (Sigma-Aldrich) were used at the indicated concentrations. Media were refreshed every 3–4 days.
Alkaline phosphatase assays
To quantify alkaline phosphatase (ALP) enzymatic activity as a measure of osteogenic differentiation, hMSCs were seeded in 96-well tissue culture plates as already described, after which cells were differentiated for 7 days in osteogenic differentiation medium. ALP enzymatic activity was quantified by measuring the formation of p-nitrophenol from p-nitrophenyl phosphate (PNPP; Sigma-Aldrich), as described previously . ALP enzymatic activity was corrected for differences in cell number, as determined by a Neutral Red assay. Cells were incubated with Neutral Red dye diluted in PBS for 1 hour at 37 °C. After washing with PBS, the dye was extracted from the cells using 0.05 M NaH2PO4 in 50 % EtOH, after which the absorbance was measured at 540 nm.
For histochemical analysis of ALP activity, cells were seeded in 48-well tissue culture plates and differentiated for 7 days in osteogenic differentiation medium. Subsequently, cells were fixed in 3.7 % formaldehyde/PBS for 10 min at 22 °C. After washing with PBS, cells were incubated for 1 hour at 37 °C in a mixture of 0.1 mg/ml naphtol AS-MX phosphate (Sigma-Aldrich), 0.5 % N,N-dimethylformamide, 2 mM MgCl2, and 0.6 mg/ml Fast Blue BB salt (Sigma-Aldrich) in 0.1 M Tris–HCl, pH 8.5.
To measure calcium deposition in the extracellular matrix, hMSCs were seeded in 24-well tissue culture plates and cultured for 13 days under osteogenic or adipogenic differentiation conditions, as indicated. Cells were subsequently washed twice with PBS after which calcium was extracted from the extracellular matrix by treatment with 150 μl of 0.5 M HCl. Calcium concentrations were measured in a colorimetric assay using o-cresolphtalein complexone as a chromogenic agent, according to the protocol provided by the manufacturer (Sigma-Aldrich).
Oil Red O staining and Triglyceride assay
To quantify adipogenic differentiation, lipid droplets were stained in mature adipocytes obtained after treatment of hMSCs for 9 days in adipogenic differentiation medium. Cells were first washed twice with PBS, fixed for 30 min with 1 % formaldehyde in PBS, and then washed once with water and twice with 60 % isopropanol. Cells were then stained for 1 hour with 0.3 % w/v Oil Red O (Sigma-Aldrich) in 60 % isopropanol. Subsequently, cells were washed once with 60 % isopropanol and twice with distilled water. For quantification of Oil Red O staining, samples were treated with 100 % isopropanol and absorbance was measured at 530 nm.
The amount of triglycerides stored in lipid droplets of mature adipocytes was quantified after treatment of hMSCs for 9 days in adipogenic differentiation medium. Cells in 96-well tissue culture plates were washed twice with PBS. Triglycerides were extracted from the lipid droplets by freezing the cells in 50 μl of a buffer containing 25 mM Tris–HCl (pH 7.5) and 1 mM EDTA, followed by addition of 40 μl tert-butanol and 10 μl methanol. Samples were heat-dried at 55 °C, after which they were resuspended in Triglycerides LiquiColor® mono reagents (HUMAN GmbH, Wiesbaden, Germany). Triglycerides were quantified by measuring the absorbance at 490 nm.
RNA isolation and real-time quantitative RT-PCR
Histone deacetylase 5 (HDAC5)-FW: 5′-ATGACAACGGGAACTTCTTTCC-3′.
Histone deacetylase 5 (HDAC5)-RV: 5′-CCATGCCACGTTCACATTGTA-3′.
Adiponectin (ADIPOQ)-FW: 5′-CCCAAAGAGGAGAGAGGAAGC-3′.
Adiponectin (ADIPOQ)-RV: 5′-GCCAGAGCAATGAGATGCAA-3′.
Alkaline phosphatase (ALPL)-FW: 5′-GATGGACAAGTTCCCCTTCGT-3′.
Alkaline phosphatase (ALPL)-RV: 5′-GGACCTGGGCATTGGTGTT-3′.
Ribosomal protein S27a (RPS27A)-FW: 5′-GTTAAGCTGGCTGTCCTGAAA-3′.
Ribosomal protein S27a (RPS27A)-RV: 5′-CATCAGAAGGGCACTCTCG-3′.
hMSCs were seeded at 4.0 × 104 cells per cm2 in six-well plates and cultured for 24 hours in the indicated differentiation media. Cells were then lysed in 250 μl of RIPA lysis buffer per well. Then 5 μl of reducing sample buffer was added to 25 μl of lysate, heated to 95 °C, and subsequently loaded onto an 8 % SDS-PAGE gel. HDAC5 was detected on blots using a goat polyclonal anti-HDAC5 antibody (G-18, 1/200 dilution) raised against the N-terminus of human HDAC5 (sc-5250; Santa Cruz Biotechnology, Dallas, TX, USA), followed by an HRP-labeled polyclonal secondary antibody. Antibodies against α-tubulin (Sigma) served as a loading control.
Gene expression microarray analysis
To identify genes that are regulated during osteogenic and adipogenic differentiation of hMSCs, a total of 54 samples (each containing 800,000 cells/ 20 cm2) were seeded in PM and grown for 24 hours. Subsequently the medium was exchanged for differentiation medium, now consisting of PM with 10−6 M DEX, 10 μg/ml insulin, 10−7 M rosiglitazone, and 50 ng/ml BMP2 (B). In addition, either 5 ng/ml TGFβ (BT), or 250 μM IBMX (BI), or 5 ng/ml TGFβ and 250 μM IBMX (BTI) were added. Samples were incubated for either 0, 1, 2, 3, or 7 days. Experiments for each group and time point were carried out as three biological replicates, while the untreated control group (time 0) consisted of six samples. RNA was isolated as already described, and hybridized onto Affymetrix HGU 133 plus 2.0 microarrays according to existing protocols .
Microarray data were analyzed with the R language for statistical computing using appropriate Bioconductor packages (http://bioconductor.org/) for reading, normalizing, and statistically evaluating the data, followed by annotation of the gene sets and integration of parallel data sources. Briefly, the analysis started with a careful quality assessment of the dataset using the automatic R pipeline AffymetrixQC , which was customized and run locally. All 54 microarrays passed the quality control and were included in the analysis, consisting of robust microarray analysis (RMA) normalization , followed by statistical analysis to find differentially expressed genes using Linear Models for Microarray Data (LIMMA) , and subsequent functional annotation and enrichment analysis using the online resource Database for Annotation, Visualization and Integrated Discovery (DAVID) [17, 18]. Finally, the list of differentially expressed genes for the contrasts of interest was crossed with the information from the DrugBank database  in order to derive the final list of candidate genes for experimental testing. All R scripts used for this analysis are available upon request. Current microarray data have been deposited in NCBI’s Gene Expression Omnibus [GEO:GSE84500] (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE84500).
Student’s t test was used for statistical comparisons. Numeric data are represented as mean ± standard deviation of triplicate experiments, unless stated otherwise.
TGFβ induces hMSCs to switch from adipogenic to osteogenic differentiation
The role of TGFβ in adipogenic and osteogenic differentiation of hMSCs is still unclear. Figure 1a also shows that adding TGFβ to adipogenic differentiation medium blocks adipogenic differentiation of hMSCs in a dose-dependent manner, both in the absence and additional presence of BMP2. Figure 1b shows that addition of TGFβ to osteogenic differentiation medium results in a similar inhibition of osteogenic differentiation, both in agreement with previous data [10, 11]. However, when hMSCs are treated with adipogenic differentiation medium (which contains a 10-fold higher concentration of DEX than osteogenic differentiation medium) in combination with BMP2, a fraction of the cells will differentiate into bone cells and a fraction into fat cells within the same well, as shown in Fig. 1c by histological staining. Addition of TGFβ under these conditions resulted is a dose-dependent increase in the number of ALP-positive bone cells, with a concomitant reduction in Oil Red O-positive fat cells. These data show that TGFβ blocks bone cell differentiation under osteogenic differentiation conditions, but enhances bone cell differentiation under adipogenic differentiation conditions. It can therefore be concluded that, under the experimental conditions used in Fig. 1c, TGFβ induces a switch from adipogenic to osteogenic differentiation.
IBMX is a critical component in the TGFβ-mediated switch in cell fate
IBMX is a phosphodiesterase inhibitor, which prevents degradation of cAMP. In order to show that enhanced cAMP levels play a role in the TGFβ-mediated switch in cell fate, we tested the effect of PGE2, a known activator of adenylate cyclase. Figure 2c shows that PGE2 is able to overcome TGFβ-induced inhibition of osteogenic differentiation of hMSCs in a dose-dependent manner, to a similar extent as IBMX.
Osteogenic differentiation requires not only matrix maturation, as indicated by alkaline phosphatase expression, but also matrix mineralization, as indicated by calcium deposition. Figure 2d shows that under osteogenic differentiation conditions BMP is required for calcium deposition by hMSCs, while both BMP and TGFβ are required under adipogenic differentiation. These data show TGFβ is able to induce fully differentiated, mineralized osteoblasts under adipogenic differentiation conditions.
IBMX suppresses HDAC5 expression
Previous studies  have shown that TGFβ-mediated inhibition of osteogenic differentiation can be overcome by trichostatin A, an inhibitor of class I and II mammalian histone deacetylases (HDACs). We confirmed this observation upon incubating hMSCs in osteogenic differentiation medium with BMP2 and TGFβ (data not shown). Kang et al.  have presented evidence that HDAC4/5 can interact with TGFβ-activated SMAD3, resulting in a complex that represses transcription of the bone marker gene osteocalcin (BGLAP). In human adipose-derived mesenchymal stem cells, trichostatin A impaired PPARG activity, at least under osteogenic differentiation conditions . Still, the role of HDACs in adipogenesis is far from clear, since HDAC inhibitors may either enhance or inhibit adipogenic differentiation .
Figure 3b shows that HDAC5 levels are similarly regulated at the protein level. Western blot analysis revealed the highest expression level in hMSCs treated with osteogenic differentiation medium containing BMP2 and TGFβ (OBT), whereas almost no protein was detected in the additional presence of IBMX (OBTI). These data indicate that under conditions whereby TGFβ enhances osteogenic differentiation, no HDAC5 is available to prevent expression of bone specific genes.
Altered gene expression during TGFβ-mediated switch in cell fate
Analysis of adipogenic differentiation inhibitors
Selected genes involved in adipogenic differentiation
Angiotensin II receptor
Bradykinin receptor B2
cAMP protein kinase
Vascular cell adhesion
Osteoporosis is a debilitating disease which affects tens of millions of people worldwide. Postmenopausal women are particularly at risk of developing this disease, since reduced estrogen activity enhances the bone-resorbing activity of osteoclasts. Additionally, aging in general results in a significant reduction in the number and function of bone-forming osteoblasts . Treatment of patients with osteoporosis may include hormone replacement therapies and the use of bisphosphonates to inhibit the activity of osteoclasts. Moreover, vitamin D3 or parathyroid hormone can be prescribed to enhance new bone formation . Recent years have seen an increasing interest in the interaction between fat and bone cells in the bone marrow. Patients suffering from osteoporosis show an increased number of adipocytes in their bone marrow, concomitant with a reduction in the pool of hMSCs that are able to differentiate into osteoblasts [2, 3]. In addition, adipokines secreted by these adipocytes may further enhance osteoclast activity , while an increasing fat content of the bone marrow may also impair the bone cell niche required for the functioning of hematopoietic stem cells .
These considerations show that there is great need for developing drugs that prevent the differentiation of hMSCs into fat cells and may thereby enhance their differentiation into bone cells. Likewise, drugs that prevent adipogenic differentiation could also be useful in the battle against obesity. In the present study we have shown that, under proper in-vitro conditions, TGFβ is able to stimulate osteogenic differentiation and prevent adipogenic differentiation within the same culture. These culture conditions require not only the presence of DEX and BMP, but also of a cAMP-enhancing stimulus such as IBMX or PGE2. Because of its strong side effects TGFβ is not suited for in-vivo applications, but our results show that the current in-vitro system can be used for identifying genes whose expression is repressed following a TGFβ-induced switch from adipogenic to osteogenic differentiation of hMSCs. By concentrating on those genes for which FDA-approved drugs are available, we have identified nine potential genes that could be tested for inhibition of adipogenic differentiation of hMSCs. Using this approach we have identified the nuclear hormone receptor PPARG, the metalloproteinase ADAMTS5, and the aldo-keto reductase AKR1B10 as potential drug targets for treatment of osteoporosis and obesity.
TGFβ is a highly pleiotropic cytokine which plays an important role in many physiological processes, as well as in cancer [28, 29]. Because of its strong inhibitory effect on the immune system, systemic treatment with TGFβ is not considered a realistic option, except for protection against autoimmune diseases . Studies in laboratory animals have shown that TGFβ treatment can result in skin fibrosis and toxicity, without displaying significant antitumor effects [31, 32]. Local injection of TGFβ into the knee joint has been shown to enhance cartilage integrity, leading to prevention of osteoarthritis . Our present study shows that the effect of TGFβ on osteogenic differentiation of hMSCs strongly depends on the culture conditions. In osteogenic differentiation medium TGFβ inhibits bone cell differentiation, while it promotes this process in adipogenic differentiation medium. Since the local conditions in the fatty bone marrow of osteoporosis patients are not well defined, it is difficult to predict whether injection of TGFβ will result in a net enhancement or a decrease of functional osteoblast cells. Furthermore, current clinical trials directed towards TGFβ are aimed at inhibiting its activity in cancer patients, since overactivation of TGFβ-induced pathways have been associated with cancer progression [28, 29].
Among the genes that were downregulated at early time points following TGFβ treatment of hMSCs in adipogenic differentiation medium is PPARG. This observation is not unexpected, since ligand-induced activation of PPARG (e.g., by rosiglitazone) is known to be essential for adipogenic differentiation. We have studied the role of PPARG in this process not by using the PPARG-specific antagonist GW9662 (see Table 1), but by omitting rosiglitazone from the culture medium, which resulted in a complete block of fat cell formation (data not shown). The promoting role of PPARG agonists in obesity is well established, but on the contrary these hormones have been shown to be clinically active as antidiabetic drugs . Because of these multiple faces of PPARG, inhibition of this nuclear hormone receptor does not seem an attractive approach for the prevention of fatty bone marrow formation in osteoporosis patients.
The second gene for which inhibitors were found to prevent fat cell differentiation is ADAMTS5, a disintegrin and metalloproteinase with thrombospondin motifs . Its main function is the cleavage of the aggrecan core protein, in which it appears more efficient than other matrix metalloproteinases (MMPs) . Recent studies on Adamts5 –/– mice, however, have indicated that ADAMTS5 may not responsible for aggrecan proteolysis, but instead regulates glucose uptake by mediating the endocytotic trafficking of LRP1 and GLUT4 . Other studies have shown that deletion of active ADAMTS5 prevents cartilage degradation in a murine model of osteoarthritis . Batimastat and the related drug Marimastat have primarily been developed as antineoplastic drugs. They prevent angiogenesis by blocking metalloproteinases including MMPs and ADAM family members . Both drugs performed poorly in clinical trials on metastatic breast cancer patients  and were therefore never marketed. Previous studies have shown that Batimastat blocks the enzymatic activity of MMP2 and MMP9, and prevents the differentiation of mouse 3T3-F442A preadipocytes into fat cells [23, 41]. Our present results show that Batimastat also prevents adipogenic differentiation of hMSCs. Long-term studies on Batimastat and Marimastat in humans are lacking and therefore it cannot be concluded whether these drugs present a real potential for the treatment of patients with osteoporosis or obesity. MMPs are not only involved in cancer progression, but also play an important role in inflammation and immunity . Multiple natural compounds have been identified which seem to block MMP activity . It will be interesting to test the effects of these compounds, some of which are used as food supplements, in relation to osteoporosis and obesity.
In this study we have made the novel observation that inhibitors of aldo-keto reductases prevent adipogenic differentiation of hMSCs. AKR1B10, which was downregulated by TGFβ in our studies, is particularly active on lipid substrates. It plays an important role in the reduction of retinaldehyde to retinol, as well as in the lipid modification of the K-Ras oncogene. In humans, this gene is particularly expressed in the intestine, adrenal gland, and liver. Inhibition of AKR1B10 prevents the outgrowth of pancreatic carcinoma cells by modulating the Ras pathway [44, 45]. Moreover, expression of AKR1B10 is considered to be a tumor marker for NSCLC  and liver tumors . To the best of our knowledge, a role of AKR1B10 in adipogenic differentiation has not yet been studied. However, an RNA-mediated knockdown study has indicated that AKR1B10 is an important regulator of fatty acid biosynthesis in human RAO-3 breast cancer cells .
AKR1B10 is structurally very similar to AKR1B1, and many inhibitors of AKR1B1 also bind AKR1B10. These include the drugs Sorbinil and Zopolrestat, which have been used in the present study. These and related FDA-approved drugs are used particularly for treatment of patients with diabetic polyneuropathy . They do so by inhibiting the metabolism of glucose by the so-called polyol pathway, which converts glucose into sorbitol. This reduced sugar accumulates in the cell and, as a result of osmotic stress, induces microvascular damage to the retina, kidney, and nerves. Although each of these drugs has its specific adverse effects, including rash, toxicity, and hypersensitivity reactions, at least some of them were well tolerated in clinical trials lasting more than a year . However, no long-term beneficial effects of these drugs were observed in diabetes patients. More recently, AKR1B1 inhibitors have also been tested as anticancer drugs . So far, no reports have been made in the literature on the effects of aldo-keto reductase inhibitors on patients with osteoporosis or obesity. Similar to observations for MMPs, multiple naturally occurring inhibitors for aldo-keto reductases have been identified, particularly from plant tissue [51, 52]. Feeding mice with the most potent of these natural AKR1B inhibitors, bisdemethoxycurcumin, has been shown to reduce particularly the incidence of intestinal cancer. Interestingly, some of these natural compounds inhibit not only AKR1B members but also MMPs. In ongoing research we are testing the effects of these natural compounds on fat cell differentiation. These studies may provide a lead towards the further development of more optimized compounds which specifically prevent adipogenic differentiation in vivo.
Using our gene expression microarray analysis we have identified novel inhibitors of adipogenic differentiation. We thereby focused on genes which were differentially expressed on days 1 and 2 following TGFβ treatment. This time frame corresponds with the upregulation of commitment genes for adipogenic differentiation such as PPARG, as shown in Fig. 5. Genes differentially expressed at later time points were particularly involved in fatty acid metabolism (data not shown). Figure 6a shows that under the experimental conditions tested inhibitors of these adipogenic commitment genes did not promote osteogenic differentiation to the same extent as TGFβ. The possibility should therefore be considered that within 24–48 hours the cells have already undergone an irreversible commitment towards either adipogenic or osteogenic differentiation. In that respect it may be interesting to study whether TGFβ added at later time points can still promote cells under adipogenic differentiation conditions to become osteoblasts. Alternatively, the possibility should be considered that under high cAMP conditions TGFβ is able to stimulate specific pathways leading to osteogenic differentiation. Our current results showing that, in the presence of cAMP-enhancing agents, TGFβ is able to prevent the formation of fat cells and promote the formation of bone cells have been obtained on commercially available mesenchymal stem cells from bone marrow of healthy human donors. Given the observation that patients suffering from osteoporosis show an increased number of adipocytes in their bone marrow, it will be interesting to carry out similar experiments with mesenchymal stem cells from such patients. This may indicate if this aberrant pattern of differentiation results from intrinsic changes in their stem cells or from an altered microenvironment in their bone marrow.
Our data show that in the presence of cAMP-enhancing agents TGFβ stimulates the ability of hMSCs to differentiate into bone cells, while impairing their ability to differentiate into fat cells. Under these conditions TGFβ treatment results in a reduced expression of genes which contribute to adipogenic differentiation, including PPARG, ADAMTS5, and AKR1B10. Since FDA-approved drugs are available for these genes, they are potential targets for treatment of patients suffering from osteoporosis or obesity.
ADM, adipogenic differentiation medium; ALP, alkaline phosphatase; BMP, bone morphogenetic protein; DAVID, Database for Annotation, Visualization and Integrated Discovery; DEX, dexamethasone; DMEM, Dulbecco’s modified Eagle’s medium; DMSO, dimethyl sulfoxide; hMSC, human mesenchymal stem cell; IBMX, 3-isobutyl-1-methylxanthine; LIMMA, Linear Models for Microarray Data; MMP, matrix metalloproteinase; ODM, osteogenic differentiation medium; PGE2, prostaglandin E2; PM, proliferation medium; PNPP, p-nitrophenyl phosphate; RMA, robust microarray analysis; TGFβ, transforming growth factor beta.
The authors would like to thank Dr Ester Piek and Ingrid de Grijs for their contributions in the initial phase of the project and Dr W. Olijve for fruitful discussions. The present studies were made possible by a grant from ERASysBio + to the LINCONET consortium (Dutch contribution obtained from ZonMw, grant 90201174). Merck Sharpe & Dohme (Oss, the Netherlands) is acknowledged for financial support in the initial phase of this project. Microarray hybridizations were carried out by the MicroArray Department at the University of Amsterdam.
EJvZ was responsible for the overall design of the study and drafted the final manuscript. ID was responsible for drafting the manuscript with respect to the microarray data and for a critical evaluation of the manuscript. JMH carried out the cell biological assays and was involved in drafting the experimental part of the manuscript. SPvdW was responsible for the design of the cell biological experiments, drafted the cell biological part of the manuscript, and made critical revisions to the final draft. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Pino AM, Rosen CJ, Rodriguez JP. In osteoporosis, differentiation of mesenchymal stem cells (MSCs) improves bone marrow adipogenesis. Biol Res. 2012;45(3):279–87.View ArticlePubMedGoogle Scholar
- Justesen J, et al. Adipocyte tissue volume in bone marrow is increased with aging and in patients with osteoporosis. Biogerontology. 2001;2(3):165–71.View ArticlePubMedGoogle Scholar
- Rosen CJ, et al. Marrow fat and the bone microenvironment: developmental, functional, and pathological implications. Crit Rev Eukaryot Gene Expr. 2009;19(2):109–24.View ArticlePubMedPubMed CentralGoogle Scholar
- Moerman EJ, et al. Aging activates adipogenic and suppresses osteogenic programs in mesenchymal marrow stroma/stem cells: the role of PPAR-gamma2 transcription factor and TGF-beta/BMP signaling pathways. Aging Cell. 2004;3(6):379–89.View ArticlePubMedPubMed CentralGoogle Scholar
- Ito H. Clinical considerations of regenerative medicine in osteoporosis. Curr Osteoporos Rep. 2014;12(2):230–4.View ArticlePubMedGoogle Scholar
- Reid IR. Short-term and long-term effects of osteoporosis therapies. Nat Rev Endocrinol. 2015;11(7):418–28.View ArticlePubMedGoogle Scholar
- Sotoca AM, Weber M, Zoelen EJJ. Gene expression regulation underlying osteo-, adipo-, and chondro-genic lineage commitment of human mesenchymal stem cells. In: Daskalaki A, editor. Medical Advancements in Aging and Regenerative Technologies: Clinical Tools and Applications. Hershey: IGI Global; 2013. p. 76–94.View ArticleGoogle Scholar
- Muruganandan S, Roman AA, Sinal CJ. Adipocyte differentiation of bone marrow-derived mesenchymal stem cells: cross talk with the osteoblastogenic program. Cell Mol Life Sci. 2009;66(2):236–53.View ArticlePubMedGoogle Scholar
- Bennett CN, et al. Regulation of osteoblastogenesis and bone mass by Wnt10b. Proc Natl Acad Sci U S A. 2005;102(9):3324–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Roelen BA, Dijke P. Controlling mesenchymal stem cell differentiation by TGFBeta family members. J Orthop Sci. 2003;8(5):740–8.View ArticlePubMedGoogle Scholar
- Rosen ED, MacDougald OA. Adipocyte differentiation from the inside out. Nat Rev Mol Cell Biol. 2006;7(12):885–96.View ArticlePubMedGoogle Scholar
- Kim EK, et al. Human mesenchymal stem cell differentiation to the osteogenic or adipogenic lineage is regulated by AMP-activated protein kinase. J Cell Physiol. 2012;227(4):1680–7.View ArticlePubMedGoogle Scholar
- Piek E, et al. Osteo-transcriptomics of human mesenchymal stem cells: accelerated gene expression and osteoblast differentiation induced by vitamin D reveals c-MYC as an enhancer of BMP2-induced osteogenesis. Bone. 2010;46(3):613–27.View ArticlePubMedGoogle Scholar
- Eijssen LM, et al. User-friendly solutions for microarray quality control and pre-processing on ArrayAnalysis.org. Nucleic Acids Res. 2013;41(Web Server issue):W71–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Irizarry RA, et al. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003;4(2):249–64.View ArticlePubMedGoogle Scholar
- Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:Article3.PubMedGoogle Scholar
- da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.View ArticleGoogle Scholar
- da Huang W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.View ArticleGoogle Scholar
- Wishart DS, et al. DrugBank: a comprehensive resource for in silico drug discovery and exploration. Nucleic Acids Res. 2006;34(Database issue):D668–72.View ArticlePubMedGoogle Scholar
- Kang JS, et al. Repression of Runx2 function by TGF-beta through recruitment of class II histone deacetylases by Smad3. Embo J. 2005;24(14):2543–55.View ArticlePubMedPubMed CentralGoogle Scholar
- Maroni P, et al. Chemical and genetic blockade of HDACs enhances osteogenic differentiation of human adipose tissue-derived stem cells by oppositely affecting osteogenic and adipogenic transcription factors. Biochem Biophys Res Commun. 2012;428(2):271–7.View ArticlePubMedGoogle Scholar
- Zhou Y, Peng J, Jiang S. Role of histone acetyltransferases and histone deacetylases in adipocyte differentiation and adipogenesis. Eur J Cell Biol. 2014;93(4):170–7.View ArticlePubMedGoogle Scholar
- Bouloumie A, et al. Adipocyte produces matrix metalloproteinases 2 and 9: involvement in adipose differentiation. Diabetes. 2001;50(9):2080–6.View ArticlePubMedGoogle Scholar
- Bermeo S, Gunaratnam K, Duque G. Fat and bone interactions. Curr Osteoporos Rep. 2014;12(2):235–42.View ArticlePubMedGoogle Scholar
- Komm BS, et al. The safety and tolerability profile of therapies for the prevention and treatment of osteoporosis in postmenopausal women. Expert Rev Clin Pharmacol. 2015;8(6):769–84.View ArticlePubMedGoogle Scholar
- Devlin MJ, Rosen CJ. The bone-fat interface: basic and clinical implications of marrow adiposity. Lancet Diabetes Endocrinol. 2015;3(2):141–7.View ArticlePubMedGoogle Scholar
- Tuljapurkar SR, et al. Changes in human bone marrow fat content associated with changes in hematopoietic stem cell numbers and cytokine levels with aging. J Anat. 2011;219(5):574–81.View ArticlePubMedPubMed CentralGoogle Scholar
- Wrzesinski SH, Wan YY, Flavell RA. Transforming growth factor-beta and the immune response: implications for anticancer therapy. Clin Cancer Res. 2007;13(18 Pt 1):5262–70.View ArticlePubMedGoogle Scholar
- Neuzillet C, et al. Targeting the TGFbeta pathway for cancer therapy. Pharmacol Ther. 2015;147:22–31.View ArticlePubMedGoogle Scholar
- Santambrogio L, et al. Studies on the mechanisms by which transforming growth factor-beta (TGF-beta) protects against allergic encephalomyelitis. Antagonism between TGF-beta and tumor necrosis factor. J Immunol. 1993;151(2):1116–27.PubMedGoogle Scholar
- Mori T, et al. Role and interaction of connective tissue growth factor with transforming growth factor-beta in persistent fibrosis: a mouse fibrosis model. J Cell Physiol. 1999;181(1):153–9.View ArticlePubMedGoogle Scholar
- Robinson SP, Rose WC. Transforming growth factor beta 1: lack of in vivo antitumor activity on A549 and Wehi 3BD+ tumors. Anticancer Res. 1992;12(1):73–9.PubMedGoogle Scholar
- Blaney Davidson EN, van der Kraan PM, van den Berg WB. TGF-beta and osteoarthritis. Osteoarthritis Cartilage. 2007;15(6):597–604.View ArticlePubMedGoogle Scholar
- Lehrke M, Lazar MA. The many faces of PPARgamma. Cell. 2005;123(6):993–9.View ArticlePubMedGoogle Scholar
- Seals DF, Courtneidge SA. The ADAMs family of metalloproteases: multidomain proteins with multiple functions. Genes Dev. 2003;17(1):7–30.View ArticlePubMedGoogle Scholar
- Durigova M, et al. MMPs are less efficient than ADAMTS5 in cleaving aggrecan core protein. Matrix Biol. 2011;30(2):145–53.View ArticlePubMedGoogle Scholar
- Gorski DJ, et al. Deletion of ADAMTS5 does not affect aggrecan or versican degradation but promotes glucose uptake and proteoglycan synthesis in murine adipose derived stromal cells. Matrix Biol. 2015;47:66–84.View ArticlePubMedGoogle Scholar
- Glasson SS, et al. Deletion of active ADAMTS5 prevents cartilage degradation in a murine model of osteoarthritis. Nature. 2005;434(7033):644–8.View ArticlePubMedGoogle Scholar
- Rothenberg ML, Nelson AR, Hande KR. New drugs on the horizon: matrix metalloproteinase inhibitors. Stem Cells. 1999;17(4):237–40.View ArticlePubMedGoogle Scholar
- Sparano JA, et al. Randomized phase III trial of marimastat versus placebo in patients with metastatic breast cancer who have responding or stable disease after first-line chemotherapy: Eastern Cooperative Oncology Group trial E2196. J Clin Oncol. 2004;22(23):4683–90.View ArticlePubMedGoogle Scholar
- Bourlier V, et al. Protease inhibitor treatments reveal specific involvement of matrix metalloproteinase-9 in human adipocyte differentiation. J Pharmacol Exp Ther. 2005;312(3):1272–9.View ArticlePubMedGoogle Scholar
- Khokha R, Murthy A, Weiss A. Metalloproteinases and their natural inhibitors in inflammation and immunity. Nat Rev Immunol. 2013;13(9):649–65.View ArticlePubMedGoogle Scholar
- Mannello F. Natural bio-drugs as matrix metalloproteinase inhibitors: new perspectives on the horizon? Recent Pat Anticancer Drug Discov. 2006;1(1):91–103.View ArticlePubMedGoogle Scholar
- Penning TM. The aldo-keto reductases (AKRs): overview. Chem Biol Interact. 2015;234:236–46.View ArticlePubMedGoogle Scholar
- Zhang W, et al. Knockdown or inhibition of aldo-keto reductase 1B10 inhibits pancreatic carcinoma growth via modulating Kras-E-cadherin pathway. Cancer Lett. 2014;355(2):273–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Fukumoto S, et al. Overexpression of the aldo-keto reductase family protein AKR1B10 is highly correlated with smokers’ non-small cell lung carcinomas. Clin Cancer Res. 2005;11(5):1776–85.View ArticlePubMedGoogle Scholar
- Liu Z, et al. Epidermal growth factor induces tumour marker AKR1B10 expression through activator protein-1 signalling in hepatocellular carcinoma cells. Biochem J. 2012;442(2):273–82.View ArticlePubMedGoogle Scholar
- Ma J, et al. Aldo-keto reductase family 1 B10 affects fatty acid synthesis by regulating the stability of acetyl-CoA carboxylase-alpha in breast cancer cells. J Biol Chem. 2008;283(6):3418–23.View ArticlePubMedGoogle Scholar
- Chalk C, Benstead TJ, Moore F. Aldose reductase inhibitors for the treatment of diabetic polyneuropathy. Cochrane Database Syst Rev. 2007;4:Cd004572.PubMedGoogle Scholar
- Liu J, Wen G, Cao D. Aldo-keto reductase family 1 member B1 inhibitors: old drugs with new perspectives. Recent Pat Anticancer Drug Discov. 2009;4(3):246–53.View ArticlePubMedGoogle Scholar
- Matsunaga T, et al. Potent and selective inhibition of the tumor marker AKR1B10 by bisdemethoxycurcumin: probing the active site of the enzyme with molecular modeling and site-directed mutagenesis. Biochem Biophys Res Commun. 2009;389(1):128–32.View ArticlePubMedGoogle Scholar
- Endo S, et al. Selective inhibition of the tumor marker AKR1B10 by antiinflammatory N-phenylanthranilic acids and glycyrrhetic acid. Biol Pharm Bull. 2010;33(5):886–90.View ArticlePubMedGoogle Scholar
- Kleiveland CR, Kassem M, Lea T. Human mesenchymal stem cell proliferation is regulated by PGE2 through differential activation of cAMP-dependent protein kinase isoforms. Exp Cell Res. 2008;314(8):1831–8.View ArticlePubMedGoogle Scholar
- Tian XY, et al. Continuous PGE2 leads to net bone loss while intermittent PGE2 leads to net bone gain in lumbar vertebral bodies of adult female rats. Bone. 2008;42(5):914–20.View ArticlePubMedGoogle Scholar