Skip to main content

Identification of a potential gene target for osteoarthritis based on bioinformatics analyses



Osteoarthritis (OA) is the most common chronic joint disease worldwide. It is characterized by pain and limited mobility in the affected joints and may even cause disability. Effective clinical options for its prevention and treatment are still unavailable. This study aimed to identify differences in gene signatures between tissue samples from OA and normal knee joints and to explore potential gene targets for OA.


Five gene datasets, namely GSE55457, GSE55235, GSE12021, GSE10575, and GSE1919, were selected from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were identified using the R programming software. The functions of these DEGs were analyzed, and a protein–protein interaction (PPI) network was constructed. Subsequently, the most relevant biomarker genes were screened using a receiver operating characteristic (ROC) curve analysis. Finally, the expression of the protein encoded by the core gene PTHLH was evaluated in clinical samples.


Eleven upregulated and 9 downregulated DEGs were shared between the five gene expression datasets. Based on the PPI network and the ROC curves of upregulated genes, PTHLH was identified as the most relevant gene for OA and was selected for further validation. Immunohistochemistry confirmed significantly higher PTHLH expression in OA tissues than in normal tissues. Moreover, similar PTHLH levels were detected in the plasma and knee synovial fluid of OA patients.


The bioinformatics analysis and preliminary experimental verification performed in this study identified PTHLH as a potential target for the treatment of OA.


Osteoarthritis (OA) is the most prevalent chronic joint disease affecting adults [1,2,3]. It ultimately causes chronic pain, restricted joint mobility, and disability [4]. The most common risk factors for OA include age, sex, prior joint injury, obesity, genetic predisposition, and mechanical factors, including malalignment and an abnormal joint shape [5,6,7,8]. The pathogenesis of OA remains unclear, and no effective method is yet available to prevent its progression [9]. Therefore, considerable research attention has been directed toward the identification of effective biomarkers and potential therapeutic targets for OA.

An increased understanding of OA has revealed it to be a total joint disease, with characteristic features such as a narrowing of the joint space, articular cartilage degradation, osteophyte formation, synovial hyperplasia, and subchondral sclerosis [10]. However, the degeneration and loss of the cartilaginous extracellular matrix (ECM) and a reduction in the chondrocyte population remain the most important characteristic changes [11]. Previous research has shown that chondrocytes are the sole type of cells present in the cartilage, where they are responsible for both the synthesis and turnover of the ECM [12]. Therefore, an understanding of the genes relevant to maintaining chondrocyte health is very important for the prevention and treatment of OA.

High-throughput sequencing technology can provide information about the expression patterns of thousands of genes in the human genome simultaneously. In recent years, the rapid development of this technology has yielded a large amount of clinical sample data. These rich datasets have enabled researchers to better understand the molecular mechanism of OA, identify disease-related gene changes, and explore potential therapeutic targets [13, 14]. However, due to heterogeneity among samples or sequencing platforms, mRNA expression data may not be consistent with different gene maps. Consequently, sequencing dataset analyses that have focused on OA have not yielded reliable results. A comprehensive bioinformatics approach would likely address the above-described shortcomings of previous work and enable the identification of the key genes involved in OA.

In this study, we selected five mRNA datasets (GSE55457, GSE55235, GSE12021, GSE10575, and GSE1919) closely related to OA from the GEO database. First, we used comprehensive bioinformatics analyses to identify DEGs between OA and normal samples. Next, we subjected these DEGs to functional enrichment and pathway analyses and constructed a PPI network, which was screened for key genes. An OA-specific biological function and pathway analysis may further our understanding of the mechanism of OA development at the molecular level. This knowledge could further promote our understanding of OA pathogenesis, in turn improving diagnosis, prognosis, and drug target identification.

Materials and methods

Data sources

