Next Article in Journal
Synthesis of Novel Arsonolipids and Development of Novel Arsonoliposome Types
Next Article in Special Issue
Insights into the Pharmacogenetics of Tacrolimus Pharmacokinetics and Pharmacodynamics
Previous Article in Journal
Gold Nanoclusters-Based NIR-II Photosensitizers with Catalase-like Activity for Boosted Photodynamic Therapy
Previous Article in Special Issue
Pharmacogenetics and Pain Treatment with a Focus on Non-Steroidal Anti-Inflammatory Drugs (NSAIDs) and Antidepressants: A Systematic Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

-Omic Approaches and Treatment Response in Rheumatoid Arthritis

1
Department of Clinical Pharmacy, San Cecilio University Hospital, Instituto de Investigación Biosanitaria de Granada (ibs.Granada), 18016 Granada, Spain
2
Institute of Parasitology and Biomedicine López-Neyra (IPBLN), Spanish National Research Council (CSIC), 18017 Granada, Spain
*
Author to whom correspondence should be addressed.
Pharmaceutics 2022, 14(8), 1648; https://doi.org/10.3390/pharmaceutics14081648
Submission received: 29 June 2022 / Revised: 22 July 2022 / Accepted: 3 August 2022 / Published: 8 August 2022
(This article belongs to the Special Issue Association Studies in Clinical Pharmacogenetics)

Abstract

:
Rheumatoid arthritis (RA) is an inflammatory disorder characterized by an aberrant activation of innate and adaptive immune cells. There are different drugs used for the management of RA, including disease-modifying antirheumatic drugs (DMARDs). However, a significant percentage of RA patients do not initially respond to DMARDs. This interindividual variation in drug response is caused by a combination of environmental, genetic and epigenetic factors. In this sense, recent -omic studies have evidenced different molecular signatures involved in this lack of response. The aim of this review is to provide an updated overview of the potential role of -omic approaches, specifically genomics, epigenomics, transcriptomics, and proteomics, to identify molecular biomarkers to predict the clinical efficacy of therapies currently used in this disorder. Despite the great effort carried out in recent years, to date, there are still no validated biomarkers of response to the drugs currently used in RA. -Omic studies have evidenced significant differences in the molecular profiles associated with treatment response for the different drugs used in RA as well as for different cell types. Therefore, global and cell type-specific -omic studies analyzing response to the complete therapeutical arsenal used in RA, including less studied therapies, such as sarilumab and JAK inhibitors, are greatly needed.

1. Introduction

Rheumatoid arthritis (RA), which affects approximately 0.5–1% of the world population, is a chronic and progressive inflammatory disorder characterized by the appearance of synovitis and severe joint destruction that causes pain and disability [1,2]. The appearance of RA is related to the subsequent activation and proliferation of both innate and adaptive immune cells, including neutrophils, T and B lymphocytes, and monocytes, as well as fibroblasts of the synovial membrane, thus leading to persistent inflammation and damage of joints and bones [3]. This aberrant cell activation is mediated by genetic and environmental factors as well as by epigenetic factors, which represent the link between genetic and environment.
Fortunately, today there are therapeutic alternatives available to reduce this excessive immune activation. According to the 2019 EULAR recommendations, different drugs that target relevant immune system molecules are currently used for the management of RA: non-steroidal anti-inflammatory drugs, glucocorticoids, immunosuppressants and disease-modifying antirheumatic drugs (DMARDs), such as methotrexate (MTX), leflunomide, sulfasalazine (SSZ) and hydroxycloroquine (HCQ), Janus kinase inhibitors (JAKi) (baricitinib, tofacitinib y upadacitinib) and biological DMARDs (TNF inhibitors (TNFi), abatacept (ABA), rituximab (RTX), tocilizumab (TCZ) and sarilumab) [4]. However, treatment response varies among patients and a significant percentage of them do not respond to first-line biological DMARDs, which results in high joint destruction and, consequently, a poor quality of life [5]. In this sense, recent evidence indicates that omics approaches are key to clarify the molecular mechanisms that influence treatment response (Figure 1) [6].
In recent years, several studies evaluating the role of -omic data in the response to RA drugs have been carried out. The aim of this review is to provide an updated overview of the potential role of omics approaches, specifically genomics, epigenomics, transcriptomics, and proteomics, to identify biomarkers to predict the clinical efficacy of therapies currently used in this disorder.

2. Methodology

2.1. Study Selection

The search was carried out in PubMed, Scopus and Google Scholar databases. The search strategy included studies published until 6 June 2022. The combinations of terms used were: “rheumatoid arthritis”, “response”, “anti-TNF”, “adalimumab”, “infliximab”, “golimumab”, “certolizumab”, “etanercept”, “abatacept “tocilizumab”, “sarilumab”, “rituximab”, “baricitinib”, “tofacitinib”, “filgotinib”, “upadacitinib”, “effectiveness”, “genomics”, “epigenomics” “methylation”, “proteomics”, “miRNA”, “proteomics”, “transcriptomics”, “single-cell”.
This review included both observational cohort, case–control, cross-sectional studies, systematic reviews and meta-analyses (written in English). Animal and experimental model studies were excluded. In addition, studies in which treatment response results could not be extracted (responders/non-responders) were also excluded. Case reports, editorials, letters to the editor and conference proceedings were excluded from the review. Finally, a total of 72 references were included in this review.

2.2. Clinical Outcome

The definitions of the different clinical outcomes that have been used in the studies evaluated are detailed below:
-
Disease Activity Score 28-joint counts (DAS28) [7]: The result is calculated by using a special calculator that includes: Tender joint count (TJC) (of 28), swollen joint count (SJC) (of 28) and global health.
-
EULAR (European Alliance of Associations for Rheumatology) response criteria [8]: This outcome classifies patients (good, moderate and non-responders) depending on the change in DAS28 and the level of disease activity reached during follow-up.
-
CDAI (Clinical Disease Activity Index) [9]: This index is calculated using TJC (of 28), SJC (of 28), and patient and physician global assessment.
-
SDAI (Simplified Disease Activity Index) [10]: Similar to the CDAI + C-Reactive Protein (CRP).
-
HAQ-DIs (Health Assessment Questionnaire–Disability Index scores) [11]: a self-reported questionnaire covering 20 items in eight domains related to measuring difficulty in performing activities of daily living.
-
ACR20 (American College of Rheumatology) [12]: The ACR20 is a composite measure defined as both improvement of 20% in the TJC and SWC, and improvement of 20% in three of the following five criteria: patient and physician global assessment, functional ability measure (HAQ), visual analog pain scale, and erythrocyte sedimentation rate or CRP.
-
ACR/EULAR remission criteria [13]: This criteria include the index SDAI and CDAI and Boolean (SWJ (of 28), TJC (of 28), patient global assessment, and CRP).
A table including the abbreviations used throughout the manuscript has been added as Appendix A.

3. Pharmacogenomics Findings in Rheumatoid Arthritis

In recent years, genome-wide association studies (GWAS), which allow analyzing up to millions of single nucleotide polymorphisms (SNPs) through the whole genome, have emerged as an essential strategy to identify genetic risk loci involved in disease susceptibility. In addition, GWAS have also been applied to the identification of genetic polymorphisms involved in response to drugs. Specifically, 11 GWAS have evaluated the role of genetic variability in drug response for the different treatment used in RA.
To date, only two GWAS have analyzed the potential role of genetic variants in the response to MTX in RA patients. In this regard, a GWAS performed in 1424 RA patients from European ancestry identified an association that almost reached genome-wide significance (p = 9.8 × 10−8) between the NRG3-rs168201 polymorphism and the change in DAS28 at 6 months [14]. However, they could not replicate this association in two independent cohorts [14]. In addition, no genomic associations were found in a GWAS performed in an Indian cohort of 457 RA patients treated with MTX [15].
The largest replication study performed to date to identify genetic variants influencing response to MTX was performed in 915 patients. A total of 25 SNPs were evaluated (14 of them selected from candidate gene studies and 11 from GWAS). Patients carrying the A-allele of the MTRR-rs1801394 polymorphism showed a decrease in DAS28 after MTX treatment [16].
Recently, a systematic review that included 35 studies (34 candidate gene studies and the GWAS published by Senapati et al.) reported six promising SNPs after multiple correction. The genetic polymorphisms ATIC-rs7563206, TYMS-rs2847153 and TYMS-rs3786362 were associated with non-response to MTX. On the other hand, the genetic variants SLC19A1-rs1051266, DHFR-rs836788 and TYMS-rs2244500 were associated with response to this drug [17].
Genetic variability influencing response to TNFi has been more extensively studied during the last years. In 2010, a candidate gene study identified an association between the rs10919563 variant, located within PTPRC, and the response to TNFi [18]. Although subsequent GWASs have not detected this association, PTPRC currently represents the most replicated genetic biomarker of response to TNFi treatment in RA patients. In this regard, in 2016 a meta-analysis involving data from four studies (n = 2158) [18,19,20,21] showed that RA patients carrying the minor allele (A) of the PTPRC SNP presented a lower response to TNFi compared with those patients carrying the PTPRC major allele [22].
Two of the first GWASs performed in RA patients to evaluate the response to TNFi drugs (ADA, ETN, IFX) failed to identify significant differences between responder and non-responder patients, likely due to the low sample size included in these studies (n = 89 and n = 196) [23,24]. Although a genetic variant (rs3794271) located at the PDE3A-SLCO1C1 locus did not reach genome-wide association with response to TNFi in the GWAS performed by Krintel et al. [24], this polymorphism reached statistical significance in a later study carried out in a Spanish cohort, when data were combined in a meta-analysis [25]. However, three subsequent replication studies, some of them performed in large cohorts of RA patients, failed to validate these results [26,27,28].
In a GWAS carried out by Plant et al. [29], four SNPs (rs12081765 and rs7305646, located in intergenic regions, and rs1532269 and rs17301249, located at PDZ2D and EYA4, respectively) were found putatively associated with good response to TNFi therapy at 6 months. Although they not reached genome-wide significance, associations were validated in two independent replication cohorts. In addition, they found that some of the SNPs associated in the discovery cohort showed opposite allelic effects in the replication cohort. This heterogeneity may be one of the reasons why the results obtained are not consistent across the different studies. A subsequent replication study performed in a cohort of 634 Spanish patients treated with ADA, ETN and IFX tried to validate these four associations [30]. After meta-analysis with previous GWAS data (n = 2.298) [24,29,31], none of these genetic variants was replicated. Only the SNP rs1532269, located at the PDZD2 gene, yielded a suggestive association (p = 0.0033) with the response to TNFi.
Another GWAS was conducted in 882 RA patients treated with ADA, ETN and IFX through the Dutch Rheumatoid Arthritis Monitoring (DREAM) registry and the database of Apotheekzorg. They also included a replication cohort of 1821 patients from four cohorts. Again, none of the analyzed polymorphisms reached the genome-wide level of significance [31]. Moreover, they could not replicate the results obtained by previous GWASs [23,24,29].
In 2013, a large GWAS analyzed the response to therapy with various TNFi drugs (ADA, ETN and IFX) in 2706 RA patients collected through an international collaboration [32], including data of two previous GWASs [23,29]. Despite the higher statistical power of this study, no association with response to TNFi therapy was identified. However, when analyzing patients treated with each drug individually, a SNP located in the CD84 gene, involved in the maturation and activation of T lymphocytes, was associated with the efficacy of ETN treatment. Furthermore, this polymorphism influenced gene expression levels, and a higher expression was associated with a better response to the drug.
Finally, other GWAS performed in European population, including a discovery cohort of 375 RA patients and a replication cohort of 245 RA patients, identified the genetic polymorphism MED15-rs113878252 as a potential biomarker of response to ETN [33]. So far, only one GWAS trying to identify genetic biomarkers of response to ADA, ETN and IFX has been performed in patients from Asian origin, specifically from Japan. The genetic variant rs284511, which is close to the MAP3K7 locus, was significantly associated with ΔDAS at 6 months in this population [34].
In recent years, several studies focused on validating genetic associations with treatment response have been published. A study of 566 RA patients of Spanish and Greek ancestries treated with IFX, ADA, and ETN analyzed 18 SNPs previously associated with TNFi response by GWAS and candidate gene studies, but they failed to replicate these associations [27]. A replication study conducted in 755 RA patients did not show an association between 12 GWAS-drawn SNPs [29,31,32] and the response to IFX, ADA and ETN [35]. However, stratification of patients according to each specific TNFi drug allowed the researchers to identify an association between the rs2378945 variant, located at the NUBPL locus, and a poor response to ETN [35]. In addition, a recent replication study of 28 GWAS-identified variants was performed in a cohort of 1361 RA patients treated with different TNFi from the REPAIR consortium and the DANBIO registry [28]. They also tried to validate the most interesting results by performing a meta-analysis with a cohort of 706 RA patients. An association between the rs7767069 SNP at the LINC02549 locus and a lower response was found. Interestingly, the T-allele of this SNP correlated with a significantly increased number of T cells (CD45RO+CD45RA+), whereas patients carrying the TT genotype showed significantly increased serum levels of CD5 and CD6, which modulate T cells and certain subsets of B cells. Moreover, they also found that patients with positive rheumatoid factor carrying the G-allele of the LRRC55 rs717117 polymorphism presented a lower response to TNFi treatment [28].
So far, only one GWAS has been performed evaluating the role of genetic variants influencing the response to CZP. A cohort of 302 RA patients from the REALISTIC trial was included in the study, and the response was evaluated at 6 and 12 months. No statistically significant differences were found [36], which may be due to the limited sample size of this study.
Finally, the role of genetics in TCZ response has also been evaluated. Specifically, a GWAS analyzing over 1600 RA patients treated with this drug identified 8 loci associated with the clinical efficacy of TCZ [37]. Notably, a subsequent replication study of this GWAS conducted in 79 RA patients treated with TCZ replicated two of these associations. In this sense, RA patients carrying the GALNT18-rs4910008 C-allele or the CD69-rs11052877 A-allele presented better clinical outcome [38].
Most of the GWAS carried out to date have been focused on genetic variants influencing the response to IFX, ADA and ETN. Despite their increasing statistical power and the subsequent replication studies performed in different cohorts, there is still no clear evidence of any biomarker that could predict the response to these drugs when evaluated together. However, more interesting results were found when GWAS data (Table 1) were stratified according to the type of drug, which indicates that within the same type of therapy, the genetic basis of the response to treatment may vary depending on the specific drug, and therefore, consideration of TNFi therapy globally may be masking association signals.

