Gene Expression Profiles of Peripheral Blood Monocytes in Osteoarthritis and Analysis of Differentially Expressed Genes

Shi, Ting; Shen, Xiongjie; Gao, Ge

doi:https://doi.org/10.1155/2019/4291689

BioMed Research International

On this page

Abstract Background Results Discussion Conclusions Methods Abbreviations Data Availability Ethical Approval Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2019 | Article ID 4291689 | https://doi.org/10.1155/2019/4291689

Gene Expression Profiles of Peripheral Blood Monocytes in Osteoarthritis and Analysis of Differentially Expressed Genes

Ting Shi,¹Xiongjie Shen,²and Ge Gao³

Academic Editor: Graziano Pesole

Received17 Apr 2019

Revised27 Jul 2019

Accepted22 Oct 2019

Published26 Nov 2019

Abstract

Background. There is little understanding of the molecular processes involved in the pathogenesis of osteoarthritis, limiting early diagnosis and effective treatment of OA. Use of genechips can provide insights into the molecular pathogenesis of diseases. In this study, determination of gene expression profiles of osteoarthritis peripheral blood mononuclear cells will allow exploration of the molecular pathogenesis of OA and find out more candidate biomarkers and potential drug targets of OA. Result. A total of 1231 DEGs were screened out including 791 upregulated DEGs and 440 downregulated DEGs. The most significant upregulated DEG was RPL38, which may inhibit chondrocyte differentiation and synthesis of the extracellular matrix. PIK3CA, PIK3CB, PIK3CD, PIK3R1, MAPK14, IL1A, JUND, FOSL2, and PPP3CA were the gene symbols of the osteoclast differentiation pathway which was the most significant pathway enriched by DEGs. However, the MAPK signaling pathway occupied the core position of all the pathways which can regulate apoptosis, cell cycle, wnt signaling pathway, p53 signaling pathway, and phosphatidylinositol signaling system. Furthermore, PI3Ks may regulate IL1A, JUND, FOSL2 and PPP3CA through the MAPK signaling pathway. Conclusion. These identified DEGs and pathways may be novel biomarkers to monitor the changes of OA and can be a potential drug target for the treatment of OA.

1. Background

Osteoarthritis (OA) is a chronic degenerative joint disease characterized by degeneration of articular cartilage, synovium inflammation, imbalance in the synthesis and catabolism of the extracellular matrix of chondrocytes, and the formation of subchondral bone and osteophytes [1]. OA is common in the elderly, especially people older than 65 [2]. It is predominant in heavily loaded joints including the knee, hip, spine, and finger joints and ultimately leads to joint dysfunction [3]. Although there are many various therapies to relieve joint pain and improve joint function, the efficacy of these treatments is limited [4]. Joint replacement surgery can only treat patients at the end stage of OA, and X-ray diagnosis is not informative without visible radiographic changes in joint tissue. There is little understanding of the molecular processes involved in the pathogenesis of OA, limiting early diagnosis and effective treatment of OA. Thus, the identification of sensitive biomarkers and the development of novel drug targets at molecular level are key goals of OA research.

Together with the Human Genome Project and the rapid development of molecular biology technology, high-throughput genechip technology has emerged, allowing the rapid and simultaneous analysis of thousands of gene loci [5]. Use of genechips can provide insights into the molecular pathogenesis of diseases. Currently, gene expression profiles of OA have mainly focused on articular cartilage, subchondral bone, and synovium [6–8], but there has not been comprehensive microarray analysis of blood monocytes. Blood is more accessible than tissue, and blood sample collection is less painful for patients, so the identification of sensitive diagnostic biomarkers of OA in the peripheral blood would be highly valuable for clinical application. Peripheral blood mononuclear cells (PBMCs) participate in the occurrence and development of osteoarthritis by promoting osteoclastogenesis and bone resorption, and inhibiting osteoclast apoptosis and interleukin 1 receptor I (IL-1RI) expression [9]. Monocytes increase the degradation of the extracellular matrix of chondrocytes by promoting the expression of matrix metalloproteinase 13 (MMP13), an enzyme that participates in the degradation of extracellular matrix proteins [10]. Monocytes promote the apoptosis of chondrocytes and ultimately lead to cartilage degeneration [11]. Therefore, determination of gene expression profiles of osteoarthritis PBMCs will allow exploration of the molecular pathogenesis of osteoarthritis and may help identify improved targets for the treatment of osteoarthritis.

In this study, gene expression profiles of osteoarthritis PBMCs were constructed by genechip technology. Differentially expressed genes (DEGs) were screened out by comparing the genechip results of osteoarthritis patients with those of normal controls. To obtain greater insights into the molecular mechanisms of OA, we applied bioinformatics analysis. Gene ontology (GO) analysis and pathway enrichment analysis were performed for DEGs on the Gene-Cloud of Biotechnology Information (GCBI) bioinformatics platform, revealing the core genes and signaling pathways in the pathogenesis of OA. In addition, the network relationships between DEGs and signaling pathways were determined by pathway relation and gene signal network analyses, revealing key players in the molecular pathogenesis of OA.

