Molecular profiling of rheumatoid arthritis patients reveals an association between innate and adaptive cell populations and response to anti-tumor necrosis factor

Patients with High Baseline Neutrophil-to-Lymphocyte Ratio Exhibit Better Response to Filgotinib as Treatment for Rheumatoid Arthritis

INTRODUCTION

High baseline neutrophil-to-lymphocyte ratio (NLR) in rheumatoid arthritis (RA) has been associated with positive responses to biologic tumor necrosis factor inhibition and negative responses to conventional synthetic disease-modifying antirheumatic drug (csDMARD) triple therapy. Datasets from three randomized clinical trials in patients with RA were used to test the hypothesis that baseline NLR is associated with improved clinical response to filgotinib in methotrexate (MTX)-naïve or MTX-experienced RA populations.

METHODS

Patients from FINCH 1 (inadequate response to MTX, MTX-IR; NCT02889796), FINCH 2 (inadequate response to biologic DMARDs; NCT02873936), and FINCH 3 (MTX-naïve; NCT02886728) were classified as baseline NLR-High or baseline NLR-Low based on a previously published cut point of 2.7. In total, 3365 patients were included across the three studies. Differences in clinical outcomes and patient-reported outcomes (PROs) were determined using linear-regression models.

RESULTS

Control-arm patients (placebo + MTX/placebo + csDMARD) classified as NLR-High exhibited worse continuous clinical and PRO responses at week 12 across clinical trials compared to NLR-Low patients. In contrast, NLR-High patients who received FIL 200 mg + MTX/csDMARD exhibited consistently better responses after 12 weeks compared to NLR-Low patients across clinical trials, clinical endpoints, and PROs. These trends were most prominent among the MTX-IR population.

CONCLUSION

The 2.7 baseline NLR cut point could be used to enrich for patients most likely to benefit from the addition of filgotinib to background MTX/csDMARD. Use of baseline NLR as part of therapeutic decision-making would not require additional diagnostics and could contribute to improved outcomes for patients with RA.

TRIAL REGISTRATION

Clinicaltrials.gov: NCT02889796; NCT02873936; NCT02886728.

Deciphering the molecular landscape of rheumatoid arthritis offers new insights into the stratified treatment for the condition

Background

For Rheumatoid Arthritis (RA), a long-term chronic illness, it is essential to identify and describe patient subtypes with comparable goal status and molecular biomarkers. This study aims to develop and validate a new subtyping scheme that integrates genome-scale transcriptomic profiles of RA peripheral blood genes, providing a fresh perspective for stratified treatments.

Methods

We utilized independent microarray datasets of RA peripheral blood mononuclear cells (PBMCs). Up-regulated differentially expressed genes (DEGs) were subjected to functional enrichment analysis. Unsupervised cluster analysis was then employed to identify RA peripheral blood gene expression-driven subtypes. We defined three distinct clustering subtypes based on the identified 404 up-regulated DEGs.

Results

Subtype A, named NE-driving, was enriched in pathways related to neutrophil activation and responses to bacteria. Subtype B, termed interferon-driving (IFN-driving), exhibited abundant B cells and showed increased expression of transcripts involved in IFN signaling and defense responses to viruses. In Subtype C, an enrichment of CD8+ T-cells was found, ultimately defining it as CD8+ T-cells-driving. The RA subtyping scheme was validated using the XGBoost machine learning algorithm. We also evaluated the therapeutic outcomes of biological disease-modifying anti-rheumatic drugs.

Conclusions

The findings provide valuable insights for deep stratification, enabling the design of molecular diagnosis and serving as a reference for stratified therapy in RA patients in the future.

Transcriptomic network analysis reveals key drivers of response to anti-TNF biologics in patients with rheumatoid arthritis

OBJECTIVE

Anti-TNF biologics have been widely used to ameliorate disease activity in patients with RA. However, a large fraction of patients show a poor response to these agents. Moreover, no clinically applicable predictive biomarkers have been established. This study aimed to identify response-associated biomarkers using longitudinal transcriptomic data in two independent RA cohorts.

