Identification of upstream regulators for synovial expression signature genes in osteoarthritis

Axenfeld-Rieger syndrome: more than meets the eye

BACKGROUND

Axenfeld-Rieger syndrome (ARS) is characterised by typical anterior segment anomalies, with or without systemic features. The discovery of causative genes identified ARS subtypes with distinct phenotypes, but our understanding is incomplete, complicated by the rarity of the condition.

METHODS

Genetic and phenotypic characterisation of the largest reported ARS cohort through comprehensive genetic and clinical data analyses.

RESULTS

128 individuals with causative variants in PITX2 or FOXC1, including 81 new cases, were investigated. Ocular anomalies showed significant overlap but with broader variability and earlier onset of glaucoma for FOXC1-related ARS. Systemic anomalies were seen in all individuals with PITX2-related ARS and the majority of those with FOXC1-related ARS. PITX2-related ARS demonstrated typical umbilical anomalies and dental microdontia/hypodontia/oligodontia, along with a novel high rate of Meckel diverticulum. FOXC1-related ARS exhibited characteristic hearing loss and congenital heart defects as well as previously unrecognised phenotypes of dental enamel hypoplasia and/or crowding, a range of skeletal and joint anomalies, hypotonia/early delay and feeding disorders with structural oesophageal anomalies in some. Brain imaging revealed highly penetrant white matter hyperintensities, colpocephaly/ventriculomegaly and frequent arachnoid cysts. The expanded phenotype of FOXC1-related ARS identified here was found to fully overlap features of De Hauwere syndrome. The results were used to generate gene-specific management plans for the two types of ARS.

CONCLUSION

Since clinical features of ARS vary significantly based on the affected gene, it is critical that families are provided with a gene-specific diagnosis, PITX2-related ARS or FOXC1-related ARS. De Hauwere syndrome is proposed to be a FOXC1opathy.

NR4A1-3 nuclear receptor activity and immune cell dysregulation in rheumatic diseases

The development and progression of immune-mediated rheumatic disease (IMRD) involves dysfunction of innate and adaptive immune cell populations leading to altered responses including inflammasome activation, dysregulated cytokine networks, increased immune cell numbers and multifaceted cell-cell communication. Several rheumatic diseases are further characterized by the presence of autoantibodies, immune complex mediated complement activation and the deficit of peripheral immune tolerance due to reduced regulatory T-lymphocyte cell function. Ultimately, in rheumatic disease the loss in cellular and tissue homeostasis culminates in the advancement of chronic inflammation. The three members of the NR4A subfamily of nuclear receptors are immediate early genes, and act as potent transcriptional responders to changes in the cellular and tissue microenvironment. Subfamily members are rapidly expressed in diseases characterized by inflammation and function to control the differentiation and activity of innate and adaptive immune cells in a cell-type and cell-context specific manner. Rheumatic disease including rheumatoid-, psoriatic-, osteo-arthritis and systemic sclerosis display altered NR4A1-3 activity in controlling immune cell migration and function, production of paracrine signaling molecules, synovial tissue hyperplasia, and regulating cartilage turn-over in vivo. Additionally, NR4A1-3 activities mediate cytokine, prostanoid and growth factor signaling to control angiogenesis, modulate the regulatory functions of mesenchymal stromal cells, alter the activation status of dendritic cells, influence the generation of peripheral myeloid and T-lymphocyte lineages and promote the maintenance of functional regulatory T-cells. Further reports uncover the potential of moderating NR4A 1-3 receptors as therapeutic targets in altering immune tolerance, pathological angiogenesis and controlling inflammation in several models of disease.

Orphan Nuclear Receptor NR4A2 Is Constitutively Expressed in Cartilage and Upregulated in Inflamed Synovium From hTNF-Alpha Transgenic Mice

