The goal of this study was to identify potential predictive biomarkers for the therapeutic effect of infliximab (IFX) in Rheumatoid arthritis (RA) and explore the potential molecular mechanism of nonresponse to IFX treatment to achieve individualized treatment of RA.
Differential gene expression between IFX responders and nonresponders in the GSE58795 and GSE78068 datasets was identified. Coexpression analysis was used to identify the modules associated with nonresponse to IFX therapy for RA, and enrichment analysis was conducted on module genes. Least absolute shrink and selection operator (LASSO) regression was used to develop a gene signature for predicting the therapeutic effect of IFX in RA, and the area under the receiver operating characteristic curve (AUC) was used to evaluate the predictive value of the signature. Correlation analysis and single-sample gene set enrichment analysis (ssGSEA) were used to explore the potential role of the hub genes. Experimental validation was conducted in synovial tissue and RA fibroblast-like synoviocytes (RA-FLSs).
A total of 46 common genes were obtained among the two datasets. The yellow-green module was identified as the key module associated with nonresponse to IFX therapy for RA. We identified a 25-gene signature in GSE78068, and the AUC for the signature was 0.831 in the internal validation set and 0.924 in the GSE58795 dataset(external validation set). 
The 25-gene signature may have potential predictive value for the therapeutic effect of IFX in RA at the beginning of IFX treatment, and autophagy may be involved in nonresponse to IFX treatment. In particular, DERL1 may be associated with the regulation of autophagy.