2. Results

2.1. Identification of Differentially Expressed Genes

Gene expression profiles of peripheral blood monocytes for OA groups and control groups were compared, revealing 1231 DEGs. Of these genes, 791 were upregulated and 440 were downregulated. We ranked the differentially expressed genes according to the value. The top thirty up- and downregulated DEGs are listed in Table 1. The lowest value of the upregulated DEGs was ribosomal protein L38 (RPL38; = 5.30 × 10⁻⁰⁵), followed by protein phosphatase 3, catalytic subunit, alpha isozyme (PPP3CA; = 7.40 × 10⁻⁰⁵). The most significant downregulated DEG was IKAROS family zinc finger 1 (IKZF1; = 1.36E − 04), followed by chromosome 7 open reading frame 43 (C7orf43; = 1.87E − 04). Cluster analysis was performed on the differentially expressed genes, and the result is represented as a heatmap (Figure 1).

2.2. GO Analysis and Pathway Enrichment Analysis of DEGs

We next performed GO analysis of the identified DEGs and found 449 biological processes that were significantly enriched. The 30 most enriched GO terms and related DEGs are shown in Table 2. We found that the DEGs were mainly enriched in small molecule metabolic process (GO:0044281, GO:0006468, and GO:0044267), immune process (GO:0045087, GO:0006955, GO:0006954, GO:0019221, and GO:0050900), cell proliferation and apoptosis (GO:0006915, GO:0008283, GO:0043065, and GO:0008219), cell cycle (GO:0000278, GO:0000090, GO:0000087, and GO:0007049), RNA splicing process (GO:0008380 and GO:0000398), and cell adhesion (GO:0030155). DEGs enriched in these biological processes included interleukin 1, alpha (IL1A), interleukin 6 receptor (IL6R), interleukin 7 receptor (IL7R), protein phosphatase 3, catalytic subunit, alpha isozyme (PPP3CA), ribosomal protein L38 (RPL38), phosphatidylinositol-4,5-bisphosphate 3-kinase (PIK3), mitogen-activated protein kinase 14 (MAPK14), tumor necrosis factor superfamily, member 10 (TNFSF10), tumor necrosis factor superfamily, member 13 (TNFSF13), activating transcription factor 2 (ATF2), and others. The DEGs mapped to 143 pathways, and the pathway enrichment analysis is partly presented in Table 3. Many of the genes with altered expression participate in osteoclast differentiation, apoptosis, focal adhesion, and cell cycle. In addition, many DEGs are parts of signaling pathways like the MAPK, PI3K-Akt, calcium, T-cell receptor, and wnt and Jak-STAT signaling pathways. Tyrosine kinase 2 (TYK2), phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit alpha (PIK3CA), phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit beta (PIK3CB), phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit delta (PIK3CD), phosphoinositide-3-kinase, regulatory subunit 1 (PIK3R1), MAPK14, IL1A, jun D proto-oncogene (JUND), FOS-like antigen 2 (FOSL2), and PPP3CA genes are related to osteoclast differentiation, and IL-1A, JUND, PPP3CA, and MAPK14 are involved in the MAPK signaling pathway.

2.3. Network Analysis

The geneSignalNetwork analysis was performed on 1231 DEGs and revealed interactions between 183 DEGs (Figure 2). The results showed that phospholipase C, beta 1 (PLCB1), PIK3CA, PIK3CB, PIK3CD, PIK3R1, v-Ki-ras2 Kirsten rat sarcoma viral oncogene homolog (KRAS), MAPK14, protein kinase, cAMP-dependent, catalytic, beta (PRKACB), phosphatidylinositol-4-phosphate 5-kinase, type I, alpha (PIP5K1A), G protein, alpha inhibiting activity polypeptide 2 (GNAI2), G protein, alpha inhibiting activity polypeptide 3 (GNAI3), calmodulin 1 (CALM1), G protein, beta polypeptide 2 (GNB2), G protein, gamma 11 (GNG11), beta-actin (ACTB), and Rac family small GTPase2 (RAC2) interact with more than 10 DEGs, making them core genes in the network. The network relationship of DEGs related to osteoclast differentiation, including PIK3CA, PIK3CB, PIK3CD, PIK3R1, MAPK14, JUND, FOSL2, IL1A, and PPP3CA, is shown in Figure 3.

Figure 2