Orphan nuclear receptor 4A2 (NR4A2/Nurr1) is a constitutively active transcription factor with potential roles in the onset and progression of inflammatory arthropathies. NR4A2 is overexpressed in synovium and cartilage from individuals with rheumatoid arthritis (RA), psoriatic arthritis, and osteoarthritis. This study documents the expression and tissue localization of NR4A2 and upstream regulator nuclear factor kappa B (NF-κB) in the human tumor necrosis factor-alpha (hTNF-α) transgenic mouse model of RA. Since TNF-α is a potent inducer of NR4A2 in vitro, we hypothesized that NR4A2 would also be upregulated and active during disease progression in this model. Expression levels of NR4A2, related receptors NR4A1 (Nur77) and 3 (NOR1), and NF-κB₁ transcripts were quantified by RT-qPCR in hTNF-α and wild-type joints at three stages of disease. The protein distribution of NR4A2 and NF-κB subunit RelA (p65) was analyzed by quantitative immunohistochemistry. Global gene expression of 88 RA-related genes was also screened and compared between groups. Consistent with previous reports on the hTNF-α model, transgenic mice exhibited significant weight loss and severely swollen paws by 19 weeks of age compared to age-matched wild-type controls. NR4A1-3 and NF-κB₁ were constitutively expressed at disease onset and in healthy joints. NF-κB₁ transcript levels increased 2-fold in hTNF-α paws with established disease (12 weeks), followed by a 2-fold increase in NR4A2 at the late disease stage (19 weeks). NR4A2 and RelA proteins were overexpressed in inflamed synovium prior to symptoms of arthritis, suggesting that gene expression changes documented in whole paws were largely driven by elevated expression in diseased synovium. Broader screening of RA-related genes by RT-qPCR identified several differentially expressed genes in hTNF-α joints including those encoding inflammatory cytokines and chemokines, matrix-degrading enzymes and inhibitors, cell surface receptors, intracellular signaling proteins and transcription factors. Consensus binding sites for NR4A receptors and NF-κB₁ were enriched in the promoters of differentially expressed genes suggesting central roles for these transcription factors in this model. This study is the first comprehensive analysis of NR4A2 in an animal model of RA and validates the hTNF-α model for testing of small molecules and genetic strategies targeting this transcription factor.

Accelerating functional gene discovery in osteoarthritis

Osteoarthritis causes debilitating pain and disability, resulting in a considerable socioeconomic burden, yet no drugs are available that prevent disease onset or progression. Here, we develop, validate and use rapid-throughput imaging techniques to identify abnormal joint phenotypes in randomly selected mutant mice generated by the International Knockout Mouse Consortium. We identify 14 genes with functional involvement in osteoarthritis pathogenesis, including the homeobox gene Pitx1, and functionally characterize 6 candidate human osteoarthritis genes in mouse models. We demonstrate sensitivity of the methods by identifying age-related degenerative joint damage in wild-type mice. Finally, we phenotype previously generated mutant mice with an osteoarthritis-associated polymorphism in the Dio2 gene by CRISPR/Cas9 genome editing and demonstrate a protective role in disease onset with public health implications. We hope this expanding resource of mutant mice will accelerate functional gene discovery in osteoarthritis and offer drug discovery opportunities for this common, incapacitating chronic disease.

Integration of transcriptome-wide association study and messenger RNA expression profile to identify genes associated with osteoarthritis

Aims

Osteoarthritis (OA) is the most prevalent joint disease. However, the specific and definitive genetic mechanisms of OA are still unclear.

Methods

Tissue-related transcriptome-wide association studies (TWAS) of hip OA and knee OA were performed utilizing the genome-wide association study (GWAS) data of hip OA and knee OA (including 2,396 hospital-diagnosed hip OA patients versus 9,593 controls, and 4,462 hospital-diagnosed knee OA patients versus 17,885 controls) and gene expression reference to skeletal muscle and blood. The OA-associated genes identified by TWAS were further compared with the differentially expressed genes detected by the messenger RNA (mRNA) expression profiles of hip OA and knee OA. Functional enrichment and annotation analysis of identified genes was performed by the DAVID and FUMAGWAS tools.

Results

