Skip to main content

Chromosome conformation signatures define predictive markers of inadequate response to methotrexate in early rheumatoid arthritis

Background

There is a pressing need in rheumatoid arthritis (RA) to identify patients who will not respond to first-line disease-modifying anti-rheumatic drugs (DMARD). We explored whether differences in genomic architecture represented by a chromosome conformation signature (CCS) in blood taken from early RA patients before methotrexate (MTX) treatment could assist in identifying non-response to DMARD and, whether there is an association between such a signature and RA specific expression quantitative trait loci (eQTL).

Methods

We looked for the presence of a CCS in blood from early RA patients commencing MTX using chromosome conformation capture by EpiSwitchâ„¢. Using blood samples from MTX responders, non-responders and healthy controls, a custom designed biomarker discovery array was refined to a 5-marker CCS that could discriminate between responders and non-responders to MTX. We cross-validated the predictive power of the CCS by generating 150 randomized groups of 59 early RA patients (30 responders and 29 non-responders) before MTX treatment. The CCS was validated using a blinded, independent cohort of 19 early RA patients (9 responders and 10 non-responders). Last, the loci of the CCS markers were mapped to RA-specific eQTL.

Results

We identified a 5-marker CCS that could identify, at baseline, responders and non-responders to MTX. The CCS consisted of binary chromosome conformations in the genomic regions of IFNAR1, IL-21R, IL-23, CXCL13 and IL-17A. When tested on a cohort of 59 RA patients, the CCS provided a negative predictive value of 90.0% for MTX response. When tested on a blinded independent validation cohort of 19 early RA patients, the signature demonstrated a true negative response rate of 86 and a 90% sensitivity for detection of non-responders to MTX. Only conformations in responders mapped to RA-specific eQTL.

Conclusions

Here we demonstrate that detection of a CCS in blood in early RA is able to predict inadequate response to MTX with a high degree of accuracy. Our results provide a proof of principle that a priori stratification of response to MTX is possible, offering a mechanism to provide alternative treatments for non-responders to MTX earlier in the course of the disease.

Background

Rheumatoid arthritis (RA) is an autoimmune inflammatory disorder driven by interacting genetic, epigenetic and environmental factors [1,2,3]. The diagnosis of RA prompts early initiation of methotrexate (MTX), the first choice of disease modifying anti-rheumatic drug (DMARD) as recommended by European League against Rheumatism (EULAR) and American College of Rheumatology (ACR) ‘treat-to-target’ strategy (target being remission or low disease activity). This approach has substantially improved outcomes in the last decade [4,5,6]. However, approximately 35–59% of patients do not achieve clinically meaningful responses after starting MTX [7]. Predicting response to MTX has been one of the main challenges in RA management for over two decades [8]. Delay of effective treatment has clinical implications as response to MTX is the most significant predictor of long-term outcome in RA [9,10,11]. The ability to identify biomarkers able to predict inadequate response has thus far proven challenging [7, 12,13,14,15]. Several attempts have been made to develop diagnostics able to predict MTX-response [16,17,18,19,20,21,22,23] (Additional file 1: Table S1) however, most studies have failed.

RA is dependent on the interaction of genetic and environmental factors [24, 25]. Epigenetic markers are closely linked to transcriptional regulation and may reflect pathogenic changes associated with disease states [26,27,28,29,30]. The first evidence suggesting that epigenetic mechanisms may play a role in autoimmune diseases came from studies performed by Richardson et al. looking at the effect of DNA methyltransferase inhibitor 5-azacytidine [31]. Other studies have also reported aberrant methylation in RA [32,33,34,35,36,37,38]. Interestingly, MTX treatment also plays a role in epigenetic regulation in RA [39,40,41]. Gene expression in mammals is regulated by non-coding elements that can impact physiology and disease; the principle mechanism of regulation is through the architectural status at both coding and non-coding genomic regions [26]. Under external perturbations, genomic regions can alter their 3-dimensional structure as a way of functional regulation of gene expression [42]. These structural changes can be measured by EpiSwitch, a high-throughput molecular technique that analyzes the spatial organization of genomic loci in a cell [43,44,45]. As multiple genomic regions contribute to phenotypic differences through changes in genomic architecture [26], this approach allows for the development of a chromosomal conformation signature (CCS) of alterations in genomic architecture between two states (disease vs. non-disease, pre-treatment vs. post-treatment). The evaluation of long range chromosome interactions has provided useful blood-based biomarkers in oncology and other non-rheumatic diseases [44,45,46,47,48]. CCS offer a stable, binary readout of cellular states and represent an emerging class of biomarkers [49]. Here, we used a chromosomal architecture based approach to predict the response to MTX by developing a blood-based classifier based on a CCS. We hypothesized that interrogation of genomic architectural changes in early RA patients would define a functional endotype able to guide clinical decision making.

Methods

Patient population