The map of DEG interaction relationship. Red nodes indicate upregulated genes, and blue nodes indicate downregulated genes. The larger the area of the node, the more important the gene in the network. The interaction relationship between genes is expressed by a connection line. The solid line without arrow indicates the interaction between DEGs is not directional, otherwise is directional. The beginning of arrow is the upstream gene, and the end of arrow is the downstream gene. The dashed line with arrow indicates the indirect interaction relationship. The solid line with flat head indicates the inhibitory interaction. The letters on the connection line are abbreviated as the type of interaction. c, compound; a, active; ind, indirect; inh, inhibition.

The pathwayRelationNetwork analysis was also performed and revealed interaction relationships between 68 pathways (Figure 4). The MAPK signaling pathway, apoptosis, pathways in cancer, cell cycle, p53 signaling pathway, calcium signaling pathway, wnt signaling pathway, adherens junction, ErbB signaling pathway, focal adhesion, regulation of actin cytoskeleton, ubiquitin-mediated proteolysis, TGF-beta signaling pathway, and pancreatic cancer interacted with more than 10 pathways, making these core pathways. The most significant pathway in this analysis was the MAPK signaling pathway and its related pathways, as shown in Figure 5. The calcium signaling pathway, cell cycle, wnt signaling pathway, TGF-beta signaling pathway, and VEGF signaling pathway were the source pathways of MAPK signaling pathway. Moreover, MAPK signaling pathway regulated apoptosis, cell cycle, wnt signaling pathway, p53 signaling pathway, and phosphatidylinositol signaling system.

3. Discussion

In the present study, the Human Genome U133A genechip from Affymetrix was used to detect the gene expression profiles of PBMCs from 19 patients with osteoarthritis and 22 healthy controls. We identified 1231 DEGs, with 791 upregulated DEGs and 440 downregulated DEGs. The top 10 most significant DEGs according to the p value were RPL38, PPP3CA, TROVE domain family, member 2 (TROVE2), PRP40 pre-mRNA processing factor 40 homolog A (PRPF40A), centrosomal protein 350 kDa (CEP350), IKAROS family zinc finger 1 (IKZF1), mitochondrial ribosomal protein S31 (MRPS31), RAB11A, member RAS oncogene family (RAB11A), zinc finger CCCH-type containing 14 (ZC3H14), and chromosome 7 open reading frame 43 (C7orf43). The most significant upregulated DEG was RPL38, which affects processes of gene expression, DNA transcription, gene translation, and cell proliferation. Normal expression of ribosomal protein L29 (RPL29) is essential for chondrocyte proliferation and skeletal growth, and high expression level of RPL29 inhibited chondrocyte terminal differentiation, keeping cells in a state of abnormal proliferation [12]. Altered expression of RPL29 also influenced the rate of extracellular matrix protein synthesis [13]. Green et al. found different expressions of ribosomal protein L10 (RPL10) during endochondral bone development, suggesting effects on cell differentiation before bone mineralization [14]. Given effects of other ribosomal proteins, we hypothesize that ribosomal protein RPL38 may participate in the pathogenesis of OA by inhibiting chondrocyte differentiation and proliferation, decreasing the synthesis of the extracellular matrix.

Osteoclasts are closely related to bone resorption and subchondral bone remodeling processes involved in the pathophysiology of OA [15]. Osteoclasts are large multinucleated cells which originate from hematopoietic precursors of the monocyte-macrophage lineage [16]. Durand et al. found PBMCs from OA patients display increased osteoclastogenesis and bone resorption [9]. Our pathway enrichment analysis of DEGs support previous findings and suggest that the osteoclast differentiation pathway contributes to OA development. PIK3CA, PIK3CB, PIK3CD, PIK3R1, MAPK14, IL1A, JUND, FOSL2, and PPP3CA were the gene symbols of osteoclast differentiation pathway, and these genes showed increased expression in PBMCs of OA patients compared to normal controls. This suggests that mononuclear cells from patients with OA have stronger osteoclast differentiation ability. PIK3CA, PIK3CB, PIK3CD, and PIK3R1 belong to the PI3Ks (phosphoinositide-3-kinases) family [17]. PI3K is not only involved in osteoclast differentiation, activation, and survival but also contributes to osteoclast-mediated bone resorption and in vivo bone homeostasis [18]. PI3K inhibitors such as wortmannin and LY294002 inhibit osteoclast chemotaxis, attachment, and spreading [19]. The PI3K isoform PIK3CB also promotes osteoclast development and bone resorption [20]. JUND and FOSL2 belong to the activator protein-1 (AP-1) transcriptional factor family, which regulates many cellular processes such as differentiation, proliferation, and apoptosis [21]. Fos is a necessary factor for the differentiation of hematopoietic precursor cells into osteoclasts. Mice deficient in Fos (Fos-/-) develop osteopetrosis and show high bone mineral density and defective osteoclast formation. Fos1 transgenic mice showed increased production of osteoclasts and active bone resorption [22]. The JUND/FOSL2 heterodimer upregulates the expression of Tcirg1 to increase osteoclastogenesis [23]. IL-1A, as a proinflammatory cytokine, can disrupt bone metastasis and stimulate osteoclast resorption activity [24]. PPP3CA, or calcineurin 1, is a serine/threonine-specific phosphatase regulated by Ca²⁺/calmodulin [25]. PPP3C has effects in tumors, cardiovascular system, immune system, and nervous system [26]. The dephosphorylation of PPP3C helps regulate cell proliferation and differentiation in the cytokine-mediated activation of immune cells [27]. Collectively, all these upregulated DEGs associated with osteoclast differentiation may promote the development of OA and have potential to serve as novel biomarkers to monitor changes of OA and are potential drug targets for the treatment of OA.