METHODS

RNA sequencing data from peripheral blood cell samples of Korean and Caucasian RA cohorts before and after initial treatment with anti-TNF biologics were analysed to assess treatment-induced expression changes that differed between highly reliable excellent responders and null responders. Weighted correlation network, immune cell composition, and key driver analyses were performed to understand response-associated transcriptomic networks and cell types and their correlation with disease activity indices.

RESULTS

In total, 305 response-associated genes showed significantly different treatment-induced expression changes between excellent and null responders. Co-expression network construction and subsequent key driver analysis revealed that 41 response-associated genes played a crucial role as key drivers of transcriptomic alteration in four response-associated networks involved in various immune pathways: type I IFN signalling, myeloid leucocyte activation, B cell activation, and NK cell/lymphocyte-mediated cytotoxicity. Transcriptomic response scores that we developed to estimate the individual-level degree of expression changes in the response-associated key driver genes were significantly correlated with the changes in clinical indices in independent patients with moderate or ambiguous response outcomes.

CONCLUSION

This study provides response-specific treatment-induced transcriptomic signatures by comparing the transcriptomic landscape between patients with excellent and null responses to anti-TNF drugs at both gene and network levels.

A proteomics study of rheumatoid arthritis patients on etanercept identifies putative biomarkers associated with clinical outcome measures

OBJECTIVES

Biologic DMARDs (bDMARDs) are widely used in patients with RA, but response to bDMARDs is heterogeneous. The objective of this work was to identify pretreatment proteomic biomarkers associated with RA clinical outcome measures in patients starting bDMARDs.

METHODS

Sequential window acquisition of all theoretical fragment ion spectra mass spectrometry (SWATH-MS) was used to generate spectral maps of sera from patients with RA before and after 3 months of treatment with the bDMARD etanercept. Protein levels were regressed against RA clinical outcome measures, i.e. 28-joint DAS (DAS28) and its subcomponents and DAS28 <2.6 (i.e. remission). The proteins with the strongest evidence for association were analysed in an independent, replication dataset. Finally, subnetwork analysis was carried out using the Disease Module Detection algorithm and biological plausibility of identified proteins was assessed by enrichment analysis.

RESULTS

A total of 180 patients with RA were included in the discovery dataset and 58 in the validation dataset from a UK-based prospective multicentre study. Ten individual proteins were found to be significantly associated with RA clinical outcome measures. The association of T-complex protein 1 subunit η with DAS28 remission was replicated in an independent cohort. Subnetwork analysis of the 10 proteins from the regression analysis identified the ontological theme, with the strongest associations being with acute phase and acute inflammatory responses.

CONCLUSION

This longitudinal study of 180 patients with RA commencing etanercept has identified several putative protein biomarkers of treatment response to this drug, one of which was replicated in an independent cohort.

The association between neutrophil-to-lymphocyte ratio and disease activity in rheumatoid arthritis

The inflammatory response is responsible for the promotion of pannus development over the joint, which is the primary factor in joint injury in rheumatoid arthritis (RA). More in-depth investigations have been conducted in recent years leading to a greater understanding of RA. Yet, it's difficult to gauge inflammation levels in RA patients. Some people who have RA do not exhibit normal symptoms, which makes it more challenging to make a diagnosis. Typical RA evaluations are subject to a few restrictions. Earlier research demonstrated that some patients continued to experience the progression of bone and joint degeneration even while in clinical remission. This progression was attributed to ongoing synovial inflammation. As a result, performing a precise evaluation of the level of inflammation is of the utmost importance. The neutrophil-to-lymphocyte ratio (NLR) has consistently been one of the most interesting novel non-specific inflammatory indicators. It is a reflection of the equilibrium between lymphocytes and neutrophils, which are inflammatory regulators and inflammatory activators, respectively. A higher NLR is linked to more severe levels of imbalance and inflammation. The aim of this study was to depict the role of NLR in RA progression and to show if NLR could predict the response to disease-modifying antirheumatic drugs (DMARDs) therapy in RA.