The SERA study is a national prospective inception cohort of patients with new RA or undifferentiated arthritis, together with age-gender matched controls [50, 51]. We used 67 SERA patients for the discovery/validation phases (plus 34 healthy controls) and an additional 19 patients from an independent cohort, who were enrolled in a blinded validation study (demographic and clinical characteristics are shown in Additional file 1: Tables S2–S4). Participants fulfilled the 2010 ACR/EULAR RA criteria at recruitment and commenced monotherapy with MTX. Importantly, patients were treated to standard of care according to SIGN Guidelines (http://www.sign.ac.uk/pdf/sign123.pdf) reflecting the true to life, population wide cohort design of SERA. Patients were treated with standard therapy, which included intra-articular or intra-muscular glucocorticoids, offered as required in the absence of low dose oral glucocorticoids. Clinical data to calculate a range of composite clinical outcomes were obtained at baseline and after 6 months of MTX treatment (Additional file 1: Tables S2–S4). DAS (Disease Activity Score) response was calculated on the basis of either DAS28 (ESR) or DAS28 (CRP) change to meet the following criteria for response. Responders (R) were defined as those achieving DAS28 remission (DAS28 < 2.6) or a good response (DAS28 improvement of > 1.2 and DAS28 ≤ 3.2) at 6 months. Non-responders (NR) were defined as patients who had no improvement in DAS28 (≤ 0.6) by 6 months. For consistency, we then calculated a CDAI (Clinical Disease Activity Index) response for each patient and used this outcome measure to ensure comparability of clinical response across all patients in the analysis [52, 53]. Baseline peripheral blood samples with EDTA were collected and centrifuged to generate a buffy layer, stored at − 80 °C. Local ethics committees approved the SERA protocol and all participants gave informed consent before enrollment.

Identification of markers by EpiSwitchâ„¢ and probe design

A pattern recognition algorithm was used to annotate the human genome for sites with the potential to form long-range chromosome conformations. The proprietary EpiSwitch pattern recognition software, based on Bayesian-modeling, provides a probabilistic score that a region is involved in long-range chromatin interactions [44, 45]. Sequences from 123 gene loci (Additional file 1: Table S5), selected based on a systematic literature review for genes that have been associated with RA [35, 54,55,56,57], were processed through the pattern recognition software to generate a list of the 13,322 chromosomal interactions to be screened for association with response to MTX in RA. Sixty 60-mer oligonucleotide probes were designed to interrogate these potential interactions and uploaded as a custom array to the Agilent SureDesign website. Importantly, each probe was present in quadruplicate on the EpiSwitchâ„¢ microarray. To subsequently evaluate a potential CCS, nested PCR (EpiSwitchâ„¢ PCR) was performed using sequence-specific oligonucleotides designed using Primer3. Oligonucleotides were tested for specificity using oligonucleotide specific BLAST.

Preparation of genomic templates

Chromatin with intact chromosome conformations from 50 µl of each blood sample was extracted using the EpiSwitch assay following the manufacturer’s instructions (Oxford BioDynamics Ltd.) [43, 45, 46]. The methods used to perform an EpiSwitch microarray and EpiSwitch PCR detection are described in detail in Supplementary Methods.

Statistical analysis

GraphPad Prism and SPSS were used for all statistical analyses of clinical data. The R statistical environment was used to analyze the CCS data. The Fisher’s exact test (for categorical variables), the t test for independent samples (for continuous normally distributed variables), and the Mann–Whitney U test (for continuous variables without normal distribution) were used to identify differences. The level of statistical significance was set at p ≥ 0.05, and all tests were 2-sided. A full description of the statistical methods used to identify the CCS can be found in Additional file 1: Additional Methods.

Mapping of CCS locations to expressed quantitative trait locus (eQTL)

The genomic locations of the five markers (IFNAR1, IL-21R, IL-23, CXCL13 and IL-17A) in the CSS were compared with the Walsh et al. [58] eQTLs data using a Window function within the Bedtools suit of genomic analysis functions [59]. Only eQTLs that were within a genomic window of 100 bp to 2 Kbp were considered when mapping to the CCS locations.

Results

Patient characteristics

Demographic and clinical characteristics was captured for the early RA patient groups (Stages I, II and III) at the time of recruitment and at 6 months following MTX therapy (Additional file 1: Table S2–S4). A total of 86 early RA patients and 34 healthy controls (HC) were included. 67 were used to develop the CCS, and an independent, blinded cohort of 19 patients used in validation (Fig. 1). The biomarker discovery pipeline comprised 3 stages (Fig. 1). Stage I, the discovery phase was performed using 8 early RA patients (4 MTX-NR and 4 MTX-R) and 4 HC. Stage II, conducted to test the CCS involved 59 early RA patients (30 MTX-NR and 29 MTX-R) and 30 HC. Stage III validation used an independent, blinded cohort of 19 early RA patients (Fig. 1).

Fig. 1
figure 1

Study design. a Samples used for biomarker discovery and validation. A subset of patient samples from the SERA Inception Cohort (86 RA patients and 34 HC) were used to discover and validate the MTX response biomarker. b Workflow for discovery and validation of the epigenetic stratifying biomarker. For Stage I and Stage II biomarker discovery and testing, clinical samples from MTX-treatment naïve patients were provided after confirmation of response by SDAI. In Stage I, an initial panel of 13,322 potential biomarkers was refined to a 5-marker chromosomal conformation signature (CCS). In Stage II, the disease specific nature of the biomarker panel was confirmed by stratification against HC and further testing was performed on 59 RA patients (29 R and 30 NR) and 30 HC. Final validation of the biomarker panel was done on an independent, blinded cohort of 19 RA patients

Initial definition of the CCS in early RA

To identify an initial set of common epigenetic profiles in early RA patients, 123 genetic loci (Additional file 1: Table S5) previously associated with RA [35, 54,55,56,57] were selected and annotated with chromosome conformation interactions predicted using the EpiSwitch in silico annotation package [44,45,46]. The in silico prediction generated 13,322 high-confidence chromosome interaction candidates, with an average of 99 per loci (99 ± 64; mean ± SD) (Additional file 1: Table S5). These candidates were used to generate a bespoke discovery CCS array (Additional file 1: Table S5) to screen baseline whole blood samples from 4 MTX-R and 4 MTX-NR (Fig. 1; Additional file 1: Table S2), and 4 HC. Age, gender makeup, DAS28, baseline CDAI and MTX doses were similar in R and NR patients (Fig. 2a, b; Additional file 1: Table S2). We identified 922 statistically significant chromosomal interactions that differentiated R, NR and HC. Of the 922 lead interactions, 420 were associated with NR, 210 with R and 159 with HC. A stepwise hierarchical clustering approach reduced the number of significant interactions from 922 to a 30-marker profile that effectively stratified patients (Additional file 1: Figure S1, Table S6).

Fig. 2
figure 2

Clinical Data for Discovery and Testing Cohorts. Discovery cohort: (a, b). a Breakdown of disease severity by CDAI scores for the discovery cohort at Baseline and 6 months. b Change in CDAI scores between responders (R) and non-responders (NR) at Baseline and 6 months following MTX therapy. Testing cohort: (c, d). c Breakdown of disease severity by CDAI scores for the testing cohort at Baseline and 6 months. d Change in CDAI scores between responders (R) and non-responders (NR) at Baseline and 6 months following MTX therapy. Significant reductions in CDAI scores were seen between R and NR at Baseline compared to R at 6 months as well as between NR and R at 6 months (****p<0.01). For CDAI scores see Additional file 1: Tables S2, S3

Refinement of the CCS

The 30 interactions identified in the initial screen were narrowed to a smaller pool using a second SERA-derived early RA patient cohort of 59 patients (30 MTX-R and 29 MTX-NR) and 30 HCs (Fig. 2c, d; Additional file 1: Table S3). Employing a stepwise approach using the EpiSwitch PCR platform (Additional file 1: Figure S1, blood samples from 24 RA patients (12 R and 12 NR) were analyzed for the interactions that best differentiated between R and NR. This resulted in a set of 10 best discriminatory interactions (Additional file 1: Table S6). These ten interactions were used to probe blood samples from an additional set of 35 RA (18 R and 17 NR). Using logistical regression on the combined data from all patients used to refine the CCS (that is 30 R and 29 NR), we identified five conformations within the IFNAR1, IL-21R, IL-23, IL-17A, and CXCL13 loci which provided the final CCS for response to MTX (Additional file 1: Figures S2–S6, Table S8). The regression coefficients and odds ratio of the logistic regression model are shown in Table 1.

Table 1 Coefficients of the logistic regression model for predicting efficacy of MTX monotherapy in baseline samples based on retrospective annotation at 6 months

The model provided a prediction probability score for R and NR, ranging from 0.0098 to 0.99 (0 = responder, 1 = non-responder). The cut-off values were set at ≤ 0.30 for R and ≥ 0.70 for NR. The score of ≤ 0.30 had a true positive rate (sensitivity) of 89% (95% confidence interval [95% CI] 72–98%), while a score of ≥ 0.70 had a true negative response rate (specificity) of 87% (95% CI 70–96%). The number of observed and predicted patients per response category is shown in Table 2. With the CCS classifier, 53 patients (90%) were correctly classified as either R or NR.

Table 2 Observed and predicted number of R and NR to MTX monotherapy at 6 months using the CCS classifier

Testing the performance of the CCS

To test the accuracy and performance robustness of the CCS classifier to discriminate non-responders to MTX in early RA prior to treatment, data from the 59 early RA patient cohort used to develop the classifier was split into 150 different training and test sets by random data re-sampling using fivefold cross-validation. While this analysis was not intended as a true validation since the classifier was tested on the same cohort used to develop it, this provided an interim measure of the classifiers ability to discriminate R and NR. The average area under the curve (AUC) of the model was 90.6% (95% CI 87–100%), with an average sensitivity for NR of 89.7% and an average specificity for R of 90.0% (Fig. 3a). To determine the predictive capability of the CCS classifier, the average model accuracy statistics were adjusted for population R/NR to MTX [60]. Using a 55% MTX response rate, the positive predictive value (PPV) was 88% while the negative predictive value (NPV) was 86.4%. When the response rate was adjusted to 60%, this decreased the PPV to 85.7% while increasing the NPV to 88.6%.

Fig. 3
figure 3

CCS performance on randomized subsets of testing cohort and validation cohort. a Receiver operating characteristic (ROC) plots for 150 runs of the EpiSwitchâ„¢ logistic classifier. Data for the 59 patient cohort was randomised 150 times using the WEKA sample randomisation function. This reorders the data prior to splitting in developing the training set, ensuring that the same starting point for the classifier is not used and allowing multiple accuracy calculations for the same data. The average area under the curve (AUC) for the 150 classifier runs was 0.90. The plot is the average ROC from the 150 test results. b Receiver operating characteristic (ROC) plots for the EpiSwitchâ„¢ logistic classifier run on the blinded validation cohort of 19 RA patients. The classifier had a sensitivity of 75.0% and specificity of 85.7% with an AUC of 0.91

Independent blinded validation

To validate the performance of the CCS, a blinded analysis was undertaken on an additional 19 early RA patients from an independent cohort selected from SERA (9 R and 10 NR) (Additional file 1: Table S4, Figure S8). This validation cohort had similar clinical characteristics to those used to develop the CCS with respect to age, gender, RA severity (as measured by CDAI, DAS28-CRP, and DAS28-ESR) and MTX dosage (Additional file 1: Tables S3, S4). The analysis of the validation cohort was done using a blinded analysis, where the actual response calls were not known until after the identification of R and NR using the CCS classifier. In this independent cohort, we obtained a PPV of 90% (95% CI 59–98%), and a NPV rate of 66.7% (95% CI 42–85%), sensitivity of 75.0% (95% CI 43–95%) and specificity of 85.7% (95% CI 42–97%). The area under the receiver operating characteristic curve (AUROC) was 0.91 (Fig. 3b). The number of observed and predicted patients per response category is shown in Table 3. The strong NPV statistics for the blinded samples confirms the accuracy and robust stratification power of the CCS for non-response to MTX prior to treatment.

Table 3 CCS classifier performance for predicting non-response to MTX in a blinded cohort of early RA patients

Biological and functional characterization of the genomic loci making up the CCS

The five genomic loci (IFNAR1, IL-21R, IL-23, IL-17A, and CXCL13) that make up the CCS fall within genes that encode for proteins involved in the immune response. Recently, eQTLs that overlap with enhancer elements have been identified in RA patients with moderate to severe disease and inadequate response to MTX [58]. eQTLs are genomic loci that contribute to variation in expression levels of mRNAs. We mapped the genomic locations of the five loci in the CCS to these previously published eQTLs [58]. The two genomic sites that corresponded to the junction of each CCS locus were mapped to eQTLs that were within 100, 200, 500, 1000 or 2000 base pairs (Additional file 1: Table S9). There was enrichment of eQTLs co-localized with the R-associated CCS loci (IFNAR1, IL-21R and IL-23). In contrast, no co-localization was observed in the NR-associated CCS loci (IL-17A and CXCL13) (Fig. 4).

Fig. 4
figure 4

Chromosome conformations in MTX responders map to RA eQTLs. Bedtools shot of CCS markers mapped within 200 bp of previously identified RA eQTLs. The overlapping CCS:eQTLs are highlighted in yellow. a The CCS regions associated with the IFNAR1 locus on chromosome 21 map to 6 eQTLs. b The CCS regions associated with the IL-21R locus on chromosome 16 map to 21 eQTLs. c The CCS regions associated with the IL-23 locus on chromosome 12 map to 4 eQTLs. d The CCS regions associated with the IL-17A locus on chromosome 6 do not map to any eQTLs. e The CCS regions associated with the CXCL13 locus on chromosome 4 do not map to any eQTLs

Discussion

To address an unmet clinical need of predicting those RA patients who will not respond to MTX, we analyzed the CCS in whole blood taken at baseline from SERA-RA patients using a well-established method of analyzing long-range chromatin interactions [43,44,45]. We identified a 5-marker panel consisting of chromosomal conformations in the genomic loci of IFNAR1, IL-21R, IL-23, IL-17A, and CXCL13 that could identify R and NR to MTX with 90% sensitivity in an independent blinded validation cohort. This could facilitate earlier access to more effective therapies, thus avoiding disease progression, unnecessary exposure to potentially toxic drugs and diminished quality of life.

Historically, identifying predictive biomarkers for MTX response has been difficult [7, 12,13,14,15]. While clinical predictors of RA disease are well established [61, 62], they do not correlate well with response to treatment at the individual level [13]. 35–59% of patients do not achieve clinically meaningful response after starting MTX [7, 63]. In the SERA study, only 30% of patients responded to MTX monotherapy [50]. While the results presented here provide a proof-of-concept and further validation is warranted, the EpiSwitch technology has several attractive features from the standpoint of a biomarker that can be used clinically [49]. First, it requires a very small amount of blood sample (typically 50 μl). Second, it utilizes an established laboratory methodology and readouts (qPCR). Last, the turnaround time is short (~ 8 h). Thus, once further validated, the approach described here faces low barriers to clinical adoption. The CCS defined a signature that suggests epigenomic control over loci with a central role in the IL-17/IL-23 axis, with the two most informative long-range chromatin interactions for predicting MTX-NR coming from IL-17A and CXCL13 loci. In agreement with our finding, earlier reports have suggested that IL-23, IFNAR1 and IL-21R are predictors of positive response in other diseases [64,65,66], while IL-17A and CXL13 are predictors of poor outcome and increased disease severity [67,68,69,70,71]. Several studies support the notion that IL-23/Th17 axis drive inflammation in chronic diseases and perhaps plays a role in the response to immunomodulating drugs [72,73,74,75]. Interestingly, no loci previously associated with MTX metabolism were implicated in our study, indicating that the response to MTX was independent of how the drug is known to be metabolized [20, 76].

Our results are consistent with previous data reporting the presence of regulatory elements (i.e. eQTL) only in the context of inflammatory diseases and not in healthy controls [77].

Recently, Walsh et al. reported that eQTL mapped from RA are particularly enriched in enhancer regions of disease related cell types such as T and B cells [58]. Previous studies have recognized that the eQTL are highly specific to different leucocyte subsets [78]. To explore the concordance of CCS data with existing genetic regulatory datasets, we mapped the genomic locations of the regions in the CCS against reported RA-specific eQTL [58]. This revealed a high level of co-localization of R-specific loci (IFNR1, IL-21R, and IL-23) to RA-specific eQTL. Notably, this level of concordance was not observed in NR loci (IL-17 and CXCL-13). The increased association of RA eQTLs with CCS regions observed in R, but not in NR, suggest dysregulation at the level of regulatory 3D genome architecture. It has been reported that eQTL are associated with expression of mRNA transcripts, with concomitant effects on protein levels [79,80,81,82]. This leaves open the possibility that the association of eQTL with CCS observed in RA-R result from phenotypic consequences due to the effects of mRNA and protein expression levels. However, the impact of eQTL on protein levels remains poorly understood.

It is now clear that the developing immune response is influenced by genetic and epigenetic factors [41]. IFNAR1, IL-21R and IL-23 loci have been reported to play a key role in the pathogenesis of RA [83,84,85] and Th17 cells are implicated in pathogenesis especially in the pre-RA phase [68, 72, 73]. Whether the CCS differences we observed represent changes impacting the pre-RA phase or are acquired during the early phase of RA is at present unclear. We hypothesize that the different genomic architecture observed for R and NR might reflect differences in epigenetic host responses to early pathogenetic events, or particular environmental exposures. We anticipate that the association observed between CCS and eQTL in RA patients may be important to understand the heterogeneity of the response observed between individuals. Furthermore, eQTLs present in active inflammatory diseases can disappear after treatment [77]. It would therefore be useful to determine whether the concordance of the CCS and eQTL observed in our study in treatment naïve RA responders will be linked only to the disease or change after MTX treatment. Moreover, it would be interesting to see whether the reported eQTL are linked to the disease state and/or the level of inflammation [68]. Our data indicate that the mapping of the QTL can reveal an altered biological status in R and NR, however further studies are needed to confirm this.

The samples used in this study came from the SERA study, a large, multi-institutional investigational program designed to identify predictive markers of RA [50]. Patients enrolled in the SERA study were carefully characterized for clinical phenotype (the majority of the patients were Caucasians, non-smokers with seropositive established disease), longitudinal follow-up of outcomes and blood samples were stored following a Standard Operating Procedure. This combination of clinical rigor and quality assurance of the inputs to the CCS are particular strengths of this study. An additional strength of the study is the approach that was used to generate the CCS. We focused on the discovery and establishment of a molecular signature using an approach informed by current biological knowledge. We evaluated a network of loci that had plausible pathophysiological relevance in RA via synovial pathogenesis studies, GWAS association and postulated MTX pharmacogenetics [35, 54,55,56,57, 86]. The step-wise selection of biomarkers used in this study, coupled with the strict separation of discovery and validation cohorts, was performed to prevent marker and model over-fitting. The robust statistical properties of the CCS classifier are another advantage of the approach presented here. Using epigenetic markers, which provide a binary readout (presence or absence) and are stable in isolated whole blood, provides high efficiency stratification. Statistical power analysis confirmed that the sample size was adequate for the development and evaluation of the signature, a critical step in biomarker development and aligned with successful development of companion diagnostics in limited size cohorts in other indications [87]. A further strength to our findings and applicability in the clinical setting is the ability to identify this molecular signature in the whole blood by using a drop of blood.

Some of the caveats associated with our findings are the relatively modest sample size and the cellular heterogeneity present in the whole blood. While the sample size for the validation cohort was not as large as previous studies seeking to identify a biomarker for MTX non-response (Additional file 1: Table S1), the use of chromosome conformations as a readout, which can generate robust signatures in smaller cohorts, provide confidence in the approach [88, 89]. The heterogeneity of cell populations is an issue that has implications for any analysis in whole blood. In RA, it is known that there is significant heterogeneity in cell populations present in the blood of RA patients [90]. While outside the scope of the current study, future studies that look at the CCS in distinct cell types in whole blood may shed greater light on the similarities and/or differences exhibited within individual populations. A final caveat of our study was limited in that it could not determine whether the observed epigenetic marks are causal or consequential (secondary to the inflammatory response). Future studies looking at larger patient sets as well as the inclusion of individuals treated with other DMARDs, including biologics, are warranted.

Conclusions

The ability to identify individuals who are unlikely to respond to a specific therapeutic agent has the potential to not only reduce side-effects to treatments that will provide no benefit but also accelerate their journey through the treatment algorithm and thus increase likelihood of a positive treatment outcome. Given the high worldwide prevalence of RA, a stratifying signature for MTX response offers direct and practical benefits. The validated classifier described here, which is based on five conditional binary biomarkers detected in the blood, offers cost-effective detection by an established ISO-certified industrial platform and the practicality that this method to stratify patients has the potential to be routinely available within clinical practice. This could be a first step toward personalized medicine in RA.

Abbreviations

GWAS:

genome wide association study

EULAR:

European League Against Rheumatism

DAS:

Disease Activity Index

ESR:

erythrocyte sedimentation rate

CRP:

C reactive protein

CDAI:

Clinical Disease Activity Index

R:

responders

NR:

non-responders

HC:

healthy controls

RA:

rheumatoid arthritis

MTX:

methotrexate

CCS:

chromatin conformation signatures

SERA:

Scottish Early Rheumatoid Arthritis

TNR:

true negative response

DMARDS:

disease modifying anti-rheumatic drug

ACR:

American College of Rheumathologists

eQTL:

expression quantitative trait loci

References

  1. Liao KP, Alfredsson L, Karlson EW. Environmental influences on risk for rheumatoid arthritis. Curr Opin Rheumatol. 2009;21:279–83.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Bottini N, Firestein GS. Epigenetics in rheumatoid arthritis: a primer for rheumatologists. Curr Rheumatol Rep. 2013;15:372–6.

    Article  PubMed  Google Scholar 

  3. McInnes IB, Schett G. The pathogenesis of rheumatoid arthritis. N Engl J Med. 2011;365:2205–19.

    Article  CAS  PubMed  Google Scholar 

  4. Smolen JS, Aletaha D, Bijlsma JWJ, Breedveld FC, Boumpas D, Burmester G, et al. Treating rheumatoid arthritis to target: recommendations of an international task force. Ann Rheum Dis. 2010;69:631–7.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Singh JA, Saag KG, Bridges SL, Akl EA, Bannuru RR, Sullivan MC, et al. 2015 American College of rheumatology guideline for the treatment of rheumatoid arthritis. Arthritis Rheumatol. 2016;68:1.

    PubMed  Google Scholar 

  6. Smolen JS, Landewé R, Breedveld FC, Buch M, Burmester G, Dougados M, et al. EULAR recommendations for the management of rheumatoid arthritis with synthetic. Ann Rheum Dis. 2014;73:492–509.

    Article  CAS  PubMed  Google Scholar 

  7. Hider SL, Silman AJ, Thomson W, Lunt M, Bunn D, Symmons DPM. Can clinical factors at presentation be used to predict outcome of treatment with methotrexate in patients with early inflammatory polyarthritis? Ann Rheum Dis. 2009;68:57–62.

    Article  CAS  PubMed  Google Scholar 

  8. O’Dell JR. Therapeutic strategies for rheumatoid arthritis. N Engl J Med. 2004;350:2591–602.

    Article  PubMed  Google Scholar 

  9. Barrera P, van der Maas A, van Ede AE, Kiemeney BA, Laan RF, van de Putte LB, et al. Drug survival, efficacy and toxicity of monotherapy with a fully human anti-tumour necrosis factor-alpha antibody compared with methotrexate in long-standing rheumatoid arthritis. Rheumatology. 2002;41:430–9.

    Article  CAS  PubMed  Google Scholar 

  10. Tm F, Lunt M, Fu B, Bunn D, Symmons D. Early treatment with and time on first disease modifying anti-rheumatic drug predicts long-term function in patients with inflammatory polyarthritis. Ann Rheum Dis. 2009;69:689–95.

    Google Scholar 

  11. Bakker MF, Jacobs JWG, Welsing PMJ, Vreugdenhil SA, van Booma-Frankfort C, Linn-Rasker SP, et al. Early clinical response to treatment predicts 5-year outcome in RA patients: follow-up results from the CAMERA study. Ann Rheum Dis. 2011;70:1099–103.

    Article  CAS  PubMed  Google Scholar 

  12. Drouin J, Haraoui B. Predictors of clinical response and radiographic progression in patients with rheumatoid arthritis treated with methotrexate monotherapy. J Rheumatol. 2010;37:1405–10.

    Article  CAS  PubMed  Google Scholar 

  13. Hider SL, Buckley C, Silman AJ, Symmons DPBI. Factors influencing response to disease modifying antirheumatic drugs in patients with rheumatoid arthritis. J Rheumatol. 2005;32:11–6.

    PubMed  Google Scholar 

  14. Scherer HU, Dörner T, Burmester GR. Patient-tailored therapy in rheumatoid arthritis: an editorial review. Curr Opin Rheumatol. 2010;22:237–45.

    Article  PubMed  Google Scholar 

  15. Ranganathan P, McLeod HL. Methotrexate pharmacogenetics: the first step toward individualized therapy in rheumatoid arthritis. Arthritis Rheum. 2006;54:1366–77.

    Article  CAS  PubMed  Google Scholar 

  16. Seitz M, Zwicker M, Villiger PM. Pretreatment cytokine profiles of peripheral blood mononuclear cells and serum from patients with rheumatoid arthritis in different american college of rheumatology response groups to methotrexate. J Rheumatol. 2003;30:28–35.

    CAS  PubMed  Google Scholar 

  17. Stuhlmüller B, Häupl T, Hernandez MM, Grützkau A, Kuban R-J, Tandon N, et al. 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–21.

    Article  PubMed  Google Scholar 

  18. Dhir V, Sandhu A, Gupta N, Dhawan V, Sharma S, Sharma A. Low serum levels of myeloid progenitor inhibitory factor-1 predict good response to methotrexate in rheumatoid arthritis. ISRN Inflamm. 2013;2013:460469.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Ally M, Hodkinson B, Meyer PWA, Musenge E, Tintinger GR, Tikly M, et al. Circulating anti-citrullinated peptide antibodies, cytokines and genotype as biomarkers of response to disease-modifying antirheumatic drug therapy in early rheumatoid arthritis. BMC Musculoskelet Disord. 2015;16:130.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Dervieux T, Greenstein N, Kremer J. Pharmacogenomic and metabolic biomarkers in the folate pathway and their association with methotrexate effects during dosage escalation in rheumatoid arthritis. Arthritis Rheum. 2006;54:3095–103.

    Article  CAS  PubMed  Google Scholar 

  21. Kumagai K, Hiyama K, Oyama T, Maeda H, Kohno N. Polymorphisms in the thymidylate synthase and methylenetetrahydrofolate reductase genes and sensitivity to the low-dose methotrexate therapy in patients with rheumatoid arthritis. Int J Mol Med. 2003;11:593–600.

    CAS  PubMed  Google Scholar 

  22. Inoue S, Hashiguchi M, Takagi K, Kawai S, Mochizuki M. Preliminary study to identify the predictive factors for the response to methotrexate therapy in patients with rheumatoid arthritis. Yakugaku Zasshi. 2009;129:843–9.

    Article  CAS  PubMed  Google Scholar 

  23. Hayashi H, Fujimaki C, Daimon T, Tsuboi S, Matsuyama T, Itoh K. Genetic polymorphisms in folate pathway enzymes as a possible marker for predicting the outcome of methotrexate therapy in Japanese patients with rheumatoid arthritis. J Clin Pharm Ther. 2009;34:355–61.

    Article  CAS  PubMed  Google Scholar 

  24. Klareskog L, Catrina AI, Paget S. Rheumatoid arthritis. Lancet. 2009;373:659–72.

    Article  CAS  PubMed  Google Scholar 

  25. Scott D, Wolfe F, Huizinga T. Rheumatoid arthritis. Lancet. 2010;376:1094–108.

    Article  PubMed  Google Scholar 

  26. Tordini F, Aldinucci M, Milanesi L, Liò P, Merelli I. The genome conformation as an integrator of multi-omic data: the example of damage spreading in cancer. Front Genet. 2016;7:1–16.

    Article  Google Scholar 

  27. Ling JQ, Hoffman AR. Epigenetics of long-range chromatin interactions. Pediatr Res. 2007;61:11R–6R.

    Article  PubMed  Google Scholar 

  28. Deng W, Blobel GA. Do chromatin loops provide epigenetic gene expression states? Curr Opin Genet Dev. 2010;20:548–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Kadauke S, Blobel GA. Chromatin loops in gene regulation. Biochim Biophys Acta Gene Regul Mech. 2009;1789:17–25.

    Article  CAS  Google Scholar 

  30. Christova R, Jones T, Wu P-J, Bolzer A, Costa-Pereira AP, Watling D, et al. P-STAT1 mediates higher-order chromatin remodelling of the human MHC in response to IFNgamma. J Cell Sci. 2007;120:3262–70.

    Article  CAS  PubMed  Google Scholar 

  31. Richardson B, Kahn L, Lovett EJ, Hudson J. Effect of an inhibitor of DNA methylation on T cells. I. 5-Azacytidine induces T4 expression on T8+ T cells. J Immunol. 1986;137:35–9.

    CAS  PubMed  Google Scholar 

  32. Glossop JR, Emes RD, Nixon NB, Haworth KE, Packham JC, Dawes PT, et al. Genome-wide DNA methylation profiling in rheumatoid arthritis identifies disease-associated methylation changes that are distinct to individual T- and B-lymphocyte populations. Epigenetics. 2014;9:1228–37.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Kennedy A, Schmidt EM, Cribbs AP, Penn H, Amjadi P, Syed K, et al. A novel upstream enhancer of FOXP3, sensitive to methylation-induced silencing, exhibits dysregulated methylation in rheumatoid arthritis Treg cells. Eur J Immunol. 2014;44:2968–78.

    Article  CAS  PubMed  Google Scholar 

  34. Nile CJ, Read RC, Akil M, Duff GW, Wilson AG. Methylation status of a single CpG site in the IL6 promoter is related to IL6 messenger RNA levels and rheumatoid arthritis. Arthritis Rheum. 2008;58:2686–93.

    Article  PubMed  Google Scholar 

  35. Nakano K, Whitaker JW, Boyle DL, Wang W, Firestein GS. DNA methylome signature in rheumatoid arthritis. Ann Rheum Dis. 2013;72:110–7.

    Article  CAS  PubMed  Google Scholar 

  36. Liu Y, Aryee MJ, Padyukov L, Fallin MD, Hesselberg E, Runarsson A, et al. Epigenome-wide association data implicate DNA methylation as an intermediary of genetic risk in rheumatoid arthritis. Nat Biotechnol. 2013;31:142–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. de la Rica L, Urquiza JM, Gómez-Cabrero D, Islam ABMMK, López-Bigas N, Tegnér J, et al. Identification of novel markers in rheumatoid arthritis through integrated analysis of DNA methylation and microRNA expression. J Autoimmun. 2013;41:6–16.

    Article  PubMed  Google Scholar 

  38. Viatte S, Plant D, Raychaudhuri S. Genetics and epigenetics of rheumatoid arthritis. Nat Rev Rheumatol. 2013;9:141–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Cribbs AP, Kennedy A, Penn H, Amjadi P, Green P, Read JE, et al. Methotrexate restores regulatory T cell function through demethylation of the FoxP3 upstream enhancer in patients with rheumatoid arthritis. Arthritis Rheumatol. 2015;67:1182–92.

    Article  CAS  PubMed  Google Scholar 

  40. Ellis JA, Munro JE, Chavez RA, Gordon L, Joo JE, Akikusa JD, et al. Genome-scale case-control analysis of CD4+ T-cell DNA methylation in juvenile idiopathic arthritis reveals potential targets involved in disease. Clin Epigenetics. 2012;4:20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Kim Y-I, Logan JW, Mason JB, Roubenoff R. DNA hypomethylation in inflammatory arthritis: reversal with methotrexate. J Lab Clin Med. 1996;128:165–72.

    Article  CAS  PubMed  Google Scholar 

  42. Boyan B, Giacomo C. Organization and function of the 3D genome. Nat Rev Genet. 2016;17:661–78.

    Article  Google Scholar 

  43. Cao F, Fang Y, Tan HK, Goh Y, Choy JYH, Koh BTH, et al. Super-enhancers and broad H3K4me3 domains form complex gene regulatory circuits involving chromatin interactions. Sci Rep. 2017;7:2186.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Bastonini E, Jeznach M, Field M, Juszczyk K, Corfield E, Dezfouli M, et al. Chromatin barcodes as biomarkers for melanoma. Pigment Cell Melanoma Res. 2014;27:788–800.

    Article  CAS  PubMed  Google Scholar 

  45. Mukhopadhyay S, Ramadass AS, Akoulitchev A, Gordon S. Formation of distinct chromatin conformation signatures epigenetically regulate macrophage activation. Int Immunopharmacol. 2014;18:7–11.

    Article  CAS  PubMed  Google Scholar 

  46. Jakub JW, Grotz TE, Jordan P, Hunter E, Pittelkow M, Ramadass A, et al. A pilot study of chromosomal aberrations and epigenetic changes in peripheral blood samples to identify patients with melanoma. Melanoma Res. 2015;25:406–11.

    Article  PubMed  Google Scholar 

  47. Watanabe T, Ishihara K, Hirosue A, Watanabe S, Hino S, Ojima H, et al. Higher-order chromatin regulation and differential gene expression in the human tumor necrosis factor/lymphotoxin locus in hepatocellular carcinoma cells. Mol Cell Biol. 2012;32:1529–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Harismendy O, Notani D, Song X, Rahim NG, Tanasa B, Heintzman N, et al. 9p21 DNA variants associated with coronary artery disease impair interferon-gamma signalling response. Nature. 2011;470:264–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Crutchley JL, Wang XQD, Ferraiuolo MA, Dostie J. Chromatin conformation signatures: ideal human disease biomarkers? Biomark Med. 2010;4:611–29.

    Article  CAS  PubMed  Google Scholar 

  50. Dale J, Paterson C, Tierney A, Ralston SH, Reid DM, Basu N, et al. The Scottish Early Rheumatoid Arthritis (SERA) Study: an inception cohort and biobank. BMC Musculoskelet Disord. 2016;17:461.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Kronisch C, McLernon DJ, Dale J, Paterson C, Ralston SH, Reid DM, et al. Brief report: predicting functional disability: one-year results from the scottish early rheumatoid arthritis inception cohort. Arthritis Rheumatol. 2016;68:1596–602.

    Article  PubMed  Google Scholar 

  52. Smolen JS, Aletaha D. Scores for all seasons: SDAI and CDAI. Clin Exp Rheumatol. 2014;32:S75–9.

    Google Scholar 

  53. Aletaha D, Martinez-Avila J, Kvien TK, Smolen JS. Definition of treatment response in rheumatoid arthritis based on the simplified and the clinical disease activity index. Ann Rheum Dis. 2012;71:1190–6.

    Article  PubMed  Google Scholar 

  54. Kooloos WM, Wessels JA, van der Straaten T, Allaart CF, Huizinga TW, Guchelaar HJ. Functional polymorphisms and methotrexate treatment outcome in recent-onset. Pharmacogenomics. 2010;11:163–75.

    Article  CAS  PubMed  Google Scholar 

  55. Stahl EA, Raychaudhuri S, Remmers EF, Xie G, Eyre S, Thomson BP, et al. Genome-wide association study meta-analysis identifies seven new rheumatoid arthritis risk loci. Nat Genet. 2010;42:508–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Eyre S, Bowes J, Diogo D, Lee A, Barton A, Martin P, et al. High-density genetic mapping identifies new susceptibility loci for rheumatoid arthritis. Nat Genet. 2012;44:1336–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Okada Y, Wu D, Trynka G, Raj T, Terao C, Ikari K, et al. Genetics of rheumatoid arthritis contributes to biology and drug discovery. Nature. 2014;506:376–81.

    Article  CAS  PubMed  Google Scholar 

  58. Walsh AM, Whitaker JW, Huang CC, Cherkas Y, Lamberth SL, Brodmerkel C, et al. Integrative genomic deconvolution of rheumatoid arthritis GWAS loci into gene and cell type associations. Genome Biol. 2016;17:79.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Rau R, Herborn G. Benefit and risk of methotrexate treatment in rheumatoid arthritis. Clin Exp Rheumatol. 2004;22:S83–94.

    CAS  PubMed  Google Scholar 

  61. Polido-Pereira J, Vieira-Sousa E, Fonseca JE. Rheumatoid arthritis: what is refractory disease and how to manage it? Autoimmun Rev. 2011;10:707–13.

    Article  PubMed  Google Scholar 

  62. Kwoh CK, Anderson LG, Greene JM, Johnson DA, O’Dell JR, Robbins ML, et al. Guidelines for the management of rheumatoid arthritis: 2002 update—American College of Rheumatology subcommittee on rheumatoid arthritis guidelines. Arthritis Rheum. 2002;46:328–46.

    Article  CAS  Google Scholar 

  63. Hazlewood GS, Barnabe C, Tomlinson G, Marshall D, Devoe DJA, Bombardier C. Methotrexate monotherapy and methotrexate combination therapy with traditional and biologic disease modifying anti-rheumatic drugs for rheumatoid arthritis: a network meta-analysis. Cochrane Database Syst Rev. 2016.

  64. Morita K, Tanaka K, Saito S, Kitamura K, Morimoto M, Sekihara H. H S. Expression of interferon receptor genes (IFNAR1 and IFNAR2 mRNA) in the liver may predict outcome after interferon therapy in patients with chronic genotype 2a or 2b hepatitis C virus infection. J Clin Gastroenterol. 1998;26:135–40.

    Article  CAS  PubMed  Google Scholar 

  65. Hsu CS, Hsu SJ, Liu WL, Chen CL, Liu CJ, Chen PJ, et al. IL-21R gene polymorphisms and serum IL-21 levels predict virological response to interferon-based therapy in Asian chronic hepatitis C patients. Antivir Ther. 2013;18:599–606.

    Article  CAS  PubMed  Google Scholar 

  66. Massirer KB, Hirata MH, Silva AEB, Ferraz MLG, Nguyen NY, Hirata RDC. Interferon-alpha receptor 1 mRNA expression in peripheral blood mononuclear cells is associated with response to interferon-alpha therapy of patients with chronic hepatitis C. Braz J Med Biol Res. 2004;37:643–7.

    Article  CAS  PubMed  Google Scholar 

  67. Kirkham BW, Lassere MN, Edmonds JP, Juhasz KM, Bird PA, Lee CS, et al. Synovial membrane cytokine expression is predictive of joint damage progression in rheumatoid arthritis: a two-year prospective study (the DAMAGE study cohort). Arthritis Rheum. 2006;54:1122–31.

    Article  CAS  PubMed  Google Scholar 

  68. Chen D-Y, Chen Y-M, Chen H-H, Hsieh C-W, Lin C-C, Lan J-L. Increasing levels of circulating Th17 cells and interleukin-17 in rheumatoid arthritis patients with an inadequate response to anti-TNF-α therapy. Arthritis Res Ther. 2011;13:R126.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Kellner H. Targeting interleukin-17 in patients with active rheumatoid arthritis: rationale and clinical potential. Ther Adv Musculoskelet Dis. 2013;5:141–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Meeuwisse CM, Van Der Linden MP, Rullmann TA, Allaart CF, Nelissen R, Huizinga TW, et al. Identification of CXCL13 as a marker for rheumatoid arthritis outcome using an in silico model of the rheumatic joint. Arthritis Rheum. 2011;63:1265–73.

    Article  CAS  PubMed  Google Scholar 

  71. Dennis G Jr, Holweg CT, Kummerfeld SK, Choy DF, Setiadi AF, Hackney JA, et al. Synovial phenotypes in rheumatoid arthritis correlate with response to biologic therapeutics. Arthritis Res Ther. 2014;16:1–18.

    Article  Google Scholar 

  72. Sherlock JP, Taylor PC, Buckley CD. The biology of IL-23 and IL-17 and their therapeutic targeting in rheumatic diseases. Curr Opin Rheumatol. 2015;27:71–5.

    Article  CAS  PubMed  Google Scholar 

  73. Toussirot E. The IL23/Th17 pathway as a therapeutic target in chronic inflammatory diseases. Inflamm Allergy Drug Targets. 2012;11:159–68.

    Article  CAS  PubMed  Google Scholar 

  74. Zhang X, Hu J, Suo L, Yang Z, Xu T, Zhang Y. IL-17 and IL23 expression as a predictor of response to infliximab treatment in Crohn’s disease. Zhonghua nei ke za zhi. 2015;54:940–4.

    CAS  PubMed  Google Scholar 

  75. Yu C, Gong X, Yang Q, Lian J, Xu K, Ruan B, et al. The serum IL-23 level predicts the response to pegylated interferon therapy in patients with chronic hepatitis B. Liver Int. 2015;35:1549–56.

    Article  CAS  PubMed  Google Scholar 

  76. Kremer JM, Galivan J, Streckfuss A, Kamen B. Methotrexate metabolism analysis in blood and liver of rheumatoid arthritis patients: association with hepatic folate deficiency and formation of polyglutamates. Arthritis Rheum. 1986;29:832–5.

    Article  CAS  PubMed  Google Scholar 

  77. Peters JE, Lyons PA, Lee JC, Richard AC, Fortune MD, Newcombe PJ, et al. Insight into genotype–phenotype associations through eQTL mapping in multiple cell types in health and immune-mediated disease. PLoS Genet. 2016;12:e1005908.

    Article  PubMed  PubMed Central  Google Scholar 

  78. Fairfax BP, Makino S, Radhakrishnan J, Plant K, Leslie S, Dilthey A, et al. Genetics of gene expression in primary immune cells identifies cell type-specific master regulators and roles of HLA alleles. Nat Genet. 2012;44:502–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  79. Battle A, Khan Z, Wang SH, Mitrano A, Ford MJ, Pritchard JK, et al. Genomic variation. Impact of regulatory variation from RNA to protein. Science. 2015;347:664–7.

    Article  CAS  PubMed  Google Scholar 

  80. Wu L, Candille SI, Choi Y, Xie D, Jiang L, Li-Pook-Than J, et al. Variation and genetic control of protein abundance in humans. Nature. 2013;499:79–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Hause RJ, Stark AL, Antao NN, Gorsic LK, Chung SH, Brown CD, et al. Identification and validation of genetic variants that influence transcription factor and cell signaling protein levels. Am J Hum Genet. 2014;95:194–208.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Garge N, Pan H, Rowland MD, Cargile BJ, Zhang X, Cooley PC, et al. Identification of quantitative trait loci underlying proteome variation in human lymphoblastoid cells. Mol Cell Proteom. 2010;9:1383–99.

    Article  CAS  Google Scholar 

  83. Fragoulis GE, Siebert S, McInnes IB. Therapeutic targeting of IL-17 and IL-23 cytokines in immune-mediated diseases. Annu Rev Med. 2016;67:337–53.

    Article  CAS  PubMed  Google Scholar 

  84. Rao DA, Gurish MF, Marshall JL, Slowikowski K, Fonseka CY, Liu Y, et al. Pathologically expanded peripheral T helper cell subset drives B cells in rheumatoid arthritis. Nature. 2017;542:110–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. de Jong T, Lubbers J, Turk S, Vosslamber S, Mantel E, Bontkes H, et al. The type I interferon signature in leukocyte subsets from peripheral blood of patients with early arthritis: a major contribution by granulocytes. Arthritis Res Ther. 2016;18:165.

    Article  PubMed  PubMed Central  Google Scholar 

  86. Firestein GS, McInnes IB. Immunopathogenesis of rheumatoid arthritis. Immunity. 2017;46:183–96.

    Article  CAS  PubMed  Google Scholar 

  87. Ray P, Le Manach Y, Riou B, Houle TT. Statistical evaluation of a biomarker. Anesthesiology. 2010;112:1023–40.

    Article  PubMed  Google Scholar 

  88. Chatterjee N, Shi J, García-Closas M. Developing and evaluating polygenic risk prediction models for stratified disease prevention. Nat Rev Genet. 2016;14210:14205–10.

    Google Scholar 

  89. Stringer S, Wray NR, Kahn RS, Derks EM. Underestimated effect sizes in GWAS: fundamental limitations of single SNP analysis for dichotomous phenotypes. PLoS ONE. 2011;6:e27964.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  90. Udalova IA, Mantovani A, Feldmann M. Macrophage heterogeneity in the context of rheumatoid arthritis. Nat Rev Rheumatol. 2016;12:472–85.

    Article  CAS  PubMed  Google Scholar 

Download references

Authors’ contributions

CC, EH, AA, IBM and CSG conceived the study and wrote the paper. EH, IBM and CSG planned and reviewed experiments and analyzed the data. AR, HW and JG performed experiments. SERA Investigators enrolled patients and provided samples and clinical data. All authors read and approved the final manuscript.

Acknowledgements

The authors would like to thank members of OBD Reference Facility Benjamin Foulkes, Chloe Bird, Emily Corfield and Matthew Salter for expedient processing of clinical samples on the EpiSwitch™ platform and Magdalena Jeznach and Willem Westra for help with preparation of the manuscript. The study employed samples from the SERA Biobank used with permission and approval of the SERA Approval Group. We gratefully acknowledge the invaluable contribution of the clinicians and operating team in SERA. We would also like to thank Prof. Raju Kucherlapati (Harvard Medical School), and Prof. Jane Mellor (Oxford Univ.), Prof. John O’Shea (National Institute of Health) and Prof. John Isaacs (New Castle Univ.) for their independent and critical review of our study.

A list of Scottish Early Rheumatoid Arthritis (SERA) inception cohort investigators is provided in Additional file 1: Additional Note.

Competing interests

CC is a former employee at Pfizer Inc. and has no other competing financial interests. ASR, JG, EH and AA are full-time employees at Oxford Biodynamics Limited and have no other competing financial interests. IBM has received consultancy fees and grant funding from Pfizer. IBM has received honorarium from AstraZeneca, BMS, Abbvie, UCB and Roche. CSG receives consultancy fees from AstraZeneca. The remaining authors report no competing financial interests.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Local ethics committees approved the SERA protocol and all participants gave informed consent before enrollment.

Funding

This work was funded by Oxford BioDynamics.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Consortia

Corresponding author

Correspondence to Claudio Carini.

Additional file

Additional file 1: Note.

Scottish Early Rheumatoid Arthritis inception cohort (SERA). Table S1. Previous studies on prediction of MTX response in RA. Table S2. Patient characteristics—discovery cohort. Table S3. Patient characteristics—test cohort. Table S4. Patient characteristics—blinded validation cohort. Additional Methods. Statistical analysis and details of assays used to generate the CCS. Table S5. Selected genes for EpiSwitch Array. Figure S1. Design for discovery and validation of epigenetic stratifying biomarker signature for SERA patients. Table S6. Stepwise marker selection. Table S7. 65 Selected genes from EpiSwitch Array analysis. Table S8. 12 Selected genes from EpiSwitch PCR. Figure S2. Graphical representation of the genomic co-ordinates of the CXCL13 CCS marker associated with non-response to MTX. Figure S3. Graphical representation of the genomic co-ordinates of the IL- 17A CCS marker associated with non-response to MTX. Figure S4. Graphical representation of the genomic co-ordinates of the IFNAR1 CCS marker associated with non-response to MTX. Figure S5. Graphical representation of the genomic co-ordinates of the IL- 21R CCS marker associated with non-response to MTX. Figure S6. Graphical representation of the genomic co-ordinates of the IL-23 CCS marker associated with non-response to MTX. Figure S7. Clinical characteristics of independent blinded validation cohort. Table S9. CCS markers associated with eQTLs.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Carini, C., Hunter, E., Scottish Early Rheumatoid Arthritis Inception cohort Investigators. et al. Chromosome conformation signatures define predictive markers of inadequate response to methotrexate in early rheumatoid arthritis. J Transl Med 16, 18 (2018). https://doi.org/10.1186/s12967-018-1387-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-018-1387-9

Keywords