Microarray gene expression profiling of osteoarthritic bone suggests altered bone remodelling, WNT and transforming growth factor-β/bone morphogenic protein signalling

Osteoarthritis (OA) is characterized by alterations to subchondral bone as well as articular cartilage. Changes to bone in OA have also been identified at sites distal to the affected joint, which include increased bone volume fraction and reduced bone mineralization. Altered bone remodelling has been proposed to underlie these bone changes in OA. To investigate the molecular basis for these changes, we performed microarray gene expression profiling of bone obtained at autopsy from individuals with no evidence of joint disease (control) and from individuals undergoing joint replacement surgery for either degenerative hip OA, or fractured neck of femur (osteoporosis [OP]). The OP sample set was included because an inverse association, with respect to bone density, has been observed between OA and the low bone density disease OP. Compugen human 19K-oligo microarray slides were used to compare the gene expression profiles of OA, control and OP bone samples. Four sets of samples were analyzed, comprising 10 OA-control female, 10 OA-control male, 10 OA-OP female and 9 OP-control female sample pairs. Print tip Lowess normalization and Bayesian statistical analyses were carried out using linear models for microarray analysis, which identified 150 differentially expressed genes in OA bone with t scores above 4. Twenty-five of these genes were then confirmed to be differentially expressed (P < 0.01) by real-time PCR analysis. A substantial number of the top-ranking differentially expressed genes identified in OA bone are known to play roles in osteoblasts, osteocytes and osteoclasts. Many of these genes are targets of either the WNT (wingless MMTV integration) signalling pathway (TWIST1, IBSP, S100A4, MMP25, RUNX2 and CD14) or the transforming growth factor (TGF)-β/bone morphogenic protein (BMP) signalling pathway (ADAMTS4, ADM, MEPE, GADD45B, COL4A1 and FST). Other differentially expressed genes included WNT (WNT5B, NHERF1, CTNNB1 and PTEN) and TGF-β/BMP (TGFB1, SMAD3, BMP5 and INHBA) signalling pathway component or modulating genes. In addition a subset of genes involved in osteoclast function (GSN, PTK9, VCAM1, ITGB2, ANXA2, GRN, PDE4A and FOXP1) was identified as being differentially expressed in OA bone between females and males. Altered expression of these sets of genes suggests altered bone remodelling and may in part explain the sex disparity observed in OA.     

Gene expression analyses of subchondral bone in early experimental osteoarthritis by microarray

Osteoarthritis (OA) is a degenerative joint disease that affects both cartilage and bone. A better understanding of the early molecular changes in subchondral bone may help elucidate the pathogenesis of OA. We used microarray technology to investigate the time course of molecular changes in the subchondral bone in the early stages of experimental osteoarthritis in a rat model. We identified 2,234 differentially expressed (DE) genes at 1 week, 1,944 at 2 weeks and 1,517 at 4 weeks post-surgery. Further analyses of the dysregulated genes indicated that the events underlying subchondral bone remodeling occurred sequentially and in a time-dependent manner at the gene expression level. Some of the identified dysregulated genes that were identified have suspected roles in bone development or remodeling; these genes include Alp, Igf1, Tgf β1, Postn, Mmp3, Tnfsf11, Acp5, Bmp5, Aspn and Ihh. The differences in the expression of these genes were confirmed by real-time PCR, and the results indicated that our microarray data accurately reflected gene expression patterns characteristic of early OA. To validate the results of our microarray analysis at the protein level, immunohistochemistry staining was used to investigate the expression of Mmp3 and Aspn protein in tissue sections. These analyses indicate that Mmp3 protein expression completely matched the results of both the microarray and real-time PCR analyses; however, Aspn protein expression was not observed to differ at any time. In summary, our study demonstrated a simple method of separation of subchondral bone sample from the knee joint of rat, which can effectively avoid bone RNA degradation. These findings also revealed the gene expression profiles of subchondral bone in the rat OA model at multiple time points post-surgery and identified important DE genes with known or suspected roles in bone development or remodeling. These genes may be novel diagnostic markers or therapeutic targets for OA.     

Boning up on glypicans—opportunities for new insights into bone biology

