Rheumatoid arthritis (RA) is therapeutically challenging due to patient heterogeneity and variability. Herein we describe a novel integration of RA synovial genome-scale transcriptomic profiling of different patient cohorts that can be used to provide predictive insights on drug responses.
A normalized compendium consisting of 256 RA synovial samples that cover an intersection of 11,769 genes from 11 datasets was build and compared with similar datasets derived from OA patients and healthy controls. Differentially expression genes (DEGs) that were identified in three independent methods were fed into functional network analysis, with subsequent grouping of the samples based on a non-negative matrix factorization method. RA-relevant pathway activation scores and four machine learning classification techniques supported the generation of a predictive model of patient treatment response. We identified 876 up-regulated DEGs including 24 known genetic risk factors and 8 drug targets.
