PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of peerjLatest ArticlesFor AuthorsEditorial BoardPeerJPeerJ
 
PeerJ. 2017; 5: e3078.
Published online 2017 March 15. doi:  10.7717/peerj.3078
PMCID: PMC5356478

Identification of dysregulated genes in rheumatoid arthritis based on bioinformatics analysis

Academic Editor: Elena Papaleo

Abstract

Background

Rheumatoid arthritis (RA) is a chronic auto-inflammatory disorder of joints. The present study aimed to identify the key genes in RA for better understanding the underlying mechanisms of RA.

Methods

The integrated analysis of expression profiling was conducted to identify differentially expressed genes (DEGs) in RA. Moreover, functional annotation, protein–protein interaction (PPI) network and transcription factor (TF) regulatory network construction were applied for exploring the potential biological roles of DEGs in RA. In addition, the expression level of identified candidate DEGs was preliminarily detected in peripheral blood cells of RA patients in the GSE17755 dataset. Quantitative real-time polymerase chain reaction (qRT-PCR) was conducted to validate the expression levels of identified DEGs in RA.

Results

A total of 378 DEGs, including 202 up- and 176 down-regulated genes, were identified in synovial tissues of RA patients compared with healthy controls. DEGs were significantly enriched in axon guidance, RNA transport and MAPK signaling pathway. RBFOX2, LCK and SERBP1 were the hub proteins in the PPI network. In the TF-target gene network, RBFOX2, POU6F1, WIPF1 and PFKFB3 had the high connectivity with TFs. The expression status of 11 candidate DEGs was detected in GSE17755, the expression levels of MAT2A and NSA2 were significantly down-regulated and CD47 had the up-regulated tendency in peripheral blood cells of patients with RA compared with healthy individuals. qRT-PCR results of MAT2A, NSA2, CD47 were compatible with our bioinformatics analyses.

Discussion

Our study might provide valuable information for exploring the pathogenesis mechanism of RA and identifying the potential biomarkers for RA diagnosis.

Keywords: Differentially expressed genes, Rheumatoid arthritis, Regulatory network, Integrated analysis, Pathogenesis

Introduction

Rheumatoid arthritis (RA) is a chronic autoinflammatory disorder of joints, which is characterized by irreversible joint damage and destruction of cartilage and bone, leading to severe disability or even premature death (Aletaha et al., 2010). It affects approximately 1% of the population worldwide, and reduces life expectancy of patients by 3–18 years. There are 25–50 new cases in a population of 100,000 worldwide each year (Gibofsky, 2014; Pincus & Callahan, 1986; Michelsen et al., 2015).

The etiology and pathophysiology of RA are not well understood (McInnes & Schett, 2011). It is widely acknowledged that RA is triggered by the complex interaction among genetic susceptibility, autoimmune disorders and dysregulated genes.

For genetic susceptibility, several studies have demonstrated that expression quantitative trait loci (eQTL) single nucleotide polymorphisms (SNPs) increase the risk of RA. A recent study indicates that HLA-DOA rs369150-A polymorphism reduces the expression of HLA-DOA, which independently associates with RA risk in the Japanese population, Asian population (Chinese and Korean) and European population (Okada et al., 2016). eQTL of GZMB identifies higher expression in the presence of the minor allele of rs8192916, which is associated with progression of joint destruction in RA (Knevel et al., 2013). For dysregulated genes, down-regulation of miR-126 inhibits PI3K/AKT signaling pathway and disrupts the imbalance between cell growth and cell death of rheumatoid arthritis synovial fibro-blasts by targeting PIK3R2 (Gao et al., 2016). Over-expression of TLR-2 and IL-6 decrease the expression of FCRL4, might contribute to the RA process (Khanzadeh et al., 2016; Eser & Sahin, 2016).

Previous studies demonstrate that early diagnosis and treatment will be beneficial to reduction of progressive joint injury of patients with RA. Currently, several biomarkers such as citrullinated proteins, rheumatoid factor, S100 proteins, matrix metalloproteinases and serum amyloid proteins have been described (Smolen, Aletaha & McInnes, 2016), but these biomarkers have limitations on sensitivity and or specificity in routine clinical practice. The early diagnosis of RA remains challenging.

In this study, we used bioinformatics methods to integrate mRNA expression datasets, which were available in the Gene Expression Omnibus database, identified DEGs in synovial tissues of patients with RA and healthy individuals, aiming to investigate the potential pathogenesis mechanism of RA and provide valuable information for the identification of potential diagnosis biomarkers in RA.

Materials and Methods

Data collection

NCBI Gene Expression Omnibus (GEO) database provides a large collection of microarray expression data (Edgar, Domrachev & Lash, 2002). “Rheumatoid Arthritis” was used as the key word to search expression profiling of RA. The inclusion criteria of datasets were set as follows: (1) whole genome expression data of RA; (2) the expression profiling of synovial tissues was generated from patients with RA and healthy individuals.