We detected 33 common genes, eight common gene ontology (GO) terms, and one common pathway for hip OA, such as calcium and integrin-binding protein 1 (CIB1) (P_TWAS = 0.025, FC_mRNA = -1.575 for skeletal muscle), adrenomedullin (ADM) (P_TWAS = 0.022, FC_mRNA = -4.644 for blood), Golgi apparatus (P_TWAS <0.001, P_mRNA = 0.012 for blood), and phosphatidylinositol 3' -kinase-protein kinase B (PI3K-Akt) signalling pathway (P_TWAS = 0.033, P_mRNA = 0.005 for blood). For knee OA, we detected 24 common genes, eight common GO terms, and two common pathways, such as histocompatibility complex, class II, DR beta 1 (HLA-DRB1) (P_TWAS = 0.040, FC_mRNA = 4.062 for skeletal muscle), Follistatin-like 1 (FSTL1) (P_TWAS = 0.048, FC_mRNA = 3.000 for blood), cytoplasm (P_TWAS < 0.001, P_mRNA = 0.005 for blood), and complement and coagulation cascades (P_TWAS = 0.017, P_mRNA = 0.001 for skeletal muscle).

Conclusion

We identified a group of OA-associated genes and pathways, providing novel clues for understanding the genetic mechanism of OA.

Cite this article:Bone Joint Res. 2020;9(3):130–138.

Sex Differences in Osteoarthritis Pathogenesis: A Comprehensive Study Based on Bioinformatics

Background

Osteoarthritis (OA) is a common disorder in the elderly. OA influences the daily life of patients and has become a worldwide health problem. It is still unclear whether the pathogenesis mechanism is different between males and females. This study investigated the differentially expressed genes (DEGs) and explored the different signaling pathways of OA between males and females.

Material/Methods

Data sets of GSE55457, GSE55584, and GSE12021 were retrieved from Gene Expression Omnibus to conduct DEGs analysis. Enrichment analysis of Kyoto Encyclopedia of Genes and Genomes pathway and Gene Ontology term was performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID) bioinformatics tool. The protein interaction network was constructed in Cytoscape 3.7.2. qRT-PCR was then performed to validate the expression of hub genes in OA patients and healthy people.

Results

In total, 4 co-upregulated and 10 co-downregulated genes were identified. We found that enriched pathways were different between males and females. BCL2L1, EEF1A1, EEF2, HNRNPD, and PABPN1 were considered as hub genes in OA pathogenesis in males, while EEF2, EEF1A1, RPL37A, FN1 were considered as hub genes in OA pathogenesis in females. Consistent with the bioinformatics analysis, the qRT-PCR analysis also showed that the gene expression of BCL2L1, HNRNPD, and PABPN1 was significantly lower in male OA patients. In contrast, EEF2, EEF1A1, and RPL37A were significantly lower in female OA patients.

Conclusions

The DEGs identified may be involved in different OA disease progression mechanisms between males and females, and they are considered as treatment targets or prognosis markers for males and females. The pathogenesis mechanism is sex-dependent.

Ingenuity pathway analysis of human facet joint tissues: Insight into facet joint osteoarthritis

Facet joint osteoarthritis (FJOA) is a common degenerative joint disorder with high prevalence in the elderly. FJOA causes lower back pain and lower extremity pain, and thus severely impacts the quality of life of affected patients. Emerging studies have focused on the histomorphological and histomorphometric changes in FJOA. However, the dynamic genetic changes in FJOA have remained to be clearly determined. In the present study, previously obtained RNA deep sequencing data were subjected to an ingenuity pathway analysis (IPA) and canonical signaling pathways of differentially expressed genes (DEGs) in FJOA were studied. The top 25 enriched canonical signaling pathways were identified and canonical signaling pathways with high absolute values of z-scores, specifically leukocyte extravasation signaling, Tec kinase signaling and osteoarthritis pathway, were investigated in detail. DEGs were further categorized by disease, biological function and toxicity (tox) function. The genetic networks between DEGs as well as hub genes in these functional networks were also investigated. It was demonstrated that C-X-C motif chemokine ligand 8, elastase, neutrophil expressed, growth factor independent 1 transcriptional repressor, Spi-1 proto-oncogene, CCAAT enhancer binding protein epsilon, GATA binding protein 1, TAL bHLH transcription factor 1, erythroid differentiation factor, minichromosome maintenance complex component 4, BTG anti-proliferation factor 2, BRCA1 DNA repair-associated, cyclin D1, chromatin assembly factor 1 subunit A, triggering receptor expressed on myeloid cells 1 and tumor protein p63 were hub genes in the top 5 IPA networks (with a score >30). The present study provides insight into the pathological processes of FJOA from a genetic perspective and may thus benefit the clinical treatment of FJOA.