The most significant pathway revealed by network analysis was the MAPK signaling pathway, with an occupied core position in the network. The MAPK signaling pathway as the source pathway can regulate apoptosis, cell cycle, wnt signaling pathway, p53 signaling pathway, and the phosphatidylinositol signaling system. MAPKs are serine/threonine protein kinases, and activating a protein kinase cascade (MAPKKK-MAPKK-MAPK) allows the transduction of extracellular signals into the nucleus. In this way, the expression of genes related to cell proliferation, differentiation, apoptosis, and stress response can be regulated [28]. The MAPK signaling pathway plays a variety of roles in OA. First, p38-MAPK promotes chondrocyte apoptosis through Fas-mediated apoptotic pathway and activation of ATF-2, caspase-3, and p53 [29]. Second, p38-MAPK induces the expression of MMP13, promotes the degradation of collagen II in the extracellular matrix of chondrocytes, and inhibits the synthesis and secretion of collagen and glycoprotein [30]. Third, p38-MAPK upregulates expression of cyclooxygenase-2 (COX2), prostaglandin E2 (PGE2), and inducible nitric oxide synthase (iNOS), and increases the synthesis of inflammatory mediators, leading to joint swelling, pain, and cartilage destruction [31]. Fourth, p38-MAPK is essential for the differentiation of osteoclast precursor cells into osteoclasts as well as osteoclast maturation and survival [32]. In our analysis of DEGs enriched in the osteoclast differentiation pathway, we found that PIK3CA, PIK3CB, PIK3CD, and PIK3R1 indirectly activate MAPK14, the prototypic member of the p38-MAPK family, through phosphorylation. MAPK14 is a key upstream gene that activates the expression of downstream genes such as IL1A, JUND, FOSL2, and PPP3CA. The network relationships between these genes revealed that PI3Ks may regulate IL1A, JUND, FOSL2, and PPP3CA through the MAPK signaling pathway. However, these regulatory relationships require further experiment to verify.

The results of this study have similarities and differences with published microarray analysis on synovial fluid, synovium membranes, subchondral bone, and cartilage. Firstly, the biological processes enriched by some DEGs are the same, such as immune response and inflammatory response (IL, TNF, TLR, and TGF-β), cell cycle (TGF-β and SMAD2), cell apoptosis (CASP and BAX). Secondly, high expression of proteinase (ADAMTS10) is same. Thirdly, some signal pathways enriched by DEGs are the same (MAPK signaling pathway, TOLL-like signaling pathway, and TGF-β signaling pathway) [6–8]. However, this study found that the main signal pathway enriched by DEGs was osteoclast differentiation, which may be related to the fact that peripheral blood mononuclear cells are precursors of osteoclasts. Only a small portion of the DEGs have the same results as previous tissue microarray because gene expression of PBMCs is affected by several processes. Attur et al. [33] showed that PBMCs can be activated by the inflammatory process of osteoarthritis when they pass through different tissues of the joint. DEGs in PBMCs may reflect the susceptibility of OA from a genetic perspective.

4. Conclusions

Gene expression profiles were determined for PBMCs in OA patients and healthy controls. Analysis revealed 1231 DEGs, with 791 upregulated DEGs and 440 downregulated DEGs. The most significant upregulated DEG was RPL38, which may inhibit chondrocyte differentiation and synthesis of the extracellular matrix. PIK3CA, PIK3CB, PIK3CD, PIK3R1, MAPK14, IL1A, JUND, FOSL2, and PPP3CA were the gene symbols of the osteoclast differentiation pathway. These genes promote the development of OA through increasing osteoclastogenesis. PI3Ks may also regulate IL1A, JUND, FOSL2, and PPP3CA through the MAPK signaling pathway. The discovery of these molecular mechanisms associated with the progression of osteoarthritis may contribute to the early diagnosis and treatment of osteoarthritis. Further work should evaluate the utility of these DEGs as biomarkers with diagnostic efficiency, as well as their potential as novel drug targets for the treatment of OA.