Peripheral blood cellular dynamics of rheumatoid arthritis treatment informs about efficacy of response to disease modifying drugs

Rheumatoid arthritis (RA) is an autoimmune disease characterized by systemic inflammation and is mediated by multiple immune cell types. In this work, we aimed to determine the relevance of changes in cell proportions in peripheral blood mononuclear cells (PBMCs) during the development of disease and following treatment. Samples from healthy blood donors, newly diagnosed RA patients, and established RA patients that had an inadequate response to MTX and were about to start tumor necrosis factor inhibitors (TNFi) treatment were collected before and after 3 months of treatment. We used in parallel a computational deconvolution approach based on RNA expression and flow cytometry to determine the relative cell-type frequencies. Cell-type frequencies from deconvolution of gene expression indicate that monocytes (both classical and non-classical) and CD4⁺ cells (T_h1 and T_h2) were increased in RA patients compared to controls, while NK cells and B cells (naïve and mature) were significantly decreased in RA patients. Treatment with MTX caused a decrease in B cells (memory and plasma cell), and a decrease in CD4 T_h cells (T_h1 and T_h17), while treatment with TNFi resulted in a significant increase in the population of B cells. Characterization of the RNA expression patterns found that most of the differentially expressed genes in RA subjects after treatment can be explained by changes in cell frequencies (98% and 74% respectively for MTX and TNFi).

Tailored therapeutic decision of rheumatoid arthritis using proteomic strategies: how to start and when to stop?

Unpredictable treatment responses have been an obstacle for the successful management of rheumatoid arthritis. Although numerous serum proteins have been proposed, there is a lack of integrative survey to compare their relevance in predicting treatment outcomes in rheumatoid arthritis. Also, little is known about their applications in various treatment stages, such as dose modification, drug switching or withdrawal. Here we present an in-depth exploration of the potential usefulness of serum proteins in clinical decision-making and unveil the spectrum of immunopathology underlying responders to different drugs. Patients with robust autoimmunity and inflammation are more responsive to biological treatments and prone to relapse during treatment de-escalation. Moreover, the concentration changes of serum proteins at the beginning of the treatments possibly assist early recognition of treatment responders. With a better understanding of the relationship between the serum proteome and treatment responses, personalized medicine in rheumatoid arthritis will be more achievable in the near future.

Synovial Inflammatory Pathways Characterize Anti-TNF-Responsive Rheumatoid Arthritis Patients

OBJECTIVE

This study was undertaken to understand the mechanistic basis of response to anti-tumor necrosis factor (anti-TNF) therapies and to determine whether transcriptomic changes in the synovium are reflected in peripheral protein markers.

METHODS

Synovial tissue from 46 rheumatoid arthritis (RA) patients was profiled with RNA sequencing before and 12 weeks after treatment with anti-TNF therapies. Pathway and gene signature analyses were performed on RNA expression profiles of synovial biopsies to identify mechanisms that could discriminate among patients with a good response, a moderate response, or no response, according to the American College of Rheumatology (ACR)/EULAR response criteria. Serum proteins encoded by synovial genes that were differentially expressed between ACR/EULAR response groups were measured in the same patients.

RESULTS

Gene signatures predicted which patients would have good responses, and pathway analysis identified elevated immune pathways, including chemokine signaling, Th1/Th2 cell differentiation, and Toll-like receptor signaling, uniquely in good responders. These inflammatory pathways were correspondingly down-modulated by anti-TNF therapy only in good responders. Based on cell signature analysis, lymphocyte, myeloid, and fibroblast cell populations were elevated in good responders relative to nonresponders, consistent with the increased inflammatory pathways. Cell signatures that decreased following anti-TNF treatment were predominately associated with lymphocytes, and fewer were associated with myeloid and fibroblast populations. Following anti-TNF treatment, and only in good responders, several peripheral inflammatory proteins decreased in a manner that was consistent with corresponding synovial gene changes.

CONCLUSION

Collectively, these data suggest that RA patients with robust responses to anti-TNF therapies are characterized at baseline by immune pathway activation, which decreases following anti-TNF treatment. Understanding mechanisms that define patient responsiveness to anti-TNF treatment may assist in development of predictive markers of patient response and earlier treatment options.

