Identification of dysregulated genes in rheumatoid arthritis based on bioinformatics analysis

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.

[1]  C. Werning.  [Rheumatoid arthritis]. , 1983, Medizinische Monatsschrift fur Pharmazeuten.

[2]  T. Pincus,et al.  Taking mortality in rheumatoid arthritis seriously--predictive markers, socioeconomic status and comorbidity. , 1986, The Journal of rheumatology.

[3]  Holger Karas,et al.  TRANSFAC: a database on transcription factors and their DNA binding sites , 1996, Nucleic Acids Res..

[4]  Y. Ko,et al.  The Ets 1 transcription factor is upregulated during inflammatory angiogenesis in rheumatoid arthritis , 2002, Journal of Molecular Medicine.

[5]  M. Gröger,et al.  Overexpression of transcription factor Ets-1 in rheumatoid arthritis synovial membrane: regulation of expression and activation by interleukin-1 and tumor necrosis factor alpha. , 2001, Arthritis and rheumatism.

[6]  Alex E. Lash,et al.  Gene Expression Omnibus: NCBI gene expression and hybridization array data repository , 2002, Nucleic Acids Res..

[7]  Minoru Kanehisa,et al.  The KEGG database. , 2002, Novartis Foundation symposium.

[8]  P. Shannon,et al.  Cytoscape: a software environment for integrated models of biomolecular interaction networks. , 2003, Genome research.

[9]  A. Vallejo,et al.  Synoviocyte-Mediated Expansion of Inflammatory T Cells in Rheumatoid Synovitis Is Dependent on CD47-Thrombospondin 1 Interaction1 , 2003, The Journal of Immunology.

[10]  Gordon K. Smyth,et al.  limma: Linear Models for Microarray Data , 2005 .

[11]  Hanlee P. Ji,et al.  The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements. , 2006, Nature biotechnology.

[12]  V. Vasiliou,et al.  Activation of a novel isoform of methionine adenosyl transferase 2A and increased S-adenosylmethionine turnover in lung epithelial cells exposed to hyperoxia. , 2006, Free radical biology & medicine.

[13]  Maqc Consortium.  The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements , 2006, Nature Biotechnology.

[14]  J. Carazo,et al.  GENECODIS: a web-based tool for finding significant concurrent annotations in gene lists , 2007, Genome Biology.

[15]  A. Reiner-Benaim.  FDR Control by the BH Procedure for Two‐Sided Correlated Tests with Implications to Gene Expression Data Analysis , 2007, Biometrical journal. Biometrische Zeitschrift.

[16]  A. Koch.  The pathogenesis of rheumatoid arthritis. , 2007, American journal of orthopedics.

[17]  J. Davis.  Bioinformatics and Computational Biology Solutions Using R and Bioconductor , 2007 .

[18]  G. Firestein,et al.  The p38 mitogen-activated protein kinase (MAPK) pathway in rheumatoid arthritis , 2007, Annals of the rheumatic diseases.

[19]  A. Silman,et al.  Rheumatoid arthritis classifi cation criteria : an American College of Rheumatology / European League Against Rheumatism collaborative initiative , 2010 .

[20]  Hongshan Zhao,et al.  NSA2, a novel nucleolus protein regulates cell proliferation and cell cycle. , 2010, Biochemical and biophysical research communications.

[21]  A. Catalano.  The Neuroimmune Semaphorin-3A Reduces Inflammation and Progression of Experimental Autoimmune Arthritis , 2010, The Journal of Immunology.

[22]  N. Nishimoto,et al.  Underexpression of mitochondrial-DNA encoded ATP synthesis-related genes and DNA repair genes in systemic lupus erythematosus , 2011, Arthritis research & therapy.

[23]  F. Rousset,et al.  RBFOX2 Is an Important Regulator of Mesenchymal Tissue-Specific Splicing in both Normal and Cancer Tissues , 2012, Molecular and Cellular Biology.

[24]  R. Weisbart,et al.  BRAF splice variants in rheumatoid arthritis synovial fibroblasts activate MAPK through CRAF. , 2013, Molecular immunology.

[25]  J. Tazi,et al.  MBNL1 and RBFOX2 cooperate to establish a splicing programme involved in pluripotent stem cell differentiation , 2013, Nature Communications.

