- Research article
- Open Access
Regulation of chondrocyte gene expression by osteogenic protein-1
Arthritis Research & Therapy volume 13, Article number: R55 (2011)
The objective of this study was to investigate which genes are regulated by osteogenic protein-1 (OP-1) in human articular chondrocytes using Affimetrix gene array, in order to understand the role of OP-1 in cartilage homeostasis.
Chondrocytes enzymatically isolated from 12 normal ankle cartilage samples were cultured in high-density monolayers and either transfected with OP-1 antisense oligonucleotide in the presence of lipofectin or treated with recombinant OP-1 (100 ng/ml) for 48 hours followed by RNA isolation. Gene expression profiles were analyzed by HG-U133A gene chips from Affimetrix. A cut-off was chosen at 1.5-fold difference from controls. Selected gene array results were verified by real-time PCR and by in vitro measures of proteoglycan synthesis and signal transduction.
OP-1 controls cartilage homeostasis on multiple levels including regulation of genes responsible for chondrocyte cytoskeleton (cyclin D, Talin1, and Cyclin M1), matrix production, and other anabolic pathways (transforming growth factor-beta (TGF-β)/ bone morphogenetic protein (BMP), insulin-like growth factor (IGF), vascular endothelial growth factor (VEGF), genes responsible for bone formation, and so on) as well as regulation of cytokines, neuromediators, and various catabolic pathways responsible for matrix degradation and cell death. In many of these cases, OP-1 modulated the expression of not only the ligands, but also their receptors, mediators of downstream signaling, kinases responsible for an activation of the pathways, binding proteins responsible for the inhibition of the pathways, and transcription factors that induce transcriptional responses.
Gene array data strongly suggest a critical role of OP-1 in human cartilage homeostasis. OP-1 regulates numerous metabolic pathways that are not only limited to its well-documented anabolic function, but also to its anti-catabolic activity. An understanding of OP-1 function in cartilage will provide strong justification for the application of OP-1 protein as a therapeutic treatment for cartilage regeneration and repair.
Cartilage degeneration is one of the features of osteoarthritis (OA). In order to identify cellular mechanisms that drive OA progression, it is necessary to understand the interplay between anabolic and catabolic processes responsible for cartilage homeostasis under physiological and pathophysiological states. Osteogenic protein-1 (OP-1) or bone morphogenetic protein-7 (BMP-7) is one of the most potent growth factors for cartilage maintenance and repair identified thus far [1, 2]. A large number of in vivo and in vitro studies have shown a high synthetic potency of human recombinant OP-1 (rhOP-1; ). In earlier work, we found that the inhibition of OP-1 gene expression by antisense oligonucleotides (ODNs) caused a significant decrease in aggrecan expression, aggrecan core protein synthesis, and proteoglycan (PG) synthesis, which resulted in the depletion of PGs from the cartilage matrix . These findings suggest that OP-1 plays a key role in maintenance of cartilage integrity and homeostasis, but further work is needed to understand the mechanisms by which OP-1 acts at the molecular level.
In the current study, we used the Affymetrix GeneChip technology to monitor OP-1 regulation of 22,000 genes from the human genome with specific emphasis on genes that are relevant to adult articular cartilage. Those included matrix proteins, anabolic and catabolic gene products, as well as their intracellular regulators and receptors. Recently, applying the same methodology differential gene expression pattern in normal and OA cartilage tissue was identified . These analyses revealed numerous interesting gene expression profiles, but per se did not allow elucidating cellular reaction patterns in response to defined extracellular stimuli. The goal of the current project was to evaluate the role OP-1 plays in regulating human articular cartilage homeostasis by using a gene array approach under conditions where endogenous OP-1 gene expression was inhibited by antisense ODNs (; OP-1AS) or OP-1 signaling was activated and/or enhanced by rhOP-1. Key microarray findings were verified by real-time PCR and additional in vitro experiments of matrix synthesis and signal transduction. We found that OP-1/BMP-7 controls numerous metabolic pathways that are not limited to its direct anabolic or anti-catabolic function, but also related to cell growth, cell proliferation, differentiation, survival, apoptosis, and death.
Materials and methods
Dulbecco's modified Eagle's medium (DMEM), fetal bovine serum (FBS), gentamicin, Ham's F-12, lipofectin, Opti-MEM, penicillin/streptomycin/fungizone (PSF), 1X Platinum Quantitative PCR SuperMix-UDG and SuperScript III reverse transcriptase with oligo (dT)12-18 were purchased from Invitrogen (Carlsbad, CA, USA). Phosphorothioate ODN was custom synthesized by Oligos Etc. (Wilsonville, OR, USA). RNeasy mini kit, QIA shredder, RNase-free DNase kit and QuantiTect Primer Assay were purchased from Qiagen (Valencia, CA, USA). Real time polymerase chain reaction (PCR) primers were custom synthesized by Integrated DNA Technologies (IDT), Coralville, IA, USA. 10,000 X SYBR Green 1 was purchased from Cambrex, Rockland, ME, USA. Recombinant human rhOP-1 was kindly provided by Stryker Biotech (Hopkinton, MA, USA).
Isolation and culture of chondrocytes
Full-thickness articular cartilage from the talus of the talocrural joint (ankle) from 12 human organ donors (age 55 to 70 years old, Collins grade 0 to 1 ) and from the femur of the tibiofemoral joint (knee) from two human organ donors (age 67 and 73 years old, Collins grade 2) was obtained from the Gift of Hope Organ and Tissue Donor Network (Elmhurst, IL, USA) with Institutional Review Board approval and appropriate consent within 24 hours of the donor's death. Knee cartilage was utilized for verification of the ankle cartilage results using real-time PCR. Chondrocytes were isolated by sequential digestion with pronase (2 mg/ml) for 60 minutes and collagenase P (0.25 mg/ml) overnight . Chondrocytes were plated in high density monolayer culture (4 × 106 cells/well in a six-well plate) and cultured for 24 hours in 50% DMEM/50% Ham's F-12 supplemented with 10% FBS, 1% PSF, and gentamicin (50 μg/ml) for attachment prior to treatment with either antisense (OP-1 AS) or recombinant OP-1 (rhOP-1). Both treatments were administered for 48 hours in the absence of serum.
Antisense ODNs were designed to be complementary to sequences in the 5'- and 3'-untranslated regions of the human OP-1 messenger RNA (mRNA) sequence (XM_030621, National Center for Biotechnology Information (NCBI)) as described . All verification experiments with appropriate negative controls (sense and scrambled probes) were performed in a previous study . For this study, the following antisense ODN was used: 5'-GGC-GAA-CGA-AAA-GGC-GAG-TGA-3' (position 237-257).
Chondrocyte cultures were divided into three experimental groups and treated for 48 hours as follows: 1) transfected with OP-1 AS in the presence of 10 μg/ml lipofectin ; 2) treated with 100 ng/ml of rhOP-1; and 3) culture control (no treatment, no serum).
Total cellular RNA was isolated using the RNeasy Mini Kit, following lysis of the cells with a Qia shredder  and included an on-column DNase digestion, according to the manufacturer's instructions (Qiagen). All samples were stored at -80°C until analyzed.
Microarray and pathway analysis
Gene expression profiles were analyzed by HG-U133A gene chips from Affimetrix (accession number: E-MTAB-571). At least 10 μg of RNA/per experimental group was required for analysis. Therefore, the RNA was pooled from donors in order to have sufficient RNA and to reduce donor-to-donor variations. Cells from all 12 donors were treated with each experimental condition. The microarray data collection was in compliance with the Minimum Information About Microarray Experiments standard . The quality of the RNA was checked by the Agilent Bioanalyzer (Agilent Technologies, Inc., Santa Clara, CA, USA), and the quality of the hybridization image was checked by the affyPLM model . To deal with the technical variation, each gene was measured by 11 different probes on the Affymetrix U133A microarray. A statistical model at the probe-level was used to identify the differentially expressed genes. To estimate the variance more efficiently with a small sample size, we utilized an empirical Bayesian correction of the linear model . Statistical significance was considered with a P-value of P < 0.001 and fold change larger than 1.5-fold between the treatment group and corresponding control. All the data analysis was conducted using the Bioconductor/R package . To interpret the biological significance of differentially expressed genes, a gene ontology analysis was conducted using DAVID/EASE .
Pathway analysis and classification by gene ontology
Regulated genes (R > 1.5-fold, P < 0.001) were used as input for both analyses. The ingenuity pathway analysis system  was used to project genes onto known biological pathways (canonical pathways). The system determines a significance value for each pathway based on an F-statistics that the input-genes occur randomly within this pathway. Grouping of genes was done by computing over-representation of regulated genes in gene ontology (GO) classes . Statistical analysis consisted of 1) analysis of differentially expressed genes under a single experimental condition in comparison to the corresponding control (up- or down-regulated in the presence of OP-1 antisense or rhOP-1); 2) analysis of differentially expressed genes when comparison is made between two treatments (OP-1 antisense and rhOP-1); and 3) gene ontology, when changes were analyzed within a family of genes according to their function (comparison was made between single treatment and control or between both treatments). Selected gene array results were verified experimentally in vitro or by real-time PCR.
Validation experiments -quantitative real time PCR
Selected gene array results were verified by real-time PCR. SuperScript III reverse transcriptase with oligo (dT)12-18 was used to transcribe 4 μg of isolated total RNA into complementary DNA (cDNA) in a total volume of 20 μl according to the manufacturer's instructions (Invitrogen). Real time PCR primer sets specific for human β-actin, GAPDH, gremlin-1, IL-6, IL-8, and LIF-1 (Table 1) were designed using the PrimerQuest program (Integrated DNA Technologies, Inc., Coralville, Iowa, USA). The specificity of the primers was verified by testing in BLAST searches . Real time PCR primer sets specific for human 18SrRNA and BMP-2 were purchased from Qiagen. Real time PCR was performed using the Smart Cycler System (Cepheid, Sunnyvale, CA, USA). Each 50 μl reaction mixture contained 1X Platinum Quantitative PCR SuperMix-UDG, 0.5X Smart Cycler additive reagent (0.1 mM Tris, pH 8.0; 0.1 mg of bovine serum albumin per ml, 75 mM trehalose, and 0.1% Tween 20), 0.5X SYBR Green 1 (vendor stock 10,000X; Cambrex, Rockland, ME), 0.2 μM each of forward and reverse primer (IDT primers) or 1X QuantiTect primers (Qiagen primers) and 1 μl cDNA (18SrRNA, β-actin, BMP-2, GAPDH, gremlin-1, IL-6, IL-8) or 2 μl cDNA (LIF-1). Cycling parameters were: preheat at 60°C for 120 seconds then 95°C for 120 seconds followed by 40 three-step cycles of 95°C for 15 seconds, various annealing temperatures and times (Table 1) and 72°C for 30 seconds. After the last amplification cycle, PCR products were analyzed by melting curve analysis in the Smart Cycler by slowly increasing the temperature to 95°C. The reactions were run in triplicate with appropriate controls (no cDNA template). The data were analyzed by using the Cepheid Smart Cycler software (version 2.0c) and reported as threshold cycle (Ct). Change in gene expression was calculated as fold change = 2-Δ(ΔCt), where Δ(ΔCt) = (Ct sample - Ct housekeeping gene) - (Ct control - Ct housekeeping gene).
Statistical analysis for real-time
PCR Data are expressed as mean +/- standard deviation. Statistical significance was assessed by the Student t-test and P-values < 0.05 were considered significant.
Microarray analysis: overview of data
GeneChip (HG-U133A) expression data from un-stimulated, rhOP-1 and OP-1AS treated chondrocytes maintained in high-density monolayer culture were generated. For the analysis of the expression data we used a three step analytical strategy: (I) processing of raw intensity values and normalization of profiles, (II) examination of expression levels of gene categories that are relevant to articular cartilage, and (III) comparison of gene expression changes between the two treatments - OP-1AS to knockdown endogenous OP-1 expression vs. addition of exogenous rhOP-1.
Analyzing the number of differentially expressed genes (fold changes of larger than 1.5 and corresponding P-values < 0.001 compared to control) after rhOP-1 or OP-1AS, we found that rhOP-1 modulated expression of 4,057 genes, while OP-1AS treatment modulated expression of only 2,618 genes respectively. More genes were down-regulated than up-regulated by either treatment: rhOP-1 down-regulated 3,365 genes vs 692 genes that were up-regulated; while OP-1AS down-regulated 2,364 genes and up-regulated only 254 genes. The functional groups of genes modulated by lack or excess of OP-1 are depicted in Figure 1. RhOP-1 primarily controlled genes responsible for molecular function, biological processes, and cellular components, while OP-1AS primarily affected genes controlling cellular processes and catalytic activity. Interestingly, either treatment up-regulated fewer functional groups than the number that were down-regulated (Figure 1). For example, rhOP-1 induced only five functional groups vs four induced by OP-1AS; while rhOP-1 down-regulated 19 functional groups vs 12 down-regulated by OP-1AS. When the results were compared between the two treatments, we found that very few gene groups with the same function were differentially regulated by both treatments (Figure 1). Groups regulated by both OP-1 conditions included the genes responsible for cellular processes (the same number of genes were up-regulated by either treatment, 100 vs 101), development, protein binding, signal transducer activity and signal transduction.
Analysis of catabolic genes: cytokines and their regulators
Previously, we showed that OP-1 was able to counteract the catabolic activity of IL-1β [16, 17] and other catabolic mediators such as fragments of cartilage matrix, fibronectin and hyaluronan [17–20]. Therefore, it was of interest to determine the effects of OP-1 on genes regulating pro-catabolic activity. Consistent with an anti-catabolic function for OP-1, a broad spectrum of genes with various pro-catabolic activities (cytokines and their regulators, matrix degrading proteinases, apoptosis-related genes, neuromediators, transcription factors, and so on) were modulated by OP-1. Multiple cytokines and chemokines, in particular members of the IL-6 family, (Figure 2), as well as their receptors and regulators of their activity (Tables 2 and 3) were found to be regulated by OP-1. Interestingly, among these mediators only members of the IL-6 family (leukemia inhibitory factor (LIF), IL-11, IL-8, and IL-6) were differentially regulated by the two treatment conditions: rhOP-1 down-regulated LIF expression by more than 15-fold, IL-11 expression by more than eight-fold, IL-8 gene by four-fold and IL-6 by two-fold, respectively (Figure 2A). Likewise, when endogenous OP-1 was inhibited by OP-1AS, expression of these four chemokines was elevated by about two-fold indicating a tight association between OP-1 levels and expression of members of the IL-6 family. Verification experiments of gene array findings included both real-time PCR analysis and in vitro metabolic tests (Figure 2). These tests confirmed that when chondrocytes in high-density monolayer cultures were treated with rhOP-1 for 48 hours, gene expression of LIF, IL-6, and IL-8 was inhibited as detected by real-time PCR, although the magnitude of changes was different from those identified by gene array (Figure 2A, B). In metabolic studies, we also found that OP-1 could overcome an inhibitory effect of IL-6 on PG synthesis in chondrocytes cultured in alginate beads (Figure 2C). In addition, our previous studies showed an ability of OP-1 to inhibit mRNA expression of IL-1, IL-6, IL-8, and other cytokines in primary and immortalized chondrocytes .
In analyzing the relationship between treatments to modulate OP-1 and the expression of genes in the IL-6 signaling pathway, we found that OP-1 not only regulates expression of the IL-6 family of cytokines but also controls expression of their receptors and downstream intracellular mediators including signal transducers and activators of transcription (STATs), mitogen activated protein (MAP) kinases, and transcription factors. This suggests OP-1 inhibits IL-6 signaling at multiple levels (Table 3). Among other genes that either regulate cytokine activity or mediate their signaling, the most affected by OP-1 were the receptors for IL-1β and tumor necrosis factor alpha (TNF-α) (see Table 2) as well as TNF-α inducible protein. Although under the experimental conditions expression of TNF-α and IL-1β genes was not influenced by OP-1, previous studies showed that injection of OP-1 into nucleus pulposus inhibited production of autocrine TNF-α and IL-1β elevated in response to injurious compression of the intervertebral discs  proving an association between OP-1 and signaling pathways of the above mentioned cytokines. In addition, several other studies have provided evidence of an ability of OP-1 to regulate either IL-1β induced responses or IL-1β downstream signaling [16–18, 22, 23].
Analysis of catabolic genes. Neuromediators
Previous studies have provided evidence that OP-1 may regulate mediators of pain-related behavior and their activation in response to injurious compression of the intervertebral disc and acute cartilage trauma [24–26]. We also reported that injection of OP-1 into nucleus pulposus down-regulated substance P expression , bradykinin and bradykinin inducible receptor β1 . Therefore, it was of interest to examine expression of neuromediators and their receptors in the present array study. After stimulation for 48 hours with rhOP-1, expression of the receptors of bradykinin and substance P was down-regulated (Table 4). Both receptors of bradykinin, constitutively expressed β2 and inducible β1, were down-regulated by the treatment with OP-1. Expression of the β1 receptor was differentially regulated under conditions of excess and lack of OP-1, that is, treatment with rhOP-1 inhibited gene expression of this receptor by 1.85-fold, while its expression was up-regulated by 1.59-fold when endogenous OP-1 expression was inhibited by antisense oligonucleotides. These results are consistent with previous data on the protein level in an in vivo model of disc herniation, where injection of OP-1 into the nucleus pulposus completely abolished bradykinin receptor β1 . Although by gene array we did not identify significant changes in the expression of bradykinin and substance P at the time point tested here, we found changes in substance P receptor and its precursor. We also found that OP-1 inhibited expression of nerve growth factor-β by almost two-fold.
Analysis of catabolic genes: Transcription factors
Besides cytokines and their receptors, OP-1 also affected gene expression of transcription factors that regulate cytokine signaling. Previously, in normal primary and immortalized chondrocytes, we found that OP-1 inhibits activation of the nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) and activator protein-1 (AP-1) transcription factors . Here, expression of a large set of transcription factors was found to be modulated by OP-1 (Table 5). In addition to common factors such as NF-κB, STAT1 and STAT6, gene array also discovered factors that repress IL-2 expression, p38 interacting protein, Runx1, and others. The majority of these transcription factors regulate directly or indirectly (as p38 interacting protein) transcriptional responses induced by various pro-inflammatory mediators (IL-1β, IL-6, matrix fragments). Others, like Runx1, are involved in the process of chondrogenesis. To further demonstrate the effect of OP-1 on activation of transcription factors, we treated cultured cells and found that OP-1 was able to at least partially inhibit activation of NF-κB in primary chondrocytes pre-treated with IL-1β or activation of Stat-1 in chondrocytes treated with IL-6 and IL-6 soluble receptor (data not shown).
Analysis of catabolic genes: Matrix degrading proteases, cathepsins, and apoptosis-related genes
Among other catabolic genes influenced by OP-1 were the matrix metalloproteinases (MMPs), cathepsins, and a number of proteases with various modes of action (Table 6). Thus, expression of membrane type-1 MMP (MMP-14) was inhibited by rhOP-1 by 1.6-fold (P < 0.001) along with tissue inhibitor of metalloproteinases (TIMP)-3 (2.06-fold, P < 0.001). At the same time, expression of MMP-2 (gelatinase A), which is activated by MMP-14 , was not affected by rhOP-1, but was down-regulated by OP-1AS (2.31-fold, (P < 0.001) as well as was MMP-9 (gelatinase B) (1.5-fold). Interestingly, the same positive correlation was found between the levels of OP-1 and expression of another TIMP, TIMP-4, which was decreased by 1.7-fold in the OP-1AS group confirming its association with MMP-2 . Parallel changes were observed in other types of proteases, such as a disintegrin and metalloproteinases (ADAM)-9, 10, and 28. Their gene expression was down-regulated under OP-1AS from 2.34 to 1.75-fold. Treatment of chondrocytes with rhOP-1 inhibited expression of ADAM-15,-19, as well as urokinase type plasminogen activator, its receptor, and transglutamianse-2. There were also some proteinases that were up-regulated by rhOP-1: ADAM-TS7, ADAM-TS12, and tissue specific plasminogen activator suggesting that perhaps these proteins are involved in anabolic/remodeling processes.
Among the proteases that were also regulated by OP-1 were cathepsins B, C, and S. So far, these lysosomal cysteine proteases have been less studied in cartilage, though cathepsin C appears to be a central coordinator for activation of many serine proteases in immune/inflammatory cells , while cathepsin B was thought to play an important role in the development of osteoarthritis . Expression of all three cathepsin genes was down-regulated under OP-1AS.
A previous study on acute impact injury in vivo  strongly suggested an anti-apoptotic effect of OP-1 in post-traumatic OA. Therefore, we expected that OP-1 may control genes involved in apoptosis-related processes. We found that rhOP-1 inhibited program cell death 8 gene (apoptosis-induced factor), Bcl-2 gene and the calpain-9 gene (Table 6). However, the key caspases that trigger and promote cell death by apoptosis were not affected. During the absence of OP-1 (antisense treatment), expression of caspases 8, 9, and 6 were inhibited and only caspase 2 was elevated (Table 6). The reason for a down-regulation of the apoptosis-related genes under conditions where OP-1 is lacking is not clear, but may be a response to help avoid cell death.
Analysis of anabolic genes: transforming growth factor-beta (TGF-β)/BMP family, their receptors and regulators of signaling
Affimetrix analysis identified a very interesting effect of OP-1 on members of the BMP/TGF-β family (Table 7). Treatment with rhOP-1 down-regulated expression of growth differentiation factor (GDF)-15, BMP-2, and Activin A, and BMP-2 inducible kinase, while inhibition of OP-1 expression up-regulated GDF-15 and Activin A. Down-regulation of BMP-2 expression in chondrocytes treated with rhOP-1 was confirmed by real-time PCR (Figure 3). Antisense reduction of OP-1 levels resulted in down-regulation of GDF-10 and TGF-α expression (Table 7). Furthermore, a correlation was also found between OP-1 and the mediators of its downstream signaling, where OP-1AS treatment inhibited expression of transcription factors, Id proteins 2 to 4 (Table 7), binding protein Gremlin (Figure 2), and MAD genes. Changes in Id genes correlated with the earlier findings from our laboratory, which demonstrated that the treatment of chondrocytes with rhOP-1 led to the elevation of Id1, Id2, and Id3 genes and proteins . Contrary to changes in the Gremlin gene, which showed a positive correlation with OP-1 levels, expression of Follistatin binding protein was inhibited by more than two-fold in chondrocytes treated with rhOP-1.
In addition, OP-1 modulated expression of the TGF-β/BMP receptors. With the exception of Activin-α RIB, which was inhibited by rhOP-1 and elevated under the lack of OP-1, expression of other receptors, Activin-α RIIB, BMPR1A, TGF-β RI, II, and III correlated positively with OP-1 expression (Table 7).
Analysis of anabolic genes: other growth factors
Previously we showed that rhOP-1 stimulated expression of insulin-like growth factor (IGF)-1 and IGF-1 receptor genes , while inhibition of OP-1 gene expression by OP-1AS down-regulated mRNA expression of these genes. We have also documented a synergistic effect of OP-1 on IGF-1 induced responses in normal and OA chondrocytes [33, 34]. Here, we confirmed an association between OP-1 and IGF-1 pathways by documenting a 1.73-fold decrease in IGF-1 receptor expression and a decrease in two IGF-1 binding proteins-5 and 7 (1.9- and 1.5-fold respectively) under OP-1AS. Furthermore, other genes within the IGF-1 signaling pathway were regulated by OP-1. Among them were PIK3R1, PRKAR2B, MAP2K2, PDE3B, and SOCS3 (Table 8).
Modulation of OP-1 levels affected mRNA expression of growth factors and some of their receptors that belong to various families, such as Nerve Growth Factor-β, Vascular Endothelial Growth Factor, Endothelial Cell Growth Factor 1 (platelet-derived), Capillary Morphogenesis Protein-1, and Fibroblast Growth Factor (FGF)-7. Their expression was inhibited by rhOP-1 from 1.93- to 1.5-fold. Contrary, the expression of the FGF-R2 and 3 receptors, and α and β receptors of Platelet-Derived Growth Factor was stimulated by rhOP-1 Table 8).
Matrix proteins and their receptors
Cartilage-specific matrix genes were found to be modulated by rhOP-1 treatment. Expression of the collagen type IX-α3 chain and cartilage oligomeric protein (COMP) was up-regulated by about 1.5-fold in chondrocytes treated with rhOP-1 (Table 9). Among proteoglycans, versican was affected the most (by about three-fold down-regulation by OP-1AS) and syndecan was differentially regulated under both rhOP-1 and OP-1AS treatments. There were a number of other matrix genes regulated by OP-1: bone sialoprotein, osteonectin, cadherins, chondroitin sulfate PG4 and dermatan sulfate PG3 (Table 9). As expected, there was a positive correlation between OP-1 and CD44 expression. Inhibition of OP-1 expression resulted in 2.34-fold reduction in CD44 expression. However, contrary to previously published data , rhOP-1 inhibited hyaluronan synthase 2 expression. A number of basement membrane proteins were modulated by OP-1: α1,2,3, and five chains of collagen type IV, laminin, versican among others. Gene expression of bamacan and laminin was inhibited by rhOP-1 and stimulated under OP-1AS. OP-1 also modulated expression of collagens that are not cartilage-specific, such as collagen type I, IV, V, VI, VIII, XIV, and XVI. Their expression was inhibited under the OP-1AS treatment (Table 9). The greatest decrease in mRNA expression was found for α1 and α2 chains of type I collagen (more than 2.6-fold).
To the best of our knowledge, this is the first report that uses microarray analysis to provide a comprehensive understanding of the role OP-1 plays in overall cartilage homeostasis. We found that OP-1 controls cartilage homeostasis on multiple levels including regulation of genes responsible for chondrocyte cytoskeleton (cyclin D, Talin1, and Cyclin M1, for example, and confirmed in ), matrix production and other anabolic pathways, as well as regulation of cytokines and various catabolic pathways responsible for matrix degradation and cell death. Importantly, in many of these cases, OP-1 modulated the expression of not only the ligands, but also their receptors, mediators of downstream signaling, kinases responsible for an activation of the pathways and transcription factors that induce transcriptional responses.
Due to high variability among human samples, only a few studies have utilized microarray analysis to test the entire human genome in primary adult articular chondrocytes [4, 37–39], and only one of Saas et al.  addressed in part the effect of BMP-7/OP-1. These analyses used the tissue from one or a maximum of two donor cartilage samples. In the present study, normal (grade 0) articular cartilage was collected from 12 donors within a similar age range. One of the limitations of the study is that we examined gene expression profiles only at one time point, after 48-hours of culture. Therefore, changes in early-response genes and late-response genes might have been missed. This could explain some results, as for example, the lack of changes in the expression of major cartilage matrix proteins. However, such an approach gave us a breath of the overall effects of OP-1 on cartilage homeostasis.
Due to the abundance of the results, we will discuss only the most relevant and those that could be explained by the current knowledge of the field. Perhaps most important was the finding that OP-1 is a unique growth factor in its capacity to display simultaneously pro-anabolic and anti-catabolic activities. It was previously shown that OP-1 stimulated expression and synthesis of collagen type II, aggrecan, hyaluronan, and CD44 [1, 2, 20, 40] as well as IGF-1, IGF-1 receptor, and responses to IGF-1 . In the current studies, we used high-density monolayers while in previous work explants or alginate beads were used with different media conditions (no serum vs serum or ITS-media). The finding that the microarray results shown here were comparable to the previous results suggest that the pro-anabolic effects of OP-1 in human articular chondrocytes are persistent. With regard to the anti-catabolic activity, the ability of OP-1 to counteract various pro-inflammatory/catabolic responses or directly inhibit expression of catabolic mediators was previously shown in primary chondrocyte cultures or in animal models of post-traumatic osteoarthritis or disc degeneration [17–19, 24, 31]. In this study, we found that OP-1 inhibits expression of IL-6 and members of the IL-6 family of chemokines as well as their receptors and signaling mediators. Furthermore, the tight association between these two classes of mediators (OP-1 and IL-6) was documented under both experimental conditions (plus or minus OP-1). Based on our new data on the role of IL-6 in acute post-traumatic responses , it is possible that OP-1 was able to protect cartilage from degenerative changes caused by acute trauma  not only due to its direct effect on matrix synthesis, but also because of its ability to inhibit IL-6, TNF-α, and the catabolic pathways induced by the fragments of the extracellular matrix: fibronectin , hyaluronan , and collagen telopeptides .
Another important effect of OP-1 may be an ability to inhibit expression of neuromediators and their receptors. Previously, an anti-pain effect of OP-1 was documented in the rat models of herniated disc or disc degeneration induced by injurious compression. In these studies, OP-1 injections reduced hyperalgesia and inhibited elevation of IL-1, TNF-α, substance P, bradykinin and their receptors in various disc tissues including spinal cord and dorsal ganglion [21, 24, 26]. In chondrocytes, it is the first report that indicates a connection between OP-1 and various neuromediators, though substance P and its NK-1 receptor were already identified in cartilage in the model of mechanical stress . Very recently, our findings received another proof in phase I placebo-controlled clinical studies on OP-1 treatment for osteoarthritis patients , in which a single injection of OP-1 reduced pain even after six months of treatment.
Interesting results were found with regard to the ability of OP-1 to regulate the TGF-β/BMP signaling pathway. The levels of OP-1 expression positively correlated with the expression of activin-like kinase (ALK)-3 or BMP-RIA, transcription factors Id proteins 2 to 4, and a binding protein Gremlin, indicating that this could be a primary route for OP-1 signaling. We also found that another binding protein, Follistatin, exhibited a negative correlation with the levels of OP-1. Thus, our results suggest a differential regulation of these two binding proteins by OP-1, which could mean that Gremlin and Follistatin perform distinct functions in the mediating BMP responses or they are involved in different stages of signaling. This point of view is supported by studies of Tardif et al.  who reported their differential expression and spatial distribution in the dog OA model. One of the most surprising results was the finding that OP-1 inhibits expression of another member of the BMP family, BMP-2, which shares the same signaling machinery and in many cases exhibits similar anabolic activities [23, 46]. This result was confirmed by real-time PCR and definitely warrants further studies in understanding the responses induced by homologous, yet very different members of the same family .
Finally, another unexpected result was the inhibition of TIMP-3 expression by rhOP-1; while previously, TGF-β has been shown to induce this inhibitor . The differences in the results could be attributed to distinct culture conditions (primary chondrocytes in our case and passaged chondrocytes as in Qureshi et al. , assessment at various time points, or distinct signaling mechanisms between TGF-β and OP-1 that are responsible for induction of TIMP-3 expression. However, on the protein level it has already been reported that OP-1 inhibits TIMP-3 protein . Furthermore, changes in TIMP-3 were parallel to changes in certain genes responsible for apoptosis, which supports the notion that in cancer cells TIMP-3 may promote cell death by apoptosis . On the other hand, TIMP-3 has been shown to inhibit aggrecanase-mediated release of glycosaminoglycans in bovine nasal cartilage . At this point, the role of TIMP-3 in human chondrocytes and its regulation by various mediators remains to be investigated.
This analysis of gene array data strongly suggests a critical role of OP-1 in human cartilage homeostasis. OP-1 regulates numerous metabolic pathways that are not only limited to its anabolic function, but also to its anti-catabolic activity. Understanding of OP-1 function in cartilage will provide strong justification for the application of OP-1 protein as therapeutic treatment for cartilage regeneration and repair.
a disintegrin and metalloproteinase
B-cell lymphoma 2
bone morphogenetic proteins
Dulbecco's modified Eagle's medium
fetal bovine serum
fibroblast growth factor
growth differentiation factor
insulin-like growth factor
leukemia inhibitory factor
- MAP kinase:
mitogen activated protein kinase
National Center for Biotechnology Information
nuclear factor kappa-light-chain-enhancer of activated B cells
polymerase chain reaction
runt-related transcription factor 1
signal transducers and activators of transcription
transforming growth factor-beta
tissue inhibitor of metalloproteinases
tumor necrosis factor-alpha
vascular endothelial growth factor.
Chubinskaya S, Kuettner KE: Regulation of osteogenic proteins by chondrocytes. Int J Biochem Cell Biol. 2003, 35: 1323-1340. 10.1016/S1357-2725(03)00035-9.
Chubinskaya S, Hurtig M, Rueger DC: (OP-1/BMP-7) in cartilage repair. In vitro and in vivo study. Bone Morphogenetic Proteins: From Local to Systemic Therapeutics. Edited by: Birkhauser VS. 2007, Berlin: Springer Verlag, 277-317.
Söder S, Hakimiyan A, Rueger DC, Kuettner KE, Aigner T, Chubinskaya S: Antisense inhibition of osteogenic protein-1 disturbs human articular cartilage integrity. Arthritis Rheum. 2005, 52: 468-478. 10.1002/art.20856.
Saas J, Haag J, Rueger DC, Chubinskaya S, Sohler F, Bartnik E, Aigner T: IL-1β, but not BMP-7 implicates a dramatic change in the gene expression pattern of human adult articular chondrocytes. Cytokines. 2006, 36: 90-99. 10.1016/j.cyto.2006.10.005.
Collins DH: Osteoarthritis. Chapter V. The Pathology of Articular and Spinal Diseases. 1949, London: Edward Arnold Publishers, 74-115.
Aydelotte MB, Kuettner KE: Differences between sub-populations of cultured bovine articular chondrocytes. I. Morphology and cartilage matrix production. Connect Tissue Res. 1988, 18: 205-222. 10.3109/03008208809016808.
Bau B, Haag J, Schmid E, Kaiser M, Gebhard PM, Aigner T: Bone morphogenetic protein-mediating receptor-associated Smads as well as common Smad are expressed in human articular chondrocytes but not up-regulated or down-regulated in osteoarthritic cartilage. J Bone Miner Res. 2002, 17: 2141-2150. 10.1359/jbmr.2002.17.12.2141.
Brazma A, Hingamp P, Quackenbush J, Sherlock G, Spellman P, Stoeckert C, Aach J, Ansorge W, Ball CA, Causton HC, Gaasterland T, Glenisson P, Holstege FC, Kim IF, Markowitz V, Matese JC, Parkinson H, Robinson A, Sarkans U, Schulze-Kremer S, Stewart J, Taylor R, Vilo J, Vingron M: Minimum information about a microarray experiment (MIAME)-toward standards for microarray data. Nat Genet. 2001, 29: 365-371. 10.1038/ng1201-365.
Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP: Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003, 31: e15-10.1093/nar/gng015.
Wettenhall JM, Smyth GK: limmaGUI: a graphical user interface for linear modeling of microarray data. Bioinformatics. 2004, 20: 3705-3706. 10.1093/bioinformatics/bth449.
Dudoit S, Gentleman RC, Quackenbush J: Open source software for the analysis of microarray data. Biotechniques. 2003, 45-51. Suppl
Hosack DA, Dennis G, Sherman BT, Lane HC, Lempicki RA: Identifying biological themes within lists of genes with EASE. Genome Biol. 2003, 4: R70-10.1186/gb-2003-4-10-r70.
Ingenuity pathway system. [http://www.ingenuity.com]
Gene Ontology analysis. [http://www.geneontology.org]
BLAST search. [http://www.ncbi.nlm.nih.gov/blast]
Chubinskaya S, Segalite D, Pikovsky D, Hakimiyan AA, Rueger DC: Effects induced by BMPS in cultures of human articular chondrocytes: comparative studies. Growth Factors. 2008, 26: 275-283. 10.1080/08977190802291733.
Im HJ, Pacione C, Chubinskaya S, Van Wijnen AJ, Sun Y, Loeser RF: Inhibitory effects of insulin-like growth factor-1 and osteogenic protein-1 on fibronectin fragment- and interleukin-1beta-stimulated matrix metalloproteinase-13 expression in human chondrocytes. J Biol Chem. 2003, 278: 25386-25394. 10.1074/jbc.M302048200.
Huch K, Wilbrink B, Flechtenmacher J, Koepp HE, Aydelotte MB, Sampath TK, Kuettner KE, Mollenhauer JA, Thonar EJ-MA: Effects of recombinant human osteogenic protein 1 on the production of proteoglycan, prostaglandin E2, and interleukin-1 receptor antagonist by human articular chondrocytes cultured in the presence of interleukin-1β. Arthritis Rheum. 1997, 40: 2157-2161. 10.1002/art.1780401209.
Koepp HE, Sampath KT, Kuettner KE, Homandberg GA: Osteogenic protein-1 (OP-1) blocks cartilage damage caused by fibronectin fragments and promotes repair by enhancing proteoglycan synthesis. Inflamm Res. 1999, 48: 199-204. 10.1007/s000110050446.
Nishida Y, Knudson CB, Knudson W: Osteogenic protein-1 inhibits matrix depletion in a hyaluronan hexasaccharide-induced model of osteoarthritis. Osteoarthritis Cart. 2004, 12: 374-382. 10.1016/j.joca.2004.01.008.
Chubinskaya S, Kawakami M, Rapoport L, Matsumoto T, Migita N, Rueger DC: Anti-catabolic effect of OP-1 in chronically compressed intervertebral discs. J Ortho Res. 2007, 25: 517-530. 10.1002/jor.20339.
Elshaier A, Hakimiyan AA, Rappoport L, Rueger DC, Chubinskaya S: Effect of interleukin-1β on osteogenic protein-1-induced signaling in adult human articular chondrocytes. Arthritis Rheum. 2009, 60: 143-154. 10.1002/art.24151.
Merrihew C, Soeder S, Rueger DC, Kuettner KE, Chubinskaya S: Differential response of endogenous osteogenic protein-1 to interleukin-1β. J Bone Joint Surg Am. 2003, 85-A: 67-74.
Kawakami M, Matsumoto T, Hashizume H, Kuribayashi K, Chubinskaya S, Yoshida M: Osteogenic protein-1 (OP-1/BMP-7) inhibits degeneration and pain-related behavior induced by chronically compressed nucleus pulposus in the rat. Spine. 2005, 30: 1933-1939. 10.1097/01.brs.0000176319.78887.64.
Willson T, Rappoport L, Wimmer M, Hakimiyan A, Shekhawat V, Pacione C, Borgia J, Hurtig M, Oegema T, Chubinskaya S: Cellular responses in acute trauma of human ankle cartilage: cell survival, catabolic cytokines and neuromediators. Proceedings of the 7th World Congress of the International Cartilage Repair Society. 2007, Warsaw, Poland, 16 (2): 148-abstract
Chubinskaya S, Kawakami M, Matsumoto T, Tsutsui S, Rappoport L, Rueger DC: Immunohistochemical evidences of the anabolic and anti-catabolic effects of epidural OP-1 injection for herniated nucleus pulposus in the rat. Trans 54th ORS. 2008, San Francisco, CA, #3-
Cowell S, Knäuper V, Stewart ML, D'Ortho MP, Stanton H, Hembry RM, López-Otín C, Reynolds JJ, Murphy G: Induction of matrix metalloproteinase activation cascades based on membrane-type 1 matrix metalloproteinase: associated activation of gelatinase A, gelatinase B and collagenase 3. Biochem J. 1998, 331 ( Pt 2): 453-458.
Baker AH, Edwards DR, Murphy G: Metalloproteinase inhibitors: biological actions and therapeutic opportunities. J Cell Sci. 2002, 115: 3719-3727. 10.1242/jcs.00063.
Hu Y, Pham CT: Dipeptidyl peptidase I regulates the development of collagen-induced arthritis. Arthritis Rheum. 2005, 52: 2553-2558. 10.1002/art.21192.
Baici A, Lang A, Zwicky R, Müntener K: Cathepsin B in osteoarthritis: uncontrolled proteolysis in the wrong place. Semin Arthritis Rheum. 2005, 34: 24-28. 10.1016/j.semarthrit.2004.03.008.
Hurtig M, Chubinskaya S, Dickey J, Rueger D: BMP-7 protects against progression of cartilage degeneration after impact injury. J Orthop Res. 2009, 27: 602-611. 10.1002/jor.20787.
Chubinskaya S, Elshaier A, Tyree E, Nicholson A, Soeder S, Loeser RF, Hakimiyan A, Rueger DC: Identification of Id proteins in human articular cartilage and the effect of growth factors on their expression. Trans 51st ORS. 2005, Washington, DC, 30: 303-
Chubinskaya S, Hakimiyan A, Pacione C, Yanke A, Rappoport L, Aigner T, Rueger D, Loeser RF: Synergistic effect of IGF-1 and OP-1 on matrix formation by normal and OA chondrocytes cultured in alginate beads. Osteoarthritis Cart. 2007, 15: 421-430. 10.1016/j.joca.2006.10.004.
Loeser RF, Pacione CA, Chubinskaya S: The combination of insulin-like growth factor 1 and osteogenic protein 1 promotes increased survival of and matrix synthesis by normal and osteoarthritic human chondrocytes. Arthritis Rheum. 2003, 48: 2188-2196. 10.1002/art.11209.
Nishida Y, Knudson CB, Eger W, Kuettner KE, Knudson WL: Osteogenic protein 1 stimulates cells-associated matrix assembly by normal human articular chondrocytes: up-regulation of hyaluronan synthase, CD44, and aggrecan. Arthritis Rheum. 2000, 43: 206-214. 10.1002/1529-0131(200001)43:1<206::AID-ANR25>3.0.CO;2-1.
Vinall RL, Lo SH, Reddi AH: Regulation of articular chondrocyte phenotype by bone morphogenetic protein 7, interleukin 1, and cellular context is dependent on the cytoskeleton. Exp Cell Res. 2002, 272: 32-44. 10.1006/excr.2001.5395.
Aigner T, McKenna L, Zien A, Fan Z, Gebhard PM, Zimmer R: Gene expression profiling of serum- and interleukin-1β-stimulated primary human adult articular chodnrocytes - A molecular analysis based on chondrocytes isolated from one donor. Cytokine. 2005, 31: 227-240. 10.1016/j.cyto.2005.04.009.
Aigner T, Zien A, Gehrsitz A, Gebhard PM, McKenna LA: Anabolic and catabolic gene expression pattern analysis in normal versus osteoarthritic cartilage using complementary DNA-array technology. Arthritis Rheum. 2001, 44: 2777-2789. 10.1002/1529-0131(200112)44:12<2777::AID-ART465>3.0.CO;2-H.
Fan Z, Chubinskaya S, Rueger DC, Bau B, Haag J, Aigner T: Regulation of anabolic and catabolic gene expression in normal and osteoarthritic adult human articular chondrocytes by osteogenic protein-1. Clin Exp Rheumatol. 2004, 22: 103-106.
Flechtenmacher J, Huch K, Thonar EJ, Mollenhauer JA, Davies SR, Schmid TM, Puhl W, Sampath TK, Aydelotte MB, Kuettner KE: Recombinant human osteogenic protein 1 is a potent stimulator of the synthesis of cartilage proteoglycans and collagens by human articular chondrocytes. Arthritis Rheum. 1996, 39: 1896-1904. 10.1002/art.1780391117.
Bajaj S, Shoemaker T, Hakimiyan AA, Rappoport L, Pascual-Garrido C, Oegema TR, Wimmer MA, Chubinskaya S: Protective effect of P188 in the model of acute trauma to human ankle cartilage: the mechanism of action. J Ortho Trauma. 2010, 24: 571-576. 10.1097/BOT.0b013e3181ec4712.
Jennings L, Wu L, King KB, Hämmerle H, Cs-Szabo G, Mollenhauer J: The effects of collagen fragments on the extracellular matrix metabolism of bovine and human chondrocytes. Connect Tissue Res. 2001, 42: 71-86. 10.3109/03008200109014250.
Millward-Sadler SJ, Mackenzie A, Wright MO, Lee HS, Elliot K, Gerrard L, Fiskerstrand CE, Salter DM, Quinn JP: Tachykinin expression in cartilage and function in human articular chondrocyte mechanotransduction. Arthritis Rheum. 2003, 48: 146-156. 10.1002/art.10711.
Hunter DJ, Pike MC, Jonas BL, Kissin E, Krop J, McAlindon T: Phase I safety and tolerability study of BMP-7 in symptomatic knee osteoarthritis. BMC Musculoskeletal Disorders. 2010, 11: 232-239. 10.1186/1471-2474-11-232.
Tardif G, Pelletier J-P, Boileau C, Martel-Pelletier J: The BMP antagonists follistatin and gremlin in normal and early osteoarthritic cartilage: an immunohistochemical study. Osteoarthritis Cart. 2009, 17: 263-270. 10.1016/j.joca.2008.06.022.
Fukui N, Zhu Y, Malloney WJ, Clohisy J, Sandell LJ: Stimulation of BMP-2 expression by pro-inflammatory cytokines IL-1 and TNF-α in normal and osteoarthritic chondrocytes. J Bone Joint Surg Am. 2003, 85-A: 59-66.
Qureshi HY, Ricci G, Zafarullah M: Smad signaling pathway is a pivotal component of tissue inhibitor of metalloproteinases-3 regulation by transforming growth factor beta in human chondrocytes. Biochim Biophys Acta. 2008, 1783: 1605-1612. 10.1016/j.bbamcr.2008.04.005.
Morris KJ, Cs-Szabo G, Cole AA: Characterization of TIMP-3 in human articular talar cartilage. Connect Tissue Res. 2010, 51 (6): 478-490. 10.3109/03008201003686958.
Gendron C, Kashiwagi M, Hughes C, Caterson B, Nagase H: TIMP-3 inhibits aggrecanase-mediated glycosaminoglycan release from cartilage explants stimulated by catabolic factors. FEBS Lett. 2003, 555: 431-436. 10.1016/S0014-5793(03)01295-X.
The authors would like to acknowledge the Gift of Hope Organ & Tissue Donor Network and donors' families as well as Dr. Margulis for tissue procurement. The work was supported by the NIH grant AR 47654, Stryker Biotech grant SC-001, and Ciba-Geigy Endowed Chair (SC). RFL was supported by NIH grant AG016697.
Stryker Biotech provided partial research support and reagents. SC served on the scientific advisory board of Stryker Biotech in 2005-2006. DCR, co-author of the manuscript, was an employee of Stryker Biotech at the time the data were collected; he has retired and no longer works at Stryker Biotech. The other authors declare that they have no competing interests.
SC was the principle investigator on the project, developed conceptual idea, wrote the research proposal, obtained research funding and IRB approval, oversaw the progress of the study and acquisition of the project related data, coordinated efforts of the participants, wrote progress reports to the funding agencies (NIH, Stryker Biotech and Ciba-Geigy Endowed Chair), and was involved in the writing of the manuscript.
LO was the research assistant in the laboratory of Dr. Susan Chubinskaya. She performed all validation experiments, starting from cell culture and real-time PCR, analyzed gene array data and prepared the first draft of the manuscript. TA provided resources and performed gene expression analysis by Affimetrix gene array. SS was a postdoctoral fellow in the laboratory of the PI, Dr. Chubinskaya. He was solely responsible for all initial experiments: design of antisense oligonucleotides, cell isolation, culture, transfection, RNA isolation, and quality control of the RNA. JB is a collaborator on the project and assisted with validation experiments: probe design and real time PCR. He also was involved in the editing of the manuscript. DCR was a senior director of R&D at Stryker Biotech. He was involved in the conceptual development of the project, its objectives, specific aims, and experimental design. He also provided all necessary tools for the project: background knowledge on OP-1 gene, OP-1 cDNA library, recombinant OP-1, and OP-1 antibodies. He also was a main consultant on the project. RFL is a long-term collaborator of the PI, who is also collaborator on this project. All IGF-related experiments were done in Dr. Loeser's laboratory as a joint effort. He was also actively involved in data analysis and production of the manuscript as a final editor.
About this article
Cite this article
Chubinskaya, S., Otten, L., Soeder, S. et al. Regulation of chondrocyte gene expression by osteogenic protein-1. Arthritis Res Ther 13, R55 (2011). https://doi.org/10.1186/ar3300
- Nucleus Pulposus
- Gene Array
- Catabolic Gene
- Cartilage Homeostasis
- Smart Cycler