Runrun Zhang1,2, Yehua Jin1,3, Cen Chang1,3, Lingxia Xu1,3, Yanqin Bian4,5, Yu Shen4,5, Yang Sun4,5, Songtao Sun6, Steven J Schrodi7, Shicheng Guo7, Dongyi He1,5. 1. Department of Rheumatology, Shanghai Guanghua Hospital, Shanghai University of Traditional Chinese Medicine, Shanghai, China. 2. Department of Rheumatology, The Second Affiliated Hospital of Shandong University of Traditional Chinese Medicine, Jinan, China. 3. Guanghua Clinical Medical College, Shanghai University of Traditional Chinese Medicine, Shanghai, China. 4. Guanghua Integrative Medicine Hospital, Shanghai, China. 5. Arthritis Institute of Integrated Traditional and Western Medicine, Shanghai Chinese Medicine Research Institute, Shanghai, China. 6. Department of Orthopedics, Shanghai Guanghua Hospital, Shanghai University of Traditional Chinese Medicine, Shanghai, China. 7. Department of Medical Genetics, School of Medicine and Public Health, University of Wisconsin-Madison, Madison, WI, United States.
Abstract
Purpose: This study aimed to provide a comprehensive understanding of the genome-wide expression patterns in the synovial tissue samples of patients with rheumatoid arthritis (RA) to investigate the potential mechanisms regulating RA occurrence and development. Methods: Transcription profiles of the synovial tissue samples from nine patients with RA and 15 patients with osteoarthritis (OA) (control) from the East Asian population were generated using RNA sequencing (RNA-seq). Gene set enrichment analysis (GSEA) was used to analyze all the detected genes and the differentially expressed genes (DEGs) were identified using DESeq. To further analyze the DEGs, the Gene Ontology (GO) functional enrichment and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed. The protein-protein interaction (PPI) network of the DEGs was constructed using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and the hub genes were identified by topology clustering with the Molecular Complex Detection (MCODE)-Cytoscape. The most important hub genes were validated using quantitative real-time PCR (qRT-PCR). Results: Of the 17,736 genes detected, 851 genes were identified as the DEGs (474 upregulated and 377 downregulated genes) using the false discovery rate (FDR) approach. GSEA revealed that the significantly enriched gene sets that positively correlated with RA were CD40 signaling overactivation, Th1 cytotoxic module, overactivation of the immune response, adaptive immune response, effective vs. memory CD8+ T cells (upregulated), and naïve vs. effective CD8+ T cells (downregulated). Biological process enrichment analysis showed that the DEGs were significantly enriched for signal transduction (P = 3.01 × 10-6), immune response (P = 1.65 × 10-24), and inflammatory response (P = 5.76 × 10-10). Molecule function enrichment analysis revealed that the DEGs were enriched in calcium ion binding (P = 1.26 × 10-5), receptor binding (P = 1.26 × 10-5), and cytokine activity (P = 2.01 × 10-3). Cellular component enrichment analysis revealed that the DEGs were significantly enriched in the plasma membrane (P = 1.91 × 10-31), an integral component of the membrane (P = 7.39 × 10-13), and extracellular region (P = 7.63 × 10-11). The KEGG pathway analysis showed that the DEGs were enriched in the cytokine-cytokine receptor interaction (P = 3.05 × 10-17), chemokine signaling (P = 3.50 × 10-7), T-cell receptor signaling (P = 5.17 × 10-4), and RA (P = 5.17 × 10-4) pathways. We confirmed that RA was correlated with the upregulation of the PPI network hub genes, such as CXCL13, CXCL6, CCR5, CXCR5, CCR2, CXCL3, and CXCL10, and the downregulation of the PPI network hub gene such as SSTR1. Conclusion: This study identified and validated the DEGs in the synovial tissue samples of patients with RA, which highlighted the activity of a subset of chemokine genes, thereby providing novel insights into the molecular mechanisms of RA pathogenesis and identifying potential diagnostic and therapeutic targets for RA.
Purpose: This study aimed to provide a comprehensive understanding of the genome-wide expression patterns in the synovial tissue samples of patients with rheumatoid arthritis (RA) to investigate the potential mechanisms regulating RA occurrence and development. Methods: Transcription profiles of the synovial tissue samples from nine patients with RA and 15 patients with osteoarthritis (OA) (control) from the East Asian population were generated using RNA sequencing (RNA-seq). Gene set enrichment analysis (GSEA) was used to analyze all the detected genes and the differentially expressed genes (DEGs) were identified using DESeq. To further analyze the DEGs, the Gene Ontology (GO) functional enrichment and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed. The protein-protein interaction (PPI) network of the DEGs was constructed using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and the hub genes were identified by topology clustering with the Molecular Complex Detection (MCODE)-Cytoscape. The most important hub genes were validated using quantitative real-time PCR (qRT-PCR). Results: Of the 17,736 genes detected, 851 genes were identified as the DEGs (474 upregulated and 377 downregulated genes) using the false discovery rate (FDR) approach. GSEA revealed that the significantly enriched gene sets that positively correlated with RA were CD40 signaling overactivation, Th1 cytotoxic module, overactivation of the immune response, adaptive immune response, effective vs. memory CD8+ T cells (upregulated), and naïve vs. effective CD8+ T cells (downregulated). Biological process enrichment analysis showed that the DEGs were significantly enriched for signal transduction (P = 3.01 × 10-6), immune response (P = 1.65 × 10-24), and inflammatory response (P = 5.76 × 10-10). Molecule function enrichment analysis revealed that the DEGs were enriched in calcium ion binding (P = 1.26 × 10-5), receptor binding (P = 1.26 × 10-5), and cytokine activity (P = 2.01 × 10-3). Cellular component enrichment analysis revealed that the DEGs were significantly enriched in the plasma membrane (P = 1.91 × 10-31), an integral component of the membrane (P = 7.39 × 10-13), and extracellular region (P = 7.63 × 10-11). The KEGG pathway analysis showed that the DEGs were enriched in the cytokine-cytokine receptor interaction (P = 3.05 × 10-17), chemokine signaling (P = 3.50 × 10-7), T-cell receptor signaling (P = 5.17 × 10-4), and RA (P = 5.17 × 10-4) pathways. We confirmed that RA was correlated with the upregulation of the PPI network hub genes, such as CXCL13, CXCL6, CCR5, CXCR5, CCR2, CXCL3, and CXCL10, and the downregulation of the PPI network hub gene such as SSTR1. Conclusion: This study identified and validated the DEGs in the synovial tissue samples of patients with RA, which highlighted the activity of a subset of chemokine genes, thereby providing novel insights into the molecular mechanisms of RA pathogenesis and identifying potential diagnostic and therapeutic targets for RA.
Authors: Paul Shannon; Andrew Markiel; Owen Ozier; Nitin S Baliga; Jonathan T Wang; Daniel Ramage; Nada Amin; Benno Schwikowski; Trey Ideker Journal: Genome Res Date: 2003-11 Impact factor: 9.043
Authors: Vamsi K Mootha; Cecilia M Lindgren; Karl-Fredrik Eriksson; Aravind Subramanian; Smita Sihag; Joseph Lehar; Pere Puigserver; Emma Carlsson; Martin Ridderstråle; Esa Laurila; Nicholas Houstis; Mark J Daly; Nick Patterson; Jill P Mesirov; Todd R Golub; Pablo Tamayo; Bruce Spiegelman; Eric S Lander; Joel N Hirschhorn; David Altshuler; Leif C Groop Journal: Nat Genet Date: 2003-07 Impact factor: 38.330
Authors: Cole Trapnell; Brian A Williams; Geo Pertea; Ali Mortazavi; Gordon Kwan; Marijke J van Baren; Steven L Salzberg; Barbara J Wold; Lior Pachter Journal: Nat Biotechnol Date: 2010-05-02 Impact factor: 54.908