Molecular Medicine Reports Molecular Medicine Reports 1791-2997 1791-3004 D.A. Spandidos 10.3892/mmr.2016.5144 mmr-13-06-4599 Articles Comparison of molecular mechanisms of rheumatoid arthritis and osteoarthritis using gene microarrays LIHONGQIANG* HAOZHENYONG* ZHAOLIQIANG LIUWEI HANYANLONG BAIYUNXING WANGJIAN Department of Orthopedics, The Harbin Fifth Hospital, Harbin, Heilongjiang 150001, P.R. China Correspondence to: Professor Jian Wang, Department of Orthopedics, The Harbin Fifth Hospital, 27 Jiankang Road, Dongli, Harbin, Heilongjiang 150001, P.R. China, E-mail: jianwangjwjw@163.com

Contributed equally

 06 2016 15 04 2016 13 6 4599 4605 02 03 2015 11 12 2015 Copyright: © Li et al. 2016 This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made.

The present study aimed to compare the molecular mechanisms of rheumatoid arthritis (RA) and osteoarthritis (OA). The microarray dataset no. GSE29746 was downloaded from Gene Expression Omnibus. After data pre-processing, differential expression analysis between the RA group and the control, as well as between the OA group and the control was performed using the LIMMA package in R and differentially expressed transcripts (DETs) with |log2fold change (FC)|>1 and P<0.01 were identified. DETs screened from each disease group were then subjected to functional annotation using DAVID. Next, DETs from each group were used to construct individual interaction networks using the BIND database, followed by sub-network mining using clusterONE. Significant functions of nodes in each sub-network were also investigated. In total, 19 and 281 DETs were screened from the RA and OA groups, respectively, with only six common DETs. DETs from the RA and OA groups were enriched in 8 and 130 gene ontology (GO) terms, respectively, with four common GO terms, of which to were associated with phospholipase C (PLC) activity. In addition, DETs screened from the OA group were enriched in immune response-associated GO terms, and those screened from the RA group were largely associated with biological processes linked with the cell cycle and chromosomes. Genes involved in PLC activity and its regulation were indicated to be altered in RA as well as in OA. Alterations in the expression of cell cycle-associated genes were indicated to be linked with the occurrence of OA, while genes participating in the immune response were involved in the occurrence of RA.

 rheumatoid arthritis osteoarthritis differentially expressed transcripts functional annotation sub-network analysis Introduction

Arthritis comprises a variety of diseases with joint pain as a common feature and includes osteoarthritis (OA), rheumatoid arthritis (RA) and psoriatic arthritis, while their underlying mechanisms are divergent (1). OA is a degenerative disease commonly manifesting with mechanical abnormalities of weight-bearing joints and hands, including knees and hips. It is characterized by loss of matrix proteoglycans, fibrillation of cartilage surface and eventual loss of collagenous matrix. Substantial studies have proven that synovial membrane inflammation, abnormal articular chondrocyte differentiation and bone remodeling contribute to the progression of OA (2–4). RA is a systemic autoimmune disease, in which the immune system targets body cells of the same organism. It features chronic inflammation of the synovium and subsequent cartilage destruction as well as bone erosion. The complement system is known to be involved in the induction and progression of inflammatory reactions in RA (5,6).

Synovial fibroblasts (SFs), the most abundant resident cell type in human synovial tissue, are thought to have an important role in the pathogenesis of chronic arthritis (7) and display marked hyperplasia in OA and RA (8). Furthermore, alterations in the expression of various genes have also been observed to be associated with the phenotypic changes in OASFs and RASFs (9,10). Previous microarray studies have confirmed a relatively high heterogeneity of the RASF phenotype (11,12). By comparison between healthy SFs, RASFs and OASFs, Del Rey et al (13) found that OASFs possessed a more homogeneous phenotype compared to RASFs. The present study subjected the microarray data from Del Rey et al (13) to a bioinformatics analysis to identify common and differential molecular mechanisms underlying the two arthritis sub-types. The transcriptional expression profiles of the OASF and RASF samples were compared with those of SFs from healthy controls and key genes were identified.

 Materials and methods Microarray data

Gene microarray dataset GSE29746 was downloaded from Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/) database (13). The data had been collected from SF cultures obtained from nine patients with RA, 11 age- and gender matched patients with OA, and 11 age- and gender-matched adult healthy donors. The platform was the Agilent-014850 Whole Human Genome Microarray 4×44 K G4112F (Agilent Technologies, Santa Clara, CA, USA).

 Microarray data pre-processing