4. Epigenomics and Treatment Response

Epigenetics is the set of inherited alterations in gene expression that are stable and do not produce any abnormality in the nucleotide sequence of DNA [39]. These variations lead to chemical alterations in DNA molecules and histones that play a crucial role in regulating gene transcription. The main epigenetic mechanisms include DNA methylation, histone modifications and regulatory non-coding RNAs. To date, studies evaluating the role of epigenetic mechanisms in the response to drugs in RA has been mainly focused on DNA methylation and microRNAs.

4.1. DNA Methylation

DNA methylation takes place in many regions of the genome and is considered one of the main mechanisms of gene expression regulation [40]. It consists of the addition of methyl groups to cytosines to form 5-methylcytosines by DNA-methyltransferases (DNMT) at four sites: repetitive sequences, “CpG islands shores”, CpG islands in the promoter region, and gene bodies throughout the genome. This methylation blocks the transcription of genes in the promoter region [41], and in general, elevated levels of 5-methylcytosine in the CpG-rich promoter region are primarily related to transcriptional repression.
Most of the studies evaluating the role of global DNA methylation in RA treatment efficacy have been focused on MTX (Table 2). In this sense, several studies have explored the impact of MTX treatment in the methylation status in different cell types, finding conflicting results. After one month of MTX treatment, significant increases in 5 mC percentage were evident in three different cell types (T and B cells and monocytes) from 19 early RA patients, most of them showing a decreased RA activity; however, no differences were observed in natural killers (NK) or polymorphonuclear leukocytes [42]. These results could not be confirmed in a cohort of 16 RA patients after 3 months of MTX, sarilumab and/or JAKi treatment. In this case, an increased DNA methylation pattern was found in lymphocytes (p = 0.033) but not in specific blood cell types (CD4+, CD8+, CD19+ and CD14+) [43]. When evaluating methylation levels of two different CD4+ T cells subsets obtained from 11 RA patients in remission after MTX treatment, 80% of the differentially methylated positions (DMPs) detected in CD4+ memory T cells showed decreased methylation levels, whereas in CD4 naïve T cells, similar percentages were found for hypo- and hypermethylated positions. Two genes, GRID2IP and PLEKHM1P1, showed decreased methylation levels in CD4+ memory T cells and in CD4+ memory T cells and naïve CD4+ T cells, respectively [44]. In addition, no statistically significant differences in the methylation status were found in whole blood leukocytes [45] and PBMCs [46] from RA patients naïve to treatment before and after 3 months of MTX treatment.
Differences in the methylome of responder and non-responder RA patients have also been evaluated in several studies. In this regard, an increased global DNA methylation level was associated with a lack of MTX response after 3 months of treatment in 181 patients with early RA naive to therapy [45]. These same authors carried out a subsequent study to identify differentially methylated regions (DMRs), in addition to DMPs, in 69 PBMC samples from RA patients before and after 3 months of treatment with MTX and corticosteroids as monotherapy or in combination with other DMARDs (SSZ and HCQ). They identified 1309 DMRs, but none of them reached genome-wide significance [46]. An earlier study analyzing the baseline DNA methylation profile of T lymphocytes from 46 patients treated with MTX in combination with other DMARDs (SSZ and HCQ) revealed two CpG sites, located near the ADAMTSL2 (hypermethylation) and BTN3A2 (hypomethylation) genes, that, when evaluated in combination, were strongly associated with response to treatment at 6 months (area under the ROC Curve (AUC): 0.85) [49]. Interestingly, both molecules play an important role in the immune system regulation [51,52]. However, the study performed by Gosselt et al. failed to validate the association between the DNA methylation level of cg14345882 (BTN3A2) and treatment response [46].
Additionally, it has also been explored whether the methylation patterns in the first weeks of MTX treatment could predict the response to this drug in 68 RA patients. At 4 weeks, 2 CpG sites (cg21040096, nearest gene RPH3AL and cg09894276, nearest gene WDR27) were associated with treatment response at 6 months [50]. The modifications in the methylation patterns identified at 4 weeks could be used as a tool to prevent a possible early failure of MTX treatment allowing drug change or escalation, if necessary.
A few studies assessing the role of methylation in clinical response to other RA therapies have also been published. Analysis of baseline methylation patterns of responder (n = 36) and non-responder patients to etanercept (ETN) (n = 36) after 3 months of treatment allowed researchers to identify five hypermethylated positions in responders. Two of the five top CpGs (cg04857395, p = 1.39 × 10−8 and cg26401028, p = 1.69 × 10−8) were located within exon 7 of the LRPAP1 gene that encodes a chaperone related to the transforming growth factor β activity [48], suggesting that this gene could be considered as a promising biomarker of ETN response. On the other hand, analysis of 26 patients treated with ETN or adalimumab (ADA) and 39 patients without biological agents failed to identify differences in the methylation profile and expression levels of DNMT1 and methyl-CpG-binding domain 2 (MBD2) between patients treated with and without TNFi [47].
Changes in the global methylation profile before and after treatment as well as in the methylation level of specific genes have been proposed as markers of response or early failure to treatment with contradictory results. Well-powered studies analyzing patients with a combined treatment of biologics and MTX are necessary to validate whether these potential response biomarkers can be extrapolated to the clinical practice.

4.2. miRNA Profiling

MicroRNAs (miRNAs) are non-coding RNAs of approximately 21 nucleotides in length. These molecules function as posttranscriptional repressors of gene expression required for appropriate cellular processes [53]. Recent findings have suggested that the epigenetic dysregulation, in particular, changes in the miRNA expression profile, could influence on the immune dysregulation observed in RA patients [54].
Different studies have evaluated the miRNA profile of rheumatoid arthritis patients treated with different drugs, identifying several miRNAs with a potential role as predictors of treatment response. The main characteristics of the studies detailed below are summarized in Table 3.
In a cohort of 95 RA patients, serum levels of miR-146a-5p, miR-125b, miR-126-3p, miR-23-3p, miR-16-5p and miR-223-3p were upregulated in patients responding to a combination of TNFi (ADA, ETN and IFX)/DMARDs at 6 months according to EULAR criteria. In addition, ROC analysis showed that increased serum levels of miR-23-3p and miR-223-3p before starting therapy were indicative of non-response with high specificity when considered together (91.5%), suggesting that they could be used as predictors of response to TNFi/DMARDs therapy [55]. Additionally, a microarray assay performed in a cohort of 108 RA patients identified a total of 59 upregulated and 78 downregulated miRNAs in PBMCs from ETN responder patients after 24 weeks of treatment. Increased levels of miR-146a-5p and decreased levels of let-7a-5p were validated by qPCR in a cohort of 92 RA patients [57]. Apart from these miRNAs, clinical outcomes, such as C-reactive protein and biologic history, were independently associated with lower clinical response. Taken together, these four biomarkers showed a high ability to predict clinical response (AUC = 0.863) [57].
Both studies found a role of miR-146a-5p in the clinical efficacy of TNFi. This miRNA, through the NF-κB pathway, is able to stimulate the release of pro-inflammatory cytokines involved in RA pathogenesis, such as TNF-α, interleukin (IL)-1b, IL-17 [64] and IL-6 [65]. In addition, an increase in miR-146a-5p serum expression after 3 months of TNFi treatment was also observed in a cohort of 13 RA patients (p = 0.033) [56]. Moreover, miR-125b was also found to be upregulated in responder patients by miRNA profiling. Similar findings have been reported in additional studies. A high baseline serum expression level of miR-125b in 32 RA patients was associated with better response to RTX after 3 months (p = 0.002) [58]. Moreover, higher miR-125b and miR-125a basal expressions were associated with better outcome in 96 active RA patients after 24 weeks of infliximab (IFX) treatment [59]. Interestingly, miR-125b plays an important role in regulating different signaling pathways involved in RA development, such as inflammation by activation of NF-κB pathway [66].
In order to identify predictor biomarkers of response to ADA + MTX, an analysis of 91 specific miRNAs was performed in 89 RA patients from the OPERA study before and after 3 and 12 months of treatment initiation. A higher pre-treatment plasma level of miR-27a-3p was significantly associated with remission at 12 months whereas increased levels were found in non-responders at 3 months post treatment. After performing two multivariate miRNA models in pre-treatment samples based on 1 (miR-19b-3p) and 10 (miR-146b-5p, -19b-3p, -27a-3p, -16-5p, -423-5p, -27b-3p, -23a-3p, -106a-5p, -29b-3p, and -17-5p) miRNAs, they found ROC curves with AUC of 67% and 84%, respectively [61].
In addition, studies analyzing expression levels of individual miRNAs have also described interesting findings regarding treatment response. The combination of low expression of miR-22 and high expression of miR-886-3p was associated with EULAR good response in 180 treatment-naïve RA patients treated with ADA [60]. miR-29b, which has a role in resistance to apoptosis, showed a decreased expression in RA patients with effective IFX therapy, but not in those treated with TCZ, which suggested that miR-29b levels may be informative with regard to immunotherapy choice [67]. In addition, serum expression levels of miRNA-5196 were significantly increased in 10 RA patients after TNFi treatment, including golimumab (GOL), ADA, and CZP. Interestingly, changes in miRNA-5196 expression positively correlated with DAS28 score. Taking this into account, miRNA-5196 could serve as a predictive biomarker of response to these drugs [62]. On the other hand, the only study assessing the predictive value of miRNAs in response to JAKi did not identify differences between RA patients in remission and not in remission after treatment, which may be due to the reduced sample size of this study (n = 16) [63].
Although there is already consensus on the influence of some miRNAs, such as miR-125b and miR-146a-5p, on the response to the different drugs used in RA, validation studies are needed to be able to use them as response biomarkers.

5. Transcriptomic Biomarkers

