- Research article
- Open Access
A genetic risk score composed of rheumatoid arthritis risk alleles, HLA-DRB1 haplotypes, and response to TNFi therapy – results from a Swedish cohort study
Arthritis Research & Therapy volume 18, Article number: 288 (2016)
To prevent debilitating and irreversible joint damage, rheumatoid arthritis (RA) is often treated with tumor necrosis factor inhibitor (TNFi), but many patients do not respond to this costly therapy. Few predictors for response are known, and it has been proposed that genetic factors which influence the development of RA may also influence disease severity and response to therapy. Several previous studies have attempted to confirm this but results remain inconclusive. We expand on previous studies by including more RA risk alleles, and maximize power by combining them into a genetic risk score.
We linked genotyped RA patients from the Epidemiological Investigation of Rheumatoid Arthritis study to the Swedish Rheumatology Quality Register, identifying patients who started a TNFi as their first biological disease-modifying anti-rheumatic drug, with a return visit within 2–8 months after treatment start (N = 867). We calculated risk scores from 76 established RA risk SNPs, and four HLA-DRB1 amino acid positions, and tested whether risk scores or individual genetic risk factors could predict the European League Against Rheumatism (EULAR) response.
We found no association between any of the risk scores or HLA-DRB1 haplotypes and EULAR response, neither overall nor stratified by anti-citrullinated protein/peptide antibody (ACPA) status. When evaluating each of the 76 SNPs, we found that the number of SNPs presenting significant associations was not higher than expected by chance (5/76 SNPs had p < 0.05 in ACPA-positive RA, 4/76 in ACPA-negative RA).
Overall, known RA risk SNPs do not predict response to TNFi, either individually or when combined into a risk score. This does not support the hypothesis that genes influencing RA onset would also influence its prognosis and treatment response.
Tumor necrosis factor inhibitor (TNFi) therapy represents an important breakthrough in the treatment of the chronic inflammatory joint disease rheumatoid arthritis (RA). Although the therapeutic utility of TNFi is well documented, the drugs are costly and patients display substantial heterogeneity in treatment response: approximately 30% of patients discontinue therapy in the first year, a decision often made due to lack of effect or adverse events .
In order to personalize treatment, substantial efforts have been made to identify factors predicting response to TNFi. Environmental or clinical factors such as smoking, gender, age, baseline disability, or presence of autoantibodies account for only a small proportion of the variability in patient response [2, 3]. Unfortunately, genome-wide association studies (GWASs) have so far not revealed any clear evidence of predictive genetic markers, except for a large collection of nominally significant markers [4–8], or a few markers which, despite having approached acceptable levels of significance, are too weak to inform any clinical decisions (such as rs6427528 on CD84 with p = 8 × 10–8 in the etanercept subsets of patients) .
Before the advent of GWASs, studies of treatment response commonly focused on individual candidate markers. This remains a viable strategy for conserving power by reducing the number of simultaneous tests, if good candidates can be found. Candidate genes in this context have mainly come in two types: those involved in TNF metabolism or mechanistic pathways, and those associated with the onset of RA. Several investigations with limited samples (less than 130 subjects) have been conducted for the –308G > A polymorphism in the TNF gene [10–18], and found that the GG genotype is associated with a better response to TNFi treatment [11, 16–18]. This finding, however, was neither supported by a following study with a larger sample  nor replicated in subsequent genome-wide interrogations [4–7, 9].
Several studies have focused instead on individual genes known to be associated with risk for RA, including HLA-DRB1, PTPN22, IL, and FCGR, but no statistically significant associations were found [2, 19–22], except for one study finding two copies of HLA-DRB1 shared epitope (SE) significantly associated with improved ACR50 response  and a recent study investigating amino acid positions, rather than the SNPs per se, finding that valine in amino acid position 11 in HLA-DRB1 (which is outside the well-described SE positions) was associated with radiological progression and response to TNFi .
GWASs have provided researchers with an expanded list of RA risk genes. In a study of the then 31 identified RA risk loci, only one SNP (rs10919563) on PTPRC/CD45 was associated with European League Against Rheumatism (EULAR) response and DAS28 changes among TNFi starters, with reasonable replications and acceptable sample sizes (ca. 1200 in each study) [25, 26]. With some possible exceptions, including the HLA region, this observation may suggest that most alleles contributing to the development of RA do not influence treatment response when the disease has been established, and would be in line with the finding that a family history of RA (a proxy of genetic liability to develop RA) does not predict RA TNFi treatment response . It could also be argued, however, that previous studies were underpowered, and only focused on a subset of the now identified genetic markers for RA.
We hypothesized that many alleles associated with the development of RA are also important in predicting TNFi response, albeit with small individual effects, and that the previous failure to show this was due to lack of power. To extend previous research, we included more RA risk alleles and maximized power by combining these multiple polymorphisms into a single parameter—a genetic risk score. We analyzed these genetic markers (all currently known RA risk SNPs tagged by the Immunochip platform, and genome-wide SNPs) and several HLA-DRB1 amino acids (positions 11, 13, 71, and 74, and haplotypes defined by these positions), grouped into scores as well as individually, to evaluate how well they predict response to TNFi therapy among patients with RA who used TNFi as their first biological disease-modifying anti-rheumatic drug (bDMARD).
We performed a cohort study in prospectively recorded data by linking all participants in the Epidemiological Investigation of Rheumatoid Arthritis (EIRA) incident case–control study, who had been genotyped with the Immunochip array, to the Swedish Rheumatology Quality Register (SRQ), identifying patients starting TNFi therapy as their first bDMARD and their response to this treatment.
Epidemiological Investigation of Rheumatoid Arthritis
EIRA is a population-based case–control study initiated in 1996. Cases were recruited from all rheumatology providers within defined areas in Sweden, within 1 year of symptom onset and initial visit to a rheumatologist. At baseline, participants completed a self-administrated questionnaire and provided blood samples for serologic (anti-citrullinated protein/peptide antibodies (ACPA)) and genetic examinations. A total of 5043 EIRA subjects (all were recruited until 2009) were available on Immunochip genotypes; after quality control, 4830 were eligible for amino acid imputation. Imputation on amino acid positions 11, 13, 71, and 74 as well as the haplotypes based on the four positions was finally successfully performed for 4726 participants (2785 patients and 1941 controls). We subsequently linked the EIRA patients with the SRQ to further identify the target patient population of the current study. All participants consented to be involved in the study.
Swedish Rheumatology Quality Register
The SRQ is a profession-driven, web-based national quality register engaging both patients and rheumatologists. This clinical register records longitudinal data entered by the patient and the treating rheumatologist at each visit. The nationwide coverage is good for patients with newly diagnosed RA (about 80%), and is excellent for RA patients treated with bDMARDs (about 90%). Of the aforementioned 2785 EIRA patients with genetic information, 2576 were registered in the SRQ (92%); 895 were registered as starting any of the five TNFi agents (etanercept, infliximab, adalimumab, certolizumab pegol, and golimumab) as their first bDMARD. Finally, a total of 867 subjects (653 ACPA-positive patients and 165 ACPA-negative patients) who had a valid visit registered in the SRQ within 7 days of starting therapy were included in the current study.
We collected genetic information on the target patient population in four forms. First, we identified RA patients with the classical HLA-DRB1 SE, genotyped as described previously . Briefly, carriers of any of the HLA-DRB1*01 (but not HLA-DRB1*0103), HLA-DRB1*04 (HLA-DRB1*0401, HLA-DRB1*0404, HLA-DRB1*0408), and HLA-DRB1*10 alleles were defined as exposed; carriers of none of these alleles were defined as unexposed.
Second, we collected data on 76 RA risk SNPs as described previously (Additional file 1: Table S1) . Briefly, based on the results of a meta-analysis of Immunochip data for 11,475 cases and 15,870 controls of European ancestry (which identified 46 RA risk SNPs), as well as a trans-ethnic meta-GWAS of 29,880 cases and 73,758 controls of European and Asian ancestry (which provided an additional 49 RA risk SNPs/loci), we managed to obtain data using our Immunochip material for all 46 SNPs from the first study, and an additional 14 SNPs from the second study; of the remaining 35 SNPs in the second study, proxy SNPs with r 2 > 0.6 were identified on 16 of them, resulting in a total of 76 SNPs. We constructed a weighted genetic risk score (GRS) by summing the alleles for each individual, weighted by the logarithm of the published odds ratio (OR) in the corresponding meta-GWAS, according to the following equation:
Third, in addition to the 76 RA risk SNPs, we further made use of the whole Immunochip and calculated a series of GRSs based on genome-wide association with being an RA patient, estimated in the full EIRA material. Those GRSs included SNPs with decreasing magnitudes of effects—from the most stringent RA associations (genome-wide significance), to the moderate associations (e.g., p < 0.0005, p < 0.005, p < 0.05), to the least associated (all Immunochip markers)—and were weighted by the logarithm of ORs derived from the corresponding association analysis.
Finally, the four amino acid positions 11 (six polymorphisms), 13 (six polymorphisms), 71 (four polymorphisms), and 74 (five polymorphisms) at HLA-DRB1, as well as the haplotypes defined by those four positions, were imputed as described previously (Additional file 2: Table S2) . Briefly, the dosages of HLA amino acids were imputed from the Immunochip data with a publicly released reference panel generate by the Type 1 Diabetes Genetics Consortium using SNP2HLA software . Long-range haplotypes across the MHC were obtained, which were used to extract the corresponding haplotypes and amino acid residues. We calculated the GRS for amino acid positions by summing up the dosages of residues for each individual (each person could have two residues at each of the four positions), weighted by the logarithm of its reported association to develop RA (measured with the OR); similarly, we calculated the GRS for haplotypes by weighting each haplotype based on its reported OR , according to the following equations:
We collected the clinical characteristics for all of the eligible RA patients at treatment start and an evaluation visit, defined as a visit within 2–8 months after starting therapy, closest to 5 months if patients had multiple follow-up visits with valid DAS28 data. Disease activity was measured with the DAS28, visual analog scale (VAS) pain, Health Assessment Questionnaire (HAQ), erythrocyte sedimentation rate (ESR), C-reactive protein (CRP), swollen joint count (SJC), tender joint count (TJC), and VAS global health. We adopted two main outcomes based on these measurements: EULAR response, categorized as good/moderate vs none; and changes in the disease activity measures, calculated by subtracting the baseline values from the values at the evaluation visit, used in a continuous fashion (e.g., ΔDAS28, ΔCRP, etc.). In addition, to test the validity of GRS, we also estimated how the distribution of the GRS differed by RA status (RA vs controls).
We analyzed the GRS of SNPs, amino acids, and haplotypes, both as linear covariates (continuous) and categorical covariates (divided into four categories according to quartiles among controls). We also directly assessed individual prediction from each SNP, amino acid, and haplotype (i.e., in its original format: 0/1/2 copies). We assessed the association between achieving good/moderate vs no EULAR response and GRS through logistic regressions. We evaluated the association between the changes of disease activity and GRS using R-square (R 2) values from linear regressions, given that delta values are virtually normally distributed even without transformation. Two post-hoc sensitivity analyses were made, where the main analysis (categorical GRS as predictor of EULAR response) was made first excluding etanercept to assess whether the prediction was different when focusing on monoclonal antibody therapies, and second restricting the time window for the evaluation visit to 2–5 months (60–150 days) to assess whether the broad time window masked a specific association to response at that time point. Finally, to illustrate the validity of the RA GRS in this dataset, its association with RA risk (comparing the RA cases with controls) was estimated by logistic regressions. All analyses were performed by adjusting for age, gender, and five principal components (PCs, to correct for population stratification ) as well as stratified by ACPA status.
Clinical presentations of RA patients
Among the 867 RA patients who started TNFi as their first bDMARD, 75% were ACPA-positive and 74% were females. At the start of TNFi treatment, we identified comparable clinical presentations between the two subsets, except for a lower ESR (24.2 vs 29.0, p = 0.0068) and a slightly higher TJC (8.9 vs 7.5, p = 0.0074) among ACPA-negative patients. At follow-up, the improvement in ACPA-positive patients and ACPA-negative patients was similar for all outcome measures (Table 1).
GRS and RA risk
We found that not all of the GRSs followed normal distributions—the risk score on position 11 was slightly skewed, and on position 13 was collapsed into three levels (Additional file 3: Table S3)—indicating a need to analyze the amino acids individually. We examined the associations between GRS and RA risk. As expected, both the linear and categorical GRS only consistently increased ACPA-positive RA risk but not ACPA-negative RA risk (Table 2 and Additional file 4: Table S4).
GRS and TNFi treatment response
In contrast, we found no association between GRS (either from the 76 SNPs, the four amino acids, or the haplotypes) and good/moderate EULAR response in RA overall or in ACPA-positive RA. Neither did the primary RA genetic risk factor, SE, predict EULAR response (Table 3). Post-hoc sensitivity analysis showed virtually identical results when excluding etanercept (Additional file 5: Table S5). Post-hoc sensitivity analysis restricting the evaluation time window to months 2–5 did not find an association with the SNP-based or haplotype-based GRS, but did find a borderline significant reduced EULAR response associated with SE alleles (any SE allele OR = 0.69 (0.49–0.97), p = 0.0347, full table in Additional file 6: Table S6). The second quartile of GRS at amino acid position 13 appeared to be significantly associated with improved EULAR response in ACPA-negative RA, possibly due to low numbers in this stratum (Additional file 7: Table S7). We then evaluated each of the 76 SNPs, each of the residues at the amino acid positions, and each of the haplotypes associated with EULAR response separately. The number of SNPs presenting significant associations was not higher than expected by chance (5/76 SNPs had p < 0.05 in ACPA-positive RA, 4/76 in ACPA-negative RA) (Additional file 8: Table S8 and Additional file 9: Figure S1). Similarly, despite the borderline significance displayed by some residues in ACPA-positive RA (glycine at position 11, serine and tyrosine at position 13) and in ACPA-negative RA (proline at position 11, arginine at position 13), none survived multiple testing corrections (Additional file 10: Table S9, 21 tests performed). In addition, none of the haplotypes, either based on four amino acids (11/13/71/74) or on three amino acids (11/71/74, 11 and 13 are in high LD), were significantly associated with good/moderate EULAR response in RA overall or ACPA-positive RA, except for a few suggestive significances in ACPA-negative RA (Table 4). We additionally performed the haplotype analyses among patients with high disease activity (baseline DAS28 > 5.1) but did not identify any significant associations (Additional file 11: Table S10). Furthermore, the sensitive analysis comparing patients with good EULAR response with those with no response (good vs no) did not reveal further evidence of associations (Additional file 12: Table S11).
We further evaluated the performance of GRS and SE in the changes of disease activity measures. As shown in Table 5 and Additional file 13: Table S12, the only consistent association with nominal significance lay in ΔHAQ in ACPA-positive RA, where the GRS together with SE explained approximately 5% of variance (range: 3.4–5.4%), although this significance would not withstand correction for the number of tests (eight tests performed per amino acid position). Moreover, amino acid positions 11 and 13 appeared to explain a moderate proportion of variance for ΔESR (ca. 11%) in ACPA-negative RA, but this was not significant after corrections for multiple tests (eight tests performed per amino acid position). In addition, we found no evidence supporting significant associations (that withstood multiple corrections) between GRS (composed of 76 SNPs) and any of the baseline clinical characteristics (Additional file 14: Table S13).
Genome-wide polygenic risk score
The null findings already presented included markers that were RA risk SNPs achieving genome-wide significance in previous studies, and weaker RA risk alleles remain to be identified, which may have an effect on treatment response. Several lines of evidence have demonstrated that currently known RA susceptibility loci only account for a small proportion of ACPA-positive RA heritability, even less in ACPA-negative RA; and that by including additional common SNPs, the proportion of RA genetic liability could be explained to a larger extent . We hypothesized that the R 2 estimates would be improved by incorporating extra SNPs. We therefore calculated several GRSs including SNPs with different magnitudes of association to RA, demonstrated their predictive capacity with regards to RA risk, and then assessed the association with response to TNFi using these scores; the results are plotted in Fig. 1. We found no visible trends that R 2 would increase as numbers of incorporated markers increased.
By linking clinical data from the Swedish Rheumatology Register to genetic data from the EIRA study, we evaluated whether known RA susceptibility genes, HLA-DRB1 amino acids, and haplotypes, either individually or integrated into a risk score, predicted treatment response to TNFi therapy. Although all were strong predictors of RA, a high genetic risk score was not associated with good EULAR response in either overall RA or stratified by ACPA status.
Our results are in line with, but also extend, the previous GWAS and candidate gene approach findings, and indicate that although treatment response in RA has been reported to be somewhat heritable , the genetic variants that influence disease onset do not necessarily influence TNFi treatment response to the same extent. This is also supported by a recent study showing that family history of RA did not predict RA TNFi treatment response , although this on its own could have been due to the information on genetic risk contained in a family history being too slight among patients who are all at high genetic risk, as evident from them having already developed the disease.
The SNP that was found to be significantly associated with both RA risk and TNFi treatment response (PTPRC/CD45 rs10919563 at chromosome 1) in a previous larger sample  could unfortunately not be replicated in our study because of lack of proxies (the closest possible SNP to the PTPRC region available in our material is 100 kbp distant). Even though weak yet nominally significant associations were indeed identified on a handful of individual genetic markers in the current study, they were neither strong enough to withstand correction for multiple testing nor close enough to be clinically informative. Unfortunately, this has also been the case in previous studies. For example, the first TNFi treatment response GWAS performed in 89 RA patients by Liu et al.  provided a reference list of 16 candidate SNPs with suggestive significance; none was replicated in a subsequent separate study with slightly larger samples (n = 151) . Plant et al.  performed a multistage GWAS in 1285 RA patients and found seven genetic loci that might influence treatment response; none survived the two additional replication attempts [4, 7]. Similarly, the majority of the markers identified by Krintel et al.  in 196 Danish RA patients did not withstand replication in another 315 Spanish subjects . However, one SNP (rs3794271) at PDE3A-SLCO1C1 reached genome-wide significance in the meta-analysis combining the Danish and Spanish cohorts. Further investigations are needed for this pharmacogenetics biomarker of interest. GWASs performed after Liu et al. attempted to both identify new predictors and replicate previous findings, neither of which succeeded. This may be due to the small sample size for GWASs, where often 2–5 million SNPs were analyzed with an average sample size less than 1000. We attempted to increase power by combining many SNPs into a single score, yet still failed to reveal any significant associations. We considered that there may still be a genetic overlap among the as yet unidentified RA alleles, but found no evidence for this because we did not observe any apparent increment of variance explained in any of the disease activity measures when a genome-wide polygenic risk score model was performed. This might indicate that RA risk alleles do not necessarily have an effect on TNFi treatment, despite most of the current treatment target genes/pathways being involved in general inflammatory response; it may be reasonable to expect more and different biological pathways via which TNFi exerts its effect. In light of the limited statistical power in the current study, however, RA risk SNPs with modest effects in TNFi treatment response cannot be ruled out. One nominally significant finding may deserve mention: rs629326 (located 23.61 kb 5′ of TAGAP, involved in T-cell activation) was fairly strongly associated with EULAR response in our material, with OR = 0.31 (0.15–0.62), and FDR-adjusted (for 76 tests) p = 0.08. Further, the HLA-DRB1 SE alleles were associated with reduced EULAR response when restricting the time window to 2–5 months (OR = 0.69 (0.49–0.97), unadjusted p = 0.0347), although this was a post-hoc analysis and the p value would not remain significant after adjusting for multiple testing. Interestingly, specific HLA-DRB1 amino acids have been associated previously with TNFi treatment response, where valine at amino acid position 11 was reported to be associated with a smaller change in Larsen score, and improved EULAR response . We unfortunately lacked data on joint erosions. Our results for treatment response, however, did not immediately support the finding that the valine-containing VKA haplotype was associated with a good EULAR response (OR in  = 1.23 (1.06–1.43)), while the OR was 1.06 (0.58–1.94) for our data among ACPA-positive RA patients. The confidence intervals overlap greatly, and we are thus unable to either refute or confirm this previous finding. We did find a nominally significant association of residue serine (OR = 1.96 (1.14–3.37)) at position 13 and tightly linked with valine at position 11. This should be interpreted with caution, however, because it was not statistically significant after correcting for the multiplicity of tests. When the haplotype analyses were restricted to patients with high disease activity RA (baseline DAS > 5.1), to maximize comparability with the UK study, the association with valine further diminished (Additional file 11: Table S10).
When addressing continuous measures of treatment response, interestingly, all of the risk scores seemed to consistently explain a small yet significant proportion (5%) of the variance in HAQ changes in ACPA-positive RA. The clinical value of such results warrants further confirmation. It has been suggested that the modest influence of genetic effects on treatment response is due to the “composite” traits of the DAS28 score that most outcomes are based upon, and which relies on information from both subjective and objective measures ; and that using a well-defined phenotype could aid in disclosing the true genetic effects . We examined each component of the DAS28 without seeing any remarkable associations, except for a small proportion of variance in ESR change explained by amino acid risk score in ACPA-negative RA, which suffered markedly from limited power.
One strength of the current study is the accuracy of both exposures and outcomes. The genetic data were genotyped and quality controlled, and the amino acid data, although imputed, presented a high concordance rate as compared with true genotyping data (a 97.3% concordance rate for two-digit SE and 95.0% for four-digit SE) . We used clinically relevant outcome measures recommended by the European guidelines collected as part of clinical practice in an unselected manner. Because no environmental confounders could in practice influence the genetic markers, we performed the analyses without including adjustment for many covariates to preserve power. One limitation is insufficient power, which could cause false positive findings, and we consider the overall lack of association to be the main conclusion of this study. Another weakness of this study is the time point at which clinical response to TNFi therapy is determined. In a highly heterogeneous disease like RA, time of clinical assessment is crucial and can substantially modify response classification. However, our register-based data could not provide evaluation at time points as exact as those usually obtained in controlled trials, but rather in a time window, reflecting the Swedish clinical practice. We performed a post-hoc sensitive analysis restricting the time window for the evaluation visit to 2–5 months, and the results remained the same.
Our negative results imply that there is no strong evidence supporting a significant role of RA risk genes in the response to TNFi treatment; and that none of the SNPs, amino acids, or haplotypes examined in our study seem to be meaningful individual predictors of TNFi therapy, although weak associations cannot be ruled out. Although this suggest that future studies of TNFi response may want to focus on other genetic risk factors, combining genetic information connected to RA onset with clinical, nongenetic predictors, or perhaps their interaction, may also potentially prove valuable.
Neovius M, Arkema EV, Olsson H, Eriksson JK, Kristensen LE, Simard JF, Askling J, Group AS. Drug survival on TNF inhibitors in patients with rheumatoid arthritis comparison of adalimumab, etanercept and infliximab. Ann Rheum Dis. 2015;74(2):354–60.
Potter C, Hyrich KL, Tracey A, Lunt M, Plant D, Symmons DP, Thomson W, Worthington J, Emery P, Morgan AW, et al. Association of rheumatoid factor and anti-cyclic citrullinated peptide positivity, but not carriage of shared epitope or PTPN22 susceptibility variants, with anti-tumour necrosis factor response in rheumatoid arthritis. Ann Rheum Dis. 2009;68(1):69–74.
Saevarsdottir S, Wedren S, Seddighzadeh M, Bengtsson C, Wesley A, Lindblad S, Askling J, Alfredsson L, Klareskog L. Patients with early rheumatoid arthritis who smoke are less likely to respond to treatment with methotrexate and tumor necrosis factor inhibitors: observations from the Epidemiological Investigation of Rheumatoid Arthritis and the Swedish Rheumatology Register cohorts. Arthritis Rheum. 2011;63(1):26–36.
Krintel SB, Palermo G, Johansen JS, Germer S, Essioux L, Benayed R, Badi L, Ostergaard M, Hetland ML. Investigation of single nucleotide polymorphisms and biological pathways associated with response to TNFalpha inhibitors in patients with rheumatoid arthritis. Pharmacogenet Genomics. 2012;22(8):577–89.
Liu C, Batliwalla F, Li W, Lee A, Roubenoff R, Beckman E, Khalili H, Damle A, Kern M, Furie R, et al. Genome-wide association scan identifies candidate polymorphisms associated with differential response to anti-TNF treatment in rheumatoid arthritis. Mol Med. 2008;14(9–10):575–81.
Plant D, Bowes J, Potter C, Hyrich KL, Morgan AW, Wilson AG, Isaacs JD, Wellcome Trust Case Control Consortium; British Society for Rheumatology Biologics Register, Barton A. Genome-wide association study of genetic predictors of anti-tumor necrosis factor treatment efficacy in rheumatoid arthritis identifies associations with polymorphisms at seven loci. Arthritis Rheum. 2011;63(3):645–53.
Umicevic Mirkov M, Cui J, Vermeulen SH, Stahl EA, Toonen EJ, Makkinje RR, Lee AT, Huizinga TW, Allaart R, Barton A, et al. Genome-wide association analysis of anti-TNF drug response in patients with rheumatoid arthritis. Ann Rheum Dis. 2013;72(8):1375–81.
Honne K, Hallgrimsdottir I, Wu C, Sebro R, Jewell NP, Sakurai T, Iwamoto M, Minota S, Jawaheer D. A longitudinal genome-wide association study of anti-tumor necrosis factor response among Japanese patients with rheumatoid arthritis. Arthritis Res Ther. 2016;18(1):12.
Cui J, Stahl EA, Saevarsdottir S, Miceli C, Diogo D, Trynka G, Raj T, Mirkov MU, Canhao H, Ikari K, et al. Genome-wide association study and gene expression analysis identifies CD84 as a predictor of response to etanercept therapy in rheumatoid arthritis. PLoS Genet. 2013;9(3):e1003394.
Cuchacovich M, Ferreira L, Aliste M, Soto L, Cuenca J, Cruzat A, Gatica H, Schiattino I, Perez C, Aguirre A, et al. Tumour necrosis factor-alpha (TNF-alpha) levels and influence of -308 TNF-alpha promoter polymorphism on the responsiveness to infliximab in patients with rheumatoid arthritis. Scand J Rheumatol. 2004;33(4):228–32.
Fonseca JE, Carvalho T, Cruz M, Nero P, Sobral M, Mourao AF, Cavaleiro J, Ligeiro D, Abreu I, Carmo-Fonseca M, et al. Polymorphism at position -308 of the tumour necrosis factor alpha gene and rheumatoid arthritis pharmacogenetics. Ann Rheum Dis. 2005;64(5):793–4.
Kang CP, Lee KW, Yoo DH, Kang C, Bae SC. The influence of a polymorphism at position -857 of the tumour necrosis factor alpha gene on clinical response to etanercept therapy in rheumatoid arthritis. Rheumatology (Oxford). 2005;44(4):547–52.
Lee YH, Rho YH, Choi SJ, Ji JD, Song GG. Association of TNF-alpha -308 G/A polymorphism with responsiveness to TNF-alpha-blockers in rheumatoid arthritis: a meta-analysis. Rheumatol Int. 2006;27(2):157–61.
Marotte H, Pallot-Prades B, Grange L, Tebib J, Gaudin P, Alexandre C, Blond JL, Cazalis MA, Mougin B, Miossec P. The shared epitope is a marker of severity associated with selection for, but not with response to, infliximab in a large rheumatoid arthritis population. Ann Rheum Dis. 2006;65(3):342–7.
Martinez A, Salido M, Bonilla G, Pascual-Salcedo D, Fernandez-Arquero M, de Miguel S, Balsa A, de la Concha EG, Fernandez-Gutierrez B. Association of the major histocompatibility complex with response to infliximab therapy in rheumatoid arthritis patients. Arthritis Rheum. 2004;50(4):1077–82.
Mugnier B, Balandraud N, Darque A, Roudier C, Roudier J, Reviron D. Polymorphism at position -308 of the tumor necrosis factor alpha gene influences outcome of infliximab therapy in rheumatoid arthritis. Arthritis Rheum. 2003;48(7):1849–52.
Padyukov L, Lampa J, Heimburger M, Ernestam S, Cederholm T, Lundkvist I, Andersson P, Hermansson Y, Harju A, Klareskog L, et al. Genetic markers for the efficacy of tumour necrosis factor blocking therapy in rheumatoid arthritis. Ann Rheum Dis. 2003;62(6):526–9.
Seitz M, Wirthmuller U, Moller B, Villiger PM. The -308 tumour necrosis factor-alpha gene polymorphism predicts therapeutic response to TNFalpha-blockers in rheumatoid arthritis and spondyloarthritis patients. Rheumatology (Oxford). 2007;46(1):93–6.
Miceli-Richard C, Comets E, Verstuyft C, Tamouza R, Loiseau P, Ravaud P, Kupper H, Becquemont L, Charron D, Mariette X. A single tumour necrosis factor haplotype influences the response to adalimumab in rheumatoid arthritis. Ann Rheum Dis. 2008;67(4):478–84.
Canete JD, Suarez B, Hernandez MV, Sanmarti R, Rego I, Celis R, Moll C, Pinto JA, Blanco FJ, Lozano F. Influence of variants of Fc gamma receptors IIA and IIIA on the American College of Rheumatology and European League Against Rheumatism responses to anti-tumour necrosis factor alpha therapy in rheumatoid arthritis. Ann Rheum Dis. 2009;68(10):1547–52.
Kastbom A, Bratt J, Ernestam S, Lampa J, Padyukov L, Soderkvist P, Skogh T. Fcgamma receptor type IIIA genotype and response to tumor necrosis factor alpha-blocking agents in patients with rheumatoid arthritis. Arthritis Rheum. 2007;56(2):448–52.
Morales-Lara MJ, Canete JD, Torres-Moreno D, Hernandez MV, Pedrero F, Celis R, Garcia-Simon MS, Conesa-Zamora P. Effects of polymorphisms in TRAILR1 and TNFR1A on the response to anti-TNF therapies in patients with rheumatoid and psoriatic arthritis. Joint Bone Spine. 2012;79(6):591–6.
Criswell LA, Lum RF, Turner KN, Woehl B, Zhu Y, Wang J, Tiwari HK, Edberg JC, Kimberly RP, Moreland LW, et al. The influence of genetic variation in the HLA-DRB1 and LTA-TNF regions on the response to treatment of early rheumatoid arthritis with methotrexate or etanercept. Arthritis Rheum. 2004;50(9):2750–6.
Viatte S, Plant D, Han B, Fu B, Yarwood A, Thomson W, Symmons DP, Worthington J, Young A, Hyrich KL, et al. Association of HLA-DRB1 haplotypes with rheumatoid arthritis severity, mortality, and treatment response. JAMA. 2015;313(16):1645–56.
Cui J, Saevarsdottir S, Thomson B, Padyukov L, van der Helm-van Mil AH, Nititham J, Hughes LB, de Vries N, Raychaudhuri S, Alfredsson L, et al. Rheumatoid arthritis risk allele PTPRC is also associated with response to anti-tumor necrosis factor alpha therapy. Arthritis Rheum. 2010;62(7):1849–61.
Plant D, Prajapati R, Hyrich KL, Morgan AW, Wilson AG, Isaacs JD, Biologics in Rheumatoid Arthritis Genetics and Genomics Study Syndicate, Barton A. Replication of association of the PTPRC gene with response to anti-tumor necrosis factor therapy in a large UK cohort. Arthritis Rheum. 2012;64(3):665–70.
Frisell T, Saevarsdottir S, Askling J. Does a family history of RA influence the clinical presentation and treatment response in RA? Ann Rheum Dis. 2016;75(6):1120–5.
Lundstrom E, Kallberg H, Smolnikova M, Ding B, Ronnelid J, Alfredsson L, Klareskog L, Padyukov L. Opposing effects of HLA-DRB1*13 alleles on the risk of developing anti-citrullinated protein antibody-positive and anti-citrullinated protein antibody-negative rheumatoid arthritis. Arthritis Rheum. 2009;60(4):924–30.
Jiang X, Frisell T, Askling J, Karlson EW, Klareskog L, Alfredsson L, Kallberg H. To what extent is the familial risk of rheumatoid arthritis explained by established rheumatoid arthritis risk factors? Arthritis Rheumatol. 2015;67(2):352–62.
Kim K, Jiang X, Cui J, Lu B, Costenbader KH, Sparks JA, Bang SY, Lee HS, Okada Y, Raychaudhuri S, et al. Interactions between amino acid-defined major histocompatibility complex class II variants and smoking in seropositive rheumatoid arthritis. Arthritis Rheumatol. 2015;67(10):2611–23.
Jia X, Han B, Onengut-Gumuscu S, Chen WM, Concannon PJ, Rich SS, Raychaudhuri S, de Bakker PI. Imputing amino acid polymorphisms in human leukocyte antigens. PLoS One. 2013;8(6):e64683.
Raychaudhuri S, Sandor C, Stahl EA, Freudenberg J, Lee HS, Jia X, Alfredsson L, Padyukov L, Klareskog L, Worthington J, et al. Five amino acids in three HLA proteins explain most of the association between MHC and seropositive rheumatoid arthritis. Nat Genet. 2012;44(3):291–6.
Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006;38(8):904–9.
Stahl EA, Wegmann D, Trynka G, Gutierrez-Achury J, Do R, Voight BF, Kraft P, Chen R, Kallberg HJ, Kurreeman FA, et al. Bayesian inference analyses of the polygenic architecture of rheumatoid arthritis. Nat Genet. 2012;44(5):483–9.
Umicevic Mirkov M, Janss L, Vermeulen SH, van de Laar MA, van Riel PL, Guchelaar HJ, Brunner HG, Albers CA, Coenen MJ. Estimation of heritability of different outcomes for genetic studies of TNFi response in patients with rheumatoid arthritis. Ann Rheum Dis. 2015;74(12):2183–7.
Suarez-Gestal M, Perez-Pampin E, Calaza M, Gomez-Reino JJ, Gonzalez A. Lack of replication of genetic predictors for the rheumatoid arthritis response to anti-TNF treatments: a prospective case-only study. Arthritis Res Ther. 2010;12(2):R72.
Acosta-Colman I, Palau N, Tornero J, Fernandez-Nebro A, Blanco F, Gonzalez-Alvaro I, Canete JD, Maymo J, Ballina J, Fernandez-Gutierrez B, et al. GWAS replication study confirms the association of PDE3A-SLCO1C1 with anti-TNF therapy response in rheumatoid arthritis. Pharmacogenomics. 2013;14(7):727–34.
This work was supported by the Swedish Foundation for Strategic Research, Swedish Research Council, ALF funding collaboration between Karolinska Institute and Stockholm County, and the IMI-funded EU project BeTheCure. XJ had full access to all data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. The funders had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Availability of data and materials
Data will not be shared.
All authors contributed to the study design and interpretation on the findings. XJ performed all statistical analyses, interpreted the findings, and drafted the manuscript. TF contributed largely to the statistical analysis, helped to interpret the results, and revised the manuscript. JA, LA, LP, and SS participated in the design and coordination of the study, as well as revising the manuscript. SV contributed to the design of the study, actively discussed all of the findings with TF and XJ, and helped to revise the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Available upon request.
Ethical approval and consent to participate
This study was conducted in accordance with Declaration of Helsinki and was approved by the Regional Ethical Review Board in Stockholm. Signed informed consent was obtained from all patients before study enrollment.
Presenting characteristics of the 76 SNPs that are known RA risk SNPs. (DOCX 27 kb)
Presenting definitions of amino acid residues and haplotypes. (DOCX 18 kb)
Presenting the distribution of risk scores calculated based on 76 SNPs, each of the four amino acid positions, and the haplotypes. (DOCX 20 kb)
Presenting associations between amino acid GRS and RA risk. (DOCX 18 kb)
Presenting GRS predicting EULAR response, excluding etanercept. (DOCX 17 kb)
Presenting GRS predicting EULAR response, at 6–8 months. (DOCX 17 kb)
Presenting associations between amino acid GRS and EULAR response. (DOCX 18 kb)
Presenting associations between each of the 76 SNPs and EULAR response. (DOCX 37 kb)
Showing a QQ plot mapping 76 SNPs predicting EULAR response. (DOCX 47 kb)
Presenting associations between individual amino acid residues and EULAR response. (DOCX 23 kb)
Presenting associations between amino acid haplotypes and EULAR response, among patients with DAS28 > 5.1. (DOCX 18 kb)
Presenting associations between amino acid haplotypes and EULAR response, comparing patients with good response with patients with no response. (DOCX 19 kb)
Presenting variance in disease activity changes by risk scores. (DOCX 24 kb)
Presenting associations between GRS SNPs and baseline clinical characteristics. (DOCX 18 kb)