Genetic dissection of canine hip dysplasia phenotypes and osteoarthritis reveals three novel loci

Background

Hip dysplasia and osteoarthritis continue to be prevalent problems in veterinary and human medicine. Canine hip dysplasia is particularly problematic as it massively affects several large-sized breeds and can cause a severe impairment of the quality of life. In Finland, the complex condition is categorized to five classes from normal to severe dysplasia, but the categorization includes several sub-traits: congruity of the joint, Norberg angle, subluxation degree of the joint, shape and depth of the acetabulum, and osteoarthritis. Hip dysplasia and osteoarthritis have been proposed to have separate genetic etiologies.

Results

Using Fédération Cynologique Internationale -standardized ventrodorsal radiographs, German shepherds were rigorously phenotyped for osteoarthritis, and for joint incongruity by Norberg angle and femoral head center position in relation to dorsal acetabular edge. The affected dogs were categorized into mild, moderate and severe dysplastic phenotypes using official hip scores. Three different genome-wide significant loci were uncovered. The strongest candidate genes for hip joint incongruity were noggin (NOG), a bone and joint developmental gene on chromosome 9, and nanos C2HC-type zinc finger 1 (NANOS1), a regulator of matrix metalloproteinase 14 (MMP14) on chromosome 28. Osteoarthritis mapped to a long intergenic region on chromosome 1, between genes encoding for NADPH oxidase 3 (NOX3), an intriguing candidate for articular cartilage degradation, and AT-rich interactive domain 1B (ARID1B) that has been previously linked to joint laxity.

Conclusions

Our findings highlight the complexity of canine hip dysplasia phenotypes. In particular, the results of this study point to the potential involvement of specific and partially distinct loci and genes or pathways in the development of incongruity, mild dysplasia, moderate-to-severe dysplasia and osteoarthritis of canine hip joints. Further studies should unravel the unique and common mechanisms for the various sub-traits.

Identification of key gene modules and transcription factors for human osteoarthritis by weighted gene co-expression network analysis

Osteoarthritis (OA) is one of the most prevalent causes of joint disease. However, the pathological mechanisms of OA have remained to be completely elucidated, and further investigation into the underlying mechanisms of OA development and the identification of novel therapeutic targets are urgently required. In the present study, the dataset GSE114007 was downloaded from the Gene Expression Omnibus database. Based on weighted gene co-expression network analysis (WGCNA) and the identification of differentially expressed genes (DEGs), the microarray data were further analyzed to identify hub genes, key transcription factors (TFs) and pivotal signaling pathways involved in the pathogenesis of OA. A total of 1,898 genes were identified to be differentially expressed between OA samples and normal samples. Based on WGCNA, the present study identified 5 hub modules closely associated with OA, and the potential key TFs for hub modules were further explored based on CisTargetX. The results demonstrated that B-Cell Lymphoma 6, Myelin Gene Expression Factor 2, Activating Transcription Factor 3, CCAAT Enhancer Binding Protein γ, Nuclear Factor Interleukin-3-Regulated, FOS Like Antigen-2, FOS-Like Antigen-1, Fos Proto-Oncogene, JunD Proto-Oncogene, Transcription Factor CP2 Like 1, RELA proto-oncogene NF-kB subunit, SRY-box transcription factor 3, V-Ets Avian Erythroblastosis Virus E26 Oncogene Homolog 2, Interferon Regulatory Factor 4 and REL proto-oncogene, NF-kB subunit were the potential key TFs. In addition, osteoclast differentiation, FoxO, MAPK and PI3K/Akt signaling pathways were revealed to be imperative for the pathogenesis of OA, as these 4 pivotal signaling pathways were observed to be tightly linked through 4 key TFs Fos Proto-Oncogene, JUN, JunD Proto-Oncogene and MYC, and 4 DEGs Vascular Endothelial Growth Factor A, Growth Arrest and DNA Damage Inducible α, Growth Arrest and DNA Damage Inducible β and Cyclin D1. The present study identified a set of potential key genes and signaling pathways, and provided an important opportunity to advance the current understanding of OA.