Bone formation is remarkable for the convergence in the activity of four major signalling pathways, the bone morphogenetic protein (BMP), fibroblast growth factor (FGF), hedgehog (HH) and wingless‐integrated (WNT) pathways. These pathways cooperate in morphogenetic, proliferative and differentiative processes that underpin the development, growth and repair of skeletal structures. They are regulated by pathway‐specific modulators and by another class of molecules, the glypicans. Glypicans are proteoglycans located on the cell surface, where they act as coreceptors to promote or inhibit signalling by ligands of the BMP, FGF, HH and WNT pathways, through protein–protein and protein–carbohydrate interactions. In this review, we discuss glypican structure, expression and function in the context of bone development and growth, with emphasis on the long bone growth plate where five of the six glypicans are expressed in overlapping patterns in the chondrogenic zone. Analyses of gene knockout models and the human conditions of Simpson–Golabi–Behmel syndrome and omodysplasia, which arise from mutations in glypican 3 (GPC3) and GPC6, respectively, highlight both subtle and striking effects of glypicans on bone growth. We draw attention to challenges and areas of opportunity, where the actions of glypicans on BMP, FGF, HH and WNT signalling might be profitably studied to help illuminate the complex interplay of signalling that drives bone growth. Copyright © 2013 John Wiley & Sons, Ltd.     

Kinetics of gene expression and bone remodelling in the clinical phase of collagen-induced arthritis

IntroductionPathological bone changes differ considerably between inflammatory arthritic diseases and most studies have focused on bone erosion. Collagen-induced arthritis (CIA) is a model for rheumatoid arthritis, which, in addition to bone erosion, demonstrates bone formation at the time of clinical manifestations. The objective of this study was to use this model to characterise the histological and molecular changes in bone remodelling, and relate these to the clinical disease development.MethodsA histological and gene expression profiling time-course study on bone remodelling in CIA was linked to onset of clinical symptoms. Global gene expression was studied with a gene chip array system.ResultsThe main histopathological changes in bone structure and inflammation occurred during the first two weeks following the onset of clinical symptoms in the joint. Hereafter, the inflammation declined and remodelling of formed bone dominated.Global gene expression profiling showed simultaneous upregulation of genes related to bone changes and inflammation in week 0 to 2 after onset of clinical disease. Furthermore, we observed time-dependent expression of genes involved in early and late osteoblast differentiation and function, which mirrored the histopathological bone changes. The differentially expressed genes belong to the bone morphogenetic pathway (BMP) and, in addition, include the osteoblast markers integrin-binding sialoprotein (Ibsp), bone gamma-carboxyglutamate protein (Bglap1), and secreted phosphoprotein 1 (Spp1). Pregnancy-associated protein A (Pappa) and periostin (Postn), differentially expressed in the early disease phase, are proposed to participate in bone formation, and we suggest that they play a role in early bone formation in the CIA model. Comparison to human genome-wide association studies (GWAS) revealed differential expression of several genes associated with human arthritis.ConclusionsIn the CIA model, bone formation in the joint starts shortly after onset of clinical symptoms, which results in bony fusion within one to two weeks. This makes it a candidate model for investigating the relationship between inflammation and bone formation in inflammatory arthritis.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-015-0531-7) contains supplementary material, which is available to authorized users.     

RNA sequence analysis identified bone morphogenetic protein-2 (BMP2) as a biomarker underlying form deprivation myopia

PurposeForm deprivation myopia (FDM) is an urgent public issue characterized by pathological changes, but the underlying mechanism remained unclear. The aim was to investigate bone morphogenetic proteins (BMPs) utilizing the pathogenesis of FDM.Material and methodsGene expression omnibus (GEO) database was used to analyze one mRNA profile (GSE89325) of FDM. Sixteen retina samples (8 FDM and 8 controls) were randomly divided into seven groups for differential gene expression analysis in R. software. The gene pathway and protein-protein interaction (PPI) analysis were performed by the DAVID and STRING databases. Cytoscape was used to draw the PPI network. The gene ontology (GO) enrichment and Kyoto encyclopedia of genes and Genomes (KEGG) analysis were determined to achieve gene annotation and visualization.ResultsA total of 18420 differentially expressed genes (DEGs) were identified associated with FDM. The only non-significant gene (BEND6) was separately analyzed between two groups. Thirteen hub genes were discovered, ACVR1, ACVR2A, ACVR2B, RGMB, BMPR2, BMPR1A, BMP2, BMPR1B, CHRD, PTH, PTH1R, PTHLH, and WNT9A. The expression alteration in FDM were mainly enriched in cytokine-cytokine, and neuroactive ligand receptor interaction pathways. BMP2 was the key gene in myopia progression.ConclusionsOf clinical perspective, our findings reveal that expression of BMP2 as an underlying mechanism of FDM, providing an insight for therapeutic interventions.     

