Differential Expression Profiles of the Transcriptome and miRNA Interactome in Synovial Fibroblasts of Rheumatoid Arthritis Revealed by Next Generation Sequencing

Abstract

Using next-generation sequencing to decipher the molecular mechanisms underlying aberrant rheumatoid arthritis synovial fibroblasts (RASF) activation, we performed transcriptome-wide RNA-seq and small RNA-seq on synovial fibroblasts from rheumatoid arthritis (RA) subject and normal donor. Differential expression of mRNA and miRNA was integrated with interaction analysis, functional annotation, regulatory network mapping and experimentally verified miRNA–target interaction data, further validated with microarray expression profiles. In this study, 3049 upregulated mRNA and 3552 downregulated mRNA, together with 50 upregulated miRNA and 35 downregulated miRNA in RASF were identified. Interaction analysis highlighted contribution of miRNA to altered transcriptome. Functional annotation revealed metabolic deregulation and oncogenic signatures of RASF. Regulatory network mapping identified downregulated FOXO1 as master transcription factor resulting in altered transcriptome of RASF. Differential expression in three miRNA and corresponding targets (hsa-miR-31-5p:WASF3, hsa-miR-132-3p:RB1, hsa-miR-29c-3p:COL1A1) were also validated. The interactions of these three miRNA–target genes were experimentally validated with past literature. Our transcriptomic and miRNA interactomic investigation identified gene signatures associated with RASF and revealed the involvement of transcription factors and miRNA in an altered transcriptome. These findings help facilitate our understanding of RA with the hope of serving as a springboard for further discoveries relating to the disease.