Identification of baseline gene expression signatures predicting therapeutic responses to three biologic agents in rheumatoid arthritis: a retrospective observational study

Seiji Nakamura, Katsuya Suzuki, Hiroshi Iijima, Yuko Hata, Chun Ren Lim, Yohei Ishizawa, Hideto Kameda, Koichi Amano, Kenichi Matsubara, Ryo Matoba, Tsutomu Takeuchi

According to EULAR recommendations, biologic DMARDs (bDMARDs) such as tumor necrosis factor inhibitor, tocilizumab (TCZ), and abatacept (ABT) are in parallel when prescribing to rheumatoid arthritis (RA) patients who have shown insufficient response to conventional synthetic DMARDs. However, most prediction studies of therapeutic response to bDMARDs using gene expression profiles were focused on a single bDMARD, and consideration of the results from the perspective of RA pathophysiology was insufficient. The aim of this study was to identify the specific molecular biological features predicting the therapeutic outcomes of three bDMARDs (infliximab [IFX], TCZ, and ABT) by studying blood gene expression signatures of patients before biologic treatment in a unified test platform. RA patients who responded inadequately to methotrexate and were later commenced on any one of IFX (n = 140), TCZ (n = 38), or ABT (n = 31) as their first biologic between May 2007 and November 2011 were enrolled. Whole-blood gene expression data were obtained before biologic administration. Patients were categorized into remission (REM) and nonremission (NON-REM) groups according to CDAI at 6 months of biologic therapy. We employed Gene Set Enrichment Analysis (GSEA) to identify functional gene sets differentially expressed between these two groups for each biologic. Then, we compiled "signature scores" for these gene sets, and the prediction performances were assessed. GSEA showed that inflammasome genes were significantly upregulated with IFX in the NON-REM group compared with the REM group. With TCZ in the REM group, B-cell-specifically expressed genes were upregulated. RNA elongation, apoptosis-related, and NK-cell-specifically expressed genes were upregulated with ABT in the NON-REM group. Logistic regression analyses showed that "signature scores" of inflammasomes, B-cell-specifically expressed, and NK-cell-specifically expressed genes were significant, independently predictive factors for treatment outcome with IFX, TCZ, and ABT, respectively. The AUCs of ROC curves of these signature scores were 0.637, 0.796, and 0.768 for IFX, TCZ, and ABT, respectively. We have identified original gene expression predictive signatures uniquely underlying the therapeutic effects of IFX, TCZ, and ABT. This is, to our knowledge, the first attempt to predict therapeutic effects of three drugs concomitantly using a unified gene expression test platform.