Dynamics of Type I and Type II Interferon Signature Determines Responsiveness to Anti-TNF Therapy in Rheumatoid Arthritis

The factors influencing long-term responses to a tumor necrosis factor inhibitor (TNFi) in rheumatoid arthritis (RA) patients currently remain unknown. Therefore, we herein conducted a multi-omics analysis of TNFi responses in a Japanese RA cohort. Blood samples were collected from 27 biological disease-modifying antirheumatic drug (DMARD)-naive RA patients at the initiation of and after three months of treatment with TNFi. Treatment responses were evaluated at one year. Differences in gene expression levels in peripheral blood mononuclear cells (PBMCs), plasma protein levels, drug concentrations, and the presence/absence of anti-drug antibodies were investigated, and a cell phenotypic analysis of PBMCs was performed using flow cytometry. After one year of treatment, thirteen patients achieved clinical remission (responders), while the others did not or switched to other biologics (non-responders). Differentially expressed genes related to treatment responses were enriched for the interferon (IFN) pathway. The expression of type I IFN signaling-related genes was higher in non-responders than in responders before and after treatment (P = 0.03, 0.005, respectively). The expression of type II IFN signaling-related genes did not significantly differ before treatment; however, it increased in non-responders and decreased in responders, with a significant difference being observed after three months of treatment (P = 1.2×10^-3). The total number of lymphocytes and C-X-C Motif Chemokine Ligand 10 (CXCL10) protein levels were associated with the type I IFN signature (P = 6.7×10^-7, 6.4×10^-3, respectively). Hepatocyte growth factor (HGF) protein levels before treatment predicted fold increases in type II IFN (P = 0.03). These IFN signature-related indices (the number of lymphocytes, CXCL10, and HGF) significantly differed between responders and non-responders (P = 0.01, 0.01, and 0.04, respectively). A single-cell analysis revealed that the type I IFN signature was more highly enriched in monocytes than in other cell types. A deconvolution analysis of bulk-RNA sequence data identified CD4+ and CD8+ T cells as the main sources of the type II IFN signature in non-responders. Collectively, the present results demonstrated that the dynamics of the type I and II IFN pathways affected long-term responses to TNFi, providing information on its biological background and potential for clinical applications.

Current status of use of high throughput nucleotide sequencing in rheumatology

OBJECTIVE

Here, we assess the usage of high throughput sequencing (HTS) in rheumatic research and the availability of public HTS data of rheumatic samples.

METHODS

We performed a semiautomated literature review on PubMed, consisting of an R-script and manual curation as well as a manual search on the Sequence Read Archive for public available HTS data.

RESULTS

Of the 699 identified articles, rheumatoid arthritis (n=182 publications, 26%), systemic lupus erythematous (n=161, 23%) and osteoarthritis (n=152, 22%) are among the rheumatic diseases with the most reported use of HTS assays. The most represented assay is RNA-Seq (n=457, 65%) for the identification of biomarkers in blood or synovial tissue. We also find, that the quality of accompanying clinical characterisation of the sequenced patients differs dramatically and we propose a minimal set of clinical data necessary to accompany rheumatological-relevant HTS data.

CONCLUSION

HTS allows the analysis of a broad spectrum of molecular features in many samples at the same time. It offers enormous potential in novel personalised diagnosis and treatment strategies for patients with rheumatic diseases. Being established in cancer research and in the field of Mendelian diseases, rheumatic diseases are about to become the third disease domain for HTS, especially the RNA-Seq assay. However, we need to start a discussion about reporting of clinical characterisation accompany rheumatological-relevant HTS data to make clinical meaningful use of this data.

Inference of an Integrative, Executable Network for Rheumatoid Arthritis Combining Data-Driven Machine Learning Approaches and a State-of-the-Art Mechanistic Disease Map