There are many studies published to date that attempt to associate transcriptomic changes with the response to the different drugs used in RA (Table 4).
Most of the studies exploring the potential role of transcriptomics in predicting clinical efficacy of RA treatments have been focused on TNFi. Prior to 2010, several genome-wide gene expression analyses assessing TNFi therapy outcome were performed, but the differentially expressed genes between responder and non-responder patients identified among studies showed low overlap. In order to validate previously reported gene expression signatures, a subsequent study linked eight previously published transcript sets [88,89,90,91,92] predicting TNFi response to the expression values of 42 RA patients treated with IFX and ADA. This approach allowed validating one of these eight predictive expression profiles. Specifically, the set of 20 genes reported by Lequerré et al. obtained the best results, with a sensitivity of 71% and a specificity of 61% for classifying RA patients [71].
In addition, expression levels of CD11c, an integrin involved in a variety of cell-matrix and cell–cell adhesion functions [93], in monocytes of RA patients treated with ADA was associated with future response [92]; however, no association with response to ADA (p = 0.33) or ETN (p = 0.13) was found in PBMCs from 75 patients [76].
Controversial results between the expression of CD39, which is primarily expressed on activated lymphoid cells, and the response to treatment have also been found in several studies. Two different studies described a lower expression of CD39 in poor responders to ADA (n = 70) [84] and MTX (n = 122) [94]. However, in a well-powered study including 2938 RA patients treated with all available TNFi drugs, a higher expression of CD39 was associated with a worse response [81]. In the first case, authors evaluated the response using the EULAR criteria, while in the last study, the clinical response was evaluated using the swollen joint count outcome. In addition to the great difference in sample size between studies, the groups of patients included were also different due to reasons of clinical practice (±MTX), which make it difficult to discern which drug is responsible for the effect.
Although no baseline differences between responders and non-responders were found in 240 RA patients treated with different TNFi, when clinical efficacy was evaluated after 14 weeks of treatment, the good responders’ group showed higher expression levels in gene co-expression modules (GCM) related to plasma, B and T cells, major histocompatibility complex, ribosomal proteins and undetermined modules and downregulation of myeloid lineage, platelets and inflammation GCM [77]. A recent study also tried to identify changes in GCMs during early TNFi treatment in two cohorts of RA patients treated with MTX (85) and ADA (70). One module was associated with ADA response and presented IFN type 1 signaling pathway genes (NFKBIE, IRF2BP2, MAPKAP-K2, IL1B and IFRD1) [87]. Moreover, correlation between GCM and response to other drugs, such as TCZ ± MTX, was studied in 60 RA patients. Network analysis within CD4+ T cells identified two GCM in the TCZ+MTX arm, four in the TCZ arm and four in the MTX arm significantly associated with sustained drug-free remission (sDFR). These modules included relevant pathways such as nuclear-transcribed mRNA catabolic processes and ribosome (TCZ + MTX arm), granulocyte migration (TCZ arm) and response to bacterium, p53 and JAK-STAT signaling (MTX arm) [80]. No differences were found in CD14+ cells.
By interrogating a synovial gene-expression dataset (GSE21537) including 62 patients treated with IFX, Dennis et al. found that good responder patients at week 16 presented higher baseline expression of a myeloid gene set (p= 0.011) and an enrichment in biological processes, such as classically activated M1 monocytes (p= 0.006) and angiogenesis (p = 0.018) [72]. Higher expression of DERL-1, a gene associated with autophagy, was associated with lack of IFX effectiveness in several cohorts of RA patients [86].
Expression of genes involved in the interferon (IFN) pathway has been widely studied in RA patients regardless of the administered drug. Patients who did not achieve IFX response showed an increased expression of the OAS1 (p = 0.033) and LGALS3BP (p = 0.041) genes [69]. Another report of 51 RA patients evidenced that upregulation of three IFN response genes negatively predicted response to RTX at week 12 (9/24 versus 20/27 [R2 = 0.17; p = 0.01]) [68]. Moreover, at week 24, eight IFN type I response genes were associated with poor response to RTX in 14 RA patients (p = 0.0074) [70]. In addition, correlation between response to TNFi therapy and expression of IFN signaling pathway genes was studied in neutrophils from patients treated with ADA, ETN and GOL. Specifically, those patients classified in the IFN-high group expression showed a better response than patients in the IFN-low group (ΔDAS28, OR: 1.4; p= 0.05) [75]. An IFN pathway downregulation after 24 weeks of RTX treatment was also observed in 68 RA responders patients [73]. When an anti-IL6 treatment was considered (TCZ), a higher expression of four type I IFN genes (IFI6, MT1G, MX2, and OASL) was associated with good response [74]. In most of the studies, a high expression of genes involved in the IFN pathway was associated with a poor response to different RA treatments. Further studies with larger sample sizes and more homogeneous in terms of treatments will provide us a greater insight into how this pathway might predict response to RA treatments.
Different signaling pathways and biomarkers have been postulated as responsible for the response to ABA in patients with RA. In a study including 19 patients, responders to MTX/ABA (n = 14) by EULAR criteria showed a significant enrichment of six signaling pathways (T cell receptor signaling, proteasome, angiogenesis, apoptosis and two mRNA processes) (Table 4) [83]. Additionally, in 45 bDMARD-naïve RA patients, a decreased type I IFN score and higher expression levels of dendritic cell-related or type I IFN-related genes (BATF2, LAMP3, CD83, CLEC4A, IDO1, IRF7, STAT1, STAT2 AND TNFSF10) were identified after treatment in responder patients by EULAR criteria, suggesting that ABA produces a reduction in IFN type I activity [82]. Upregulated expression of genes related to elongation, arrest and recovery (OR: 6.85; p = 0.03309) as well as NK-cell-related genes (OR: 6.46; p = 0.00388) have been associated with the lack of ABA effectiveness by CDAI [78].
Recently, the analysis of synovial tissue of 50 RA patients evidenced that a higher baseline expression of genes involved in the myeloid leukocyte and T cell activation pathways was related to a better response after treatment with DMARDs (MTX, ADA, ABA, RTX, TCZ) [85], suggesting that a high baseline immune activation may predict the response to this therapy.
On the other hand, single-cell sequencing technology is a novel strategy that is being used to have a better understanding of the cell subpopulations involved in immune-mediated pathologies [95]. In the case of RA, several studies have characterized the role of different cell types in disease pathogenesis at the single cell level [96], which could lead to the identification of new surface molecules potentially useful as therapeutic targets. In this sense, the molecules involved in the development of the disease could be more selectively modulated. Unfortunately, to date, there is no single-cell study evaluating the response to treatment in patients with RA. It is to be expected that the application of this strategy to the identification of molecular profiles of response to drugs represents a great advance for personalized medicine.
A significant number of candidate genes, GCM and pathways have been proposed as biomarkers of treatment response in patients with RA (TRAF6, CD11c, CD39, CHI3L1, EPPK1, CDC20, CXCR2, MPO, TNFAIP6, MYADM, FCGR2B, RFX2, IRF8 and TAF1, FOXO4, TAF11, IFN pathway and dendritic cell-related genes) but the differences in the characteristics of the studies published so far, in terms of type of tissue analyzed, time of response evaluation, sample size, therapy administered, patient characteristics, and analysis techniques, make it very difficult to reproduce the results obtained by the different studies.

6. Identification of Response Biomarkers by Proteomics