The following gene expression profile datasets were downloaded from the GEO database: GSE55457, GSE55235, GSE12021, GSE10575, and GSE1919. The five GEO gene datasets contained 37 OA samples and 36 control samples (the corresponding information is shown in Table 1); all rheumatoid arthritis (RA) samples were excluded. The profiles of GSE55457, GSE55235, and GSE12021 were based on the GPL96 [HG-U133A] Affymetrix Human Genome U133A Array. The profile of GSE10575 was based on the GPL570 [HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array. The profile of GSE1919 was based on the GPL91 [HG_U95A] Affymetrix Human Genome U95A Array.

Table 1 Sample statistics of the five GEO gene datasets

Identification of DEGs

The limma package for the R platform was used to screen for DEGs between OA and normal knee joints, including synovial tissue samples (OA, 35; normal, 34) and cell samples (healthy chondrocytes, 2; OA osteocytes, 2). A P value < 0.05 and |log2FC| > 1 were used as the screening criteria for both upregulated and downregulated DEGs. Venn diagrams of genes shared between and unique to the five datasets were drawn using the R package. We also standardized the log2FC values of the genes shared between the five datasets and used the R package pheatmap to produce a cluster analysis heatmap of these shared genes.

GO enrichment and KEGG pathway analyses

Next, Gene Ontology (GO, and Kyoto Encyclopedia of Genes and Genomes (KEGG, enrichment analyses were performed using the R package clusterProfiler to detect the functions and pathways affected by the DEGs. The R package ggplot2 was used to plot the results for which KEGG annotation entries existed in at least four datasets, and for which GO subclass annotation entries existed in at least three datasets according to the annotation statistics contained in each dataset.

PPI network construction

To construct the PPI network of shared DEGs, we used the Search Tool for the Retrieval of Interacting Genes (STRING), an online database. Using the Cytoscape3.6.0 software, the gene with the highest degree of connectivity to the surrounding genes was selected based on the PPI network.

Collection of clinical samples and culture of chondrocytes

Twelve samples of knee cartilage were obtained from human patients who had been diagnosed with primary knee OA (grades III–IV) and underwent total knee arthroplasty. Six normal knee cartilage specimens were collected from patients with osteosarcoma who underwent segmental resection and artificial tumor knee prosthesis replacement or trauma amputation surgery. The clinical sample information is shown in Table 2. The study was conducted in accordance with the Declaration of Helsinki (2000 revision), and the protocol was approved by the Ethics Committee of the Second Xiangya Hospital, Central South University in Changsha, China (No. 2019106). The primary culture of human chondrocytes was the same as reported previously [15]. Cultured chondrocytes from no later than the first passage was used for the experiments.

Table 2 Patient information of knee cartilage and chondrocytes samples

Immunohistochemistry (IHC) and immunocytochemistry (ICC)

Six OA samples were randomly selected from 12 OA cartilage samples and verified with 6 normal cartilage samples. Cartilage tissues were fixed in 4% paraformaldehyde for 48 h and then decalcified in 10% ethylenediaminetetraacetic acid for approximately 4 weeks. Each tissue was embedded in paraffin and sectioned at a thickness of 6 μm. The experimental process was performed as previously described [16]. Additionally, 6 OA chondrocytes samples were randomly selected from 12 OA chondrocytes samples and verified with 6 normal chondrocytes samples. Human primary chondrocytes were seeded into chamber slides and then fixed with 4% paraformaldehyde for 10 min, permeabilized with 0.1% Triton X-100 in phosphate-buffered saline (PBS) for 15 min, and treated with 1% bovine serum albumin–PBS for 30 min. All fixed tissues and cells were incubated with a rabbit anti-PTHLH antibody (Santa Cruz Biotechnology, USA) for 2 h, followed by a peroxidase-conjugated anti-goat IgG secondary antibody for 30 min. The nuclei were counterstained with hematoxylin. All procedures were performed at room temperature.

Enzyme-linked immunosorbent assay (ELISA)

The blood samples of 20 normal healthy people came from the Health Management Center of the Second Xiangya Hospital, and their physical examination results showed that they were healthy. The blood and synovial fluid samples of 23 OA patients came from the Department of Orthopedics of the Second Xiangya Hospital. Their physical examination results showed that there was no parathyroid-related disease. Blood and synovial fluid samples were collected, centrifuged, and stored immediately at − 80 °C until analysis. The clinical sample information is shown in Table 3. The PTHLH concentrations in the plasma and synovial fluid were determined using a commercial ELISA (USCN Life Science & Technology Co Ltd., Wuhan, China) according to the manufacturer’s instructions.

Table 3 Patient information of blood samples and synovial fluid samples

Statistical analysis

R version 3.5.3 was used for the statistical analyses of the GEO datasets. In addition, the GraphPad Prism Software, version 7.0 was used for the statistical analysis of experimental data. The results are expressed as means ± standard deviations (SDs). A receiver operating characteristic (ROC) curve analysis was performed to identify possibly optimal models and discard suboptimal ones independently of (and prior to specifying) the cost context or the class distribution. Unpaired t test was used to analyze the statistical differences in age and BMI between OA group and control group. In addition, we used D’Agostino-Pearson, Shapiro-Wilk, and Kolmogorov–Smirnov normality test to test the normal distribution of data. Statistical analysis of immunohistochemical data showed that the data are normal distribution, although the sample size is relatively small. However, the statistical analysis of ELISA data showed that the data are not normal distribution. Therefore, we choose Mann-Whitney U test for statistical analysis. P < 0.05 was considered to indicate a significant difference between the groups.


Identification of 20 DEGs shared between five GEO profiles

We first analyzed the DEGs in the five GEO gene datasets and identified 1334, 1521, 1384, 3393, and 732 DEGs, including 309, 970, 304, 2253, and 403 upregulated genes and 1025, 551, 1080, 1140, and 329 downregulated genes, in GSE55457, GSE55235, GSE12021, GSE10575, and GSE1919, respectively. Twenty DEGs, including 11 upregulated and 9 downregulated genes, were shared between all five datasets, as shown by Venn diagrams in Fig. 1 a and b. In addition, we performed a cluster analysis of these 20 common DEGs in each dataset (Fig. 1c–g). As shown in the figure, downregulated and upregulated genes are represented in pink and blue, respectively. The OA and normal groups are depicted in purple and green, respectively, while red and blue coloring indicates high and low expression levels, respectively.

Fig. 1

DEGs screened from five GEO datasets. a Eleven upregulated DEGs and b nine downregulated were shared between the five GEO datasets. cg Cluster analysis heatmaps of the 20 shared DEGs in each dataset

GO functional enrichment analysis

Next, a GO functional enrichment analysis was performed using the R package clusterProfiler to identify the functions of the DEGs, which were classified into the groups of biological processes (BP), cellular component (CC), and molecular function (MF) (Fig. 2). Notably, the common genes were significantly associated with several GO terms. Particularly, the most strongly enriched BP terms included “leukocyte migration,” “response to steroid hormone,” “ossification,” and “extracellular structure organization.” The most strongly enriched MF terms were “extracellular matrix structural constituent” and “glycosaminoglycan binding.” The most common CC terms were “cell leading edge,” “collagen-containing extracellular matrix,” and “extracellular matrix.” These functions are known to be closely related to OA.

Fig. 2

GO enrichment analyses of DEGs. The image indicates GO subclass annotations existing in at least three datasets (dark blue columns of bars indicate). The shared genes classified into the molecular function (MF), biological processes (BP), and cellular component (CC) categories are summarized

KEGG pathway enrichment analysis

A KEGG pathway enrichment analysis revealed that the shared genes were mainly enriched in OA-related signaling pathways, including the PI3K-AKT, MAPK, and NF-kappa B signaling pathways, as well as in osteoclast differentiation [17,18,19] (Fig. 3).

Fig. 3

Results of the KEGG pathway enrichment analyses of DEGs. KEGG annotation entries existed in at least four datasets (dark blue columns of bars indicate) according to the annotation statistics

PPI network analysis

To reveal the interactive relationships between proteins encoded by the DEGs, we constructed a PPI network of the 20 shared DEGs using STRING and Cytoscape software (Fig. 4a). Here, the top 10 hub genes were HGF, PTHLH, CDH2, KLF4, WIF1, STC1, SLC7A5, IGLL5, STMN2, and GAP43. Four of the top five genes were found to be upregulated (HGF, PTHLH, CDH2, WIF1) and were thus selected for further analysis. Subsequently, we identified the upregulated genes (HGF, PTHLH, and CDH2) with higher scores in the sub-network analysis and selected the best potential target gene using a ROC curve. As shown in Fig. 4 b and c, in the GSE 55457 data set, HGF (AUC = 0.87) has the largest AUC area. However, there are many studies on the relationship between HGF and osteoarthritis, and the results are controversial or even contradictory [20,21,22]. Therefore, we finally decided to choose PTHLH (AUC = 0.91) in the GSE55235. These findings suggest that PTHLH is the most relevant biomarker for OA.

Fig. 4

Construction of a PPI network using the DEGs. a Schematic of the PPI network constructed from the DEGs. The scores are indicated by colors, with the top 5 DEGs indicated in red. The ROC-AUCs for b GSE55457 and c GSE 55235

IHC and ICC validation of the gain of PTHLH expression

To confirm the high expression of PTHLH in relation with OA, we subjected samples of normal and OA cartilage to IHC. Notably, the PTHLH expression was significantly higher in OA cartilage than in normal cartilage (P < 0.001; Fig. 5a). We also assessed PTHLH expression in normal and OA primary chondrocytes and found higher expression in the latter (P < 0.01; Fig. 5b).

Fig. 5

Increased PTHLH expression in b OA cartilage and b OA chondrocytes. The numbers of PTHLH (+) cells in cartilage tissues and primary chondrocytes subjected to IHC are presented as means ± SDs in the lower panel. Student’s t test. a ***P < 0.001, n = 6 (scale bars = 100 μm; magnification: × 10). b **P < 0.01, n = 6 (scale bars = 100 μm; magnification: × 20)

PTHLH concentrations in the plasma and synovial fluid

The healthy human knee cavity contains only a small amount of synovial fluid, which is almost impossible to obtain. In contrast, OA patients exhibit a significantly increased amount of synovial fluid. Thus, we used ELISA to evaluate the PTHLH levels in the synovial fluid and plasma of OA patients. Notably, no significant differences were observed in the PTHLH levels between these sample types (P > 0.05; Fig. 6a). We also measured the PTHLH levels in the plasma samples of healthy people and found significantly higher levels in them compared with those in the plasma samples of OA patients (P < 0.001; Fig. 6b).

Fig. 6

PTHLH levels in the plasma and synovial fluid. a PTHLH levels in the plasma and synovial fluid samples of OA patients, Mann-Whitney U test, P > 0.05. b PTHLH levels in the plasma samples of healthy and OA patients, Mann-Whitney U test, ***P < 0.001


OA is one of the leading causes of disability worldwide and thus places heavy burdens on the healthcare systems and global economy [23, 24]. Although the early diagnosis and treatment of OA are important, the mechanisms underlying OA pathogenesis remain unclear. It is essential to understand OA pathogenesis as this knowledge will facilitate diagnosis, prognosis, and drug target discovery [25,26,27]. High-throughput sequencing and microarray technology have been widely used to predict potential diagnostic and therapeutic targets for OA. In this study, we identified 20 genes that were differentially expressed between OA and normal samples and common to all five selected GEO mRNA expression datasets. Through GO and KEGG enrichment analyses, we determined that the processes affected by these genes were consistent with the findings of previous studies in which these pathways were identified as contributors to the pathological process of OA [28,29,30].

We also constructed PPI networks containing these DEGs. Among the top 10 hub genes, namely HGF, PTHLH, CDH2, KLF4, WIF1, STC1, SLC7A5, IGLL5, STMN2, and GAP43, previous studies identified HGF as a multifunctional contributor to OA. On the one hand, HGF can promote synovial and vascular proliferation, leading to inflammation and osteophyte formation [20, 22, 31]; on the other hand, it supports bone and cartilage repair [21]. Other reports have described CDH2 polymorphisms and KLF4 as potentially associated with the development of OA [32, 33]. PTHLH has been identified primarily as a factor that promotes chondrocyte proliferation and inhibits chondrocyte hypertrophy and terminal differentiation [34, 35]. Few reports have described correlations of the remaining hub genes with OA.

Our bioinformatics analysis results and existing literature evidence led us to speculate that PTHLH is a diagnostic and predictive biomarker of OA. Our IHC analysis to validate the aberrant gain-of-expression of PTHLH in OA cartilage tissue confirmed the significantly increased PTHLH levels in OA cartilage tissues and chondrocytes relative to those in normal control samples. However, an ELISA-based experiment revealed a significantly higher PTHLH level in the plasma of healthy people than in the plasma of OA patients. We evaluated why the PTHLH level is higher in OA knee cartilage or chondrocytes than in their normal counterparts and why the plasma PTHLH level is lower in OA patients than in normal people and identified the following possible reasons. First, under normal circumstances, PTHLH (the protein encoded by PTHLH, also known as PTHrP) plays predominantly paracrine and/or autocrine roles [36]. PTHLH exerts paracrine actions that affect physiological homeostasis in many tissues and processes, including hair follicles, cartilage, vascular smooth muscle, bone, pancreas, mammary gland development, and tooth eruption [35]. The secretion of PTHLH in some tissues and organs, such as the vascular smooth muscle, was shown to decrease with aging [37], and most patients with OA tend to be older [5]. Second, chondrocytes and synovial tissues can produce PTLHL directly [38, 39]. Therefore, in the pathological state of OA, chondrocytes and synovial tissues may secrete more PTHLH to maintain their own cell activity (Fig. 7).

Fig. 7

Flowchart of the study

Our study had some limitations. Firstly, there were statistical differences in age and BMI between the OA group and the control group (Table 2 and Table 3). At present, there are few studies on BMI and PTHLH levels. We are not sure whether BMI will affect PTHLH level. We emphasize that the biggest weakness of the study is that we do not know if OA itself causes the differences or the differences of age and BMI in the groups. Therefore, we designed a study that can obviate this issue. In future experiments, we will collect clinical samples in groups according to different age stages and different BMI intervals, and expand the sample size as much as possible. Secondly, although our analysis included five public datasets, the sample size was relatively small. In addition, the sample size used to verify the PTHLH expression data was also relatively small. Further experimental studies with larger sample sizes are needed to confirm our results.


In summary, we used bioinformatics analysis to identify 20 DEGs between normal and OA cartilage tissues and identified PTHLH as the main upregulated gene in the resulting PPI network. Our findings indicate that PTHLH may serve as a potential therapeutic target.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.





Differentially expressed genes


Protein–protein interaction


Receiver operating characteristic


Extracellular matrix


Gene Ontology


Kyoto Encyclopedia of Genes and Genomes


Search Tool for the Retrieval of Interacting Genes


Biological processes


Cellular component


Molecular function


  1. 1.

    Martel-Pelletier J, Barr AJ, Cicuttini FM, Conaghan PG, Cooper C, Goldring MB, et al. Osteoarthritis. Nat Rev Dis Primers. 2016;2:16072.

    Article  Google Scholar 

  2. 2.

    Snoeker B, Turkiewicz A, Magnusson K, Frobell R, Yu D, Peat G, Englund M: Risk of knee osteoarthritis after different types of knee injuries in young adults: a population-based cohort study. Br J Sports Med 2019.

  3. 3.

    Michael JW, Schluter-Brust KU, Eysel P. The epidemiology, etiology, diagnosis, and treatment of osteoarthritis of the knee. Dtsch Arztebl Int. 2010;107:152–62.

    PubMed  PubMed Central  Google Scholar 

  4. 4.

    O'Neill TW, Felson DT. Mechanisms of osteoarthritis (OA) pain. Curr Osteoporos Rep. 2018;16:611–6.

    Article  Google Scholar 

  5. 5.

    Palazzo C, Nguyen C, Lefevre-Colau MM, Rannou F, Poiraudeau S. Risk factors and burden of osteoarthritis. Ann Phys Rehabil Med. 2016;59:134–8.

    Article  Google Scholar 

  6. 6.

    Zengini E, Finan C, Wilkinson JM. The genetic epidemiological landscape of hip and knee osteoarthritis: where are we now and where are we going? J Rheumatol. 2016;43:260–6.

    CAS  Article  Google Scholar 

  7. 7.

    O'Neill TW, McCabe PS, McBeth J. Update on the epidemiology, risk factors and disease outcomes of osteoarthritis. Best Pract Res Clin Rheumatol. 2018;32:312–26.

    Article  Google Scholar 

  8. 8.

    Li YS, Xiao WF, Luo W. Cellular aging towards osteoarthritis. Mech Ageing Dev. 2017;162:80–4.

    Article  Google Scholar 

  9. 9.

    Runhaar J, Zhang Y. Can we prevent OA? Epidemiology and public health insights and implications. Rheumatology (Oxford). 2018;57:iv3–9.

    Article  Google Scholar 

  10. 10.

    Chen D, Shen J, Zhao W, Wang T, Han L, Hamilton JL, et al. Osteoarthritis: toward a comprehensive understanding of pathological mechanism. Bone Res. 2017;5:16044.

    CAS  Article  Google Scholar 

  11. 11.

    Gao Y, Liu S, Huang J, Guo W, Chen J, Zhang L, et al. The ECM-cell interaction of cartilage extracellular matrix on chondrocytes. Biomed Res Int. 2014;2014:648459.

    PubMed  PubMed Central  Google Scholar 

  12. 12.

    Loeser RF. Aging and osteoarthritis: the role of chondrocyte senescence and aging changes in the cartilage matrix. Osteoarthr Cartil. 2009;17:971–9.

    CAS  Article  Google Scholar 

  13. 13.

    Liu CC, Lee HC, Peng YS, Tseng AH, Wu JL, Tsai WY, Wong CS, Su LJ: Transcriptome analysis reveals novel genes associated with cartilage degeneration in posttraumatic osteoarthritis progression. Cartilage 2019:1947603519847744.

  14. 14.

    Park R, Ji JD. Unique gene expression profile in osteoarthritis synovium compared with cartilage: analysis of publicly accessible microarray datasets. Rheumatol Int. 2016;36:819–27.

    CAS  Article  Google Scholar 

  15. 15.

    Duan ZX, Tu C, Liu Q, Li SQ, Li YH, Xie P, et al. Adiponectin receptor agonist AdipoRon attenuates calcification of osteoarthritis chondrocytes by promoting autophagy. J Cell Biochem. 2020;121:3333–44.

    CAS  Article  Google Scholar 

  16. 16.

    Li Y, Xiao W, Sun M, Deng Z, Zeng C, Li H, et al. The expression of osteopontin and Wnt5a in articular cartilage of patients with knee osteoarthritis and its correlation with disease severity. Biomed Res Int. 2016;2016:9561058.

    PubMed  PubMed Central  Google Scholar 

  17. 17.

    Xie L, Xie H, Chen C, Tao Z, Zhang C, Cai L. Inhibiting the PI3K/AKT/NF-kappaB signal pathway with nobiletin for attenuating the development of osteoarthritis: in vitro and in vivo studies. Food Funct. 2019;10:2161–75.

    CAS  Article  Google Scholar 

  18. 18.

    Wu Z, Lu H, Yao J, Zhang X, Huang Y, Ma S, et al. GABARAP promotes bone marrow mesenchymal stem cells-based the osteoarthritis cartilage regeneration through the inhibition of PI3K/AKT/mTOR signaling pathway. J Cell Physiol. 2019;234:21014–26.

    CAS  Article  Google Scholar 

  19. 19.

    Wongwichai T, Teeyakasem P, Pruksakorn D, Kongtawelert P, Pothacharoen P. Anthocyanins and metabolites from purple rice inhibit IL-1beta-induced matrix metalloproteinases expression in human articular chondrocytes through the NF-kappaB and ERK/MAPK pathway. Biomed Pharmacother. 2019;112:108610.

    CAS  Article  Google Scholar 

  20. 20.

    Dankbar B, Neugebauer K, Wunrau C, Tibesku CO, Skwara A, Pap T, et al. Hepatocyte growth factor induction of macrophage chemoattractant protein-1 and osteophyte-inducing factors in osteoarthritis. J Orthop Res. 2007;25:569–77.

    CAS  Article  Google Scholar 

  21. 21.

    Adamopoulos IE, Athanasou NA. Hepatocyte growth factor in normal and diseased bone and joint tissues. Curr Rheumatol Rev. 2006;2:1–7.

    CAS  PubMed  Google Scholar 

  22. 22.

    Nagashima M, Hasegawa J, Kato K, Yamazaki J, Nishigai K, Ishiwata T, et al. Hepatocyte growth factor (HGF), HGF activator, and c-met in synovial tissues in rheumatoid arthritis and osteoarthritis. J Rheumatol. 2001;28:1772–8.

    CAS  PubMed  Google Scholar 

  23. 23.

    Cross M, Smith E, Hoy D, Nolte S, Ackerman I, Fransen M, et al. The global burden of hip and knee osteoarthritis: estimates from the global burden of disease 2010 study. Ann Rheum Dis. 2014;73:1323–30.

    Article  Google Scholar 

  24. 24.

    Fransen M, Bridgett L, March L, Hoy D, Penserga E, Brooks P. The epidemiology of osteoarthritis in Asia. Int J Rheum Dis. 2011;14:113–21.

    Article  Google Scholar 

  25. 25.

    Deng Z, Jia Y, Liu H, He M, Yang Y, Xiao W, et al. RhoA/ROCK pathway: implication in osteoarthritis and therapeutic targets. Am J Transl Res. 2019;11:5324–31.

    CAS  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Li H, George DM, Jaarsma RL, Mao X. Metabolic syndrome and components exacerbate osteoarthritis symptoms of pain, depression and reduced knee function. Ann Transl Med. 2016;4:133.

    CAS  Article  Google Scholar 

  27. 27.

    Tu C, He J, Wu B, Wang W, Li Z. An extensive review regarding the adipokines in the pathogenesis and progression of osteoarthritis. Cytokine. 2019;113:1–12.

    CAS  Article  Google Scholar 

  28. 28.

    Wang W, Li J, Li F, Peng J, Xu M, Shangguan Y, et al. Scutellarin suppresses cartilage destruction in osteoarthritis mouse model by inhibiting the NF-kappaB and PI3K/AKT signaling pathways. Int Immunopharmacol. 2019;77:105928.

    Article  Google Scholar 

  29. 29.

    Liu F, Li L, Lu W, Ding Z, Huang W, Li YT, et al. Zhanghui, yin Z: Scutellarin ameliorates cartilage degeneration in osteoarthritis by inhibiting the Wnt/beta-catenin and MAPK signaling pathways. Int Immunopharmacol. 2019;105954.

  30. 30.

    Hu SL, Wang K, Shi YF, Shao ZX, Zhang CX, Sheng KW, et al. Downregulating Akt/NF-kappaB signaling and its antioxidant activity with Loureirin a for alleviating the progression of osteoarthritis: in vitro and vivo studies. Int Immunopharmacol. 2019;105953.

  31. 31.

    Lin YM, Huang YL, Fong YC, Tsai CH, Chou MC, Tang CH. Hepatocyte growth factor increases vascular endothelial growth factor-a production in human synovial fibroblasts through c-met receptor pathway. PLoS One. 2012;7:e50924.

    CAS  Article  Google Scholar 

  32. 32.

    Shang H, Hao Y, Hu W, Hu X, Jin Q. CDH2 gene rs11564299 polymorphism is a risk factor for knee osteoarthritis in a Chinese population: a case-control study. J Orthop Surg Res. 2019;14:208.

    Article  Google Scholar 

  33. 33.

    Fan L, Li M, Cao FY, Zeng ZW, Li XB, Ma C, et al. Astragalus polysaccharide ameliorates lipopolysaccharide-induced cell injury in ATDC5 cells via miR-92a/KLF4 mediation. Biomed Pharmacother. 2019;118:109180.

    CAS  Article  Google Scholar 

  34. 34.

    Reynard LN, Loughlin J. Insights from human genetic studies into the pathways involved in osteoarthritis. Nat Rev Rheumatol. 2013;9:573–83.

    CAS  Article  Google Scholar 

  35. 35.

    Martin TJ. Parathyroid hormone-related protein, its regulation of cartilage and bone development, and role in treating bone diseases. Physiol Rev. 2016;96:831–71.

    CAS  Article  Google Scholar 

  36. 36.

    Wysolmerski JJ. Parathyroid hormone-related protein: an update. J Clin Endocrinol Metab. 2012;97:2947–56.

    CAS  Article  Google Scholar 

  37. 37.

    Jono S, Nishizawa Y, Shioi A, Morii H. Parathyroid hormone-related peptide as a local regulator of vascular calcification. Its inhibitory action on in vitro calcification by bovine vascular smooth muscle cells. Arterioscler Thromb Vasc Biol. 1997;17:1135–42.

    CAS  Article  Google Scholar 

  38. 38.

    Kohno H, Shigeno C, Kasai R, Akiyama H, Iida H, Tsuboyama T, et al. Synovial fluids from patients with osteoarthritis and rheumatoid arthritis contain high levels of parathyroid hormone-related peptide. J Bone Miner Res. 1997;12:847–54.

    CAS  Article  Google Scholar 

  39. 39.

    Pelosi M, Lazzarano S, Thoms BL, Murphy CL. Parathyroid hormone-related protein is induced by hypoxia and promotes expression of the differentiated phenotype of human articular chondrocytes. Clin Sci (Lond). 2013;125:461–70.

    CAS  Article  Google Scholar 

Download references


We sincerely thank the researchers for providing their GEO databases information online; it is our pleasure to acknowledge their contributions.


This work was supported by the Fundamental Research Funds for the Central Universities of Central South University (NO. 2017zzts231), the National Natural Science Foundation of China (No. 81372180, 81902745), the Hunan Provincial Science and Technology Association key-point program (No.2017DK2013), and Natural Science Foundation of Hunan Province, China (2018JJ3716).

Author information




Zhi-xi Duan and Li Zhi-hong participated in the design of the study; Zhi-xi Duan and Yu-sheng Li performed the experiments; Chao Tu, Lin Qi, Peng Xie, and Yi-han Li were responsible for data collection and analysis; Zhi-xi Duan and Yu-sheng Li wrote the manuscript. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Zhi-hong Li.

Ethics declarations

Ethics approval and consent to participate

This experimental study was approved by the Ethics Committee of the Second Xiangya Hospital, Central South University in Changsha, China.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Duan, Z., Li, Y., Tu, C. et al. Identification of a potential gene target for osteoarthritis based on bioinformatics analyses. J Orthop Surg Res 15, 228 (2020).

Download citation


  • Osteoarthritis
  • GEO database
  • Differentially expressed genes
  • Protein–protein interaction network