Rheumatoid arthritis (RA) is a multifactorial, complex autoimmune disease that involves various genetic, environmental, and epigenetic factors. Systems biology approaches provide the means to study complex diseases by integrating different layers of biological information. Combining multiple data types can help compensate for missing or conflicting information and limit the possibility of false positives. In this work, we aim to unravel mechanisms governing the regulation of key transcription factors in RA and derive patient-specific models to gain more insights into the disease heterogeneity and the response to treatment. We first use publicly available transcriptomic datasets (peripheral blood) relative to RA and machine learning to create an RA-specific transcription factor (TF) co-regulatory network. The TF cooperativity network is subsequently enriched in signalling cascades and upstream regulators using a state-of-the-art, RA-specific molecular map. Then, the integrative network is used as a template to analyse patients' data regarding their response to anti-TNF treatment and identify master regulators and upstream cascades affected by the treatment. Finally, we use the Boolean formalism to simulate in silico subparts of the integrated network and identify combinations and conditions that can switch on or off the identified TFs, mimicking the effects of single and combined perturbations.

Application of NGS Technology in Understanding the Pathology of Autoimmune Diseases

NGS technologies have transformed clinical diagnostics and broadly used from neonatal emergencies to adult conditions where the diagnosis cannot be made based on clinical symptoms. Autoimmune diseases reveal complicate molecular background and traditional methods could not fully capture them. Certainly, NGS technologies meet the needs of modern exploratory research, diagnostic and pharmacotherapy. Therefore, the main purpose of this review was to briefly present the application of NGS technology used in recent years in the understanding of autoimmune diseases paying particular attention to autoimmune connective tissue diseases. The main issues are presented in four parts: (a) panels, whole-genome and -exome sequencing (WGS and WES) in diagnostic, (b) Human leukocyte antigens (HLA) as a diagnostic tool, (c) RNAseq, (d) microRNA and (f) microbiome. Although all these areas of research are extensive, it seems that epigenetic impact on the development of systemic autoimmune diseases will set trends for future studies on this area.

Early Prediction of Clinical Response to Anti-TNF Treatment using Multi-omics and Machine Learning in Rheumatoid Arthritis

Objectives

Advances in immunotherapy by blocking TNF have remarkably improved treatment outcomes for Rheumatoid arthritis (RA) patients. Although treatment specifically targets TNF, the downstream mechanisms of immune suppression are not completely understood. The aim of this study was to detect biomarkers and expression signatures of treatment response to TNF inhibition.

Methods

Peripheral blood mononuclear cells (PBMCs) from 39 female patients were collected before anti-TNF treatment initiation (day 0) and after 3 months. The study cohort included patients previously treated with MTX who failed to respond adequately. Response to treatment was defined based on the EULAR criteria and classified 23 patients as responders and 16 as non-responders. We investigated differences in gene expression in PBMCs, the proportion of cell types and cell phenotypes in peripheral blood using flow cytometry and the level of proteins in plasma. Finally, we used machine learning models to predict non-response to anti-TNF treatment.

Results

The gene expression analysis in baseline samples revealed notably higher expression of the gene EPPK1 in future responders. We detected the suppression of genes and proteins following treatment, including suppressed expression of the T cell inhibitor gene CHI3L1 and its protein YKL-40. The gene expression results were replicated in an independent cohort. Finally, machine learning models mainly based on transcriptomic data showed high predictive utility in classifying non-response to anti-TNF treatment in RA.

Conclusions

Our integrative multi-omics analyses identified new biomarkers for the prediction of response, found pathways influenced by treatment and suggested new predictive models of anti-TNF treatment in RA patients.

A Molecular Signature Response Classifier to Predict Inadequate Response to Tumor Necrosis Factor-α Inhibitors: The NETWORK-004 Prospective Observational Study

Introduction

Timely matching of patients to beneficial targeted therapy is an unmet need in rheumatoid arthritis (RA). A molecular signature response classifier (MSRC) that predicts which patients with RA are unlikely to respond to tumor necrosis factor-α inhibitor (TNFi) therapy would have wide clinical utility.

Methods