5. Methods

5.1. Research Subjects

The research subjects consisted of 19 OA patients and 22 healthy controls. All the research subjects were collected from Department of Rheumatology and Immunology, the second Xiangya Hospital of Central South University. The symptomatic primary OA patients who met the inclusion criteria were diagnosed by American College of Rheumatology(ACR) standard diagnostic criteria [34], with the Kellgren and Lawrence scoring system rating over grade 1 [35]. The Kellgren and Lawrence scoring system is divided into five grades: grade 0 means normal radiograph; grade 1 means suspicious pathology; grade 2 means small osteophytes and possible joint space narrowing; grade 3 means moderate osteophytes, definite joint space narrowing, and some of subchondral bone sclerosis; and grade 4 means large osteophytes, severe joint space narrowing, and obvious subchondral bone sclerosis. Within the 19 OA patients analyzed, there were 3 OA patients of grade 1, 4 OA patients of grade 2, 6 OA patients of grade 3, and 6 OA patients of grade 4. Exclusion criteria included the diagnosis of RA, osteoporosis, ankylosing spondylitis, neoplastic disease, or any other inflammatory autoimmune diseases. The healthy controls were collected from the same period of physical examination without clinical and radiographic osteoarthritis according to ACR criteria. The age, height, weight, and sex of OA patients and healthy controls were comparable. The research process followed the ethical standards of human experiment in Xiangya Medical College of Central South University and obtained the consent of the research subjects.

5.2. Monocyte Extraction and RNA Isolation

We collected fresh peripheral blood from each subject and separated circulating mononuclear cells (MNCs) with UNI-SEP separation tube (Novamed INC, Illinois, USA) according to manufacturer’s recommendations. UNI-SEP products are sterile plastic centrifuge tubes containing a solution of 5.6% polysucrose and 9.6% sodium metrizoate. Then, we used monocyte negative isolation kit (Dynal Biotech ASA, Oslo, Norway) to isolate monocytes from MNCs. The antibody mix contains mouse IgG antibodies for CD2, CD7, CD16, CD19, CD56 and CD235a bound the surface of depletion dynabeads, and we discarded beads with unwanted cells and collected untouched monocytes. Then, the monocytes were lysed in TRIzol Reagent (Life Technologies, Bleiswijk, The Netherlands) for 5 min, chloroform was added, and then centrifuged (15 min at 12 000). Transferring the upper aqueous phase to another RNase-free EP tube and adding isopropanol, the tube was placed stably for 10 min and then centrifuged (15 min at 12 000). The sediments (RNA) were washed with 75% ethanol and dissolved in DEPC-treated water. Quality and quantity of the isolated RNA were measured by analyzing samples using the 2100 Bioanalyzer (Agilent Technologies, Amstelveen, The Netherlands) and Nanodrop spectrophotometer.

5.3. Microarray Analysis

Total RNA was reverse transcribed to synthesis first-strand cDNA by using a T7-Oligao (dT) Promoter Primer, and then second-strand cDNA was synthesized and purified. Next, in in vitro transcription reaction, second-strand cDNA serves as a template performed for complementary RNA (cRNA) amplification and biotin labeling and was cleaned up and fragmented. All the operations were following the instructions of MessageAmp™ Premier RNA Amplification kit (Invitrogen, USA). Subsequently, the biotinylated cRNAs were hybridized to Affymetrix Human Genome U133A genechip overnight at 45°C with rocking. Genechip was then washed and stained according to the protocol of GeneChip® Hybridization, Wash, and Stain kit (Thermo Fisher Scientific, USA). Finally, Affymetrix Gene Array Scanner G2500A scanned the chip, and the scanned images were analyzed by MAS5.0 software.

5.4. Screening of DEGs