EGR1 promotes the cartilage degeneration and hypertrophy by activating the Krüppel-like factor 5 and β-catenin signaling

Osteoarthritis is one of the most common orthopedic diseases in elderly people who have lost their mobility. In this study，we observed abnormally high EGR1 expression in the articular cartilage of patients with osteoarthritis. We also found significantly high EGR1 expression in the articular cartilage of mice with destabilized medial meniscus (DMM)-induced osteoarthritis and 20-month-old mice. In vitro experiments indicated that IL-1β could significantly enhance EGR1 expression in primary mouse chondrocytes. EGR1 over-expression in chondrocytes using adenovirus could inhibit COl2A1 expression and enhance MMP9 and MMP13 expression. And silencing EGR1, using RNAi, had the opposite effects. Moreover, EGR1 over-expression accelerated chondrocyte hypertrophy in vitro, and EGR1 knockdown reversed this effect. We then explored the underlying mechanism. EGR1 over-expression increased Kruppel-Like Factor 5 (KLF5) protein level without influencing its synthesis. Enhanced EGR1 expression induced its integration with KLF5, leading to suppressed ubiquitination of KLF5. Moreover, EGR1 prompted β-catenin nuclear transportation to control chondrocyte hypertrophy. Ectopic expression of EGR1 in articular cartilage aggravated the degradation of the cartilage matrix in vivo. The EGR1 inhibitor, ML264, protected chondrocytes from IL-1β-mediated cartilage matrix degradation in vitro and DMM-induced osteoarthritis in vivo. Above all, we demonstrate the effect and mechanisms of EGR1 on osteoarthritis and provide evidence that the ML264 might be a potential drug for treating osteoarthritis in the future.

Identification of Differentially Expressed Genes under the Regulation of Transcription Factors in Osteosarcoma

The present study was to investigate and identify the differentially expressed genes (DEGs) in the transcriptional regulatory network of osteosarcoma (OS). The gene expression dataset from Gene Expression Omnibus (GEO) datasets was downloaded. DEGs were identified and their functional annotation was also conducted. In addition, differentially expressed transcription factors (TFs) and the regulatory genes were identified. The electronic validation was used to verify the expression of selected genes. The integrated analysis led to 932 DEGs. The results of functional annotation indicated that these DEGs significantly enriched in the p53 signaling pathway, Jak-STAT signaling pathway and Wnt signaling pathway. ZNF354C, NFIC, NFATC2, SP2, FOXO3, EGR1, ZEB1, RREB1, EGR2 and SRF were covered by most TFs. The expression levels of NFIC and EGR2 in electronic validation were compatible with our bio-informatics result. In conclusion, the deregulation of these genes may provide valuable information in understanding the underlying molecular mechanism in the OS.

Analysis of gene expression and methylation datasets identified ADAMTS9, FKBP5, and PFKBF3 as biomarkers for osteoarthritis

BACKGROUND

Osteoarthritis (OA) is a kind of chronic osteoarthropathy and degenerative joint disease. Epigenetic regulation in the gene expression dynamics has become increasingly important in OA. We performed a combined analysis of two types of microarray datasets (gene expression and DNA methylation) to identify methylation-based key biomarkers to provide a better understanding of molecular biological mechanisms of OA.