The protein–protein interaction map specific to the rheumatoid arthritis pathophysiology and gene expression data in blood patient samples was used to discover a molecular signature of non-response to TNFi therapy. Inadequate response predictions were validated in blood samples from the CERTAIN cohort and a multicenter blinded prospective observational clinical study (NETWORK-004) among 391 targeted therapy-naïve and 113 TNFi-exposed patient samples. The primary endpoint evaluated the ability of the MSRC to identify patients who inadequately responded to TNFi therapy at 6 months according to ACR50. Additional endpoints evaluated the prediction of inadequate response at 3 and 6 months by ACR70, DAS28-CRP, and CDAI.

Results

The 23-feature molecular signature considers pathways upstream and downstream of TNFα involvement in RA pathophysiology. Predictive performance was consistent between the CERTAIN cohort and NETWORK-004 study. The NETWORK-004 study met primary and secondary endpoints. A molecular signature of non-response was detected in 45% of targeted therapy-naïve patients. The MSRC had an area under the curve (AUC) of 0.64 and patients were unlikely to adequately respond to TNFi therapy according to ACR50 at 6 months with an odds ratio of 4.1 (95% confidence interval 2.0–8.3, p value 0.0001). Odds ratios (3.4–8.8) were significant (p value < 0.01) for additional endpoints at 3 and 6 months, with AUC values up to 0.74. Among TNFi-exposed patients, the MSRC had an AUC of up to 0.83 and was associated with significant odds ratios of 3.3–26.6 by ACR, DAS28-CRP, and CDAI metrics.

Conclusion

The MSRC stratifies patients according to likelihood of inadequate response to TNFi therapy and provides patient-specific data to guide therapy choice in RA for targeted therapy-naïve and TNFi-exposed patients.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40744-021-00330-y.

A blood-based molecular signature response classifier (MSRC) integrating next-generation RNA sequencing data with clinical features predicts the likelihood that a patient with rheumatoid arthritis will have an inadequate response to TNFi therapy. Treatment selection guided by test results, with likely inadequate responders appropriately redirected to a different therapy, could improve response rates to TNFi therapies, generate healthcare cost savings, and increase rheumatologists’ confidence in prescribing decisions and altered treatment choices. The MSRC described in this study predicts the likelihood of inadequate response to TNFi therapies among targeted therapy-naïve and TNFi-exposed patients in a multicenter, 24-week blinded prospective clinical study: NETWORK-004. Patients with a molecular signature of non-response are less likely to have an adequate response to TNFi therapies than those patients lacking the signature according to ACR50, ACR70, CDAI, and DAS28-CRP with significant odds ratios of 3.4–8.8 for targeted therapy-naïve patients and 3.3–26.6 for TNFi-exposed patients. This MSRC provides a solution to the long-standing need for precision medicine tools to predict drug response in rheumatoid arthritis—a heterogeneous and progressive disease with an abundance of therapeutic options. These data validate the performance of the MSRC in a blinded prospective clinical study of targeted therapy-naïve and TNFi therapy-exposed patients.

Novel approaches to develop biomarkers predicting treatment responses to TNF-blockers

Introduction: Chronic inflammatory diseases (CIDs) cause significant morbidity and are a considerable burden for the patients in terms of pain, impaired function, and diminished quality of life. Important progress in CID treatment has been obtained with biological therapies, such as tumor-necrosis-factor blockers. However, more than a third of the patients fail to respond to these inhibitors and are exposed to the side effects of treatment, without the benefits. Therefore, there is a strong interest in developing tools to predict response of patients to biologics. Areas covered: The authors searched PubMed for recent studies on biomarkers for disease assessment and prediction of therapeutic responses, focusing on the effect of TNF blockers on immune responses in spondyloarthritis (SpA), and other CID, in particular rheumatoid arthritis and inflammatory bowel disease. Conclusions will be drawn about the possible development of predictive biomarkers for response to treatment. Expert opinion: No validated biomarker is currently available to predict treatment response in CID. New insight could be generated through the development of new bioinformatic modeling approaches to combine multidimensional biomarkers that explain the different genetic, immunological and environmental determinants of therapeutic responses.