Large-scale gene expression profiles, differentially represented in osteoarthritic synovium of the knee joint using cDNA microarray technology

Osteoarthritis (OA) is one of the most common age-related chronic disorders of articular cartilage, joints and bone tissue. Diagnosis of OA commonly depends on clinical and radiographic findings. However, changes in cartilage associated with the early stage of OA cannot be detected using radiographs, because significant cartilage degeneration must occur before radiographic findings show alterations of the appearance of cartilage. To identify new biomarkers of OA, we analysed gene expression profiles of synovium from 43 patients with OA, ten patients with rheumatoid arthritis (RA), and eight non-OA/non-RA patients using a novel cDNA microarray chip. We identified 21 genes with simultaneous significant differences in expression between OA and non-OA/non-RA groups and between OA and RA groups. Linear discriminant analysis showed that the three groups could be well separated using those 21 genes. Statistical analysis also revealed that several of the 21 genes were associated with disease progression and clinical presentation. The graphical modelling method indicated that some of the 21 genes are significantly associated with a particular clinical presentation, suggesting biological relationships among those genes. This is the first report of the use of cDNA microarray technology to create large-scale gene expression profiles differentially expressed in situ in OA synovium of the knee joint.     

BMP inhibition stimulates WNT-dependent generation of chondrogenic mesoderm from embryonic stem cells

WNT and bone morphogenetic protein (BMP) signaling are known to stimulate hemogenesis from pluripotent embryonic stem (ES) cells. However, osteochondrogenic mesoderm was generated effectively when BMP signaling is kept to a low level, while WNT signaling was strongly activated. When mesoderm specification from ES cells was exogenous factor dependent, WNT3a addition supported the generation of cardiomyogenic cells expressing lateral plate/extraembryonic mesoderm genes, and this process involved endogenous BMP activities. Exogenous BMP4 showed a similar effect that depended on endogenous WNT activities. However, neither factor induced robust chondrogenic activity. In support, ES cell differentiation in the presence of either WNT3a or BMP4 was associated with elevated levels of both Bmp and Wnt mRNAs, which appeared to provide sufficient levels of active BMPs and WNTs to promote the nonchondrogenic mesoderm specification. The osteochondrogenic mesoderm expressed PDGFRα, which also expressed genes that mark somite and rostral presomitic mesoderm. A strong WNT signaling was required for generating the mesodermal progeny, while approximately 50- to 100-fold lower concentration of WNT3a was sufficient for specifying axial mes(end)oderm. Thus, depending on the dose and cofactor (BMP), WNT signaling stimulates the generation of different biological activities and specification of different types of mesodermal progeny from ES cells.     

Decrease in expression of bone morphogenetic proteins 4 and 5 in synovial tissue of patients with osteoarthritis and rheumatoid arthritis

Bone morphogenetic proteins (BMPs) have been identified as important morphogens with pleiotropic functions in regulating the development, homeostasis and repair of various tissues. The aim of this study was to characterize the expression of BMPs in synovial tissues under normal and arthritic conditions. Synovial tissue from normal donors (ND) and from patients with osteoarthritis (OA) and rheumatoid arthritis (RA) were analyzed for BMP expression by using microarray hybridization. Differential expression of BMP-4 and BMP-5 was validated by semiquantitative RT-PCR, in situ hybridization and immunohistochemistry. Activity of arthritis was determined by routine parameters for systemic inflammation, by histological scoring of synovitis and by semiquantitative RT-PCR of IL-1β, TNF-α, stromelysin and collagenase I in synovial tissue. Expression of BMP-4 and BMP-5 mRNA was found to be significantly decreased in synovial tissue of patients with RA in comparison with ND by microarray analysis (p < 0.0083 and p < 0.0091). Validation by PCR confirmed these data in RA (p < 0.002) and also revealed a significant decrease in BMP-4 and BMP-5 expression in OA compared with ND (p < 0.015). Furthermore, histomorphological distribution of both morphogens as determined by in situ hybridization and immunohistochemistry showed a dominance in the lining layer of normal tissues, whereas chronically inflamed tissue from patients with RA revealed BMP expression mainly scattered across deeper layers. In OA, these changes were less pronounced with variable distribution of BMPs in the lining and sublining layer. BMP-4 and BMP-5 are expressed in normal synovial tissue and were found decreased in OA and RA. This may suggest a role of distinct BMPs in joint homeostasis that is disturbed in inflammatory and degenerative joint diseases. In comparison with previous reports, these data underline the complex impact of these factors on homeostasis and remodeling in joint physiology and pathology.     