[26]  Shelly C. Lu,et al.  Polyamine and methionine adenosyltransferase 2A crosstalk in human colon and liver cancer. , 2013, Experimental cell research.

[27]  Wei-Qing Chen,et al.  Associations of polymorphisms in the genes of FGFR2, FGF1, and RBFOX2 with breast cancer risk by estrogen/progesterone receptor status , 2013, Molecular carcinogenesis.

[28]  M. Hua,et al.  A novel human TINP1 gene promotes cell proliferation through inhibition of p53 and p21 expression. , 2013, Oncology reports.

[29]  Haitao Zhu,et al.  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. , 2013, Acta histochemica.

[30]  Z. Cai,et al.  Expression of methionine adenosyltransferase 2A in renal cell carcinomas and potential mechanism for kidney carcinogenesis , 2014, BMC Cancer.

[31]  A. Zhernakova,et al.  A genetic variant in granzyme B is associated with progression of joint destruction in rheumatoid arthritis. , 2013, Arthritis and rheumatism.

[32]  Shui-Tein Chen,et al.  MicroRNA-21-3p, a Berberine-Induced miRNA, Directly Down-Regulates Human Methionine Adenosyltransferases 2A and 2B and Inhibits Hepatoma Cell Growth , 2013, PloS one.

[33]  A. Gibofsky.  Epidemiology, pathophysiology, and diagnosis of rheumatoid arthritis: A Synopsis. , 2014, The American journal of managed care.

[34]  G. Christofori,et al.  The RNA-binding protein Rbfox2: an essential regulator of EMT-driven alternative splicing and a mediator of cellular invasion , 2014, Oncogene.

[35]  João C. Setubal,et al.  Genetic resources for advanced biofuel production described with the Gene Ontology , 2014, Front. Microbiol..

[36]  Z. Zhan,et al.  Increased phosphorylation of ezrin is associated with the migration and invasion of fibroblast-like synoviocytes from patients with rheumatoid arthritis. , 2014, Rheumatology.

[37]  G. Haugeberg,et al.  A Comparison of Disease Burden in Rheumatoid Arthritis, Psoriatic Arthritis and Axial Spondyloarthritis , 2015, PloS one.

[38]  Christie S. Chang,et al.  The BioGRID interaction database: 2015 update , 2014, Nucleic Acids Res..

[39]  Stephan J Sanders,et al.  De novo mutations in congenital heart disease with neurodevelopmental and other congenital anomalies , 2015, Science.

[40]  P. Gaffney,et al.  Lupus Risk Variant Increases pSTAT1 Binding and Decreases ETS1 Expression. , 2015, American journal of human genetics.

[41]  D. Roberts,et al.  CD47-dependent regulation of H₂S biosynthesis and signaling in T cells. , 2015, Methods in enzymology.

[42]  R. Kong,et al.  microRNA-126 targeting PIK3R2 promotes rheumatoid arthritis synovial fibro-blasts proliferation and resistance to apoptosis by regulating PI3K/AKT pathway. , 2016, Experimental and molecular pathology.

[43]  B. Eser,et al.  Evaluation of tool-like receptor-2 and 4 and interleukin-6 gene expressions in Turkish rheumatoid arthritis patients , 2016, Clinical Rheumatology.

[44]  X. Montalban,et al.  Differential expression of sema3A and sema7A in a murine model of multiple sclerosis: Implications for a therapeutic design. , 2016, Clinical immunology.

[45]  Simon W. Jones,et al.  Review: Long Noncoding RNAs in the Regulation of Inflammatory Pathways in Rheumatoid Arthritis and Osteoarthritis , 2016, Arthritis & rheumatology.

[46]  Jia Cao.  The Common Prescription Patterns Based on the Hierarchical Clustering of Herb-Pairs Efficacies , 2016, Evidence-based complementary and alternative medicine : eCAM.

[47]  M. Kanai,et al.  Contribution of a Non-classical HLA Gene, HLA-DOA, to the Risk of Rheumatoid Arthritis. , 2016, American journal of human genetics.

[48]  A. Hosseini,et al.  Investigation of the human FCRL1, 2, and 4 gene expressions in patients with rheumatoid arthritis , 2016, Rheumatology International.