The gene expression profile data were extracted using the Linear Models for Microarray Data (LIMMA) package in R, followed by normexp background correction (14) and subsequent quantile normalization (15). The 95th quantile of all negative controls in each chip was calculated, and only probes with expression values larger than this value in all the samples were retained. According to the annotation platform, the values of probes corresponding to the same transcript were averaged and then defined as the final expression value of a transcript.

 Screening of differentially expressed transcripts (DETs)

Volcano plots were drawn to display the differential expression profiles of transcripts for each disease group and the control group. Differential expression analysis of transcripts between the RA group and the control group was performed using the t-test with the LIMMA package, as well as between the OA group and the control group (16). Transcripts with |log2fold change (FC)|>1 and P≤0.01 were screened as DETs. DETs from the two disease groups were compared, and a Venn diagram was used to display the result. A heatmap was also used to exhibit the overall expression profiles of DETs across all samples.

 Functional annotation of DETs

The screened DETs were submitted to the online tool Database for Annotation, Visualization and Integrated Discovery (DAVID; http://david.abcc.ncifcrf.gov/summary.jsp) to perform functional annotation based on the gene ontology (GO) database using Fisher's exact test (17).

 Analysis of interaction networks of DETs

DETs screened from each disease group were submitted to the Search Tool for the Retrieval of Interacting Genes/Proteins (http://www.stringdb.org/) to construct protein-protein interaction (PPI) networks based on the Biomolecular Interaction Network Database (18). The sub-networks within each network were then detected by clustering analysis using clusterONE (19), and genes involved in each cluster were then subjected to functional enrichment analysis with DAVID (P<0.05 as cut off value).

 Results Screening of differentially expressed transcripts

The differential expression profiles between each disease group and the control group were displayed in Volcano plots (Fig. 1). In total, 19 DETs were screened from the RA group, with |log2FC| values of 1–2.8, among which eight genes were downregulated and 11 genes were upregulated. Furthermore, 281 DETs were screened from the OA group, with |log2FC| values of 1.0–3.3, among which 113 genes were downregulated and 168 were upregulated. As illustrated by the Venn diagram (Fig. 2), the OA group contained a greater number of DETs than the RA group, while the two groups had six DETs in common (Table I). In addition, the heatmaps indicated that the disease samples may be separated from the control samples using the identified DETs (Fig. 3).

 Functional annotation of DETs from each disease group

According to functional annotation, DETs from the RA and OA groups were enriched in 8 and 130 GO terms, respectively. Certain DETs from the RA group (e.g. PLCH2, PLCB1, NPTX1 and MASP1) and from the OA group (e.g. F2RL2, PLCB4, PLCH2, PLCB1 and PLCXD1) were commonly enriched in four GO terms, including two molecular function terms associated with phospholipase C activity (Table II). In addition to the four common GO terms, DETs screened from the RA patients were also associated with immune response (e.g. MASP1, IL27RA and ITGB2), while those from the OA group were predominantly associated with the cell cycle and chromosomes, including NEK2, TTK, PTTG2, MAF, CENPN and SGOL2.

 Analysis of PPI networks of DETs

The interaction networks of DETs from the RA group and the OA group are shown in Figs. 4 and 5, respectively. It was discovered that the PPI network based on DETs from the OA group was more complex than that of the RA group. The interaction network of the RA group contained one sub-network, while that of the OA group had four sub-networks. DEGs including PSMB9, MAF, HTR2B and HTR2A were the hubs of corresponding PPI sub-networks of DEGs from the OA group, while PLCB1 was a hub of the PPI sub-network of DEGs from the RA group (Table III).

 Discussion

The present study revealed that in OA patients, the number of DEGs was higher compared to that in RA patients, which was consistent with the findings of Del Rey et al (13). This previous study identified 2,050 DEGs, several of which were also identified in the screening performed in the present study, such as ITGB2, PIP4K2C and NRXN2, although with different magnitudes of differential expression. Notably, the present study aimed to unravel the mechanisms underlying RA and OA from the perspective of PPI by building PPI networks, the constuction of which was not conducted in the previous study by Del Rey et al (13).

Two DEGs, PLCH2 and PLCB1, screened from the RA patients as well as the OA patients, were enriched in two PLC (phospholipase C) activity-associated GO terms, which may infer that the two arthritis types share certain common mechanisms regarding phospholipase C activity. PLCH2 and PLCB1 encode two members of the phosphoinositide-specific PLC superfamily, PLC-eta 2 and -beta 1, respectively, and PLC catalyzes the hydrolysis of phosphatidylinositol 4,5-bisphosphate into inositol 1,4,5-trisphosphate and 1,2-diacylglycerol) (20). Previous studies have reported a notable elevation of the pro-inflammatory enzyme PLA2 (phospholipase A2), another phospholipase type catalyzing the hydrolysis of membrane glycerophospholipids to release arachidonic acid and lysophospholipids in synovial fluids and sera of RA patients, and its expression has been proven to positively correlate with the disease activity in RA (21,22). Vignon et al (23) observed that PLA2 activity in RA and OA patients was similar, implying that pathological changes mediated by PLA2 are common in RA and OA. However, to the best of our knowledge, PLC activity has not been previously reported in RA or OA. As the two phospholipase types hydrolyze phospholipids at different sites, PLC may presumably also have a role in RA and OA, which requires further experimental validation.

Furthermore, in the present study, two PLC-regulating genes, HTR2B and HTR2A, which encode two members of the 5-hydroxytryptamine 2 receptor family that binds to the neurotransmitter serotonin, were also observed to be specifically differentially expressed in OA patients. These receptors activate PLC to initiate PLC-mediated signal transduction pathways (24). Of note, HTR2A was significantly upregulated, while HTR2B was significantly downregulated, implying their different roles in regulating PLC. The two genes were observed to be enriched in the GO biological process terms phosphoinositide-mediated signaling and second-messenger-mediated signaling as well as in the cellular component term plasma membrane-associated processes. This suggests that alterations in PLC-associated biological functions may be the predominant aberrations in OA patients.

In addition, DEGs screened from OA samples were also specifically and predominantly enriched in GO terms associated with the cell cycle and chromosomes. Among the abundance of DEGs screened from the OA samples, PSMB9 and MAF were the key hubs of the sub-PPI network. PSMB9, which was downregulated in OA, encodes a member of the proteasome B-type family, and proteasomes cleave peptides via an adenosine triphosphate/ubiquitin-dependent pathway. This result was in accordance with previous studies, as the ubiquitin-proteasome pathway has been implicated in the pathogenesis of OA (25), and Rollin et al (26) has reported that PSMB9 is aberrantly expressed in chondrocytes of OA patients. As the present study indicated that PSMB9 was mainly enriched in the GO biological process terms protein metabolic processes and cell cycle, it may be hypothesized that the downregulation of PSMB9 expression induces OA by disturbing normal protein metabolism and the cell cycle. MAF, a homolog of the avian musculoaponeurotic fibrosarcoma oncogene V-Maf, encodes a DNA-binding, leucine zipper-containing transcription factor. The increase in its expression in chondrocytes from patients with OA has been validated by Li et al (27), and it has been suggested that the upregulation of c-MAF expression may alter the phenotype of chondrocytes via regulating corresponding target genes (28), or interacting with other genes associated with chondrogenic differentiation (29). As in the present study, c-MAF was mainly enriched in chromosome-associated GO cellular component terms, it is presumed that it causes abnormal chondrogenic differentiation via inducing chromosomal abnormalities in OA patients.

Compared to those of OA, DETs screened from RA samples were exclusively enriched in immune response-associated GO biological process terms in addition to the four common GO terms, implying that the immune system has a particularly critical role in the occurrence of RA, which is consistent with the fact that RA is a systemic autoimmune disease. MASP1 encodes a serine protease that functions as a component of the lectin pathway of complement activation (30). The upregulation of MASP1 expression, which was also observed by Rioja et al (31), indicated enhanced activation of the lectin pathway as well as the activation of the complement system in RA. Enhanced complement activation has been indicated to be potentially associated with the occurrence and/or augmentation of inflammation in RA (5). Activated synovial fibroblasts and upregulation of the expression of various adhesion molecules that mediate their attachment to the cartilage have been confirmed in rheumatoid arthritis (32). Among these adhesion molecules, the expression of integrin α group proteins, including VLA-3, -4 and -5 has been indicated to increase most significantly (33,34). ITGB2 encodes integrin β2, the expression of which was observed to be upregulated in the present study, which was also reported in RASFs by Del Rey et al (35); it is therefore indicated that ITGB2 is linked with RA.

In conclusion, the present study indicated that genes involved in PLC activity, including PLCH2 and PLCB1, and its regulation, including HTR2A and HTR2B, are aberrantly expressed in RA as well as in OA. Alterations in the expression of genes associated with the cell cycle, including PSMB9 and MAF, were indicated to be linked with OA, while genes participating in the immune response were linked with RA.

 References RoivainenASöderströmKOPiriläLAroHKortekangasPMerilahti-PaloRYli-JamaTToivanenAToivanenPOncoprotein expression in human synovial tissue: An immunohistochemical study of different types of arthritisBrit J Rheumatol35933942199610.1093/rheumatology/35.10.933 PulligOWeselohGRonnebergerDKäkönenSMSwobodaBChondrocyte differentiation in human osteoarthritis: Expression of osteocalcin in normal and osteoarthritic cartilage and boneCalcif Tissue Int67230240200010.1007/s00223000110810954778 DrissiHZuscikMRosierRO'KeefeRTranscriptional regulation of chondrocyte maturation: Potential involvement of transcription factors in OA pathogenesisMol Aspects Med26169179200510.1016/j.mam.2005.01.00315811433 KwonDRParkGYIntra-articular injections for the treatment of osteoarthritis: Focus on the clinical use of several regimensOsteoarthritis-Diagnosis, Treatment and SurgeryChenQInTechRijeka671002012 NakagawaKSakiyamaHTsuchidaTYamaguchiKToyoguchiTMasudaRMoriyaHComplement C1 s activation in degenerating articular cartilage of rheumatoid arthritis patients: Immunohistochemical studies with an active form specific antibodyAnn Rheum Dis58175181199910.1136/ard.58.3.175103649161752845 ChenMDahaMRKallenbergCGThe complement system in systemic autoimmune diseaseJ Autoimmun34J276J286201010.1016/j.jaut.2009.11.014 BartokBFiresteinGSFibroblast-like synoviocytes: Key effector cells in rheumatoid arthritisImmunol Rev233233255201010.1111/j.0105-2896.2009.00859.x201930032913689 IzquierdoECañeteJDCelisRDel ReyMJUsateguiAMarsalSSanmartíRCriadoGPablosJLSynovial fibroblast hyperplasia in rheumatoid arthritis: Clinicopathologic correlations and partial reversal by anti-tumor necrosis factor therapyArthritis Rheum6325752583201110.1002/art.3043321547893 StanczykJOspeltCKarouzakisEFilerARazaKKollingCGayRBuckleyCDTakPPGaySKyburzDAltered expression of microRNA-203 in rheumatoid arthritis synovial fibroblasts and its role in fibroblast activationArthritis Rheum63373381201110.1002/art.30115212799943116142 KarouzakisEGayREGaySNeidhartMEpigenetic control in rheumatoid arthritis synovial fibroblastsNat Rev Rheumatol5266272200910.1038/nrrheum.2009.5519412193 Van BaarsenLGWijbrandtsCATimmerTCVan Der Pouw KraanTCTakPPVerweijCLSynovial tissue heterogeneity in rheumatoid arthritis in relation to disease activity and biomarkers in peripheral bloodArthritis Rheum6216021607201010.1002/art.2741520178127 KasperkovitzPVTimmerTCSmeetsTJVerbeetNLTakPPvan BaarsenLGBaltusBHuizingaTWPietermanEFeroMFibroblast-like synoviocytes derived from patients with rheumatoid arthritis show the imprint of synovial tissue heterogeneity: Evidence of a link between an increased myofibroblast like phenotype and high-inflammation synovitisArthritis Rheum52430441200510.1002/art.2081115692990 Del ReyMJUsateguiAIzquierdoECañeteJDBlancoFJCriadoGPablosJLTranscriptome analysis reveals specific changes in osteoarthritis synovial fibroblastsAnn Rheum Dis71275280201210.1136/annrheumdis-2011-200281 RitchieMESilverJOshlackAHolmesMDiyagamaDHollowayASmythGKA comparison of background correction methods for two-colour microarraysBioinformatics2327002707200710.1093/bioinformatics/btm41217720982 SmythGKSpeedTNormalization of cDNA microarray dataMethods31265273200310.1016/S1046-2023(03)00155-514597310 SmythGKLimma: Linear models for microarray dataBioinformatics and computational biology solutions using R and BioconductorGentlemanRCareyVJHuberWIrizarryRADudoitSSpringerNew York397420200510.1007/0-387-29362-0_23 Huang daWShermanBTLempickiRASystematic and integrative analysis of large gene lists using DAVID bioinformatics resourcesNat Protoc44457200910.1038/nprot.2008.21119131956 FranceschiniASzklarczykDFrankildSKuhnMSimonovicMRothALinJMinguezPBorkPvon MeringCJensenLJSTRING v9. 1: Protein-protein interaction networks, with increased coverage and integrationNucleic Acids Res41Database issueD808D815201310.1093/nar/gks10943531103 NepuszTYuHPaccanaroADetecting overlapping protein complexes in protein-protein interaction networksNat Methods9471472201210.1038/nmeth.1938224264913543700 McLaughlinSWangJGambhirAMurrayDPIP(2) and proteins: Interactions, organization, and information flowAnnu Rev Biophys Biomol Struct31151175200210.1146/annurev.biophys.31.082901.13425911988466 PruzanskiWKeystoneECSternbyBBombardierCSnowKMVadasPSerum phospholipase A2 correlates with disease activity in rheumatoid arthritisJ Rheumatol151351135519883199395 LinMFarewellVVadasPBookmanAAKeystoneECPruzanskiWSecretory phospholipase A2 as an index of disease activity in rheumatoid arthritis. Prospective double blind study of 212 patientsJ Rheumatol231162116619968823686 VignonEBalblancJCMathieuPLouisotPRichardMMetalloprotease activity, phospholipase A2 activity and cytokine concentration in osteoarthritis synovial fluidsOsteoarthritis Cartilage1115120199310.1016/S1063-4584(05)80026-38886087 BergKAMaayaniSClarkeWPInteractions between effectors linked to serotonin receptorsAnn NY Acad Sci861111120199810.1111/j.1749-6632.1998.tb10181.x JangBCSungSHParkJGParkJWBaeJHShinDHParkGYHanSBSuhSIGlucosamine hydrochloride specifically inhibits COX-2 by preventing COX-2 N-glycosylation and by increasing COX-2 protein turnover in a proteasome-dependent mannerJ Biol Chem2822762227632200710.1074/jbc.M61077820017635918 RollinRMarcoFCamafeitaECalvoELópez DuránLJoverJLópezJAFernández-GutiérrezBDifferential proteome of bone marrow mesenchymal stem cells from osteoarthritis patientsOsteoarthritis Cartilage16929935200810.1016/j.joca.2007.12.00618222713 LiTXiaoJWuZQiuGOver-expression of c-maf by chondrocytes in osteoarthritisJ Int Med Res37129135200910.1177/14732300090370011519215682 MacleanHEKimJIGlimcherMJWangJKronenbergHMGlimcherLHAbsence of transcription factor c-maf causes abnormal terminal differentiation of hypertrophic chondrocytes during endochondral bone developmentDev Biol2625163200310.1016/S0012-1606(03)00324-514512017 HuangWLuNEberspaecherHDe CrombruggheBA new long form of c-Maf cooperates with Sox9 to activate the type II collagen geneJ Biol Chem2775066850675200210.1074/jbc.M20654420012381733 GoeldnerISkareTBoldtABNassFRMessias-ReasonIJUtiyamaSRAssociation of MASP-2 levels and MASP2 gene polymorphisms with rheumatoid arthritis in patients and their relativesPLoS One9e90979201410.1371/journal.pone.0090979246325983954616 RiojaIClaytonCLGrahamSJLifePFDicksonMCGene expression profiles in the rat streptococcal cell wall-induced arthritis model identified using microarray analysisArthritis Res Ther7R101R117200510.1186/ar1458156421301064886 NikkariLAhoHYli-JamaTLarjavaHJalkanenMHeinoJExpression of integrin family of cell adhesion receptors in rheumatoid synovium. Alpha 6 integrin subunit in normal and hyperplastic synovial lining cell layerAm J Pathol14210191027199376827571886885 RinaldiNSchwarz-EywillMWeisDLeppelmann-JansenPLukoschekMKeilholzUBarthTFIncreased expression of integrins on fibroblast-like synoviocytes from rheumatoid arthritis in vitro correlates with enhanced binding to extracellular matrix proteinsAnn Rheum Dis564551199710.1136/ard.56.1.4590591411752256 IshikawaHHirataSAndohYKuboHNakagawaNNishibayashiYMizunoKAn immunohistochemical and immunoelectron microscopic study of adhesion molecules in synovial pannus formation in rheumatoid arthritisRheumatol Int165360199610.1007/BF018164368853226 Del ReyMJIzquierdoEUsateguiAGonzaloEBlancoFJAcquadroFPablosJLThe transcriptional response of normal and rheumatoid arthritis synovial fibroblasts to hypoxiaArthritis Rheum6235843594201010.1002/art.2775020848564

Volcano plots of expression values of all the transcripts. RA, synovial fibroblasts from rheumatoid arthritis group; OA, synovial fibroblasts from the osteoarthritis group; HSF, healthy synovial fibroblasts; DEG, differentially expressed gene.

Venn diagram of differentially expressed transcripts between the rheumatoid arthritis group and the osteoarthritis group. RA, synovial fibroblasts from rheumatoid arthritis group; OA, synovial fibroblasts from the osteoarthritis group.

Clustering analysis of samples based on the screened differentially expressed transcripts. (A) Clustering analysis was based on the expression levels of 19 transcripts differentially expressed between the RA group and the healthy controls. (B) Clustering analysis based on the expression levels of 281 transcripts differentially expressed between the synovial fibroblasts in the OA group and the healthy controls. Each row represents a single gene and each column represents a tissue sample. The colored bar at the top indicates the group: Gray, HSF; blue, synovial fibroblasts from RA group; pink, synovial fibroblasts from OA group. Expression levels were indicated in a heat-map style with dark red indicating high expression and dark green low expression. RA, rheumatoid arthritis; OA, osteoarthritis; HSF, healthy synovial fibroblasts.

Interaction network of differentially expressed transcripts screened from synovial fibroblasts of patients with rheumatoid arthritis. Green squares represent downregulated DEGs and red squares indicate upregulated DEGs. DEG, differentially expressed gene.

Interaction network of DEGs screened from synovial fibroblasts from patients in osteoarthritis group. Green squares represent downregulated DEGs and red squares indicate upregulated DEGs. DEG, differentially expressed gene.

Common differentially expressed transcripts between the RA group and the OA group.

Transcript Gene symbol RA group
OA group
logFC P-value logFC P-value
NM_182734 PLCB1 −1.18252 7.74×10−5 −1.03053 7.02×10−3
NM_000211 ITGB2 2.81476 4.41×10−4 3.31223 3.08×10−5
NM_001853 COL9A3 −1.12924 4.07×10−3 −1.03848 5.18×10−3
NM_139125 MASP1 −1.31698 4.33×10−3 −1.27800 3.61×10−3
NM_002522 NPTX1 −1.69186 5.62×10−3 −1.48981 9.55×10−3
NM_014638 PLCH2 1.37193 7.12×10−3 1.29122 7.54×10−3
A_32_P4882 NA −1.62833 8.70×10−3 −2.66329 3.58×10−5

RA, the rheumatoid arthritis; OA, the osteoarthritis; FC, fold change.

Common GO terms between the RA group and the OA group.

GO term and function Differentially expressed genes
RA group OA group
GO:0004435 - Phosphoinositide phospholipase C activity PLCH2, PLCB1 F2RL2, PLCB4, PLCH2, PLCB1
GO:0004629 - Phospholipase C activity PLCH2, PLCB1 F2RL2, PLCB4, PLCH2, PLCB1, PLCXD1
GO:0005509 - Calcium ion binding NPTX1, MASP1, PLCH2, PLCB1 F10, MASP1, NRXN2, PCDHB3, PCDHB2, VWCE, COLEC12, GRIN3A, C9ORF140, SLIT2, NPTX1, CDH15, CLGN, PLCB4, HPSE, PLCH2, PCDHB16, PLCB1, THBS1, DTNA
GO:0007267 – Cell-cell signaling NPTX1, WISP1, ITGB2 FGF18, ACHE, NTF3, PCDHB3, DLGAP5, EFNB2, MME, PCDHB2, ITGB2, RIMS1, NPTX1, PCDHB16, DTNA, HTR2A

RA, rheumatoid arthritis; OA, the osteoarthritis; GO, gene ontology.

Identified sub-networks and functional annotation.

Variable OA group
RA group
Cluster1 Cluster2 Cluster3 Cluster4 Cluster1
Nodes (n) 4 4 3 6 3
Density 1 0.667 0.667 0.6 1
Quality 1 1 1 0.529 1
P-value 0.008 0.011 0.026 0.046 0.016
Genes ABCD3, PEX19, ABCD2, ABCD1 BACH1, MAF, MAFB, ATF4 HTR2A, MPDZ, HTR2B NCOA3, ESR1, NCOA1, PSMB9, GRIP1, ESR2 PLCB1, PARD3, PARD6A
Top significant GO terms Peroxisomal membrane Sequence-specific DNA binding Serotonin binding Steroid hormone receptor signaling pathway Occluding junction

RA, rheumatoid arthritis; OA, osteoarthritis.