The image signals of the chip were transformed into the original data of gene expression level by MAS5.0 software. OA group and control group were set up on the GCBI working platform (https://www.gcbi.com.cn/, GMINIX Informatics Ltd. Co, Shanghai, P.R. China). The original data CEL files were imported into the two groups. Using SAM (significance analysis of microarray) R program [36], DEGs between OA group and control group were screened according to the screening criteria of , , and fold change > 1.2(Q value is the false discovery rate, which represents the expected proportion of false-positive genes in differentially expressed genes. value is adjusted by the Q value. When Q < 0.05, it can assure lower false discovery rate and higher reliability of results). Hierarchical clustering (a heatmap) was performed with these screened DEGs to identify whether they can make an obvious distinction between OA samples and control samples.

5.5. GO Analysis and Pathway Enrichment Analysis of DEGs

We also performed GO (gene ontology) analysis (, FDR, false discovery rate < 0.05) and pathway enrichment analysis () for DEGs on the GCBI working platform. In GO analysis, biological processes, cell components, and molecular processes are mainly involved. In the screening of DEGs, more information is provided by biological processes. Therefore, in GO analysis on the GCBI platform, biological processes should be selected as the “analysis type” in parameter setting.

5.6. Network Analysis

For DEGs and its associated signaling pathways, we did geneSignalNetwork and pathwayRelationNetwork analysis on the GCBI platform, so the interaction relationship of DEGs and signaling pathways was obtained.

Abbreviations

OA:	Osteoarthritis
PBMCs:	Peripheral blood mononuclear cells
IL-1RI:	Interleukin 1 receptor type I
MMP13:	Matrix metallopeptidase 13
DEGs:	Differentially expressed genes
GO:	Gene ontology
GCBI:	Gene-Cloud of Biotechnology Information
ACR:	American College of Rheumatology
RA:	Rheumatoid arthritis
RPL38:	Ribosomal protein L38
PPP3CA:	Protein phosphatase 3, catalytic subunit, alpha isozyme
IKZF1:	IKAROS family zinc finger 1
C7orf43:	Chromosome 7 open reading frame 43
IL1A:	Interleukin 1, alpha
IL6R:	Interleukin 6 receptor
IL7R:	Interleukin 7 receptor
PIK3:	Phosphoinositide-3-kinase
MAPK14:	Mitogen-activated protein kinase 14
TNFSF10:	Tumor necrosis factor (ligand) superfamily, member 10
TNFSF13:	Tumor necrosis factor (ligand) superfamily, member 13
ATF2:	Activating transcription factor 2
TYK2:	Tyrosine kinase 2
PIK3CA:	Phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit alpha
PIK3CB:	Phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit beta
PIK3CD:	Phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit delta
PIK3R1:	Phosphoinositide-3-kinase, regulatory subunit 1 (alpha);
MAPK14:	Mitogen-activated protein kinase 14
JUND:	Jun D proto-oncogene
FOSL2:	FOS-like antigen 2.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Ethical Approval

This research process followed the ethical standards of human experiment in Xiangya Medical College of Central South University and obtained the consent of the research subjects.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

TS collected samples of research subjects and completed all microarray-related experiments. XJS performed the data analyses. GG and TS helped perform the analysis with constructive discussions. TS and XJS wrote the manuscript. All authors read and approved the final manuscript.

Acknowledgments

The authors would like to acknowledge the online analysis platform GCBI for providing us technical assistance. This work was supported by the Hunan Provincial Natural Science Fund (grant no. 2016JJ4053).

References

G. Musumeci, M. A. Szychlinska, and A. Mobasheri, “Age-related degeneration of articular cartilage in the pathogenesis of osteoarthritis: molecular markers of senescent chondrocytes,” Histology and Histopathology, vol. 30, no. 1, pp. 1–12, 2015.
View at: Google Scholar
B. Xia, D. Chen, J. Zhang, S. Hu, H. Jin, and P. Tong, “Osteoarthritis pathogenesis: a review of molecular mechanisms,” Calcified Tissue International, vol. 95, no. 6, pp. 495–505, 2014.
View at: Publisher Site | Google Scholar
M. J. Thomas, G. Peat, T. Rathod et al., “The epidemiology of symptomatic midfoot osteoarthritis in community-dwelling older adults: cross-sectional findings from the clinical assessment study of the foot,” Arthritis Research & Therapy, vol. 17, no. 178, pp. 1–11, 2015.
View at: Publisher Site | Google Scholar
F. Russo, M. D’Este, G. Vadala et al., “Platelet rich plasma and hyaluronic acid blend for the treatment of osteoarthritis: rheological and biological evaluation,” PLoS One, vol. 11, no. 6, Article ID e015704, 2016.
View at: Publisher Site | Google Scholar
S. Warner and A. Valdes, “The genetics of osteoarthritis: a review,” Journal of Functional Morphology and Kinesiology, vol. 1, no. 1, pp. 140–153, 2016.
View at: Publisher Site | Google Scholar
W. Wang, Y. Liu, J. Hao et al., “Comparative analysis of gene expression profiles of hip articular cartilage between non-traumatic necrosis and osteoarthritis,” Gene, vol. 591, no. 1, pp. 43–47, 2016.
View at: Publisher Site | Google Scholar
I. Prasadam, J. Batra, S. Perry, W. Gu, R. Crawford, and Y. Xiao, “Systematic identification, characterization and target gene analysis of microRNAs involved in osteoarthritis subchondral bone pathogenesis,” Calcified Tissue International, vol. 99, no. 1, pp. 43–55, 2016.
View at: Publisher Site | Google Scholar
R. Park and J. D. Ji, “Unique gene expression profile in osteoarthritis synovium compared with cartilage: analysis of publicly accessible microarray datasets,” Rheumatology International, vol. 36, no. 6, pp. 819–827, 2016.
View at: Publisher Site | Google Scholar
M. Durand, S. V. Komarova, A. Bhargava et al., “Monocytes from patients with osteoarthritis display increased osteoclastogenesis and bone resorption: the in vitro osteoclast differentiation in arthritis study,” Arthritis & Rheumatism, vol. 65, no. 1, pp. 148–158, 2013.
View at: Publisher Site | Google Scholar
P. Castrogiovanni, M. D. Rosa, C. Guglielmino et al., “Moderate physical activity as a prevention method for knee osteoarthritis and the role of synoviocytes as biological key,” International Journal of Molecular Sciences, vol. 20, no. 3, 2019.
View at: Publisher Site | Google Scholar
H. Raghu, C. M. Lepus, Q. Wang et al., “CCL2/CCR2, but not CCL5/CCR5, mediates monocyte recruitment, inflammation and cartilage destruction in osteoarthritis,” Annals of the Rheumatic Diseases, vol. 76, no. 5, pp. 914–922, 2017.
View at: Publisher Site | Google Scholar
S. A. Miller, A. J. Brown, M. C. Farach-Carson, and C. B. Kirn-Safran, “HIP/RPL29 down-regulation accompanies terminal chondrocyte differentiation,” Differentiation, vol. 71, no. 6, pp. 322–336, 2003.
View at: Publisher Site | Google Scholar
S. S. Jiang, C.-H. Chen, K.-Y. Tseng et al., “Gene expression profiling suggests a pathological role of human bone marrow-derived mesenchymal stem cells in aging-related skeletal diseases,” Aging, vol. 3, no. 7, pp. 672–684, 2011.
View at: Publisher Site | Google Scholar
H. Green, A. E. Canfield, M. C. Hillarby et al., “The Ribosomal protein QM is expressed differentially during vertebrate endochondral bone development,” Journal of Bone and Mineral Research, vol. 15, no. 6, pp. 1066–1075, 2000.
View at: Publisher Site | Google Scholar
J. Geurts, A. Patel, M. T. Hirschmann et al., “Elevated marrow inflammatory cells and osteoclasts in subchondral osteosclerosis in human knee osteoarthritis,” Journal of Orthopaedic Research, vol. 34, no. 2, pp. 262–269, 2016.
View at: Publisher Site | Google Scholar
D. S. Amarasekara, H. Yun, S. Kim, N. Lee, H. Kim, and J. Rho, “Regulation of osteoclast differentiation by cytokine networks,” Immune Network, vol. 18, no. 1, p. e8, 2018.
View at: Publisher Site | Google Scholar
D. A. Fruman, H. Chiu, B. D. Hopkins, S. Bagrodia, L. C. Cantley, and R. T. Abraham, “The PI3K pathway in human disease,” Cell, vol. 170, no. 4, pp. 605–635, 2017.
View at: Publisher Site | Google Scholar
W. Xuan, X. Feng, C. Qian et al., “Osteoclast differentiation gene expression profiling reveals chemokine CCL4 mediates RANKL-induced osteoclast migration and invasion via PI3K pathway,” Cell Biochemistry and Function, vol. 35, no. 3, pp. 171–177, 2017.
View at: Publisher Site | Google Scholar
S. E. Lee, K. M. Woo, S. Y. Kim et al., “The phosphatidylinositol 3-kinase, p38, and extracellular signal-regulated kinase pathways are involved in osteoclast differentiation,” Bone, vol. 30, no. 1, pp. 71–77, 2002.
View at: Publisher Site | Google Scholar
D. Győri, D. Csete, S. Benkő et al., “The phosphoinositide 3-kinase isoform PI3Kβ regulates osteoclast-mediated bone resorption in humans and mice,” Arthritis & Rheumatology, vol. 66, no. 8, pp. 2210–2221, 2014.
View at: Publisher Site | Google Scholar
A. Papoudou-Bai, E. Hatzimichael, A. Barbouti, and P. Kanavaros, “Expression patterns of the activator protein-1 (AP-1) family members in lymphoid neoplasms,” Clinical and Experimental Medicine, vol. 17, no. 3, pp. 291–304, 2017.
View at: Publisher Site | Google Scholar
K. Matsuo, J. M. Owens, M. Tonko, C. Elliott, T. J. Chambers, and E. F. Wagner, “Fosl1 is a transcriptional target of c-Fos during osteoclast differentiation,” Nature Genetics, vol. 24, no. 2, pp. 184–187, 2000.
View at: Publisher Site | Google Scholar
G. E. Beranger, D. Momier, J.-M. Guigonis, M. Samson, G. F. Carle, and J.-C. Scimeca, “Differential binding of poly(ADP-Ribose) polymerase-1 and JunD/Fra2 accounts for RANKL-induced Tcirg1 gene expression during osteoclastogenesis,” Journal of Bone and Mineral Research, vol. 22, no. 7, pp. 975–983, 2007.
View at: Publisher Site | Google Scholar
S. Wijekoon, E. C. Bwalya, J. Fang, S. Kim, K. Hosoya, and M. Okumura, “Chronological differential effects of pro-inflammatory cytokines on RANKL-induced osteoclast differentiation of canine bone marrow-derived macrophages,” Journal of Veterinary Medical Science, vol. 79, no. 12, pp. 2030–2035, 2017.
View at: Publisher Site | Google Scholar
T. Mizuguchi, M. Nakashima, M. Kato et al., “Loss-of-function and gain-of-function mutations in PPP3CA cause two distinct disorders,” Human Molecular Genetics, vol. 27, no. 8, pp. 1421–1433, 2018.
View at: Publisher Site | Google Scholar
P. N. Gabrovska, R. A. Smith, L. M. Haupt, and L. R. Griffiths, “Investigation of two Wnt signalling pathway single nucleotide polymorphisms in a breast cancer-affected Australian population,” Twin Research and Human Genetics, vol. 14, no. 6, pp. 562–567, 2011.
View at: Publisher Site | Google Scholar
C. Hundhausen, C. Boesch-Saadatmandi, N. Matzner et al., “Ochratoxin a lowers mRNA levels of genes encoding for key proteins of liver cell metabolism,” Cancer Genomics Proteomics, vol. 5, no. 6, pp. 319–332, 2008.
View at: Google Scholar
Y. Sun, W.-Z. Liu, T. Liu, X. Feng, N. Yang, and H.-F. Zhou, “Signaling pathway of MAPK/ERK in cell proliferation, differentiation, migration, senescence and apoptosis,” Journal of Receptors and Signal Transduction, vol. 35, no. 6, pp. 600–604, 2015.
View at: Publisher Site | Google Scholar
C. Zhang, L. Yu, Y. Zhou, Q. Zhao, and S.-Q. Liu, “Chitosan oligosaccharides inhibit IL-1β-induced chondrocyte apoptosis via the P38 MAPK signaling pathway,” Glycoconjugate Journal, vol. 33, no. 5, pp. 735–744, 2016.
View at: Publisher Site | Google Scholar
H. Lim and H. P. Kim, “Matrix metalloproteinase-13 expression in IL-1β-treated chondrocytes by activation of the p38 MAPK/c-Fos/AP-1 and JAK/STAT pathways,” Archives of Pharmacal Research, vol. 34, no. 1, pp. 109–117, 2011.
View at: Publisher Site | Google Scholar
H. Joos, W. Albrecht, S. Laufer, and R. E. Brenner, “Differential effects of p38MAP kinase inhibitors on the expression of inflammation-associated genes in primary, interleukin-1β-stimulated human chondrocytes,” British Journal of Pharmacology, vol. 160, no. 5, pp. 1252–1262, 2010.
View at: Publisher Site | Google Scholar
K. Lee, Y. H. Chung, H. Ahn, H. Kim, J. Rho, and D. Jeong, “Selective regulation of MAPK signaling mediates RANKL-dependent osteoclast differentiation,” International Journal of Biological Sciences, vol. 12, no. 2, pp. 235–245, 2016.
View at: Publisher Site | Google Scholar
M. Attur, I. Belitskaya-Lévy, C. Oh et al., “Increased interleukin-1β gene expression in peripheral blood leukocytes is associated with increased pain and predicts risk for progression of symptomatic knee osteoarthritis,” Arthritis & Rheumatism, vol. 63, no. 7, pp. 1908–1917, 2011.
View at: Publisher Site | Google Scholar
A. K. Singh, M. Kalaivani, A. Krishnan, P. K. Aggarwal, and S. K. Gupta, “Prevalence of osteoarthritis of knee among elderly persons in urban slums using american college of rheumatology (ACR) criteria,” Journal of Clinical and Diagnostic Research, vol. 8, no. 9, pp. 9–11, 2014.
View at: Publisher Site | Google Scholar
J. H. Kellgren and J. S. Lawrence, “Radiological assessment of osteo-arthrosis,” Annals of the Rheumatic Diseases, vol. 16, no. 4, pp. 494–502, 1957.
View at: Publisher Site | Google Scholar
V. G. Tusher, R. Tibshirani, and G. Chu, “Significance analysis of microarrays applied to the ionizing radiation response,” Proceedings of the National Academy of Sciences, vol. 98, no. 9, pp. 5116–5121, 2001.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2019 Ting Shi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1346

Downloads

1435

Citations