Activation of the pro-migratory bone morphogenetic protein receptor 1B gene in human MDA-MB-468 triple-negative breast cancer cells that over-express CYP2J2

Secondary metastases are the leading cause of mortality in patients with breast cancer. Cytochrome P450 (CYP) 2J2 (CYP2J2) is upregulated in many human tumors and generates epoxyeicosanoids from arachidonic acid that promote tumorigenesis and metastasis, but at present there is little information on the genes that mediate these actions. In this study MDA-MB-468 breast cancer cells were stably transfected with CYP2J2 (MDA-2J2 cells) and Affymetrix microarray profiling was undertaken. We identified 182 genes that were differentially expressed in MDA-2J2 cells relative to control (MDA-CTL) cells (log[fold of control] ≥2). From gene ontology pathway analysis bone morphogenetic protein (BMP) receptor 1B (BMPR1B) emerged as an important upregulated gene in MDA-2J2 cells. Addition of the BMPR1B ligand BMP2 stimulated the migration of MDA-2J2 cells, but not MDA-CTL cells, from 3D-matrigel droplets. Migration of MDA-2J2 cells was prevented by the BMPR antagonist dorsomorphin. These findings indicate that over-expression of CYP2J2 in MDA-MB-468-derived breast cancer cells activates BMPR1B expression that may contribute to increased migration. Targeting BMPR1B may be a novel approach to inhibit the metastatic activity of breast cancers that contain high levels of CYP2J2.     

Large-scale gene expression profiles,differentially represented in osteoarthritic synovium of the knee joint using cDNA microarray technology

Osteoarthritis (OA) is one of the most common age-related chronic disorders of articular cartilage, joints and bone tissue. Diagnosis of OA commonly depends on clinical and radiographic findings. However, changes in cartilage associated with the early stage of OA cannot be detected using radiographs, because significant cartilage degeneration must occur before radiographic findings show alterations of the appearance of cartilage. To identify new biomarkers of OA, we analysed gene expression profiles of synovium from 43 patients with OA, ten patients with rheumatoid arthritis (RA), and eight non-OA/non-RA patients using a novel cDNA microarray chip. We identified 21 genes with simultaneous significant differences in expression between OA and non-OA/non-RA groups and between OA and RA groups. Linear discriminant analysis showed that the three groups could be well separated using those 21 genes. Statistical analysis also revealed that several of the 21 genes were associated with disease progression and clinical presentation. The graphical modelling method indicated that some of the 21 genes are significantly associated with a particular clinical presentation, suggesting biological relationships among those genes. This is the first report of the use of cDNA microarray technology to create large-scale gene expression profiles differentially expressed in situ in OA synovium of the knee joint.     

Up-regulation of the Wnt, estrogen receptor, insulin-like growth factor-I, and bone morphogenetic protein pathways in C57BL/6J osteoblasts as opposed to C3H/HeJ osteoblasts in part contributes to the differential anabolic response to fluid shear

C57BL/6J (B6), but not C3H/HeJ (C3H), mice responded to mechanical loading with an increase in bone formation. A 30-min steady fluid shear of 20 dynes/cm(2) increased [(3)H]thymidine incorporation and alkaline phosphatase activity and up-regulated the expression of early mechanoresponsive genes (integrin beta1 (Igtb1) and cyclooxygenase-2 (Cox-2)) in B6 but not C3H osteoblasts, indicating that the differential mechanosensitivity was intrinsic to osteoblasts. In-house microarray analysis with 5,500 gene fragments revealed that the expression of 669 genes in B6 osteoblasts and 474 genes in C3H osteoblasts was altered 4 h after the fluid shear. Several genes associated with the insulin-like growth factor (IGF)-I, the estrogen receptor (ER), the bone morphogenetic protein (BMP)/transforming growth factor-beta, and Wnt pathways were differentially up-regulated in B6 osteoblasts. In vitro mechanical loading also led to up-regulation of these genes in the bones of B6 but not C3H mice. Pretreatment of B6 osteoblasts with inhibitors of the Wnt pathway (endostatin), the BMP pathway (Noggin), or the ER pathway (ICI182780) blocked the fluid shear-induced proliferation. Inhibition of integrin and Cox-2 activation by echistatin and indomethacin, respectively, each blocked the fluid shear-induced up-regulation of genes associated with these four pathways. In summary, up-regulation of the IGF-I, ER, BMP, and Wnt pathways is involved in mechanotransduction. These four pathways are downstream to the early mechanoresponsive genes, i.e. Igtb1 and Cox-2. In conclusion, differential up-regulation of these anabolic pathways may in part contribute to the good and poor response, respectively, in the B6 and C3H mice to mechanical loading.     

Free fatty acid receptor 4-β-arrestin 2 pathway mediates the effects of different classes of unsaturated fatty acids in osteoclasts and osteoblasts

Bone is a dynamic tissue that is constantly remodelled by bone resorbing osteoclasts and bone forming osteoblasts, respectively. A breakdown in the remodelling process underlies several bone diseases such as osteoporosis. Unsaturated fatty acids (UFAs) have been shown to have beneficial effects on bone health. However, the mechanism of action of UFAs in bone remains unclear. Free fatty acid receptor 4 (FFAR4) is expressed in bone cells and preferentially binds ω−3 and ω−7 UFAs. Therefore, we sought to determine if FFAR4 influenced the action of different classes of UFAs in bone cells. FFAR4 and potential signalling pathways, β-arrestin 2 (βarr2) and G_αq, were silenced in RAW264.7 murine macrophages (pre-osteoclasts) and MC3T3-E1 murine pre-osteoblasts. Cell differentiation, activation of signalling pathways and expression of regulatory genes were evaluated. The ω−3 UFAs, docosahexaenoic acid (DHA) and eicosapentaenoic acid (EPA), and the ω−7 UFA, palmitoleic acid (PLA), were shown to require the FFAR4/βarr2 signalling pathway to inhibit osteoclast differentiation in RAW264.7 murine macrophages. The ω−6 UFA, arachidonic acid, and the ω−9 UFA, oleic acid (OA), were shown to inhibit osteoclast formation but did not use FFAR4. DHA, EPA, PLA and OA enhanced osteoblast signalling through the FFAR4/βarr2 signalling axis. This study reveals that FFAR4/βarr2 signalling may mediate the bone protective effects of different classes of UFAs in osteoclasts and osteoblasts.     

Up‐regulated CST5 inhibits bone resorption and activation of osteoclasts in rat models of osteoporosis via suppression of the NF‐κB pathway

Here, we aim at exploring the effect of CST5 on bone resorption and activation of osteoclasts in osteoporosis (OP) rats through the NF‐κB pathway. Microarray analysis was used to screen the OP‐related differentially expressed genes. Osteoporosis was induced in rats by intragastric retinoic acid administration. The serum levels of tartrate‐resistant acid phosphatase (TRAP), bone alkaline phosphatase (BALP) and osteocalcin (OC) and the expression of CD61 on the surface of osteoclasts were examined. The number of osteoclasts and the number and area of resorption pits were detected. Besides, the pathological changes and bone mineral density in bone tissues of rats were assessed. Also, the relationship between CST5 and the NF‐κB pathway was identified through determining the expression of CST5, RANKL, RANK, OPG, p65 and IKB. Poorly expressed CST5 was indicated to affect the OP. CST5 elevation and inhibition of the NF‐κB pathway decreased serum levels of TRAP, BALP and OC and expression of CD61 in vivo and in vitro. In OP rats, CST5 overexpression increased trabecular bones and bone mineral density of bone tissues, but decreased trabecular separation, fat within the bone marrow cavities and the number of osteoclasts through inhibiting the NF‐κB pathway. In vivo experiments showed that CST5 elevation inhibited growth in number and area of osteoclastic resorption pits and restrained osteoclastic bone absorption by inhibiting the NF‐κB pathway. In summary, overexpression of CST5 suppresses the activation and bone resorption of osteoclasts by inhibiting the activation of the NF‐κB pathway.