Total of nine mRNA expression datasets of synovial tissues of RA were incorporated into our study. Four datasets were based on GPL570 [HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array, four datasets were based on GPL96 [HG-U133A] Affymetrix Human Genome U133A Array, and one dataset was based on GPL91 [HG_U95A] Affymetrix Human Genome U95A Array.

Data preprocessing and Identification of DEGs

The raw expression datasets were downloaded and preprocessed by log2 transformation in R language (Shi et al., 2006). The Linear Models “limma” package in R language was used to analyze the microarray datasets (Smyth, 2005). Differentially expressed genes were identified in patients with RA compared to healthy individuals. The false discovery rate (FDR) (Reiner-Benaim, 2007) was utilized for multiple testing corrections by using the Benjamini and Hochberg method. FDR < 0.05 was set as the threshold of DEGs.

Heatmap analysis

To assess the similarity of gene expression patterns between two samples, two-way hierarchical clustering analysis (Cao, 2016) was performed. The heatmap was constructed using ‘pheatmap’ package in R language.

Functional annotation of DEGs

The significantly dysregulated genes between patients with RA and healthy individuals were filtered as DEGs. In order to understand the biological roles of DEGs, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and Gene Ontology (GO) terms were enriched by the online tool GeneCoDis3 (Kanehisa, 2002; Carmona-Saez et al., 2007; Torto-Alalibo et al., 2014). The threshold of GO function and KEGG pathway enrichment were FDR < 0.05.

Protein–protein interaction network

In order to explore the interaction between the top 10 up- and down-regulated DEGs in synovial tissues of patients with RA, BioGRID, a database of known protein interactions, was used to predict the PPI association among DEGs and the PPI interaction network was visualized by Cytoscape (Shannon et al., 2003; Chatr-Aryamontri et al., 2015).

Construction of transcription factors regulatory network

TRANSFAC (Wingender et al., 1996) is a database covering details of TFs and their DNA binding sites. In our analysis, the TFs in the human genome and the motifs of genomic binding sites were downloaded from the TRANSFAC. Moreover, the 2 KB sequence in the upstream promoter region of top 10 up- and down-regulated DEGs was downloaded from UCSC (http://www.genome.ucsc.edu/cgi-bin/hgTables). Target sites of potential TFs were then distinguished. The TFs and its target genes were used for constructing regulatory network, and the network was visualized by Cytoscape software (http://cytoscape.org/) (Shannon et al., 2003). Nodes represented TFs or target genes and lines represented association between TFs and target genes.

The DEGs were analyzed in the GSE17755 datasets

We obtained a number of key DEGs in RA synovial tissues based on the bioinformatics analyses. In order to analyze whether those DEGs were dysregulated in peripheral blood cells of patients with RA, the expression levels of representative DEGs were preliminarily validated in the expression profiling of peripheral blood cells of patients with RA from the publicly available GEO dataset. “Rheumatoid arthritis” was used as the key word to search the expression profiling in the GEO database. The mRNA expression profiling datasets were included based on the following criteria: (1) the dataset was generated from the peripheral blood cells of patients with RA; (2) The expression profiling of peripheral blood cells of RA patients and healthy individuals were available in the dataset; (3) the RA patients had no treatment record before peripheral blood cells collection. Finally, GSE17755 and GSE15573 were left. A total of 18 patients with RA and 15 healthy individuals were included in the GSE15573 datasets, which had the smaller sample size compared to GSE17755. There were 112 patients with RA and 45 healthy individuals in GSE17755 dataset. Based on aforementioned information, GSE17755 was incorporated into ours study (Lee et al., 2011). Box-plot analysis was performed to preliminarily detect the expression levels of DEGs in peripheral blood cells of RA patients and health individuals.

Quantitative real-time polymerase chain reaction (qRT-PCR)

In order to validate the expression levels of candidate genes in RA, 10 peripheral blood samples were collected from RA patients, who were diagnosed at the Affiliated Zhongshan Hospital. Moreover, 10 matched peripheral blood samples were collected from healthy individuals. This work was approved by the Ethics Committee of the Affiliated Zhonghan Hospital and informed written consent was obtained from all patients. The research complied with the principles of the Declaration of Helsinki.

Total RNA of peripheral blood samples were extracted by using Trizol (Invitrogen, Carlbad, CA, USA) according to the manufacture instructions. FastQuant cDNA Synthesis Kit (Tiangen, Beijing, China) was used to synthesize the cDNA of mRNA. qRT-PCR reactions were performed by using SuperReal PreMix Plus SYBR Green Kit (Tiangen, Beijing, China) on Applied Biosystems 7500 (Applied Biosystems, Foster City, CA, USA). GAPDH were used as internal control for mRNA detection. The PCR primers used in our study were shown as follows: GAPDH forward primer 5′ GGA GCG AGA TCC CTC CAA AAT 3′, reverse primer 5′ GGC TGT TGT CAT ACT TCT CAT GG 3′; MAT2A forward primer 5′ ATG AAC GGA CAG CTC AAC GG 3′, reverse primer 5′ CCA GCA AGA AGG ATC ATT CCA G 3′; NSA2 forward primer 5′ CAC CGT AAA CGC TAT GGA TAC C 3′, reverse primer 5′ GCT TAG CCT TCA GAC CAA TCA TT 3′; CD47 forward primer 5′ TCC GGT GGT ATG GAT GAG AAA 3′, reverse primer 5′ ACC AAG GCC AGT AGC ATT CTT 3′. At least triple experiments were subjected to qRT-PCR verification. The relative expression of candidate genes was calculated by using the ΔCT equation methods. Mean ± standard deviation and independent-samples t-test was used in the statistical analysis. P < 0.05 was considered as significant difference. * indicated P < 0.05; ** indicated P < 0.01.

Results

Differentially expressed genes in RA

Nine mRNA expression profiles of RA synovial tissues, including 168 patients with RA and 41 healthy individuals, were incorporated into our study (Table 1). In our analysis, a total of 378 DEGs including 202 up- and 176 down-regulated genes were identified in patients with RA compared to healthy individuals. LCK and MAT2A were the most significantly up- and down-regulated genes in RA, respectively (Table 2). The expression pattern of the top 200 DEGs in RA patients and healthy individuals was shown in Fig. 1, which indicated that the mRNA expression profiling of synovial tissues of RA patients was different from healthy individuals.

Table 1
The expression profiling of patients with RA and normal controls.
Table 2
The top 10 up- and down-regulated DEGs in RA.
Figure 1
Hierarchical clustering analysis based on the expression profile of the top 200 discriminatory DEGs between synovial tissues of RA and normal controls.

GO classification of DEGs in RA

DEGs in patients with RA were performed to GO annotation for investigating the biological roles. The threshold of the GO terms was set as FDR < 0.001. As Table 3 shown, RNA splicing (GO: 0007585), signal transduction (GO: 0007165) and gene expression (GO: 0010467) were the highly significant enrichment of the GO biological process. Furthermore, protein binding (GO: 0005515), nucleotide binding (GO: 0000166) and ATP binding (GO: 0005524) were the highest enrichment of the GO molecular function.

Table 3
GO annotation of dysregulated DEGs in RA.

Pathway enrichment analysis

KEGG analysis was used to understand the signaling pathway enrichment of DEGs in RA. The threshold was FDR < 0.05. As Table 4 shown, the significant enrichment of pathways focused on axon guidance (hsa04360), RNA transport (hsa03013), protein processing in endoplasmic reticulum (hsa04141) and MAPK signaling pathway (hsa04010).

Table 4
KEGG enrichment of dysregulated DEGs in RA.

Protein-protein interaction network

The PPI network between the top 10 up- and down-regulated DEGs in RA were constructed by Cytoscape based on the BioGRID database (Fig. 2). The network consisted of 419 nodes and 480 edges. In the network, the nodes with high degree were defined as hub proteins. The hub proteins in the PPI network including RBFOX2 (degree = 90), LCK (degree = 78), SERBP1 (degree = 45), TCOF1 (degree = 36) and TARS (degree = 32), as Fig. 2 shown. A total of 402 predicted proteins interaction with top 10 up- and down-regulated were overlapped with 202 up- and 168 down-regulated DEGs, respectively. Lastly, 15 DEGs interacted with top 10 up- and down-regulated DEGs were obtained and their interactions were shown in Table S2. RBFOX2 interacted with 6 dysregulated DEGs including DDX5, ATN1, ESR2, RPA2, RRPA3 and ATXN2; LCK interacted with WASL; SERBP1 interacted with PIAS1; TCOF1 interacted with RPA2, RPA3, TOP1 and NOP56.

Figure 2
The protein-protein interaction network of the top 10 up- and down-regulated DEGs in synovial tissues of RA.

Construction of transcription factors regulatory network

In our study, TFs regulatory network under the regulation of TFs and target genes were constructed. Based on the TRANSFAC database, total of 42 proteins were identified as TFs of top 10 up- and down-regulated DEGs in RA. The identified TFs and target genes were applied to construct the TFs-target genes network by Cytoscape software. As shown in Fig. 3, the network consisted of 61 nodes and 142 edges. The DEGs, such as RBFOX2 (degree = 16), POU6F1 (degree = 11), WIPF1 (degree = 10), PFKFB3 (degree = 9) and EIF5B (degree = 9) had the high connectivity with TFs.

Figure 3
The transcription factors-target genes regulatory network of top 10 up and down-regulated DEGs in synovial tissues of RA.

The expression levels of DEGs were analyzed in the GSE17755 datasets

GSE17755 covered the expression profiling of peripheral blood cells from 112 RA patients and 45 health individuals. The descriptive statistics (spread and distribution) of DEGs expression in RA patients compared to healthy individuals were depicted by using box-plot analysis, which was visually illustrated by median and inter-quartile range. As Figs. 4F, ,4G4G and and4K4K shown, the expression levels of MAT2A and NSA2 were significantly down-regulated and CD47 had the up-regulated tendency in peripheral blood cells of patients with RA, which were consistent with our integrated analysis in RA synovial tissues. However, the expression of LCK, SERBP1, WIPF1, PFKFB3, RBFOX2, ATF6B, POU6F1 and CA5A in GSE17755 was incompatible with our integrated analysis in RA synovial tissues (Figs. 4A4E, ,4H4H4J).

Figure 4
The expression levels of DEGs was analyzed in the GSE17755 dataset.

qRT-PCR validation of the expression levels of candidate genes in RA

qRT-PCR was subjected to validate the expression levels of dysregulated genes in peripheral blood samples of 10 RA patients and 10 healthy individuals. Three candidate DEGs including MAT2A, NSA2 and CD47 were used to apply for qRT-PCR verification. As Figs. 5A and and5B5B shown, MAT2A (P < 0.01) and NSA2 (P < 0.05) were obviously down-regulated in peripheral blood samples of RA patients compared with healthy individuals. However, the expression level of CD47 (P > 0.05) in peripheral blood samples was not significant difference between RA patients and healthy individuals, but its expression displayed the up-regulated tendency in peripheral blood samples of RA patients.

Figure 5
qRT-PCR validation of candidate DEGs in peripheral blood samples of RA patients and healthy individuals.

Discussion

Rheumatoid arthritis, a chronic inflammatory disease of joint, is triggered by the complex interaction between genetic susceptibility and dysregulated genes. The clinical manifestations of RA patients present pannus and the abnormal proliferative synovial tissue. Early diagnosis would improve the prognosis and the life quality of RA patients. However, the pathogenesis mechanism of RA is ambiguous and there are no feasible biomarkers for early diagnosis in clinical practice. In our study, bioinformatics analysis-based identification of key genes and pathways in RA were conducted, aiming to provide valuable ground work for further investigation of pathogenesis mechanism of RA.

RBFOX2 belongs to binding protein, is thought to be a key regulator of alternative exon splicing in the nervous system and other cell types (Venables et al., 2013a). It was one of top 10 up-regulated genes in RA synovial tissues (Table 2). In the TFs regulatory network, RBFOX2 had the highest connectivity with TFs, was regulated by 16 TFs. RBFOX2 is involved in various diseases and biological processes including congenital heart disease, breast cancer, epithelial-mesenchymal transition and pluripotent stem cell differentiation (Homsy et al., 2015; Braeutigam et al., 2014; Cen et al., 2013; Venables et al., 2013b). The mutations of RBFOX2 are associated with congenital heart disease and neuro developmental disabilities (Homsy et al., 2015). During the epithelial-mesenchymal transition, RBFOX2 is an essential regulator to mediate cellular invasion (Braeutigam et al., 2014). To our knowledge, our study first reported the aberrant expression of RBFOX2 in RA, however, the biological function of RBFOX2 in RA needs to be further elucidated.

CD47 encodes CD47 molecule, a membrane protein, is involved in the increase of intracellular calcium concentration and may play an essential role in membrane transport and signal transduction (Kaur et al., 2015). CD47 was significantly down-regulated in RA (Table 2). In the TFs regulatory network, CD47 was targeted by 6 TFs including v-Myb, SOX-9, HNF-1, COMP-1, E2F and Elk-1. In RA, the interaction between CD47 and thrombospondin-1 could trigger T cell infiltration and expansion in the rheumatoid synovium, and perpetuates the inflammatory process in the rheumatoid joint (Vallejo et al., 2003). CD47 might play essential roles in the progression of RA. In our study, the expression level of CD47 was increased in synovial tissues and peripheral blood cells of patients with RA (Fig. 4F).

MAT2A encodes methionine adenosyltransferase 2A, is an enzyme that catalyzes the production of S-adenosylmethionine (AdoMet) from methionine and ATP, and is a key methyl donor in cellular processes (Panayiotidis et al., 2006). In our study, MAT2A was the most significantly down-regulated DEG in synovial tissues of RA (Table 2). The previous studies display that MAT2A is associated with uncontrolled cell proliferation in cancer. The protein expression level of MAT2A is decreased in renal cell carcinoma compared to normal tissues (Wang et al., 2014), while MAT2A is over-expressed in various gastrointestinal cancers, such as gastric cancer, colon cancer and liver cancer (Zhang et al., 2013; Tomasi et al., 2013). miR-21-3p targets down-regulation of MAT2A and inhibits cell growth in hepatoma (Lo, Tsai & Chen, 2013). Based on aforementioned information, down-regulated MAT2A in synovial tissues of RA might implicate in cell proliferation, invasiveness of fibroblast-like synoviocytes, which results in erosion of bone and cartilage in RA. The expression of MAT2A in synovial tissues and FLS needs to be validated through quantitative real-time polymerase chain reaction in a larger sample size of patients with RA and healthy individuals, in addition, the biological roles of MAT2A in FLS needs to be further investigated through in vivo and in vitro experiment. MAT2A had the lower expression in both of peripheral blood cells and synovial tissues of RA patients compared to healthy individuals (Fig. 4G).

NSA2 encodes NSA2 (Nop seven-associated 2), ribosome biogenesis homolog, locates in the nucleolus of the cell, and is involved in cell cycle regulation and proliferation (Zhang et al., 2010). It had the lower expression level in synovial tissues and peripheral blood cells of RA patients compared to healthy individuals (Table 2, Fig. 4K). Over-expression of the NSA2 protein promotes cell growth and regulates the G1/S transition in the cell cycle in different cell lines including HeLa (human cervical carcinoma), HEK293T (human embryonic kidney) and A549 (human epithelial lung adenocarcinoma) (Zhang et al., 2010; Li et al., 2013). The association between dysregulated NSA2 and pathogenesis of RA has not been reported. The functions of NSA2 in initiation and progression of RA need to be elucidated in the future study.

Dysregulated DEGs were identified in synovial tissues of patients with RA, and the mRNA expression level of representative DEGs was preliminarily detected in peripheral blood cells of patients with RA based on GSE17755 dataset. As Fig. 4 shown, LCK, SERBP1, WIPF1, PFKFB3 and RBFOX2 were down-regulated in peripheral blood cells of RA, but those genes were up-regulated in synovial tissues of RA; POU6F1 and CA5A were down-regulated in peripheral blood cells of RA, but they were up-regulated in synovial tissues of RA. The contradictory expression status of DEGs in synovial tissues and peripheral blood cells might attribute to different characteristics of expression profiling between synovial tissues and peripheral blood cells in patients with RA. In our bioinformatics analyses, the expression levels of MAT2 and NSA2 were obviously down-regulated both in synovial tissues and peripheral blood samples of RA patients; and CD47 expression was up-regulated both in synovial tissues and with the up-regulated tendency in peripheral blood samples of RA patients. The expression status of MAT2, NSA2 and CD47 in peripheral blood samples of RA patients were validate through qRT-PCR. The qRT-PCR results indicated MAT2A and NSA2 were significantly down-regulated in peripheral blood samples of RA patients and CD47 had the up-regulated tendency in peripheral blood samples of RA patients compared with healthy individuals (Fig. 5). In general, qRT-PCR results of MAT2A, NSA2 and CD47 were compatible with our bioinformatics analyses. The associations between a majority of identitied TFs in the TF-gene regulatory network and RA remain unknown. The roles of a number of identified TFs in other autoimmue disorders are explored. c-Ets-1 belongs to the ETS family of transcript factors and regulates the expression of various immune-related genes. In our work, it was predicted to bind with the promoter of SAR1A, SERBP and TARS. Over-expression of c-Ets-1 in synovial membrane regulates inflammatory angiogenesis in RA through being activated by interleukin-1 and tumor necrosis factor alpha (Wernert et al., 2002; Redlich et al., 2001). A eQTL analysis indicates that the allele of rs6590330 in ETS1 is associated with decreased ETS1 expression and increases the risk of systemic lupus erythematosus by enhancing the binding of pSTAT1 (Lu et al., 2015). The regulatory significance of identified TFs in development and progression of RA needs to be profoundly investigated in our future work.

The dysregulated genes in RA were significantly enriched in KEGG pathways including axon guidance, RNA transport and MAPK signaling pathway. Axon guidance represents a key stage in the formation of neuronal network. Axons are guided by a variety of guidance factors, such as netrins, ephrins, slits, and semaphorins. Catalano reports that semaphoring-3A (SEMA3A), a number of semaphorin family, attenuates inflammation and progression of collagen-induced arthritis (Catalano, 2010). SEMA3A is decreased in the immune system of experimental autoimmune encephalomyelitis murine model of multiple sclerosis, which is one of autoimmunity diseases (Gutierrez-Franco et al., 2016). In our study, down-regulated SEMA6A and up-regulated SEMA4D were significantly enriched in axon guidance pathway. THOC2 (interacted with NR5A1 in PPI network) and EIF58 (one of the top 10 up-regulated DEGs in RA) were significantly enriched in RNA transport pathways. The mitogen-activated protein kinase (MAPK) cascade is involved in various cellular functions, including cell proliferation, differentiation and migration. A series of articles depict that MAPK signaling pathway associates with RA (Pearson & Jones, 2016; Weisbart et al., 2013; Schett, Zwerina & Firestein, 2008). Increased expression of phosphorylation ezrin may contribute to migration and invasion of fibroblast-like synoviocytes in RA, which are mediated by rho kinase and the p38 MAPK pathway phosphorylation (Xiao et al., 2014). BRAF splice variants activate MAPK through CRAF, increase expression of MT1-MMP, and enhance fibroblast invasion of collagen in RA. Based on the abovementioned information, axon guidance, RNA transport and MAPK signaling pathway might contribute to the progression of RA.

In our study, RBFOX2 and those dysregulated DEGs interacted with it including DDX5, ATN1, ESR2, RPA2, RPA3 and ATXN2 were significantly enriched in six significant enrichment of biological process and molecular function, such as RNA splicing, signaling transduction, gene expression, protein binding , nucleotide binding and ATP binding (Table S1). NPM1 (interacted with NAS2) was significantly enriched in signal transduction and protein binding (Table S1). Our finding indicated that the protein-protein interactions between dysregulated DEGs might play essential roles in progression and development of RA through those biological processes and molecular functions.

There are limitations in our study. Firstly, further biological function studies should be performed to investigate the roles of these DEGs and pathways in the progression of RA. Secondly, whether those dysregulated genes in RA have the diagnostic value in RA diagnosis needs to be validated through large cohorts of RA patients and healthy individuals.

Conclusions

In conclusion, we identified 378 DEGs in synovial tissues of RA patients compared to healthy individuals. The top 10 up- and down-regulated genes were selected to establish the PPI network and TF-target genes network. The expression levels of candidate DEGs was preliminarily detected in peripheral blood cells of RA according to GSE17755 dataset. Our study indicated that dysregulated genes in RA might contribute to the progression of RA by regulating biological process, molecular function and signaling pathways based on complex gene regulatory network interaction, such as PPI network and TF-genes regulatory network. Our findings may pave the road for exploring the underlying mechanisms of RA and identifying the potential biomarkers for RA diagnosis.

Supplemental Information

10.7717/peerj.3078/supp-1

Table S1

The top 10 up-/down-regulated DEGs and those DEGs interacted with top 10 up-/down-regulated DEGs enriched in GO terms:

DEGs, differentially expressed genes; FDR, false discovery rate; NA, not available.

10.7717/peerj.3078/supp-2

Table S2

The DEGs interacted with the top 10 up- and down-regulated DEGs in RA:

NA: not applicable; DEGs: differentially expressed genes; RA: rheumatoid arthritis. The predicted proteins of TRD, POU6F1 and CA5A were not available in the BioGRID database and were not shown in Table S2.

Funding Statement

The authors received no funding for this work.

Additional Information and Declarations

Competing Interests

Haiwei Du is an employee of Beijing Medintell Biomed Co., Ltd. Ruihu Hao and Haiwei Du equally contribute to the work and are co-first authors.

Author Contributions

Ruihu Hao and Haiwei Du performed the experiments, analyzed the data, contributed reagents/materials/analysis tools, wrote the paper, prepared figures and/or tables, reviewed drafts of the paper.

Lin Guo conceived and designed the experiments, wrote the paper, prepared figures and/or tables, reviewed drafts of the paper.

Fengde Tian and Ning An performed the experiments, reviewed drafts of the paper.

Tiejun Yang and Changcheng Wang analyzed the data, reviewed drafts of the paper.

Bo Wang and Zihao Zhou contributed reagents/materials/analysis tools, reviewed drafts of the paper.

Data Availability

The following information was supplied regarding data availability:

NCBI Gene Expression Omnibus database: GSE77298, GSE48780, GSE55235; GSE55457, GSE55584, GSE36700, GSE24742, GSE12021, GSE1919, GSE17755.

References

Aletaha et al. (2010) Aletaha D, Neogi T, Silman AJ, Funovits J, Felson DT, Bingham 3rd CO, Birnbaum NS, Burmester GR, Bykerk VP, Cohen MD, Combe B, Costenbader KH, Dougados M, Emery P, Ferraccioli G, Hazes JM, Hobbs K, Huizinga TW, Kavanaugh A, Kay J, Kvien TK, Laing T, Mease P, Menard HA, Moreland LW, Naden RL, Pincus T, Smolen JS, Stanislawska-Biernat E, Symmons D, Tak PP, Upchurch KS, Vencovsky J, Wolfe F, Hawker G. 2010 rheumatoid arthritis classification criteria: an American College of Rheumatology/European league against rheumatism collaborative initiative. Annals of the Rheumatic Diseases. 2010;69:1580–1588. doi: 10.1136/ard.2010.138461. [PubMed] [Cross Ref]
Braeutigam et al. (2014) Braeutigam C, Rago L, Rolke A, Waldmeier L, Christofori G, Winter J. The RNA-binding protein Rbfox2: an essential regulator of EMT-driven alternative splicing and a mediator of cellular invasion. Oncogene. 2014;33:1082–1092. doi: 10.1038/onc.2013.50. [PubMed] [Cross Ref]
Cao (2016) Cao J. The common prescription patterns based on the hierarchical clustering of herb-pairs efficacies. Evidence-Based Complementary and Alternative Medicine. 2016;2016:6373270. doi: 10.1155/2016/6373270. [PMC free article] [PubMed] [Cross Ref]
Carmona-Saez et al. (2007) Carmona-Saez P, Chagoyen M, Tirado F, Carazo JM, Pascual-Montano A. GENECODIS: a web-based tool for finding significant concurrent annotations in gene lists. Genome Biology. 2007;8 doi: 10.1186/gb-2007-8-1-r3. Article gb-2007-8-1-r3. [PMC free article] [PubMed] [Cross Ref]
Catalano (2010) Catalano A. The neuroimmune semaphorin-3A reduces inflammation and progression of experimental autoimmune arthritis. Journal of Immunology. 2010;185:6373–6383. doi: 10.4049/jimmunol.0903527. [PubMed] [Cross Ref]
Cen et al. (2013) Cen YL, Qi ML, Li HG, Su Y, Chen LJ, Lin Y, Chen WQ, Xie XM, Tang LY, Ren ZF. Associations of polymorphisms in the genes of FGFR2, FGF1, and RBFOX2 with breast cancer risk by estrogen/progesterone receptor status. Molecular Carcinogenesis. 2013;52(Suppl 1):E52–59. doi: 10.1002/mc.21979. [PubMed] [Cross Ref]
Chatr-Aryamontri et al. (2015) Chatr-Aryamontri A, Breitkreutz B-J, Oughtred R, Boucher L, Heinicke S, Chen D, Stark C, Breitkreutz A, Kolas N, O’Donnell L. The BioGRID interaction database: 2015 update. Nucleic Acids Research. 2015;43:D470–D478. doi: 10.1093/nar/gku1204. [PMC free article] [PubMed] [Cross Ref]
Edgar, Domrachev & Lash (2002) Edgar R, Domrachev M, Lash AE. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Research. 2002;30:207–210. doi: 10.1093/nar/30.1.207. [PMC free article] [PubMed] [Cross Ref]
Eser & Sahin (2016) Eser B, Sahin N. Evaluation of tool-like receptor-2 and 4 and interleukin-6 gene expressions in Turkish rheumatoid arthritis patients. Clinical Rheumatology. 2016;35:2693–2697. doi: 10.1007/s10067-016-3282-1. [PubMed] [Cross Ref]
Gao et al. (2016) Gao J, Zhou XL, Kong RN, Ji LM, He LL, Zhao DB. microRNA-126 targeting PIK3R2 promotes rheumatoid arthritis synovial fibro-blasts proliferation and resistance to apoptosis by regulating PI3K/AKT pathway. Experimental and Molecular Pathology. 2016;100:192–198. doi: 10.1016/j.yexmp.2015.12.015. [PubMed] [Cross Ref]
Gibofsky (2014) Gibofsky A. Epidemiology, pathophysiology, and diagnosis of rheumatoid arthritis: a synopsis. American Journal of Managed Care. 2014;20:S128–S135. [PubMed]
Gutierrez-Franco et al. (2016) Gutierrez-Franco A, Costa C, Eixarch H, Castillo M, Medina-Rodriguez EM, Bribian A, De Castro F, Montalban X, Espejo C. Differential expression of sema3A and sema7A in a murine model of multiple sclerosis: implications for a therapeutic design. Clinical Immunology. 2016;163:22–33. doi: 10.1016/j.clim.2015.12.005. [PubMed] [Cross Ref]
Homsy et al. (2015) Homsy J, Zaidi S, Shen Y, Ware JS, Samocha KE, Karczewski KJ, DePalma SR, McKean D, Wakimoto H, Gorham J, Jin SC, Deanfield J, Giardini A, Porter GA, Kim Jr R, Bilguvar K, Lopez-Giraldez F, Tikhonova I, Mane S, Romano-Adesman A, Qi H, Vardarajan B, Ma L, Daly M, Roberts AE, Russell MW, Mital S, Newburger JW, Gaynor JW, Breitbart RE, Iossifov I, Ronemus M, Sanders SJ, Kaltman JR, Seidman JG, Brueckner M, Gelb BD, Goldmuntz E, Lifton RP, Seidman CE, Chung WK. De novo mutations in congenital heart disease with neurodevelopmental and other congenital anomalies. Science. 2015;350(6265):1262–1266. doi: 10.1126/science.aac9396. [PMC free article] [PubMed] [Cross Ref]
Kanehisa (2002) Kanehisa M. The KEGG database. Novartis Foundation Symposium. 2002;247:91–101. doi: 10.1002/0470857897.ch8. [PubMed] [Cross Ref]
Kaur et al. (2015) Kaur S, Schwartz AL, Miller TW, Roberts DD. CD47-dependent regulation of H(2)S biosynthesis and signaling in T cells. Methods in Enzymology. 2015;555:145–168. doi: 10.1016/bs.mie.2014.11.023. [PubMed] [Cross Ref]
Khanzadeh et al. (2016) Khanzadeh A, Habibagahi Z, Hosseini A, Amirghofran Z. Investigation of the human FCRL1, 2, and 4 gene expressions in patients with rheumatoid arthritis. Rheumatology International. 2016;36:1149–1156. doi: 10.1007/s00296-016-3495-2. [PubMed] [Cross Ref]
Knevel et al. (2013) Knevel R, Krabben A, Wilson AG, Brouwer E, Leijsma MK, Lindqvist E, De Rooy DP, Daha NA, Van der Linden MP, Tsonaka S, Zhernakova A, Westra HJ, Franke L, Houwing-Duistermaat JJ, Toes RE, Huizinga TW, Saxne T, Van der helm-van Mil AH. A genetic variant in granzyme B is associated with progression of joint destruction in rheumatoid arthritis. Arthtitis and Rheumatism. 2013;65:582–589. doi: 10.1002/art.37808. [PubMed] [Cross Ref]
Lee et al. (2011) Lee HM, Sugino H, Aoki C, Nishimoto N. Underexpression of mitochondrial-DNA encoded ATP synthesis-related genes and DNA repair genes in systemic lupus erythematosus. Arthritis Research & Therapy. 2011;13 doi: 10.1186/ar3317. Article 3317. [PMC free article] [PubMed] [Cross Ref]
Li et al. (2013) Li W, Song AP, Zhao F, Hu YM, Hua M. A novel human TINP1 gene promotes cell proliferation through inhibition of p53 and p21 expression. Oncology Reports. 2013;30:1848–1852. doi: 10.3892/or.2013.2647. [PubMed] [Cross Ref]
Lo, Tsai & Chen (2013) Lo TF, Tsai WC, Chen ST. MicroRNA-21-3p, a berberine-induced miRNA, directly down-regulates human methionine adenosyltransferases 2A and 2B and inhibits hepatoma cell growth. PLOS ONE. 2013;8(9):e75628. doi: 10.1371/journal.pone.0075628. [PMC free article] [PubMed] [Cross Ref]
Lu et al. (2015) Lu X, Zoller EE, Weirauch MT, Wu Z, Namjou B, Williams AH, Ziegler JT, Comeau ME, Marion MC, Glenn SB, Adler A, Shen N, Nath SK, Stevens AM, Freedman BI, Tsao BP, Jacob CO, Kamen DL, Brown EE, Gilkeson GS, Alarcon GS, Reveille JD, Anaya JM, James JA, Sivils KL, Criswell LA, Vila LM, Alarcon-Riquelme ME, Petri M, Scofield RH, Kimberly RP, Ramsey-Goldman R, Joo YB, Choi J, Bae SC, Boackle SA, Graham DC, Vyse TJ, Guthridge JM, Gaffney PM, Langefeld CD, Kelly JA, Greis KD, Kaufman KM, Harley JB, Kottyan LC. Lupus risk variant increases pSTAT1 binding and decreases ETS1 expression. American Journal of Human Genetics. 2015;96:731–739. doi: 10.1016/j.ajhg.2015.03.002. [PubMed] [Cross Ref]
McInnes & Schett (2011) McInnes IB, Schett G. The pathogenesis of rheumatoid arthritis. New England Journal of Medicine. 2011;365:2205–2219. doi: 10.1056/NEJMra1004965. [PubMed] [Cross Ref]
Michelsen et al. (2015) Michelsen B, Fiane R, Diamantopoulos AP, Soldal DM, Hansen IJ, Sokka T, Kavanaugh A, Haugeberg G. A comparison of disease burden in rheumatoid arthritis, psoriatic arthritis and axial spondyloarthritis. PLOS ONE. 2015;10(4):e0123582 doi: 10.1371/journal.pone.0123582. [PMC free article] [PubMed] [Cross Ref]
Okada et al. (2016) Okada Y, Suzuki A, Ikari K, Terao C, Kochi Y, Ohmura K, Higasa K, Akiyama M, Ashikawa K, Kanai M, Hirata J, Suita N, Teo YY, Xu H, Bae SC, Takahashi A, Momozawa Y, Matsuda K, Momohara S, Taniguchi A, Yamada R, Mimori T, Kubo M, Brown MA, Raychaudhuri S, Matsuda F, Yamanaka H, Kamatani Y, Yamamoto K. Contribution of a non-classical HLA gene, HLA-DOA, to the risk of rheumatoid arthritis. American Journal of Human Genetics. 2016;99:366–374. doi: 10.1016/j.ajhg.2016.06.019. [PubMed] [Cross Ref]
Panayiotidis et al. (2006) Panayiotidis MI, Stabler SP, Ahmad A, Pappa A, Legros Jr LH, Hernandez-Saavedra D, Schneider BK, Allen RH, Vasiliou V, McCord JM, Kotb M, White CW. Activation of a novel isoform of methionine adenosyl transferase 2A and increased S-adenosylmethionine turnover in lung epithelial cells exposed to hyperoxia. Free Radical Biology and Medicine. 2006;40:348–358. doi: 10.1016/j.freeradbiomed.2005.09.004. [PubMed] [Cross Ref]
Pearson & Jones (2016) Pearson MJ, Jones SW. Review: long non-coding RNAs in the regulation of inflammatory pathways in rheumatoid arthritis and osteoarthritis. Arthritis & Rheumatology. 2016;68:2575–2583. doi: 10.1002/art.39759. [PMC free article] [PubMed] [Cross Ref]
Pincus & Callahan (1986) Pincus T, Callahan LF. Taking mortality in rheumatoid arthritis seriously—predictive markers, socioeconomic status and comorbidity. Journal of Rheumatology. 1986;13:841–845. [PubMed]
Redlich et al. (2001) Redlich K, Kiener HP, Schett G, Tohidast-Akrad M, Selzer E, Radda I, Stummvoll GH, Steiner CW, Groger M, Bitzan P, Zenz P, Smolen JS, Steiner G. Overexpression of transcription factor Ets-1 in rheumatoid arthritis synovial membrane: regulation of expression and activation by interleukin-1 and tumor necrosis factor alpha. Arthtitis and Rheumatism. 2001;44:266–274. doi: 10.1002/1529-0131(200102)44:2<266::AID-ANR43>3.0.CO;2-G. [PubMed] [Cross Ref]
Reiner-Benaim (2007) Reiner-Benaim A. FDR control by the BH procedure for two-sided correlated tests with implications to gene expression data analysis. Biometrical Journal. 2007;49:107–126. doi: 10.1002/bimj.200510313. [PubMed] [Cross Ref]
Schett, Zwerina & Firestein (2008) Schett G, Zwerina J, Firestein G. The p38 mitogen-activated protein kinase (MAPK) pathway in rheumatoid arthritis. Annals of the Rheumatic Diseases. 2008;67:909–916. doi: 10.1136/ard.2007.074278. [PMC free article] [PubMed] [Cross Ref]
Shannon et al. (2003) Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Research. 2003;13:2498–2504. doi: 10.1101/gr.1239303. [PubMed] [Cross Ref]
Shi et al. (2006) Shi L, Reid LH, Jones WD, Shippy R, Warrington JA, Baker SC, Collins PJ, De Longueville F, Kawasaki ES, Lee KY. The microarray quality control (MAQC) project shows inter-and intraplatform reproducibility of gene expression measurements. Nature Biotechnology. 2006;24:1151–1161. doi: 10.1038/nbt1239. [PMC free article] [PubMed] [Cross Ref]
Smolen, Aletaha & McInnes (2016) Smolen JS, Aletaha D, McInnes IB. Rheumatoid arthritis. Lancet. 2016;388:2023–2038. doi: 10.1016/S0140-6736(16)30173-8. [PubMed] [Cross Ref]
Smyth (2005) Smyth GK. Bioinformatics and computational biology solutions using R and bioconductor. New York: Springer; 2005. Limma: linear models for microarray data; pp. 397–420.
Tomasi et al. (2013) Tomasi ML, Ryoo M, Skay A, Tomasi I, Giordano P, Mato JM, Lu SC. Polyamine and methionine adenosyltransferase 2A crosstalk in human colon and liver cancer. Experimental Cell Research. 2013;319:1902–1911. doi: 10.1016/j.yexcr.2013.04.005. [PMC free article] [PubMed] [Cross Ref]
Torto-Alalibo et al. (2014) Torto-Alalibo T, Purwantini E, Lomax J, Setubal JC, Mukhopadhyay B, Tyler BM. Genetic resources for advanced biofuel production described with the Gene Ontology. Frontiers in Microbiology. 2014;5 doi: 10.3389/fmicb.2014.00528. Article 528. [PMC free article] [PubMed] [Cross Ref]
Vallejo et al. (2003) Vallejo AN, Yang H, Klimiuk PA, Weyand CM, Goronzy JJ. Synoviocyte-mediated expansion of inflammatory T cells in rheumatoid synovitis is dependent on CD47-thrombospondin 1 interaction. Journal of Immunology. 2003;171:1732–1740. doi: 10.4049/jimmunol.171.4.1732. [PubMed] [Cross Ref]
Venables et al. (2013a) Venables JP, Brosseau JP, Gadea G, Klinck R, Prinos P, Beaulieu JF, Lapointe E, Durand M, Thibault P, Tremblay K, Rousset F, Tazi J, Abou Elela S, Chabot B. RBFOX2 is an important regulator of mesenchymal tissue-specific splicing in both normal and cancer tissues. Molecular and Cellular Biology. 2013a;33:396–405. doi: 10.1128/MCB.01174-12. [PMC free article] [PubMed] [Cross Ref]
Venables et al. (2013b) Venables JP, Lapasset L, Gadea G, Fort P, Klinck R, Irimia M, Vignal E, Thibault P, Prinos P, Chabot B, Abou Elela S, Roux P, Lemaitre JM, Tazi J. MBNL1 and RBFOX2 cooperate to establish a splicing programme involved in pluripotent stem cell differentiation. Nature Communications. 2013b;4 doi: 10.1038/ncomms3480. Article 3480. [PubMed] [Cross Ref]
Wang et al. (2014) Wang X, Guo X, Yu W, Li C, Gui Y, Cai Z. Expression of methionine adenosyltransferase 2A in renal cell carcinomas and potential mechanism for kidney carcinogenesis. BMC Cancer. 2014;14:196. doi: 10.1186/1471-2407-14-196. [PMC free article] [PubMed] [Cross Ref]
Weisbart et al. (2013) Weisbart RH, Chan G, Li E, Farmani N, Heinze E, Rubell A, Nishimura RN, Colburn K. BRAF splice variants in rheumatoid arthritis synovial fibroblasts activate MAPK through CRAF. Molecular Immunology. 2013;55:247–252. doi: 10.1016/j.molimm.2013.02.001. [PubMed] [Cross Ref]
Wernert et al. (2002) Wernert N, Justen HP, Rothe M, Behrens P, Dreschers S, Neuhaus T, Florin A, Sachinidis A, Vetter H, Ko Y. The Ets 1 transcription factor is upregulated during inflammatory angiogenesis in rheumatoid arthritis. Journal of Molecular Medicine. 2002;80:258–266. doi: 10.1007/s00109-001-0316-0. [PubMed] [Cross Ref]
Wingender et al. (1996) Wingender E, Dietze P, Karas H, Knuppel R. TRANSFAC: a database on transcription factors and their DNA binding sites. Nucleic Acids Research. 1996;24:238–241. doi: 10.1093/nar/24.1.238. [PMC free article] [PubMed] [Cross Ref]
Xiao et al. (2014) Xiao Y, Sun M, Zhan Z, Ye Y, Huang M, Zou Y, Liang L, Yang X, Xu H. Increased phosphorylation of ezrin is associated with the migration and invasion of fibroblast-like synoviocytes from patients with rheumatoid arthritis. Rheumatology. 2014;53:1291–1300. doi: 10.1093/rheumatology/keu013. [PubMed] [Cross Ref]
Zhang et al. (2010) Zhang H, Ma X, Shi T, Song Q, Zhao H, Ma D. NSA2, a novel nucleolus protein regulates cell proliferation and cell cycle. Biochemical and Biophysical Research Communications. 2010;391:651–658. doi: 10.1016/j.bbrc.2009.11.114. [PubMed] [Cross Ref]
Zhang et al. (2013) Zhang T, Zheng Z, Liu Y, Zhang J, Zhao Y, Liu Y, Zhu H, Zhao G, Liang J, Li Q, Xu H. Overexpression of methionine adenosyltransferase II alpha (MAT2A) in gastric cancer and induction of cell cycle arrest and apoptosis in SGC-7901 cells by shRNA-mediated silencing of MAT2A gene. Acta Histochemica. 2013;115:48–55. doi: 10.1016/j.acthis.2012.03.006. [PubMed] [Cross Ref]

Articles from PeerJ are provided here courtesy of PeerJ, Inc