METHODS

We obtained two expression profiling datasets (GSE55235, GSE55457) and one DNA methylation profiling data set (GSE63695) from the Gene Expression Omnibus. First, differentially expressed genes (DEGs) between patients with OA and controls were identified using the Limma package in R(v3.4.4). Then, function enrichment analysis of DEGs was performed using a DAVID database. For DNA methylation datasets, ChAMP methylation analysis package was used to identify differential methylation genes (DMGs). Finally, a comprehensive analysis of DEGs and DMGs was conducted to identify genes that exhibited differential expression and methylation simultaneously.

RESULTS

We identified 112 DEGs and 2,896 DMGs in patients with OA compared with controls. Functional analysis of DEGs obtained that inflammatory responses, immune responses, and positive regulation of apoptosis, tumor necrosis factor (TNF) signaling pathway, and osteoclast differentiation may be involved in the pathogenesis of OA. Cross-analysis revealed 26 genes that exhibited differential expression and methylation in OA. Among them, ADAMTS9, FKBP5, and PFKBF3 were identified as valuable methylation-based biomarkers for OA.

CONCLUSION

In summary, our study identified different molecular features between patients with OA and controls. This may provide new clues for clarifying the pathogenetic mechanisms of OA.

Integrated microarray analysis provided novel insights to the pathogenesis of glaucoma

Glaucoma is characterized as a visual field defect, which is the second most common cause of blindness. The present study performed an integrated analysis of microarray studies of glaucoma derived from Gene Expression Omnibus (GEO). Following the identification of the differentially expressed genes (DEGs) in glaucoma compared with normal control (NC) tissues, the functional annotation, glaucoma-specific protein-protein interaction (PPI) network and transcriptional regulatory network constructions were performed. The acute intraocular pressure (IOP) elevation rat models were established and reverse transcription-quantitative polymerase chain reaction (RT-qPCR) was performed for DEGs expression confirmation. Three datasets were downloaded from GEO. A total of 97 DEGs, 82 upregulated and 15 downregulated were identified in glaucoma compared with NC groups with false discovery rate <0.05. Response to virus and immune response were two significantly enriched GO terms in glaucoma. Valine, leucine and isoleucine degradation was a significantly enriched pathway of DEGs in glaucoma. According to the PPI network, HDAC1, HBN, UBR4 and PDK1 were hub proteins in glaucoma. FOXD3, HNF-4 and AP-1 were the three transcription factors (TFs) derived from top 10 TFs which covered the majority of downstream DEGs in glaucoma. Based on the RT-qPCR results, the expression levels of 3 DEGs, raftlin, lipid raft linker 1 (RFTN1), PBX homeobox 1 (PBX1), HDAC1 were significantly upregulated and the expression of GEM was significantly downregulated in acute IOP elevation rat model at the first and fifth day. These four DEGs had the same expression pattern with our integrated analysis. Therefore, the current study concluded that 6 DEGs, including HEPH, SELENBP1, RFTN1, ID1, HDAC-1 and PBX1 and three TFs, including FOXD3, HNF-4 and AP-1 may be involved with the pathogenesis of glaucoma. The findings of the current study may improve diagnosis and drug design for glaucoma.

E2F1 and TFDP1 Regulate PITX1 Expression in Normal and Osteoarthritic Articular Chondrocytes