Evidence obtained in recent years suggests that proteome profiling also represents a useful tool for identifying response biomarkers. Table 5 summarized the more recent proteomics studies associated with the effectiveness to biological therapy in RA.
A serum proteome analysis was carried out in seven female RA patients treated with TCZ. Proteins expression levels were measured baseline and previous to the first and second dose of TCZ. A moderate/good response after two doses of TCZ was associated with an increased expression of 10 proteins, including apolipoproteins A-I, A-II, C-I, and C-II, Retinol binding protein 4, Selectin-L, Superoxide dismutase 3 or MCAM/MUC18/CD146). A decreased expression of seven proteins was observed after 2 months of TCZ treatment, such as pregnancy zone protein, a1-acid glycoprotein, C-reactive protein, haptoglobin, and the serine protease inhibitor clade A (a1-antitrypsin, leucine-rich a2-glycoprotein [97].
In addition to anti-citrullinated peptide antibodies (ACPA), RA patients can present other autoantibodies. In this sense, a study has analyzed the potential role of 376 different autoantibodies in treatment response by performing a multiplex bead-based assay in 286 RA patients treated with ADA or MTX [99]. By this approach, they identified two autoantibodies (citrullinated HNRNPA1 and citrullinated vimentin) significantly associated with treatment outcome by EULAR criteria at 3/6 months. In addition, the analysis of the ACPA patient subgroup identified an association of the presence of citrullinated CPSF6 with poor response at 3/6 months. Nevertheless, ACPA seropositivity is still the best marker to predict response to ADA/MTX.
Serum protein profiles at weeks 0 and 14 were investigated in 20 RA patients treated with IFX after one month of MTX and leflunomide. A total of 5 from the 13 differentially expressed proteins that overlapped (fibrinogen beta chain, haptoglobin, testicular tissue protein Li 70, C-reactive protein and serotransferrin) were validated, but only the serotransferrin protein was significant after verification by parallel reaction monitoring. Indeed, this protein involved in the hypoxia-inducible factor-1 pathway and ferroptosis, was upregulated in the responder patients (n = 15) and downregulated in the non-responder patients (n = 5) treated with IFX, MTX and leflunomide by EULAR criteria after 14 weeks of therapy, indicating the possible role of serotransferrin in resistance to these medications [100]. Although haptoglobin could not be validated as a response biomarker in this study, high levels of haptoglobin were associated with a good response in a cohort of 50 RA patients baseline and after 6 months of ETN treatment [98].

7. Multi-Omic Approaches for Response Prediction

Although individual -omic approaches have been useful in the identification of molecular biomarkers of response prediction, using integrative approaches that combine different -omic layers result essential to better understand how different mechanisms act in a complementary way to modify treatment response. Several studies that jointly evaluate the effect of the genome, methylome, transcriptome and proteome on the response to RA treatment have been published in recent years. The main characteristics of the studies detailed below are summarized in Table 6.
The effect on gene expression and protein levels of three drugs with different mechanisms of action has been analyzed in 45 patients with RA and 35 controls. At week 24, 600 transcripts were differentially expressed (FDR < 0.05) in patients treated with IFX or TCZ but not with MTX, and this effect was greater for TCZ than for IFX. In addition, in the case of TCZ, most changes in gene expression were in the direction toward the healthy state. In addition, when they studied the influence of TCZ and IFX on molecular remission at 14 weeks, TCZ and IFX, but not MTX, were able to achieve molecular remission at the protein level and, again, this effect was greater for TCZ. Although it would be necessary to carry out additional studies to confirm these results, it seems that TCZ has a more powerful effect than IFX and normalizes the molecular profiles of RA patients at transcriptome and protein levels. Interestingly, these drugs also produced transcriptional changes mainly in genes that were expressed at high or low levels in neutrophils, suggesting that the neutrophil signature was normalized by the drug treatments [102].
By integrating transcriptomic and genomic data, Aterido et al. identified a gene signature associated with TNFi response. The analysis of transcriptomic data from the synovium of 11 RA patients yielded an association between 13 GCMs and the response to TNFi at week 14. Notably, two of these GCMs were also associated with the response to ADA and IFX at the genetic level in GWAS data from a Spanish cohort (348 patients). In addition, the ADA-associated module, which was significantly enriched for genes involved in the nucleotide metabolism and epigenetic marks from CD4+ regulatory T cells, was validated in an independent GWAS set (2706 patients) (p < 0.05) [101], displaying the relevant role of these cells in mediating the response to TNFi [101].
Recently, a study evaluated gene expression and protein levels in PBMCs from 39 female RA patients and developed machine learning models to predict treatment response [104]. Analysis of baseline gene expression levels identified 192 differentially expressed genes between future responders and non-responders. Specifically, the genes EPPK1, BCL6-AS1 and CDC20 showed the highest differences between both subgroups. Some changes were also showed during treatment; in this regard, CXCR2, MPO, TNFAIP6 and MYADM were downregulated by treatment and FCGR2B appeared to be upregulated. Furthermore, some proteins, as CRP, IL-6, MMP-1, SAA, TNF-RI, VEGF, YKL-40, MIP-1 beta and MIG, were significantly suppressed in responders (FDR < 0.05) during anti-TNF treatment. In addition, a suppression of gene expression levels of CHI3L1 was evident upon TNFi treatment. Interestingly, this gene encodes YKL-40, one of the proteins that showed decreased levels after treatment [104]. Finally, they developed machine learning models that showed high prediction capacity in classifying non-responders RA patients before TNFi treatment, especially the model based on transcriptomic data (AUC = 0.81).
Differentially expressed genes and DMPs were found between responders and non-responders in PBMcs, monocytes and CD4+ T cells of 80 RA patients treated with ADA or ETN. In PBMCs, large differences were found in the gene expression profiles of responder and non-responder patients treated with ADA or ETN, including some genes involved in DNA nucleotide binding, specifically RFX2, IRF8 and TAF1 for the ADA cohort, and FOXO4 and TAF11 for the ETN cohort. TRAF6 involved in TNF receptor signaling, was only differentially expressed between responders and non-responders in the ETN cohort [103]. Interestingly, a low percentage of these differentially expressed genes overlap between both treatments. DMPs were strongly hypermethylated in ETN responders (76%) compared to ADA responders (46%). In addition, genes associated with the TNF signaling pathway were dysregulated in CD4+ T cells but not in monocytes of patients treated with ADA, which may be due to a clearer molecular TNF signaling signature associated with CD4+ T cells. CD4+ T cells from ETN responders showed upregulated genes in the FoxO signaling pathway and downregulated genes in the NOD-like receptor and JAK/STAT signaling pathways. With these results, the authors performed different machine learning models trying to predict the response to these drugs. The final model based on the best expression and methylation models predicted that approximately 30% of patients will not respond to ADA or ETN [103].
The DNA methylation profile in PBMCs and six immune system cells was evaluated recently in a discovery cohort of 62 RA and a validation cohort of 60 patients treated with different anti-TNF. At baseline after validation, 11 CpG sites from monocytes, 3 from NK, 2 from CD4+ T cells and 2 from neutrophils were associated with response. At week 12, 38 CpGs replicated in NK cells, 24 in neutrophils, 21 in B cells, 19 in monocytes, 13 in CD8+ T cells and 6 in CD4+ T cells were also associated with response. The cell-type deconvolution approach identified CpG sites in CD4+ T cells, NK cells, neutrophils and monocytes that were significantly associated with the response to TNFi. They found that from the 99 pathways modified at the epigenetic level in responders, 73 (73.7%) were also significantly altered at the transcriptomic level compared to 11 pathways statistically modified at the transcriptomic level from the 57 differentially methylated pathways in non-responders [105].

8. Conclusions

In contrast to studies focused on analyzing specific candidate molecules, -omic strategies allow the evaluation of a large number of potential biomarkers, thus representing a powerful tool to identify molecular signatures predicting treatment response. In this regard, the application of genomic, epigenomic, transcriptomic, and proteomic approaches to the discovery of response biomarkers has yielded relevant findings in RA.
However, despite the great effort carried out in recent years to clarify the role of epigenetic modifications, genes and proteins in the ineffectiveness of RA treatments, to date, there are still no validated biomarkers of response to the drugs currently used in this pathology. In this sense, the inconsistent results observed appear to be largely due to the heterogeneity across studies with respect to different variables that directly influence treatment outcome, such as patient phenotype, timing of sample collection, type of outcome, or response criteria. Therefore, it results essential to perform well-powered validation studies with homogenized conditions in order to be able to identify reliable treatment response biomarkers. In addition, -omic studies have evidenced significant differences in the molecular profiles associated with treatment response for the different drugs used in RA, as in the case of anti-TNF drugs, and also for different cell types. This highlights the importance of analyzing each drug independently as well as in individual cell types.
Therefore, global, single-cell and cell type-specific studies analyzing response to the complete therapeutical arsenal used in RA, including both biological and targeted synthetic therapies (sarilumab and JAK inhibitors), are greatly needed.
In addition, due to the complexity of the molecular mechanisms involved in treatment response, studies integrating different omics layers are essential to unravel the molecular network that determines clinical response, as well as to develop predictive models that allow the identification of RA patients with a greater probability of responding efficiently to treatment, thus leading to better clinical management of these patients.

Author Contributions

Conceptualization, A.M.-P. and A.M.; writing—original draft preparation, A.M.-P.; writing—review and editing, A.M. and J.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported through a Contract “Rio Hortega” [CM20/00141] (Instituto de Salud Carlos III) and by Redes de Investigación Cooperativa Orientadas a Resultados en Salud (RICORS) (RD21/0002/0039) from Instituto de Salud Carlos III.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Abbreviations.
Table A1. Abbreviations.
ABA (Abatacept)JAKi (Janus kinase inhibitors)
ACPA (Anti-Citrullinated Peptide Antibodies)IFN (Interferon)
ACR (American College of Rheumatology)IFX (Infliximab)
ADA (Adalimumab)MTX (Methotrexate)
CDAI (Clinical Disease Activity Index)NK (Natural Killers)
CPZ (Certolizumab pegol)RDM (Region Differentially Methylated)
CRP (C-Reactive Protein)RA (Rheumatoid Arthritis)
DAS28 (Disease Activity Score 28-joint counts) RTX (Rituximab)
DMARDs (Disease-Modifying Antirheumatic Drugs)SDAI (Simplified Disease Activity Index)
DMRs (Differentially Methylated Regions)TCZ (Tocilizumab)
DNMT (DNA-Methyltransferases)TNFi (TNF inhibitors)
ETN (Etanercept) Treg (Regulatory T cell)
EULAR (European Alliance of Associations for Rheumatology)PBMC (Peripheral Blood Mononuclear Cell)
GCM (Gene Coexpression Modules)PDM (Positions Differentially Methylated)
GWAS (Genome-Wide Association Study)RTX (Rituximab)
GOL (Golimumab)sDFR (Sustained Drug-Free Remission)
HAQ-DIs (Health Assessment Questionnaire–Disability Index scores)SNP (Single Nucleotide Polymorphisms)
HCQ (Hydroxicloroquine)SSZ (Sulfasalazine)

References

  1. Karami, J.; Masoumi, M.; Khorramdelazad, H.; Bashiri, H.; Darvishi, P.; Sereshki, H.A.; Shekarabi, M.; Sahebkar, A. Role of Autophagy in the Pathogenesis of Rheumatoid Arthritis: Latest Evidence and Therapeutic Approaches. Life Sci. 2020, 254, 117734. [Google Scholar] [CrossRef]
  2. Smolen, J.S.; Aletaha, D.; McInnes, I.B. Rheumatoid Arthritis. Lancet 2016, 388, 2023–2038. [Google Scholar] [CrossRef]
  3. Edilova, M.I.; Akram, A.; Abdul-Sater, A.A. Innate Immunity Drives Pathogenesis of Rheumatoid Arthritis. Biomed. J. 2021, 44, 172–182. [Google Scholar] [CrossRef]
  4. Smolen, J.S.; Landewé, R.B.M.; Bijlsma, J.W.J.; Burmester, G.R.; Dougados, M.; Kerschbaumer, A.; McInnes, I.B.; Sepriano, A.; Van Vollenhoven, R.F.; De Wit, M.; et al. EULAR Recommendations for the Management of Rheumatoid Arthritis with Synthetic and Biological Disease-Modifying Antirheumatic Drugs: 2019 Update. Ann. Rheum. Dis. 2020, 79, S685–S699. [Google Scholar] [CrossRef] [Green Version]
  5. Callaghan, C.A.; Boyter, A.C.; Mullen, A.B.; McRorie, E.R. Biological Therapy for Rheumatoid Arthritis: Is Personalised Medicine Possible? Eur. J. Hosp. Pharm. Sci. Pract. 2014, 21, 229–237. [Google Scholar] [CrossRef]
  6. Ota, M.; Fujio, K. Multi-Omics Approach to Precision Medicine for Immune-Mediated Diseases. Inflamm. Regen. 2021, 41, 23. [Google Scholar] [CrossRef] [PubMed]
  7. Prevoo, M.L.; van’t Hof, M.A.; Kuper, H.H.; van Leeuwen, M.A.; van de Putte, L.B.; van Riel, P.L. Modified Disease Activity Scores That Include Twenty-Eight-Joint Counts. Development and Validation in a Prospective Longitudinal Study of Patients with Rheumatoid Arthritis. Arthritis Rheum. 1995, 38, 44–48. [Google Scholar] [CrossRef] [Green Version]
  8. van Gestel, A.M.; Prevoo, M.L.; van’t Hof, M.A.; van Rijswijk, M.H.; van de Putte, L.B.; van Riel, P.L. Development and Validation of the European League Against Rheumatism Response Criteria for Rheumatoid Arthritis. Comparison with the Preliminary American College of Rheumatology and the World Health Organization/International League Against Rheumatism Criteria. Arthritis Rheum. 1996, 39, 34–40. [Google Scholar] [CrossRef]
  9. Aletaha, D.; Nell, V.P.K.; Stamm, T.; Uffmann, M.; Pflugbeil, S.; Machold, K.; Smolen, J.S. Acute Phase Reactants Add Little to Composite Disease Activity Indices for Rheumatoid Arthritis: Validation of a Clinical Activity Score. Arthritis Res. Ther. 2005, 7, R796–R806. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Smolen, J.S.; Breedveld, F.C.; Schiff, M.H.; Kalden, J.R.; Emery, P.; Eberl, G.; van Riel, P.L.; Tugwell, P. A Simplified Disease Activity Index for Rheumatoid Arthritis for Use in Clinical Practice. Rheumatology 2003, 42, 244–257. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Fries, J.F.; Spitz, P.; Kraines, R.G.; Holman, H.R. Measurement of Patient Outcome in Arthritis. Arthritis Rheum. 1980, 23, 137–145. [Google Scholar] [CrossRef] [PubMed]
  12. Felson, D.T.; Anderson, J.J.; Boers, M.; Bombardier, C.; Furst, D.; Goldsmith, C.; Katz, L.M.; Lightfoot, R.J.; Paulus, H.; Strand, V. American College of Rheumatology. Preliminary Definition of Improvement in Rheumatoid Arthritis. Arthritis Rheum. 1995, 38, 727–735. [Google Scholar] [CrossRef]
  13. Felson, D.T.; Smolen, J.S.; Wells, G.; Zhang, B.; van Tuyl, L.H.D.; Funovits, J.; Aletaha, D.; Allaart, C.F.; Bathon, J.; Bombardieri, S.; et al. American College of Rheumatology/European League against Rheumatism Provisional Definition of Remission in Rheumatoid Arthritis for Clinical Trials. Ann. Rheum. Dis. 2011, 70, 404–413. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Taylor, J.C.; Bongartz, T.; Massey, J.; Mifsud, B.; Spiliopoulou, A.; Scott, I.C.; Wang, J.; Morgan, M.; Plant, D.; Colombo, M.; et al. Genome-Wide Association Study of Response to Methotrexate in Early Rheumatoid Arthritis Patients. Pharm. J. 2018, 18, 528–538. [Google Scholar] [CrossRef]
  15. Senapati, S.; Singh, S.; Das, M.; Kumar, A.; Gupta, R.; Kumar, U.; Jain, S.; Juyal, R.C.; Thelma, B.K. Genome-Wide Analysis of Methotrexate Pharmacogenomics in Rheumatoid Arthritis Shows Multiple Novel Risk Variants and Leads for TYMS Regulation. Pharmacogenet. Genom. 2014, 24, 211–219. [Google Scholar] [CrossRef]
  16. López-Rodríguez, R.; Ferreiro-Iglesias, A.; Lima, A.; Bernardes, M.; Pawlik, A.; Paradowska-Gorycka, A.; Świerkot, J.; Slezak, R.; Dolžan, V.; González-Álvaro, I.; et al. Replication Study of Polymorphisms Associated with Response to Methotrexate in Patients with Rheumatoid Arthritis. Sci. Rep. 2018, 8, 7342. [Google Scholar] [CrossRef]
  17. Eektimmerman, F.; Swen, J.J.; Madhar, M.B.; Allaart, C.F.; Guchelaar, H.-J. Predictive Genetic Biomarkers for the Efficacy of Methotrexate in Rheumatoid Arthritis: A Systematic Review. Pharm. J. 2020, 20, 159–168. [Google Scholar] [CrossRef]
  18. Cui, J.; Saevarsdottir, S.; Thomson, B.; Padyukov, L.; Van Der Helm-van Mil, A.H.M.; Nititham, J.; Hughes, L.B.; De Vries, N.; Raychaudhuri, S.; Alfredsson, L.; et al. Rheumatoid Arthritis Risk Allele PTPRC Is Also Associated with Response to Anti-Tumor Necrosis Factor α Therapy. Arthritis Rheum. 2010, 62, 1849–1861. [Google Scholar] [CrossRef] [PubMed]
  19. Ferreiro-Iglesias, A.; Montes, A.; Perez-Pampin, E.; Cañete, J.D.; Raya, E.; Magro-Checa, C.; Vasilopoulos, Y.; Sarafidou, T.; Caliz, R.; Ferrer, M.A.; et al. Replication of PTPRC as Genetic Biomarker of Response to TNF Inhibitors in Patients with Rheumatoid Arthritis. Pharm. J. 2016, 16, 137–140. [Google Scholar] [CrossRef] [PubMed]
  20. Zervou, M.I.; Myrthianou, E.; Flouri, I.; Plant, D.; Chlouverakis, G.; Castro-Giner, F.; Rapsomaniki, P.; Barton, A.; Boumpas, D.T.; Sidiropoulos, P.; et al. Lack of Association of Variants Previously Associated with Anti-TNF Medication Response in Rheumatoid Arthritis Patients: Results from a Homogeneous Greek Population. PLoS ONE 2013, 8, e74375. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Plant, D.; Prajapati, R.; Hyrich, K.L.; Morgan, A.W.; Wilson, A.G.; Isaacs, J.D.; Barton, A. Replication of Association of the PTPRC Gene with Response to Anti-Tumor Necrosis Factor Therapy in a Large UK Cohort. Arthritis Rheum. 2012, 64, 665–670. [Google Scholar] [CrossRef] [Green Version]
  22. Lee, Y.H.; Bae, S.-C. Associations between PTPRC Rs10919563 A/G and FCGR2A R131H Polymorphisms and Responsiveness to TNF Blockers in Rheumatoid Arthritis: A Meta-Analysis. Rheumatol. Int. 2016, 36, 837–844. [Google Scholar] [CrossRef]
  23. Liu, C.; Batliwalla, F.; Li, W.; Lee, A.; Roubenoff, R.; Beckman, E.; Khalili, H.; Damle, A.; Kern, M.; Furie, R.; et al. Genome-Wide Association Scan Identifies Candidate Polymorphisms Associated with Differential Response to Anti-TNF Treatment in Rheumatoid Arthritis. Mol. Med. 2008, 14, 575–581. [Google Scholar] [CrossRef]
  24. Krintel, S.B.; Palermo, G.; Johansen, J.S.; Germer, S.; Essioux, L.; Benayed, R.; Badi, L.; Ostergaard, M.; Hetland, M.L. Investigation of Single Nucleotide Polymorphisms and Biological Pathways Associated with Response to TNFα Inhibitors in Patients with Rheumatoid Arthritis. Pharm. Genom. 2012, 22, 577–589. [Google Scholar] [CrossRef]
  25. Acosta-Colman, I.; Palau, N.; Tornero, J.; Fernández-Nebro, A.; Blanco, F.; González-Alvaro, I.; Cañete, J.D.; Maymó, J.; Ballina, J.; Fernández-Gutiérrez, B.; et al. GWAS Replication Study Confirms the Association of PDE3A-SLCO1C1 with Anti-TNF Therapy Response in Rheumatoid Arthritis. Pharmacogenomics 2013, 14, 727–734. [Google Scholar] [CrossRef]
  26. Smith, S.L.; Plant, D.; Lee, X.H.; Massey, J.; Hyrich, K.; Morgan, A.W.; Wilson, A.G.; Isaacs, J.; Barton, A. Previously Reported PDE3A-SLCO1C1 Genetic Variant Does Not Correlate with Anti-TNF Response in a Large UK Rheumatoid Arthritis Cohort. Pharmacogenomics 2016, 17, 715–720. [Google Scholar] [CrossRef] [Green Version]
  27. Lopez-Rodriguez, R.; Perez-Pampin, E.; Marquez, A.; Blanco, F.J.; Joven, B.; Carreira, P.; Ferrer, M.A.; Caliz, R.; Valor, L.; Narvaez, J.; et al. Validation Study of Genetic Biomarkers of Response to TNF Inhibitors in Rheumatoid Arthritis. PLoS ONE 2018, 13, e0196793. [Google Scholar] [CrossRef] [Green Version]
  28. Sánchez-Maldonado, J.M.; Cáliz, R.; López-Nevot, M.Á.; Cabrera-Serrano, A.J.; Moñiz-Díez, A.; Canhão, H.; Ter Horst, R.; Quartuccio, L.; Sorensen, S.B.; Glintborg, B.; et al. Validation of GWAS-Identified Variants for Anti-TNF Drug Response in Rheumatoid Arthritis: A Meta-Analysis of Two Large Cohorts. Front. Immunol. 2021, 12, 672255. [Google Scholar] [CrossRef]
  29. Plant, D.; Bowes, J.; Potter, C.; Hyrich, K.L.; Morgan, A.W.; Wilson, A.G.; Isaacs, J.D.; Barton, A. Genome-Wide Association Study of Genetic Predictors of Anti-Tumor Necrosis Factor Treatment Efficacy in Rheumatoid Arthritis Identifies Associations with Polymorphisms at Seven Loci. Arthritis Rheum. 2011, 63, 645–653. [Google Scholar] [CrossRef] [Green Version]
  30. Márquez, A.; Ferreiro-Iglesias, A.; Dávila-Fajardo, C.L.; Montes, A.; Pascual-Salcedo, D.; Perez-Pampin, E.; Moreno-Ramos, M.J.; García-Portales, R.; Navarro, F.; Moreira, V.; et al. Lack of Validation of Genetic Variants Associated with Anti-Tumor Necrosis Factor Therapy Response in Rheumatoid Arthritis: A Genome-Wide Association Study Replication and Meta-Analysis. Arthritis Res. Ther. 2014, 16, R66. [Google Scholar] [CrossRef] [Green Version]
  31. Mirkov, M.U.; Cui, J.; Vermeulen, S.H.; Stahl, E.A.; Toonen, E.J.M.; Makkinje, R.R.; Lee, A.T.; Huizinga, T.W.J.; Allaart, R.; Barton, A.; et al. Genome-Wide Association Analysis of Anti-TNF Drug Response in Patients with Rheumatoid Arthritis. Ann. Rheum. Dis. 2013, 72, 1375–1381. [Google Scholar] [CrossRef] [PubMed]
  32. Cui, J.; Stahl, E.A.; Saevarsdottir, S.; Miceli, C.; Diogo, D.; Trynka, G.; Raj, T.; Mirkov, M.U.; Canhao, H.; Ikari, K.; et al. Genome-Wide Association Study and Gene Expression Analysis Identifies CD84 as a Predictor of Response to Etanercept Therapy in Rheumatoid Arthritis. PLoS Genet. 2013, 9, e1003394. [Google Scholar] [CrossRef] [Green Version]
  33. Julià, A.; Fernandez-Nebro, A.; Blanco, F.; Ortiz, A.; Cañete, J.D.; Maymó, J.; Alperi-López, M.; Fernández-Gutierrez, B.; Olivè, A.; Corominas, H.; et al. A Genome-Wide Association Study Identifies a New Locus Associated with the Response to Anti-TNF Therapy in Rheumatoid Arthritis. Pharm. J. 2016, 16, 147–150. [Google Scholar] [CrossRef]
  34. Honne, K.; Hallgrímsdóttir, I.; Wu, C.; Sebro, R.; Jewell, N.P.; Sakurai, T.; Iwamoto, M.; Minota, S.; Jawaheer, D. A Longitudinal Genome-Wide Association Study of Anti-Tumor Necrosis Factor Response among Japanese Patients with Rheumatoid Arthritis. Arthritis Res. Ther. 2016, 18, 12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Ferreiro-Iglesias, A.; Montes, A.; Perez-Pampin, E.; Cañete, J.D.; Raya, E.; Magro-Checa, C.; Vasilopoulos, Y.; Caliz, R.; Ferrer, M.A.; Joven, B.; et al. Evaluation of 12 GWAS-Drawn SNPs as Biomarkers of Rheumatoid Arthritis Response to TNF Inhibitors. A Potential SNP Association with Response to Etanercept. PLoS ONE 2019, 14, e0213073. [Google Scholar] [CrossRef]
  36. White, I.R.; Kleinstein, S.E.; Praet, C.; Chamberlain, C.; McHale, D.; Maia, J.M.; Xie, P.; Goldstein, D.B.; Urban, T.J.; Shea, P.R. A Genome-Wide Screen for Variants Influencing Certolizumab Pegol Response in a Moderate to Severe Rheumatoid Arthritis Population. PLoS ONE 2022, 17, e0261165. [Google Scholar] [CrossRef] [PubMed]
  37. Wang, J.; Bansal, A.T.; Martin, M.; Germer, S.; Benayed, R.; Essioux, L.; Lee, J.S.; Begovich, A.; Hemmings, A.; Kenwright, A.; et al. Genome-Wide Association Analysis Implicates the Involvement of Eight Loci with Response to Tocilizumab for the Treatment of Rheumatoid Arthritis. Pharm. J. 2013, 13, 235–241. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Maldonado-Montoro, M.; Cañadas-Garre, M.; González-Utrilla, A.; Plaza-Plaza, J.C.; Calleja-Hernández, M.Á. Genetic and Clinical Biomarkers of Tocilizumab Response in Patients with Rheumatoid Arthritis. Pharmacol. Res. 2016, 111, 264–271. [Google Scholar] [CrossRef]
  39. Cavalli, G.; Heard, E. Advances in Epigenetics Link Genetics to the Environment and Disease. Nature 2019, 571, 489–499. [Google Scholar] [CrossRef] [Green Version]
  40. Guo, S.; Xu, L.; Chang, C.; Zhang, R.; Jin, Y.; He, D. Epigenetic Regulation Mediated by Methylation in the Pathogenesis and Precision Medicine of Rheumatoid Arthritis. Front. Genet. 2020, 11, 811. [Google Scholar] [CrossRef]
  41. Ahmadi, M.; Gharibi, T.; Dolati, S.; Rostamzadeh, D.; Aslani, S.; Baradaran, B.; Younesi, V.; Yousefi, M. Epigenetic Modifications and Epigenetic Based Medication Implementations of Autoimmune Diseases. Biomed. Pharmacother. 2017, 87, 596–608. [Google Scholar] [CrossRef]
  42. de Andres, M.C.; Perez-Pampin, E.; Calaza, M.; Santaclara, F.J.; Ortea, I.; Gomez-Reino, J.J.; Gonzalez, A. Assessment of Global DNA Methylation in Peripheral Blood Cell Subpopulations of Early Rheumatoid Arthritis before and after Methotrexate. Arthritis Res. Ther. 2015, 17, 233. [Google Scholar] [CrossRef] [Green Version]
  43. Liebold, I.; Grützkau, A.; Göckeritz, A.; Gerl, V.; Lindquist, R.; Feist, E.; Zänker, M.; Häupl, T.; Poddubnyy, D.; Zernicke, J.; et al. Peripheral Blood Mononuclear Cells Are Hypomethylated in Active Rheumatoid Arthritis and Methylation Correlates with Disease Activity. Rheumatology 2021, 60, 1984–1995. [Google Scholar] [CrossRef]
  44. Guderud, K.; Sunde, L.H.; Flåm, S.T.; Mæhlen, M.T.; Mjaavatten, M.D.; Lillegraven, S.; Aga, A.B.; Evenrød, I.M.; Norli, E.S.; Andreassen, B.K.; et al. Rheumatoid Arthritis Patients, Both Newly Diagnosed and Methotrexate Treated, Show More DNA Methylation Differences in CD4+ Memory Than in CD4+ Naïve T Cells. Front. Immunol. 2020, 11, 194. [Google Scholar] [CrossRef]
  45. Gosselt, H.R.; Van Zelst, B.D.; De Rotte, M.C.F.J.; Hazes, J.M.W.; De Jonge, R.; Heil, S.G. Higher Baseline Global Leukocyte DNA Methylation Is Associated with MTX Non-Response in Early RA Patients. Arthritis Res. Ther. 2019, 21, 157. [Google Scholar] [CrossRef] [Green Version]
  46. Gosselt, H.R.; Vallerga, C.L.; Mandaviya, P.R.; Lubberts, E.; Hazes, J.M.W.; De Jonge, R.; Heil, S.G. Epigenome Wide Association Study of Response to Methotrexate in Early Rheumatoid Arthritis Patients. PLoS ONE 2021, 16, e0247709. [Google Scholar] [CrossRef]
  47. Liu, C.C.; Fang, T.J.; Ou, T.T.; Wu, C.C.; Li, R.N.; Lin, Y.C.; Lin, C.H.; Tsai, W.C.; Liu, H.W.; Yen, J.H. Global DNA Methylation, DNMT1, and MBD2 in Patients with Rheumatoid Arthritis. Immunol. Lett. 2011, 135, 96–99. [Google Scholar] [CrossRef]
  48. Plant, D.; Webster, A.; Nair, N.; Oliver, J.; Smith, S.L.; Eyre, S.; Hyrich, K.L.; Wilson, A.G.; Morgan, A.W.; Isaacs, J.D.; et al. Differential Methylation as a Biomarker of Response to Etanercept in Patients with Rheumatoid Arthritis. Arthritis Rheumatol. 2016, 68, 1353–1360. [Google Scholar] [CrossRef]
  49. Glossop, J.R.; Nixon, N.B.; Emes, R.D.; Sim, J.; Packham, J.C.; Mattey, D.L.; Farrell, W.E.; Fryer, A.A. DNA Methylation at Diagnosis Is Associated with Response to Disease-Modifying Drugs in Early Rheumatoid Arthritis. Epigenomics 2017, 9, 419–428. [Google Scholar] [CrossRef] [Green Version]
  50. Nair, N.; Plant, D.; Plant, D.; Verstappen, S.M.; Verstappen, S.M.; Isaacs, J.D.; Morgan, A.W.; Hyrich, K.L.; Hyrich, K.L.; Hyrich, K.L.; et al. Differential DNA Methylation Correlates with Response to Methotrexate in Rheumatoid Arthritis. Rheumatology 2020, 59, 1364–1371. [Google Scholar] [CrossRef] [Green Version]
  51. Le Goff, C.; Morice-Picard, F.; Dagoneau, N.; Wang, L.W.; Perrot, C.; Crow, Y.J.; Bauer, F.; Flori, E.; Prost-Squarcioni, C.; Krakow, D.; et al. ADAMTSL2 Mutations in Geleophysic Dysplasia Demonstrate a Role for ADAMTS-like Proteins in TGF-Beta Bioavailability Regulation. Nat. Genet. 2008, 40, 1119–1123. [Google Scholar] [CrossRef] [Green Version]
  52. Smith, I.A.; Knezevic, B.R.; Ammann, J.U.; Rhodes, D.A.; Aw, D.; Palmer, D.B.; Mather, I.H.; Trowsdale, J. BTN1A1, the Mammary Gland Butyrophilin, and BTN2A2 Are Both Inhibitors of T Cell Activation. J. Immunol. 2010, 184, 3514–3525. [Google Scholar] [CrossRef]
  53. Stanczyk, J.; Ospelt, C.; Karouzakis, E.; Filer, A.; Raza, K.; Kolling, C.; Gay, R.; Buckley, C.D.; Tak, P.P.; Gay, S.; et al. Altered Expression of MicroRNA-203 in Rheumatoid Arthritis Synovial Fibroblasts and Its Role in Fibroblast Activation. Arthritis Rheum. 2011, 63, 373–381. [Google Scholar] [CrossRef] [Green Version]
  54. Huang, R.-Y.; Wu, J.-Q.; Liu, Z.-H.; Sun, S.-L. MicroRNAs in Rheumatoid Arthritis: What Is the Latest with Regards to Diagnostics? Expert Rev. Mol. Diagn. 2019, 19, 363–366. [Google Scholar] [CrossRef]
  55. Castro-Villegas, C.; Pérez-Sánchez, C.; Escudero, A.; Filipescu, I.; Verdu, M.; Ruiz-Limón, P.; Aguirre, M.A.; Jiménez-Gomez, Y.; Font, P.; Rodriguez-Ariza, A.; et al. Circulating MiRNAs as Potential Biomarkers of Therapy Effectiveness in Rheumatoid Arthritis Patients Treated with Anti-TNFα. Arthritis Res. Ther. 2015, 17, 49. [Google Scholar] [CrossRef] [Green Version]
  56. Bogunia-Kubik, K.; Wysoczańska, B.; Piątek, D.; Iwaszko, M.; Ciechomska, M.; Świerkot, J. Significance of Polymorphism and Expression of MiR-146a and NFkB1 Genetic Variants in Patients with Rheumatoid Arthritis. Arch. Immunol. Ther. Exp. 2016, 64, 131–136. [Google Scholar] [CrossRef] [Green Version]
  57. Liu, Y.; Han, Y.; Qu, H.; Fang, J.; Ye, M.; Yin, W. Correlation of MicroRNA Expression Profile with Clinical Response to Tumor Necrosis Factor Inhibitor in Treating Rheumatoid Arthritis Patients: A Prospective Cohort Study. J. Clin. Lab. Anal. 2019, 33, e22953. [Google Scholar] [CrossRef] [Green Version]
  58. Duroux-Richard, I.; Pers, Y.M.; Fabre, S.; Ammari, M.; Baeten, D.; Cartron, G.; Touitou, I.; Jorgensen, C.; Apparailly, F. Circulating MiRNA-125b Is a Potential Biomarker Predicting Response to Rituximab in Rheumatoid Arthritis. Mediat. Inflamm. 2014, 2014, 342524. [Google Scholar] [CrossRef]
  59. Cheng, P.; Wang, J. The Potential of Circulating MicroRNA-125a and MicroRNA-125b as Markers for Inflammation and Clinical Response to Infliximab in Rheumatoid Arthritis Patients. J. Clin. Lab. Anal. 2020, 34, e23329. [Google Scholar] [CrossRef] [Green Version]
  60. Krintel, S.B.; Dehlendorff, C.; Hetland, M.L.; Hørslev-Petersen, K.; Andersen, K.K.; Junker, P.; Pødenphant, J.; Ellingsen, T.; Ahlquist, P.; Lindegaard, H.M.; et al. Prediction of Treatment Response to Adalimumab: A Double-Blind Placebo-Controlled Study of Circulating MicroRNA in Patients with Early Rheumatoid Arthritis. Pharm. J. 2016, 16, 141–146. [Google Scholar] [CrossRef]
  61. Sode, J.; Krintel, S.B.; Carlsen, A.L.; Hetland, M.L.; Johansen, J.S.; HØrslev-Petersen, K.; Stengaard-Pedersen, K.; Ellingsen, T.; Burton, M.; Junker, P.; et al. Plasma MicroRNA Profiles in Patients with Early Rheumatoid Arthritis Responding to Adalimumab plus Methotrexate vs. Methotrexate Alone: A Placebo-Controlled Clinical Trial. J. Rheumatol. 2018, 45, 53–61. [Google Scholar] [CrossRef] [PubMed]
  62. Ciechomska, M.; Bonek, K.; Merdas, M.; Zarecki, P.; Swierkot, J.; Gluszko, P.; Bogunia-Kubik, K.; Maslinski, W. Changes in MiRNA-5196 Expression as a Potential Biomarker of Anti-TNF-α Therapy in Rheumatoid Arthritis and Ankylosing Spondylitis Patients. Arch. Immunol. Ther. Exp. 2018, 66, 389–397. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Fernández-Ruiz, J.C.; Ramos-Remus, C.; Sánchez-Corona, J.; Castillo-Ortiz, J.D.; Castañeda-Sánchez, J.J.; Bastian, Y.; Romo-García, M.F.; Ochoa-González, F.; Monsivais-Urenda, A.E.; González-Amaro, R.; et al. Analysis of MiRNA Expression in Patients with Rheumatoid Arthritis during Remission and Relapse after a 5-Year Trial of Tofacitinib Treatment. Int. Immunopharmacol. 2018, 63, 35–42. [Google Scholar] [CrossRef]
  64. Taganov, K.D.; Boldin, M.P.; Chang, K.-J.; Baltimore, D. NF-KappaB-Dependent Induction of MicroRNA MiR-146, an Inhibitor Targeted to Signaling Proteins of Innate Immune Responses. Proc. Natl. Acad. Sci. USA 2006, 103, 12481–12486. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Liu, T.; Zhang, L.; Joo, D.; Sun, S.C. NF-ΚB Signaling in Inflammation. Signal Transduct. Target. Ther. 2017, 2, 17023. [Google Scholar] [CrossRef] [Green Version]
  66. Zhang, B.; Wang, L.-S.; Zhou, Y.-H. Elevated MicroRNA-125b Promotes Inflammation in Rheumatoid Arthritis by Activation of NF-ΚB Pathway. Biomed. Pharmacother. 2017, 93, 1151–1157. [Google Scholar] [CrossRef]
  67. Ren, B.; Liu, J.; Wu, K.; Zhang, J.; Lv, Y.; Wang, S.; Liu, L.; Liu, D. TNF-α-Elicited MiR-29b Potentiates Resistance to Apoptosis in Peripheral Blood Monocytes from Patients with Rheumatoid Arthritis. Apoptosis 2019, 24, 892–904. [Google Scholar] [CrossRef]
  68. Thurlings, R.M.; Boumans, M.; Tekstra, J.; van Roon, J.A.; Vos, K.; van Westing, D.M.; van Baarsen, L.G.; Bos, C.; Kirou, K.A.; Gerlag, D.M.; et al. Relationship between the Type I Interferon Signature and the Response to Rituximab in Rheumatoid Arthritis Patients. Arthritis Rheum. 2010, 62, 3607–3614. [Google Scholar] [CrossRef]
  69. van Baarsen, L.G.M.; Wijbrandts, C.A.; Rustenburg, F.; Cantaert, T.; van der Pouw Kraan, T.C.T.M.; Baeten, D.L.; Dijkmans, B.A.C.; Tak, P.P.; Verweij, C.L. Regulation of IFN Response Gene Activity during Infliximab Treatment in Rheumatoid Arthritis Is Associated with Clinical Response to Treatment. Arthritis Res. Ther. 2010, 12, R11. [Google Scholar] [CrossRef] [Green Version]
  70. Raterman, H.G.; Vosslamber, S.; de Ridder, S.; Nurmohamed, M.T.; Lems, W.F.; Boers, M.; van de Wiel, M.; Dijkmans, B.A.C.; Verweij, C.L.; Voskuyl, A.E. The Interferon Type I Signature towards Prediction of Non-Response to Rituximab in Rheumatoid Arthritis Patients. Arthritis Res. Ther. 2012, 14, R95. [Google Scholar] [CrossRef] [Green Version]
  71. Toonen, E.J.M.; Gilissen, C.; Franke, B.; Kievit, W.; Eijsbouts, A.M.; den Broeder, A.A.; van Reijmersdal, S.V.; Veltman, J.A.; Scheffer, H.; Radstake, T.R.D.J.; et al. Validation Study of Existing Gene Expression Signatures for Anti-TNF Treatment in Patients with Rheumatoid Arthritis. PLoS ONE 2012, 7, e0033199. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Dennis, G.; Holweg, C.T.J.; Kummerfeld, S.K.; Choy, D.F.; Setiadi, A.F.; Hackney, J.A.; Haverty, P.M.; Gilbert, H.; Lin, W.Y.; Diehl, L.; et al. Synovial Phenotypes in Rheumatoid Arthritis Correlate with Response to Biologic Therapeutics. Arthritis Res. Ther. 2014, 16, R90. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Sellam, J.; Marion-Thore, S.; Dumont, F.; Jacques, S.; Garchon, H.J.; Rouanet, S.; Taoufik, Y.; Hendel-Chavez, H.; Sibilia, J.; Tebib, J.; et al. Use of Whole-Blood Transcriptomic Profiling to Highlight Several Pathophysiologic Pathways Associated with Response to Rituximab in Patients with Rheumatoid Arthritis: Data from a Randomized, Controlled, Open-Label Trial. Arthritis Rheumatol. 2014, 66, 2015–2025. [Google Scholar] [CrossRef] [Green Version]
  74. Sanayama, Y.; Ikeda, K.; Saito, Y.; Kagami, S.I.; Yamagata, M.; Furuta, S.; Kashiwakuma, D.; Iwamoto, I.; Umibe, T.; Nawata, Y.; et al. Prediction of Therapeutic Responses to Tocilizumab in Patients with Rheumatoid Arthritis: Biomarkers Identified by Analysis of Gene Expression in Peripheral Blood Mononuclear Cells Using Genome-Wide Dna Microarray. Arthritis Rheumatol. 2014, 66, 1421–1431. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Wright, H.L.; Thomas, H.B.; Moots, R.J.; Edwards, S.W. Interferon Gene Expression Signature in Rheumatoid Arthritis Neutrophils Correlates with a Good Response to TNFi Therapy. Rheumatology 2015, 54, 188–193. [Google Scholar] [CrossRef] [Green Version]
  76. Smith, S.L.; Eyre, S.; Yarwood, A.; Hyrich, K.; Morgan, A.W.; Wilson, A.G.; Isaacs, J.; Plant, D.; Barton, A. Biologics in Rheumatoid Arthritis Genetics Investigating CD11c Expression as a Potential Genomic Biomarker of Response to TNF Inhibitor Biologics in Whole Blood Rheumatoid Arthritis Samples. Arthritis Res. Ther. 2015, 17, 359. [Google Scholar] [CrossRef] [Green Version]
  77. Oswald, M.; Curran, M.E.; Lamberth, S.L.; Townsend, R.M.; Hamilton, J.D.; Chernoff, D.N.; Carulli, J.; Townsend, M.J.; Weinblatt, M.E.; Kern, M.; et al. Modular Analysis of Peripheral Blood Gene Expression in Rheumatoid Arthritis Captures Reproducible Gene Expression Changes in Tumor Necrosis Factor Responders. Arthritis Rheumatol. 2015, 67, 344–351. [Google Scholar] [CrossRef] [Green Version]
  78. Nakamura, S.; Suzuki, K.; Iijima, H.; Hata, Y.; Lim, C.R.; Ishizawa, Y.; Kameda, H.; Amano, K.; Matsubara, K.; Matoba, R.; et al. Identification of Baseline Gene Expression Signatures Predicting Therapeutic Responses to Three Biologic Agents in Rheumatoid Arthritis: A Retrospective Observational Study. Arthritis Res. Ther. 2016, 18, 159. [Google Scholar] [CrossRef] [Green Version]
  79. Muskardin, T.W.; Vashisht, P.; Dorschner, J.M.; Jensen, M.A.; Chrabot, B.S.; Kern, M.; Curtis, J.R.; Danila, M.I.; Cofield, S.S.; Shadick, N.; et al. Increased Pretreatment Serum IFN-β/α Ratio Predicts Non-Response to Tumour Necrosis Factor α Inhibition in Rheumatoid Arthritis. Ann. Rheum. Dis. 2016, 75, 1757–1762. [Google Scholar] [CrossRef]
  80. Teitsma, X.M.; Jacobs, J.W.G.; Mokry, M.; Borm, M.E.A.; Pethö-Schramm, A.; van Laar, J.M.; Bijlsma, J.W.J.; Lafeber, F.P.J. Identification of Differential Co-Expressed Gene Networks in Early Rheumatoid Arthritis Achieving Sustained Drug-Free Remission after Treatment with a Tocilizumab-Based or Methotrexate-Based Strategy. Arthritis Res. Ther. 2017, 19, 170. [Google Scholar] [CrossRef] [Green Version]
  81. Spiliopoulou, A.; Colombo, M.; Plant, D.; Nair, N.; Cui, J.; Coenen, M.J.H.; Ikari, K.; Yamanaka, H.; Saevarsdottir, S.; Padyukov, L.; et al. Association of Response to TNF Inhibitors in Rheumatoid Arthritis with Quantitative Trait Loci for CD40 and CD39. Ann. Rheum. Dis. 2019, 78, 1055–1061. [Google Scholar] [CrossRef] [Green Version]
  82. Yokoyama-Kokuryo, W.; Yamazaki, H.; Takeuchi, T.; Amano, K.; Kikuchi, J.; Kondo, T.; Nakamura, S.; Sakai, R.; Hirano, F.; Nanki, T.; et al. Identification of Molecules Associated with Response to Abatacept in Patients with Rheumatoid Arthritis. Arthritis Res. Ther. 2020, 22, 46. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Derambure, C.; Dzangue-Tchoupou, G.; D’Agostino, M.A.; Lequerre, T.; Vittecoq, O. Gene Expression Regulated by Abatacept Associated with Methotrexate and Correlation with Disease Activity in Rheumatoid Arthritis. PLoS ONE 2020, 15, e0237143. [Google Scholar] [CrossRef] [PubMed]
  84. Oliver, J.; Nair, N.; Orozco, G.; Smith, S.; Hyrich, K.L.; Morgan, A.; Isaacs, J.; Wilson, A.G.; Barton, A.; Plant, D. Correction to: Transcriptome-Wide Study of TNF-Inhibitor Therapy in Rheumatoid Arthritis Reveals Early Signature of Successful Treatment. Arthritis Res. Ther. 2021, 23, 80. [Google Scholar] [CrossRef] [PubMed]
  85. Triaille, C.; Durez, P.; Sokolova, T.; Tilman, G.; de Bellefon, L.M.; Galant, C.; Coulie, P.; Lauwerys, B.R.; Limaye, N. Common Transcriptomic Effects of Abatacept and Other DMARDs on Rheumatoid Arthritis Synovial Tissue. Front. Immunol. 2021, 12, 724895. [Google Scholar] [CrossRef]
  86. Cai, Y.; Xu, K.; Aihaiti, Y.; Li, Z.; Yuan, Q.; Xu, J.; Zheng, H.; Yang, M.; Wang, B.; Yang, Y.; et al. Derlin-1, as a Potential Early Predictive Biomarker for Nonresponse to Infliximab Treatment in Rheumatoid Arthritis, Is Related to Autophagy. Front. Immunol. 2022, 12, 795912. [Google Scholar] [CrossRef]
  87. Sutcliffe, M.; Nair, N.; Oliver, J.; Morgan, A.W.; Isaacs, J.D.; Wilson, A.G.; Verstappen, S.M.M.; Viatte, S.; Hyrich, K.L.; Morris, A.P.; et al. Pre-Defined Gene Co-Expression Modules in Rheumatoid Arthritis Transition towards Molecular Health Following Anti-TNF Therapy. Rheumatology 2022, 61 (Suppl. S1), keac133-186. [Google Scholar] [CrossRef]
  88. Sekiguchi, N.; Kawauchi, S.; Furuya, T.; Inaba, N.; Matsuda, K.; Ando, S.; Ogasawara, M.; Aburatani, H.; Kameda, H.; Amano, K.; et al. Messenger Ribonucleic Acid Expression Profile in Peripheral Blood Cells from RA Patients Following Treatment with an Anti-TNF-Alpha Monoclonal Antibody, Infliximab. Rheumatology 2008, 47, 780–788. [Google Scholar] [CrossRef] [Green Version]
  89. Lequerré, T.; Gauthier-Jauneau, A.-C.; Bansard, C.; Derambure, C.; Hiron, M.; Vittecoq, O.; Daveau, M.; Mejjad, O.; Daragon, A.; Tron, F.; et al. Gene Profiling in White Blood Cells Predicts Infliximab Responsiveness in Rheumatoid Arthritis. Arthritis Res. Ther. 2006, 8, R105. [Google Scholar] [CrossRef] [Green Version]
  90. Tanino, M.; Matoba, R.; Nakamura, S.; Kameda, H.; Amano, K.; Okayama, T.; Nagasawa, H.; Suzuki, K.; Matsubara, K.; Takeuchi, T. Prediction of Efficacy of Anti-TNF Biologic Agent, Infliximab, for Rheumatoid Arthritis Patients Using a Comprehensive Transcriptome Analysis of White Blood Cells. Biochem. Biophys. Res. Commun. 2009, 387, 261–265. [Google Scholar] [CrossRef]
  91. Julià, A.; Erra, A.; Palacio, C.; Tomas, C.; Sans, X.; Barceló, P.; Marsal, S. An Eight-Gene Blood Expression Profile Predicts the Response to Infliximab in Rheumatoid Arthritis. PLoS ONE 2009, 4, e7556. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  92. Stuhlmüller, B.; Häupl, T.; Hernandez, M.M.; Grützkau, A.; Kuban, R.J.; Tandon, N.; Voss, J.W.; Salfeld, J.; Kinne, R.W.; Burmester, G.R. CD11c as a Transcriptional Biomarker to Predict Response to Anti-TNF Monotherapy with Adalimumab in Patients with Rheumatoid Arthritis. Clin. Pharmacol. Ther. 2010, 87, 311–321. [Google Scholar] [CrossRef]
  93. Mazzone, R.; Zwergel, C.; Artico, M.; Taurone, S.; Ralli, M.; Greco, A.; Mai, A. The Emerging Role of Epigenetics in Human Autoimmune Disorders. Clin. Epigenetics 2019, 11, 34. [Google Scholar] [CrossRef] [Green Version]
  94. Peres, R.S.; Liew, F.Y.; Talbot, J.; Carregaro, V.; Oliveira, R.D.; Almeida, S.L.; França, R.F.O.; Donate, P.B.; Pinto, L.G.; Ferreira, F.I.S.; et al. Low Expression of CD39 on Regulatory T Cells as a Biomarker for Resistance to Methotrexate Therapy in Rheumatoid Arthritis. Proc. Natl. Acad. Sci. USA 2015, 112, 2509–2514. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  95. Papalexi, E.; Satija, R. Single-Cell RNA Sequencing to Explore Immune Cell Heterogeneity. Nat. Rev. Immunol. 2018, 18, 35–45. [Google Scholar] [CrossRef]
  96. Cheng, L.; Wang, Y.; Wu, R.; Ding, T.; Xue, H.; Gao, C.; Li, X.; Wang, C. New Insights From Single-Cell Sequencing Data: Synovial Fibroblasts and Synovial Macrophages in Rheumatoid Arthritis. Front. Immunol. 2021, 12, 709178. [Google Scholar] [CrossRef]
  97. Yanagida, M.; Kawasaki, M.; Fujishiro, M.; Miura, M.; Ikeda, K.; Nozawa, K.; Kaneko, H.; Morimoto, S.; Takasaki, Y.; Ogawa, H.; et al. Serum Proteome Analysis in Patients with Rheumatoid Arthritis Receiving Therapy with Tocilizumab: An Anti-Interleukin-6 Receptor Antibody. Biomed. Res. Int. 2013, 2013, 607137. [Google Scholar] [CrossRef]
  98. Blaschke, S.; Rinke, K.; Maring, M.; Flad, T.; Patschan, S.; Jahn, O.; Mueller, C.A.; Mueller, G.A.; Dihazi, H. Haptoglobin-A1, -A2, Vitamin D-Binding Protein and Apolipoprotein C-III as Predictors of Etanercept Drug Response in Rheumatoid Arthritis. Arthritis Res. Ther. 2015, 17, 45. [Google Scholar] [CrossRef] [Green Version]
  99. Ling, S.F.; Nair, N.; Verstappen, S.M.M.; Barton, A.; Zucht, H.D.; Budde, P.; Schulz-Knappe, P.; Plant, D. Proteomic Analysis to Define Predictors of Treatment Response to Adalimumab or Methotrexate in Rheumatoid Arthritis Patients. Pharm. J. 2020, 20, 516–523. [Google Scholar] [CrossRef]
  100. Chen, J.; Li, S.; Ge, Y.; Kang, J.; Liao, J.; Du, J.; Tian, J.; Xie, X.; Li, F. ITRAQ and PRM-Based Proteomic Analysis Provides New Insights into Mechanisms of Response to Triple Therapy in Patients with Rheumatoid Arthritis. J. Inflamm. Res. 2021, 14, 6993–7006. [Google Scholar] [CrossRef]
  101. Aterido, A.; Cañete, J.D.; Tornero, J.; Blanco, F.; Fernández-Gutierrez, B.; Pérez, C.; Alperi-López, M.; Olivè, A.; Corominas, H.; Martínez-Taboada, V.; et al. A Combined Transcriptomic and Genomic Analysis Identifies a Gene Signature Associated with the Response to Anti-TNF Therapy in Rheumatoid Arthritis. Front. Immunol. 2019, 10, 1459. [Google Scholar] [CrossRef]
  102. Tasaki, S.; Suzuki, K.; Kassai, Y.; Takeshita, M.; Murota, A.; Kondo, Y.; Ando, T.; Nakayama, Y.; Okuzono, Y.; Takiguchi, M.; et al. Multi-Omics Monitoring of Drug Response in Rheumatoid Arthritis in Pursuit of Molecular Remission. Nat. Commun. 2018, 9, 2755. [Google Scholar] [CrossRef] [Green Version]
  103. Tao, W.; Concepcion, A.N.; Vianen, M.; Marijnissen, A.C.A.; Lafeber, F.P.G.J.; Radstake, T.R.D.J.; Pandit, A. Multiomics and Machine Learning Accurately Predict Clinical Response to Adalimumab and Etanercept Therapy in Patients With Rheumatoid Arthritis. Arthritis Rheumatol. 2021, 73, 212–222. [Google Scholar] [CrossRef]
  104. Yoosuf, N.; Maciejewski, M.; Ziemek, D.; Jelinsky, S.A.; Folkersen, L.; Müller, M.; Sahlström, P.; Vivar, N.; Catrina, A.; Berg, L.; et al. Early Prediction of Clinical Response to Anti-TNF Treatment Using Multi-Omics and Machine Learning in Rheumatoid Arthritis. Rheumatology 2022, 61, 1680–1689. [Google Scholar] [CrossRef]
  105. Julià, A.; Gómez, A.; López-Lasanta, M.; Blanco, F.; Erra, A.; Fernández-Nebro, A.; Mas, A.J.; Pérez-García, C.; Vivar, M.L.G.; Sánchez-Fernández, S.; et al. Longitudinal Analysis of Blood DNA Methylation Identifies Mechanisms of Response to Tumor Necrosis Factor Inhibitor Therapy in Rheumatoid Arthritis. EBioMedicine 2022, 80, 104053. [Google Scholar] [CrossRef]
Figure 1. -Omic approaches for a personalized therapy.
Figure 1. -Omic approaches for a personalized therapy.
Pharmaceutics 14 01648 g001
Table 1. GWAS and replication studies and response to DMARDs in rheumatoid arthritis patients.
Table 1. GWAS and replication studies and response to DMARDs in rheumatoid arthritis patients.
DrugClinical OutcomeDiscovery Cohort (n)Replication Cohort (n)Main Associationsp-ValueReplicationReference
MTX∆DAS28
at 6 months
1424429
177
NRG3 (rs168201)9.8 × 10−8-[14]
MTX∆DAS28
at 6 months
457-ARL14|PPM1L (rs7624766)3.9 × 10−7-[15]
ADA, ETN, IFX∆DAS28/EULAR
at 3 months
89-MAFB (rs6028945)2 × 10−7No [28,33][23]
IFNK (rs7046653)5 × 10−7-
ADA, ETN, IFX∆DAS28/EULAR
at 3 months
196-NR2F2 (rs10520789)6 × 10−7No [28][24]
PDE3A-SLCO1C1 (rs3794271)3.5 × 10−6Yes [25]
No [26,27,28]
ADA, ETN, IFX∆DAS28
at 6 months
566774EYA (rs17301249)6 × 10−5No [30][29]
Intergenic region (rs12081765)7 × 10−4No [30]
PDZD2 (rs1532269)7 × 10−4Yes [30]
Intergenic region (rs7305646)1 × 10−4No [30]
ADA, ETN, IFX∆DAS28
at 3 months
8821821Intergenic region (rs4411591)5 × 10−5Yes [35][31]
ADA, ETN, IFX∆DAS28/EULAR
from 3 to 12 months
2706290CD84 (rs6427528) 8 × 10−8Yes [35][32]
ADA, ETN, IFX∆DAS28 at 12 weeks375245MED15 (rs113878252) 1.2 × 10−8-[33]
ADA, ETN, IFX∆DAS28
at 6 months
487-rs2845112.5 × 10−8No [27][34]
CZPΔACR20, ΔDAS28 at week 6 and ΔDAS28 at week 12302-rs122873155.7 × 10−8-[36]
rs353550831.5 × 10−7
TCZ ∆DAS28 at 4 months
ACR20 at 6 months
1157526CD69 (rs11052877)4 × 10−3Yes [38][37]
ACR (American College of Rheumatology); ADA (Adalimumab); CZP (Certolizumab pegol); DAS (Disease Activity Score); ETN (Etanercept); IFX (Infliximab); MTX (Methotrexate); TCZ (Tocilizumab).
Table 2. Methylation studies and response to DMARDs in rheumatoid arthritis patients.
Table 2. Methylation studies and response to DMARDs in rheumatoid arthritis patients.
StudyPatients (n)DrugSampleTime SampleOutcome
Liu et al., 2011 [47]65ETN or ADA (26)Peripheral bloodBaselineDAS28
De Andrés et al., 2015 [42]19 early RA patientsMTX (12 GR, 1 MR, 2 NR) /17 controlsT, B, NK, monocytes and polymorphonuclear leukocytes from whole bloodBaseline/
after 1 month
DAS28 at 6 months
Plant et al., 2016 [48]72ETN (36 GR/36 PR)Whole bloodBaselineDAS28 at 3 months
Glossop et al., 2017 [49]46MTX, SSZ and HCQ (35 GR/11 NR)Whole bloodBaselineEULAR criteria at 6 months
Gosselt et al., 2019 [45]181MTX or MTX + SSZ + HCQ + corticosteroids (140 MR/GR and 41 NR)Whole blood leukocytesBaseline and at 3 monthsDAS28
Liebold et al., 2021 [43]16 RA
17 controls
MTX, sarilumab,
Janus kinase inhibitors (8 GR-MR/8 NR)
Peripheral blood and CD4+, CD8+, CD14+ and CD19+Baseline/3 monthsDAS28-ESR DAS28-CRP
Guderud et al., 2020 [44]72MTX (36 GR + 36 PR)Whole bloodBaseline and 4 weeks after MTXEULAR criteria at 6 months
Nair et al., 2020 [50]68MTX (34 GR + 34 PR)Whole bloodBaseline and 4 weeks after MTXDAS28 at 6 months
Gosselt et al., 2021 [46]69MTX or MTX + SSZ + HCQ + corticosteroidsWhole bloodBaselineDAS28 at 3 months
ADA (Adalimumab); CRP (C-Reactive Protein); DAS28 (Disease Activity Score in 28 joints); ETN (Etanercept); ESR (Erythrocyte Sedimentation Rate); GR (Good Response); HCQ (Hydroxicloroquine); MTX (Methotrexate); MR (Moderate Response); NR (Non-Responders); PR (Poor Response); SJC (Swollen Joint Count); SSZ (Sulphasalazine); TJC (Tender Joint Count).
Table 3. MiRNA studies and response to DMARDs in rheumatoid arthritis patients.
Table 3. MiRNA studies and response to DMARDs in rheumatoid arthritis patients.
StudyPatients (n)Drug (n Patients)Response (n)SampleTime SampleOutcome
Castro Villegas et al., 2015 [55]Study cohort (10); Replication cohort (85)ADA (15), ETN (25) and IFX (55)GR (85) NR (10)SerumBaseline and at 6 monthsEULAR criteria at 6 months
Bogunia-Kubik et al., 2016 [56]13anti-TNF-αNot specifiedSerumBefore and after 3 months of TNFiEULAR criteria 3 months
Liu et al., 2019 [57]Study cohort (16); Replication cohort (92)ETN8 GR; 8 NR 60 GR; 32 NRPBMCsBaselineEULAR criteria at week 24
Duroux-Richard et al., 2014 [58]32RTX16 GR; 16 NRBlood (16) and serum samples (32)BaselineEULAR criteria at 3 months
Cheng et al., 2020 [59]96IFX69 GR; 27 NRPeripheral blood samplesBaseline, 4, 12 and 24 weeksEULAR criteria at week 24
Krintel et al., 2015 [60]180ADA (89) or ADA ± MTX (91) EULAR criteria
Sode et al., 2018 [61]89ADA + MTX (89)ADA + MTX: 40 GR; 46 NRPlasmaBaseline and at 3 monthsACR/EULAR remission at 3 and 12 months
Ciechomska et al., 2018 [62]10ETN (7) ADA (3)Not specifiedSerumBaseline and after TNFiDAS28
Fernandez-Ruiz et al., 2018 [63]16Tofacitinib10 Remission; 6 No remissionBloodAt the first month after the last dose of tofacitinibRemission ((DAS28) <2.6 and no swollen joints)
ACR (American College of Rheumatology); ADA (Adalimumab); DAS (Disease Activity Score); ETN (Etanercept); GR (Good Response); IFX (Infliximab); MTX (Methotrexate); NR (Non-Responders); PBMC (Peripheral Blood Mononuclear Cell); RTX (Rituximab).
Table 4. Studies that Explore Transcriptomic Biomarkers and Response to DMARDs in RA Patients.
Table 4. Studies that Explore Transcriptomic Biomarkers and Response to DMARDs in RA Patients.
StudyPatients (n)Drug (n Patients)ResponseSampleTime SampleOutcome
Thurlings et al., 2010 [68]51RTXNot specifiedPBMCBaselineEULAR criteria at weeks 12 and 24
Van Baarsen et al., 2010 [69]33IFX(12 GR and 6 PR)Whole bloodBefore/after 1 monthDAS, tender joint counts and HAQ-DIs criteria at week 16
Raterman et al., 2012 [70]14RTX8 GR; 6 NRWhole bloodBaselineEULAR at week 24
Toonen et al., 2012 [71]42IFX (27) or ADA (15)(18 GR and 24 NR)Whole bloodBaselineEULAR criteria at week 14
Glynn Dennis et al., 2014 [72]GSE21537 dataset (62)IFXNot specifiedSynovialBaselineEULAR at week 16
Sellam et al., 2014 [73]68RTX44 GR; 24 NRPBMCsBaseline and 24 weeksEULAR at week 24
Sanayama et al., 2014 [74]40 + 20TCZGR 29 NR 8
GR 15 NR 5
PBMCBaseline, 3 and 6 monthsphysician’s global assessment and CDAI at 6 months
Wright et al., 2015 [75]20ADA (13), ETN (5), GOL (2)5 GR; 13 MR; 2 NRNeutrophilsBaselineDAS28 at week 12
Smith et al., 2015 [76]75ADA (25) ETN (50)ADA (16 GR, 9 NR) ETN (25 GR, 25 NR)Whole bloodBaselineEULAR criteria at month 3
Oswald et al., 2015 [77]240ABCoN (IFX 20, ETN 21, ADA 9)
GO-FURTHER (GOL 72) BATTER-UP (IFX 23, ETN 31, GOL 9, ADA 41, CZP 14)
ABConN (GR 35, NR 15)
GO-FURTHER (GR 66, NR 6)
BATTER-UP (GR 79, NR 39)
Whole bloodBaseline/after 14 weeksEULAR at 14 weeks
Nakamura et al., 2016 [78]209IFX (140), TCZ (38), or ABA (31)IFX (30% REM), TCZ (21.1% REM), ABA (22.6% REM)Whole bloodBaselineCDAI at 6 months
Wampler Muskardin et al., 2016 [79]Test cohort:32 (ABCoN) Validation cohort: 92 (TETRAD registry)IFX (19), ADA (37), ETN (60), GOL (2), CZP (6)Test cohort: 13 NR and 19 GR
Validation cohort: 44 NR, 30 MR and 18 GR
Serum sampleBaselineEULAR at 14 weeks EULAR at 12 weeks
Teitsma et al., 2017 [80]60MTX + TCZ (19) MTX + TCZ (24) MTX + TCZ (17)14 sDFR 5 control 13 sDFR 11 controls 10 sDFR 7 controlsWhole bloodBaselinesDFR
Sipiliopoulou et al., 2019 [81]2938 (BRAGGSS, DREAM, EIRA, ReAct, WTCCC, Other cohorts)IFX (792), ADA (1255), ETN (721), GOL (17), CZP (34)Not specifiedWhole bloodBaselineESR and SJC baseline and between 3–6 months after treatment
Yokoyama-Kokuryo et al., 2020 [82]45ABA ± MTX27 GR; 8 MR/NRWhole bloodBaseline and 6 monthsEULAR at 6 months
Derambure et al., 2020 [83]19ABA + MTX14 GR; 5 NRWhole bloodBaseline and 6 monthsDAS28-CRP at 6 months
Oliver et al., 2021 [84]70ADA50 GR; 20 NRWhole bloodBaseline and 3 monthsEULAR at 3 months
Triaille et al., 2021 [85]50MTX, ADA, ABA, RTX, TCZNot specifiedSynovial tissueBaseline and after 16 weeksEULAR at 16 weeks
Cai et al., 2022 [86]Test cohorts: GSE58795, GSE78068
Validation cohorts: GSE77298, GSE55457, and GSE89408 datasets
IFXGSE58795 36 GR; 23 NR GSE78068 42 GR; 98 NR GSE77298: 16 RA GSE55457: 13 RA GSE89408: 152 RAWhole blood SynoviumBaselineESR and CRP
Sutcliffe et al., 2022 [87]155 RAMS (MTX) BRAGGSS cohort (ADA)MTX (85) or ADA (70)42 GR; 43 NR 50 GR; 20 NRWhole bloodBaseline and at 4 weeks Baseline and at 3 monthsEULAR criteria after 3 months
EULAR criteria after 6 months
ABA (Abatacept); ADA (Adalimumab); CDAI (Clinical Disease Activity Index); CRP (C-Reactive Protein); CZP (Certolizumab pegol); DAS (Disease Activity Score); ESR (Erythrocyte Sedimentation Rate); ETN (Etanercept); GOL (Golimumab); GR (Good Response); HAQ-DIs (Health Assessment Questionnaire–Disability Index scores); IFX (Infliximab); MTX (Methotrexate); NR (Non-Responders); MR (Moderate Response); REM (Remission); RTX (Rituximab); NON-REM (Non-Remission); sDFR (sustained Drug-Free Remission); SJC (Swollen Joint Count); TCZ (Tocilizumab).
Table 5. Summary of the proteomics studies associated with the effectiveness to biological therapy in RA.
Table 5. Summary of the proteomics studies associated with the effectiveness to biological therapy in RA.
StudyPatients (n)Drug (n Patients)ResponseSampleTime SampleOutcome
Yanagida et al., 2013 [97]7TCZ7 MR or GRSerumBaseline, 4 and 8 weeksDAS28 (Baseline, 4 and 8 weeks)
Blaschke et al., 2015 [98]50ETN31 GR, 19 NRSerumBaseline/after 12 and 24 weeksEULAR criteria at 6 months
Ling et al., 2020 [99]286BRAGGSS cohort- ADA (150) RAMS cohort- MTX (136)ADA: 58 GR, 58 MR and 34 PR. MTX: 59 GR, 2 MR and 75 PRSerumBaselineEULAR criteria at 3 months (BRAGGSS)
EULAR criteria at 6 months (RAMS)
Chen et al., 2021 [100]20IFX + MTX + Leflunomide5 NR, 15 GRSerumBaseline and after 14 weeksEULAR criteria after 14 weeks
ADA (Adalimumab); DAS (Disease Activity Score); ETN (Etanercept); GR (Good Responders); IFX (Infliximab); MTX (methotrexate); MR (Moderate Responders); NR (Non-Responders); TCZ (Tocilizumab).
Table 6. Multi-omic studies and response to DMARDs in rheumatoid arthritis patients.
Table 6. Multi-omic studies and response to DMARDs in rheumatoid arthritis patients.
StudyPatients (n)Drug (n Patients)ResponseSampleTime SampleOutcomeOmics
Aterido et al., 2019 [101]11ADAGR 5 NR 3 MR 3Synovial biopsiesbaselineEULAR criteria at week 14Transcriptomic and genomic
Tasaki et al., 2018 [102]34 RA 35 controlsMTX (21), TCZ (13) and IFX (18)IR: MTX (11), TCZ (3) and IFX (8) GR: MTX (10), TCZ (10) and IFX (10)26 cell types from whole bloodBaseline, 4, 8, 12 and 24 weeksEULAR criteria at week 24 (DAS28-ESR)Transcriptomics and proteomics
Tao et al., 2021 [103]80ETN (38) or ADA (42)ADA (20 GR/MR; 18 NR) ETN (19GR/MR; 23 NR)CD14+, CD4+ from whole bloodBaselineEULAR criteria at 6 monthsTranscriptomics and epigenomics
Yoosuf et al., 2022 [104]39 femaleIFX (16), ADA (11), ETN (8), GOL (2), CZP (2)23 (GR/MR); 16 (NR)PBMCs from whole bloodBaseline and 3 monthsEULAR criteriaTranscriptomics and proteomics
Julià et al., 2022 [105]Discovery cohort (62)
Validation cohort (60)
ADA (5), CZP (10), ETN (34), GOL (12), IFX (1) ADA (7), CZP (13), ETN (31), GOL (9), IFX (0)Week 0: GR (50); NR (12); Week 12: GR (44); NR (7) Week 0: GR (49); NR (10); Week 12: GR (48); NR (10)Whole blood and neutrophils, macrophages, CD4+ T, CD8+ T, B and NK cellsBaseline and at week 12EULAR criteria at week 12Transcriptomics and epigenomics
ADA (Adalimumab); CZP (Certolizumab pegol); DAS (Disease Activity Score); ESR (Erythrocyte Sedimentation Rate); ETN (Etanercept); GOL (Golimumab); GR (Good Response); IFX (Infliximab); IR (Inadequate Response); MR (Moderate Response) MTX (Methotrexate); NR (Non-Responders); PBMC (Peripheral Blood Mononuclear Cell); TCZ (tocilizumab).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Madrid-Paredes, A.; Martín, J.; Márquez, A. -Omic Approaches and Treatment Response in Rheumatoid Arthritis. Pharmaceutics 2022, 14, 1648. https://doi.org/10.3390/pharmaceutics14081648

AMA Style

Madrid-Paredes A, Martín J, Márquez A. -Omic Approaches and Treatment Response in Rheumatoid Arthritis. Pharmaceutics. 2022; 14(8):1648. https://doi.org/10.3390/pharmaceutics14081648

Chicago/Turabian Style

Madrid-Paredes, Adela, Javier Martín, and Ana Márquez. 2022. "-Omic Approaches and Treatment Response in Rheumatoid Arthritis" Pharmaceutics 14, no. 8: 1648. https://doi.org/10.3390/pharmaceutics14081648

APA Style

Madrid-Paredes, A., Martín, J., & Márquez, A. (2022). -Omic Approaches and Treatment Response in Rheumatoid Arthritis. Pharmaceutics, 14(8), 1648. https://doi.org/10.3390/pharmaceutics14081648

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop