- Research
- Open access
- Published:
Chromosome conformation signatures define predictive markers of inadequate response to methotrexate in early rheumatoid arthritis
Journal of Translational Medicine volume 16, Article number: 18 (2018)
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).
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).
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.
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.
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%.
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.
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).
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
Liao KP, Alfredsson L, Karlson EW. Environmental influences on risk for rheumatoid arthritis. Curr Opin Rheumatol. 2009;21:279–83.
Bottini N, Firestein GS. Epigenetics in rheumatoid arthritis: a primer for rheumatologists. Curr Rheumatol Rep. 2013;15:372–6.
McInnes IB, Schett G. The pathogenesis of rheumatoid arthritis. N Engl J Med. 2011;365:2205–19.
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.
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.
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.
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.
O’Dell JR. Therapeutic strategies for rheumatoid arthritis. N Engl J Med. 2004;350:2591–602.
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.
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.
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.
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.
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.
Scherer HU, Dörner T, Burmester GR. Patient-tailored therapy in rheumatoid arthritis: an editorial review. Curr Opin Rheumatol. 2010;22:237–45.
Ranganathan P, McLeod HL. Methotrexate pharmacogenetics: the first step toward individualized therapy in rheumatoid arthritis. Arthritis Rheum. 2006;54:1366–77.
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.
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.
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.
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.
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.
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.
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.
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.
Klareskog L, Catrina AI, Paget S. Rheumatoid arthritis. Lancet. 2009;373:659–72.
Scott D, Wolfe F, Huizinga T. Rheumatoid arthritis. Lancet. 2010;376:1094–108.
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.
Ling JQ, Hoffman AR. Epigenetics of long-range chromatin interactions. Pediatr Res. 2007;61:11R–6R.
Deng W, Blobel GA. Do chromatin loops provide epigenetic gene expression states? Curr Opin Genet Dev. 2010;20:548–54.
Kadauke S, Blobel GA. Chromatin loops in gene regulation. Biochim Biophys Acta Gene Regul Mech. 2009;1789:17–25.
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.
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.
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.
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.
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.
Nakano K, Whitaker JW, Boyle DL, Wang W, Firestein GS. DNA methylome signature in rheumatoid arthritis. Ann Rheum Dis. 2013;72:110–7.
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.
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.
Viatte S, Plant D, Raychaudhuri S. Genetics and epigenetics of rheumatoid arthritis. Nat Rev Rheumatol. 2013;9:141–53.
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.
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.
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.
Boyan B, Giacomo C. Organization and function of the 3D genome. Nat Rev Genet. 2016;17:661–78.
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.
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.
Mukhopadhyay S, Ramadass AS, Akoulitchev A, Gordon S. Formation of distinct chromatin conformation signatures epigenetically regulate macrophage activation. Int Immunopharmacol. 2014;18:7–11.
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.
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.
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.
Crutchley JL, Wang XQD, Ferraiuolo MA, Dostie J. Chromatin conformation signatures: ideal human disease biomarkers? Biomark Med. 2010;4:611–29.
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.
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.
Smolen JS, Aletaha D. Scores for all seasons: SDAI and CDAI. Clin Exp Rheumatol. 2014;32:S75–9.
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.
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.
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.
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.
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.
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.
Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–2.
Rau R, Herborn G. Benefit and risk of methotrexate treatment in rheumatoid arthritis. Clin Exp Rheumatol. 2004;22:S83–94.
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.
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.
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.
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.
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.
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.
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.
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.
Kellner H. Targeting interleukin-17 in patients with active rheumatoid arthritis: rationale and clinical potential. Ther Adv Musculoskelet Dis. 2013;5:141–52.
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.
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.
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.
Toussirot E. The IL23/Th17 pathway as a therapeutic target in chronic inflammatory diseases. Inflamm Allergy Drug Targets. 2012;11:159–68.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Firestein GS, McInnes IB. Immunopathogenesis of rheumatoid arthritis. Immunity. 2017;46:183–96.
Ray P, Le Manach Y, Riou B, Houle TT. Statistical evaluation of a biomarker. Anesthesiology. 2010;112:1023–40.
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.
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.
Udalova IA, Mantovani A, Feldmann M. Macrophage heterogeneity in the context of rheumatoid arthritis. Nat Rev Rheumatol. 2016;12:472–85.
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
Consortia
Corresponding author
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.
About this article
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
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12967-018-1387-9