We previously reported a loss-of-PITX1 expression in patients suffering of knee/hip osteoarthritis (OA). Search for the mechanism underlying this event led us to discover that PITX1 repression was triggered by the aberrant nuclear accumulation of Prohibitin (PHB1), an E2F1 co-repressor, in OA articular chondrocytes. In the current study, we assessed in details the involvement of E2F transcription factors in regulating PITX1 expression. We also analyzed other genes that are similarly regulated by E2F in regard to osteoarthritis. The transcriptional regulation of the PITX1 promoter by E2F1 was analyzed with the luciferase reporter assay, and chromatin immunoprecipitation assays, which confirmed direct E2F1-PITX1 interactions. The probable binding sites for E2F1 in the PITX1 promoter were identified by DNA pulldown experiments. In silico and in vitro analyses show that the PITX1 proximal promoter region contains 2 specific sequences that are bound by E2F1. Overexpression of E2F1 enhances PITX1 promoter activity and mRNA transcription. In primary control and osteoarthritis chondrocytes, real time RT-PCR was used to measure the mRNA expression levels of candidate genes under E2F1 transcriptional control. Transcription Factor Dp-1 (TFDP1) knockdown experiments confirmed that the E2F1-TFDP1 complex regulates PITX1. Knockdown of TFDP1, an E2F1 dimerization partner, inhibits the activating effect of E2F1 and reduces both PITX1 promoter activity and mRNA transcription. Real time RT-PCR results reveal reduced expression of TFDP1 and a similar downregulation of their targets PITX1, BRCA1, CDKN1A, and RAD51 in mid-stage OA chondrocytes. Collectively, our data define a previously uncharacterized role for E2F1 and TFDP1 in the transcriptional regulation of PITX1 in articular chondrocytes. Additional E2F1 targets may be affected in OA pathogenesis.

Silencing of microRNA-138-5p promotes IL-1β-induced cartilage degradation in human chondrocytes by targeting FOXC1: miR-138 promotes cartilage degradation

Objectives

Osteoarthritis (OA) is characterised by articular cartilage degradation. MicroRNAs (miRNAs) have been identified in the development of OA. The purpose of our study was to explore the functional role and underlying mechanism of miR-138-5p in interleukin-1 beta (IL-1β)-induced extracellular matrix (ECM) degradation of OA cartilage.

Materials and Methods

Human articular cartilage was obtained from patients with and without OA, and chondrocytes were isolated and stimulated by IL-1β. The expression levels of miR-138-5p in cartilage and chondrocytes were both determined. After transfection with miR-138-5p mimics, allele-specific oligonucleotide (ASO)-miR-138-5p, or their negative controls, the messenger RNA (mRNA) levels of aggrecan (ACAN), collagen type II and alpha 1 (COL2A1), the protein levels of glycosaminoglycans (GAGs), and both the mRNA and protein levels of matrix metalloproteinase (MMP)-13 were evaluated. Luciferase reporter assay, quantitative real-time polymerase chain reaction (qRT-PCR), and Western blot were performed to explore whether Forkhead Box C1 (FOCX1) was a target of miR-138-5p. Further, we co-transfected OA chondrocytes with miR-138-5p mimics and pcDNA3.1 (+)-FOXC1 and then stimulated with IL-1β to determine whether miR-138-5p-mediated IL-1β-induced cartilage matrix degradation resulted from targeting FOXC1.

Results

MiR-138-5p was significantly increased in OA cartilage and in chondrocytes in response to IL-1β-stimulation. Overexpression of miR-138-5p significantly increased the IL-1β-induced downregulation of COL2A1, ACAN, and GAGs, and increased the IL-1β-induced over expression of MMP-13.We found that FOXC1 is directly regulated by miR-138-5p. Additionally, co-transfection with miR-138-5p mimics and pcDNA3.1 (+)-FOXC1 resulted in higher levels of COL2A1, ACAN, and GAGs, but lower levels of MMP-13.

Conclusion

miR-138-5p promotes IL-1β-induced cartilage degradation in human chondrocytes, possibly by targeting FOXC1.

Cite this article: Y. Yuan, G. Q. Zhang, W. Chai,M. Ni, C. Xu, J. Y. Chen. Silencing of microRNA-138-5p promotes IL-1β-induced cartilage degradation in human chondrocytes by targeting FOXC1: miR-138 promotes cartilage degradation. Bone Joint Res 2016;5:523–530. DOI: 10.1302/2046-3758.510.BJR-2016-0074.R2.