- Research article
- Open Access
Rheumatoid arthritis response to treatment across IgG1 allotype – anti-TNF incompatibility: a case-only study
Arthritis Research & Therapyvolume 17, Article number: 63 (2015)
We have hypothesized that incompatibility between the G1m genotype of the patient and the G1m1 and G1m17 allotypes carried by infliximab (INX) and adalimumab (ADM) could decrease the efficacy of these anti-tumor necrosis factor (anti-TNF) antibodies in the treatment of rheumatoid arthritis (RA).
The G1m genotypes were analyzed in three collections of patients with RA totaling 1037 subjects. The first, used for discovery, comprised 215 Spanish patients. The second and third were successively used for replication. They included 429 British and Greek patients and 393 Spanish and British patients, respectively. Two outcomes were considered: change in the Disease Activity Score in 28 joint (ΔDAS28) and the European League Against Rheumatism (EULAR) response criteria.
An association between less response to INX and incompatibility of the G1m1,17 allotype was found in the discovery collection at 6 months of treatment (P = 0.03). This association was confirmed in the replications (P = 0.02 and 0.08, respectively) leading to a global association (P = 0.001) that involved a mean difference in ΔDAS28 of 0.4 units between compatible and incompatible patients (2.3 ± 1.5 in compatible patients vs. 1.9 ± 1.5 in incompatible patients) and an increase in responders and decrease in non-responders according to the EULAR criteria (P = 0.03). A similar association was suggested for patients treated with ADM in the discovery collection, but it was not supported by replication.
Our results suggest that G1m1,17 allotypes are associated with response to INX and could aid improved therapeutic targeting in RA.
Advances in the treatment of rheumatoid arthritis (RA) including anti-tumor necrosis factor (anti-TNF) monoclonal antibodies have led to successful control of the disease in many patients . However, it is still necessary to change the initial drug because of poor efficacy in a significant fraction of them. This trial and error approach increases the burden of RA and could consume the early months after RA onset when a better long-term prognosis could be obtained . Confronted with this problem rheumatologists have sought predictive biomarkers to orient drug choice [3,4]. A few have already been found as autoantibody seronegativity that identifies patients with poor response to rituximab, but more are needed [3-6].
An area of recent progress has been awareness of the importance of the blood levels of the biologics at the trough between two treatment doses and of anti-drug antibodies [7-13]. These antibodies could work in two ways to decrease the drug’s efficacy: neutralizing the biologic and increasing its clearance. They are present in most patients showing infusion reactions to infliximab (INX) [12,14,15].
A related area of concern for biologics that bear the fragment crystallizable (Fc) of immunoglobulin G (IgG) is the possibility of inducing anti-allotype antibodies or T cell reactions in incompatible patients [16-20]. The allotypes are protein polymorphisms in the IgGs that are able to induce antibodies when injected in incompatible subjects (those not bearing them). There are several human immunoglobulin allotypes, but for the biologics used to treat RA the most relevant are in the heavy chain of IgG1 (Figure 1). INX and adalimumab (ADM), two of the most commonly used anti-TNF monoclonal antibodies, bear the G1m1 and G1m17 allotypes [16,17]. At a genetic level, these two allotypes are in perfect linkage disequilibrium (LD) in Europeans (r2 = 1) and, therefore, are referred as the G1m1,17 allotype but the anti-allotype antibodies are directed either against the G1m1 epitope or against the G1m17 epitope . These allotypes are missing in a large fraction of Europeans, who are susceptible to mounting an anti-allotype response when exposed to INX or ADM. Detection of the anti-allotype antibodies targeting biologics has been rare, but this could be due to the technical challenges presented by these antibodies [16,17,19,22]. In addition, exposition to the G1m1 allotype in incompatible subjects induces T-cell responses directed against a different peptide in the IgG1 molecule . The G1m1,17 allotypes could also influence treatment response by modifying the strength of antibody immune responses, as has been shown for several antigens, for rheumatoid factor (RF), and for anti-ADM antibodies [20,23-25]. These differences could be mediated by the allotypes themselves or by other variants in the immunoglobulin heavy chain locus (IGH) through LD with the allotypes [26,27].
The IGH locus displays a high level of structural complexity, with highly homologous genes (Figure 1) arising from multiple segmental duplication events [26,27]. Consequently, this locus is not adequately covered on genome-wide SNP chips [23,28]. Specifically, none of the SNPs coding for G1m allotypes or their tagging SNPs was included in the genome-wide association studies (GWAS) analyzing the response of patients with RA to treatment with anti-TNF drugs [29-32].
Here, we show that incompatibility at the G1m allotype was associated with less response to treatment with INX. This association was supported by three sets of patient samples, one used for discovery and the other two for replication. The difference between G1m compatible and incompatible patients was small, but of possible utility because it is similar to that observed between seronegative and seropositive patients in the response to rituximab, which is used in clinical practice, and because it was stronger in some patient subgroups. A similar association was found for the patients treated with ADM in the discovery set, but it was not reproduced in the first replication set, which included a considerably larger number of ADM treated patients.
The three sets of patients with RA included in the study are shown in Figure 2. The discovery set included 215 patients of Spanish ancestry with RA who were recruited in six Spanish hospitals. Evaluations included the disease activity score in 28 joints (DAS28), which was available at the start of treatment with INX or ADM and after three, six and twelve months. The first replication set consisted of 429 patients treated with INX or ADM, 384 from the Biologics in Rheumatoid Arthritis Genetics and Genomics Study Syndicate (BRAGGSS), all of them with UK white ancestry, and 45 patients from two Greek hospitals, all of them of Greek ancestry (Figure 2). The DAS28 of these patients was evaluated at the start and after six months of treatment with INX or ADM. The second replication set of 393 patients with RA, all treated with INX, was subsequently collected in twelve Spanish hospitals (nine new hospitals plus new patients from three of the hospitals contributing to the discovery set, n = 234 Spanish with white ancestry) and by BRAGGSS (n = 159 UK whites). Analyses were limited to the patients with valid genotypes (95.3, 94.2 and 98.2% in the three patient sets, respectively): 205 from the discovery set (183 with three months follow-up, 186 with six months follow-up and 150 with twelve months follow-up); 404 from the first replication set (all with six month follow-up); and 386 from the second replication set (97 with three months of follow-up, 362 with six months of follow-up and 169 with twelve months of follow-up). All were biologic-naïve at treatment start. The indication of treatment, choice of drug and control of evolution were performed independently of this study. All patients provided their written informed consent. Collection of samples was approved by the local ethics committees (listed in the Appendix) and the study was approved by the Comité Ético de Investigación Clínica de Galicia (registry numbers 2009/173, 2011/162 and 2013/156).
G1m allotype genotyping
The G1m allotypes were studied at the DNA level [17,18]. We have genotyped two nsSNPs encoding them: rs1071803 (which A allele codes for the G1m17 allotype, whereas the G allele codes for the G1m3 allotype) and rs11621259 (which C allele codes for one of the two amino-acids comprising the G1m1 allotype, whereas the A allele codes for a protein that does not induce antibodies, namely a null allotype). Due to paralog sequences in other IGH genes, there are no tests able to discriminate with confidence between G/A and A/A genotypes at rs1071803 and between C/A and C/C genotypes at rs11621259 (Figure 1). A carrier analysis was therefore performed that distinguished carriers and non-carriers of the G1m1,17 allotypes. Genotypes were obtained by PCR amplification followed by single-base extension with the SNaPshot Multiplex Kit (Applied Biosystems, Foster City, California). Samples with different genotypes were sequenced to assess the accuracy of results. The primers used are listed in Additional file 1: Table S1.
The Statistica 7.0 (Statsoft) software was used throughout. Appropriate analyses were performed to evaluate differences in clinical characteristics stratified according to the specific anti-TNF (Table 1). Thus, 2 × 2 tables with Chi-square tests were used for comparing dichotomous clinical characteristics, 2 × 3 tables for variables with three levels as the EULAR response criteria, while t-tests were used to evaluate differences between quantitative clinical characteristics as ΔDAS28, baseline C-reactive protein (CRP) or baseline erythrocyte sedimentation rate (ESR). Response to treatment was considered as ΔDAS28 or according to the EULAR criteria . A generalized linear model for ΔDAS28 and a logistic regression model for the EULAR criteria (confronting good responders with non-responders; moderate responders were not considered) were fitted. Unadjusted and adjusted analyses were conducted. Included covariates are indicated in each analysis. Samples with missing data for any of the variables implicated in each specific analysis were excluded. Meta-analysis of the beta coefficients corresponding to the regression models for the three sets of patients was done according to a fixed effects model using inverse-variance weights as implemented in the meta library of the R project .
Association in the discovery patient set
Patients with RA in the discovery collection showed a median of six years between diagnosis and start of the treatment with INX or ADM (Table 1). They had clinical characteristics of severe disease: 85.3% were seropositive, 84.9% had developed joint erosions and they had an active disease (baseline DAS28 = 5.9 ± 1.2) after treatment with a mean of 2.5 disease-modifying antirheumatic drugs (DMARDs). The assessed treatment was with the first biologic administered and it was most often INX (73.7%). The patients treated with INX showed a lower decrease in DAS28 at three and six months than the patients treated with ADM, indicating the need to consider the drug as an important variable in subsequent analyses. There were no other significant differences between the patients treated with the two anti-TNF drugs.
Results of the two nsSNPs for the G1m1,17 allotypes were fully concordant confirming the perfect LD between them. Slightly less than half the patients (49.3%) did not carry the G1m1,17 allotypes and, therefore, were incompatible with the INX and ADM allotypes. A significant association was found in unadjusted analysis between lower response (ΔDAS28) and incompatibility at the G1m1,17 allotypes after six months of treatment (Table 2). This association remained significant after adjusting for gender, baseline DAS28, rheumatoid factor (RF) and anti-TNF (frequency data provided in Additional file 1: Table S2). Stratified analysis by anti-TNF showed significant association in the two strata, but only with INX in the adjusted analysis (Table 2). Association in patients treated with ADM was dubious because it disappeared after adjusting for the covariates, but the effect size was equal to that shown in patients treated with INX (Beta = 0.16). Therefore, the lack of association could be attributed to the lower number of patients treated with ADM than with INX and reduced statistical power (Table 2). No association between ΔDAS28 and G1m1,17 carrier status was found either at three or twelve months of treatment, but nominal differences were in the same direction as that observed at six months (data not shown). No significant association was detected using the EULAR response criteria at any time of follow-up (Additional file 1: Table S3).
First replication set
New patient samples from the UK and from Greece were obtained to replicate the previous findings. These patients showed some clinical features that were different from the discovery set of patients in baseline DAS28, baseline Health Assessment Questionnaire (HAQ), erosions, time since disease diagnosis and smoking (Tables 1 and 3). However, these patients also showed an established disease (median of 10 years from diagnosis to the start of the first anti-TNF treatment), showing high activity (baseline DAS28 = 6.5 ± 1.1), frequent erosions (58.5%) and seropositivity (79.7%). There were some differences between patients treated with INX and those treated with ADM in this replication set: the INX treated patients were younger at diagnosis, showed more joint erosions and responded with less decrease in DAS28 than the patients treated with ADM. Therefore, it was also necessary to consider the anti-TNF in the analysis.
More than half of the patients in this set (56.4%) were carriers of the G1m1,17 allotypes. Non-carrier status was associated with a lower response as evaluated with ΔDAS28, but only in the patients treated with INX (n = 199; Table 4). The association was significant, both in the unadjusted and in the adjusted analyses, replicating the finding of the discovery set. No association was found in the patients treated with ADM (Table 4). Analysis of the secondary outcome, EULAR response criteria, also showed significant association between incompatible allotypes and worse response to INX (P = 0.03; Additional file 1: Table S4), but not to ADM.
Second replication set
Subsequently, a new set of patients with RA treated with INX was collected in Spanish hospitals and in the BRAGGSS. These patients also showed active disease with high baseline DAS28 (6.1 ± 1.2) and HAQ (1.9 (1.4 to 2.3)), a long evolution and common erosions (77.6%), but with a lower frequency of seropositivity than the two previous patient sets (Table 3). Response to treatment at six months was more similar to the discovery collection, with more balanced proportions of the three response classes, than to the first replication set.
Patients in this second replication set also showed an association between lower response to treatment with INX at six months and incompatibility at the G1m1,17 allotype (Table 4, lower rows), but the difference was only significant before adjusting for covariates (P = 0.03, versus P = 0.08 after adjustment by baseline DAS28 and gender). Comparison of responder with non-responder patients according to the EULAR criteria did not show significant differences (Additional file 1: Table S4).
Results from the patients with RA in the three previous sets were combined by simple pooling and by meta-analysis. Pooled analysis showed a significantly worse response at six months in the patients treated with INX who were carriers of the incompatible allotypes (Table 5). Very similar results were obtained by meta-analysis: significant association between higher ΔDAS28 and G1m1,17 compatibility (Beta = 0.12, 95% confidence interval (C.I.) = 0.05 to 0.19, P = 0.0005) with no heterogeneity between the three patient sets (I2 = 0%, 95% C.I. = 0.0 to 81.4%). In addition, analysis according to the EULAR criteria showed that the compatible patients were more commonly responders and less commonly non-responders to INX than the incompatible patients (30.1% versus 21.5% responders and 16.9% versus 22.4% non-responders; P = 0.03). On the contrary, combined analysis of response to ADM at six months (Table 5), or to any of the two anti-TNF at three or twelve months (not shown) did not show significant differences between G1m1,17 carriers and non-carriers.
Analysis of ΔDAS28 in function of the G1m1,17 status was repeated adjusting for all characteristics that were different between the discovery and replication sets (baseline DAS28, gender, baseline HAQ, age at diagnosis, years from diagnosis to treatment, erosive arthritis, anti-citrullinated peptide antibodies (ACPA) and smoking; Additional file 1: Table S5). Results were very similar to those already shown in spite of the exclusion of patients lacking data for any of the new covariates (Beta = 0.13, P = 0.01 for the 304 patients treated with INX who had complete information). We also checked the effect of stratification by the clinical characteristics taken individually. There was a strong effect of the ACPA status. A stronger association was observed in the ACPA negative patients receiving INX (n = 132, Beta = 0.26, P = 0.00092) than in the whole set of patients, whereas it was reduced and non-significant in the ACPA positive subgroup (n = 342, Beta = 0.06, P = 0.19). In addition, stratification by ACPA status of the patients receiving INX also showed association at 12 months of follow-up in the ACPA negative patients (n = 81, Beta = 0.22, P = 0.022 versus n = 164, Beta = −0.08, P = 0.3 in the ACPA positive patients). No effect was found at three months. Also, the association was stronger in the patients who initiated treatment at a younger age than the median of the whole set of patients (56 years) than in the group that started treatment older than 56 years (Beta = 0.20, P = 0.000040 in the 369 patients starting INX at ≤56 years versus Beta = −0.005, P = 0.9 in the 313 older patients) and in those in whom treatment was delayed less since diagnosis of RA (Beta = 0.20, P = 0.000054 in the 315 patients treated with INX ≤8 years after diagnosis, versus Beta = 0.02, P = 0.7 in the 359 with >8 years of delay). These two variables, age or treatment initiation and delay since diagnosis were weakly correlated (r = 0.22), and, consequently, the effect of one was not explained by the effect of the other. In addition, we checked the effect of considering the concomitant use of methotrexate, which was the DMARD used by 85.1% of the 883 patients with this information. The results changed very little, both when considering all the patients, or separately the patients treated with INX (n = 661, Beta = 0.11, P = 0.0016). The association showed a trend to be stronger in the patients taking methotrexate (Beta = 0.12), but the low number of patients in the group not taking this DMARD made the results inconclusive. Similarly, no differences in the effect of the Gm1,17 allotype were detected with the other clinical variables analyzed (smoking habit, RF status and age at diagnosis).
Our results showed a consistent association between the G1m allotypes and response to treatment with INX at six months of follow-up. A decreased ΔDAS28 was observed in the patients with incompatible allotypes in the discovery patients and in the two replication sets. In addition, fewer responders and more non-responders were present among the G1m1,17 incompatible patients than among the compatible patients. The effect of G1m incompatibility was weak in the whole set of patients, but stronger in some patient subgroups. Therefore, we think that it could be of clinical utility if it is confirmed in further studies.
We consider the association of response to INX as consistently reproduced in spite of the P = 0.08 obtained in the second replication, which is over the pre-specified significance threshold. This opinion is supported by two arguments. First, this study hypothesized a lower response in the incompatible patients, which is the result that has been found. Therefore, single-tailed tests could have been used in place of two-tailed tests . Single-tailed tests would lead to more significant differences in the three patient sets including P = 0.04 in the second replication. Second, some authors have recommended applying a lower significance threshold for replication than for discovery to compensate for the overestimated effect sizes commonly obtained in the first study showing significant association [36,37]. In any case, confirmation by independent researchers is needed for validation of this biomarker. Unfortunately, we cannot find support in the published GWAS of response to anti-TNF because they did not include the allotype-specific nsSNPs and showed low coverage of the IGH locus [23,28-32].
The requirement of replication is still more important for the subgroup analyses that we only report here as hypothesis generating exploratory results. This caution is necessary because subgroup analyses are especially prone to false positive results . Therefore, the stronger associations found in ACPA negative patients than in ACPA positive patients, in patients younger than the median age at treatment initiation and in those patients with a short delay in treatment since diagnosis should be taken with caution. These results, if confirmed, could reflect the differences in immune responsiveness between the subgroups of patients.
In effect, our results are concordant with the initial hypothesis involving induction of immune responses against the incompatible allotype of INX. These immune responses could include anti-allotype antibodies [16,17,19] or a T cell response independent of anti-allotype antibodies . The anti-allotype antibodies would target the INX molecule leading to its accelerated clearance. Anti-allotype antibodies have been found in polytransfused subjects and in multiparous women  who were exposed to blood containing allotype incompatible immunoglobulins. However, these antibodies have not been found in some studies that have searched for them after administration of IgG1-based drugs [17,20,39]. These negative results have been attributed to the technical difficulty of the assays [16,17,19,22,39]. In this respect, the authors who were successful in finding anti-allotype antibodies in patients treated with IgG1-based drugs had used very sensitive assays, which included few washes to avoid dissociation of low affinity antibodies and acid dissociation of pre-existing immunocomplexes . Our approach, at the DNA level instead of the anti-allotype antibodies, has avoided these technical difficulties, has permitted us to use patients lacking sera at the relevant point in the treatment course, and allows for other forms of involvement of the allotypes beyond the induction of antibodies. One of these alternative mechanisms is a different type of immune response against the G1m1 allotype that involves CD4 T cell activation and production of cytokines . The T cell antigen does not include the allotype but a nearby peptide that becomes accessible to antigen processing and presentation because the G1m1 allotype introduces an asparaginyl endopeptidase cleavage site . None of the two alternatives, antibodies or T cell responses, could be tested in our patients because of unavailability of the relevant samples.
Alternative mechanisms for the association that do not involve anti-drug immune reactions could also be invoked. They are suggested by differences in the strength of antibody responses against autoantigens, pathogens and cancer antigens associated with the G1m allotypes [20,23-25]. These associations could be due to the allotypes themselves, because they modify the interaction with Fc-gamma receptors , or to other variants in the immunoglobulin heavy chain locus that is very polymorphic, displays extensive LD and includes the sequences of all immunoglobulin heavy chains [26,27].
The associations were observed in the three sets of patients at six months with ΔDAS28, but only in the first replication and in the combined analysis with the EULAR response. This discordance between the outcomes is congruent with the lower sensitivity to change of the EULAR response, which was considered as a dichotomous variable, than of ΔDAS28, which is a continuous variable [41,42]. Another aspect of our study deserving comment is that results at three or twelve months of follow-up were non-significant (although the same direction of change was observed). These negative results could be ascribed to less powerful analysis due to fewer patients with data available at these times than at six months (Figure 2), or to latency in the development of an immune response preventing association at three months, or to change of treatment in patients showing poor response at six months. In an attempt to compensate by the latter confounding factor, we included the 33 patients that we knew had interrupted treatment after six months due to inefficacy in the EULAR response analysis at twelve months (as non-responders), but this analysis did not show association (not shown). In addition, it is possible that differences in the composition of the patient sets at the different times of evaluation contribute to the discrepant results. This alternative is suggested by the demonstration of association at twelve months of follow-up in the subset of ACPA negative patients, which was also the subset with the strongest association at six months of follow-up. At present, we cannot discriminate between the discussed possibilities.
Differences between the discovery and the replication sets included origin within Europe, baseline disease activity, time of evolution, seropositivity, erosions and smoking. Reproducibility of the association across these differences suggests that it is robust and that could be confirmed in additional studies. This replication is a requirement before defining the applicability of our findings. However, we can already comment on the effect size of the INX association we have observed. The mean difference in ΔDAS28 (which is larger than the attributable effect according to linear regression analysis probably because it includes effects attributable to other covariates) between G1m1,17 compatible and incompatible patients after six months in INX was 0.4. This is below 0.6, the change considered significant in the EULAR criteria  but similar to the difference in ΔDAS28 between seropositive and seronegative RA patients treated with rituximab at six months (0.3 for RF and 0.6 for ACPA according to a meta-analysis of patient registers ; and 0.35 for seropositive patients from a meta-analysis of clinical trials ), which is considered sufficient to help in the choice of treatment. Therefore, it is possible that the G1m1,17 allotypes could become useful, after further confirmation and analysis, to inform drug choice in RA treatment. This could be especially true in some subgroups of patients as suggested by the strong association in ACPA negative patients and in young patients and those initiating treatment early.
In contrast with INX, patients treated with ADM did not show consistent association between response and G1m1,17 compatibility. The total number of ADM treated patients was not enough to exclude with sufficient power an effect as was observed with INX (643 patients will be needed to exclude with 0.8 power an effect size of slope = 0.11). However, it is possible that G1m incompatibility does not affect response to ADM because there is already precedent of a study showing that response of RA patients to ADM is not affected by their allotype . In addition, several differences between INX and ADM make more likely the development of immune responses against the first. They include the larger immunogenicity of other parts of the INX molecule, which is a mouse/human chimera, than of ADM, which is a human antibody; the different routes of administration and much wider fluctuations of INX than of ADM, and the development of neutralizing anti-idiotypic antibodies against ADM that could reflect the immunodominance of the ADM paratope over other epitopes . The large fluctuations of INX concentration are due to its intravenous administration every eight weeks in contrast with the subcutaneous administration every two weeks of ADM .
Other biologics used for the treatment of RA bear the G1m1 or G1m17 allotypes and it will be interesting to study if the response to them is affected by the allotype of the patients. Rituximab and tocilizumab bear the G1m1,17 allotypes, whereas abatacept has the G1m1 allotype and certolizumab pegol the G1m17 allotype due to the incomplete presence of the IgG1 constant region in their molecules [17,45-47]. Etanercept, in turn, cannot induce any anti-allotype responses because it lacks the G1m3/G1m17 allotypes and bears the null allotype at G1m1 . We were unable to find information about the allotypes of golimumab (although, it is described as bearing the INX IgG1 heavy chain ).
The G1m1,17 allotypes of patients with RA have shown an association with response to INX treatment at six months of follow-up. This association was replicated in independent patient sets. The patients who were incompatible with the allotypes in the INX molecule showed a poorer response than the compatible ones. This pattern of response is consistent with the development of immune reactivity against INX. The difference in response between allotype compatible and incompatible patients was small, but of possible utility because it was of similar magnitude to the difference in response to rituximab between seronegative and seropositive patients, and because it was larger in some subgroups of patients.
anti-citrullinated peptide antibodies
Biologics in Rheumatoid Arthritis Genetics and Genomics Study Syndicate
disease activity score in 28 joints
disease-modifying antirheumatic drugs
erythrocyte sedimentation rate
European League against Rheumatism
genetic marker in constant region of immunoglobulin gamma heavy chain of subclass1
genome-wide association study
Health Assessment Questionnaire
immunoglobulin G subclass 1
immunoglobulin heavy chain
immunoglobulin heavy constant gamma
non synonymous single nucleotide polymorphism
polymerase chain reaction
single nucleotide polymorphism
tumor necrosis factor
Smolen JS, Aletaha D, Bijlsma JW, 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.
Finckh A, Liang MH, van Herckenrode CM, de Pablo P. Long-term impact of early treatment on radiographic progression in rheumatoid arthritis: a meta-analysis. Arthritis Rheum. 2006;55:864–72.
Emery P, Dorner T. Optimising treatment in rheumatoid arthritis: a review of potential biological markers of response. Ann Rheum Dis. 2011;70:2063–70.
Isaacs JD, Ferraccioli G. The need for personalised medicine for rheumatoid arthritis. Ann Rheum Dis. 2011;70:4–7.
Chatzidionysiou K, Lie E, Nasonov E, Lukina G, Hetland ML, Tarp U, et al. Highest clinical effectiveness of rituximab in autoantibody-positive patients with rheumatoid arthritis and in those for whom no more than one previous TNF antagonist has failed: pooled data from 10 European registries. Ann Rheum Dis. 2011;70:1575–80.
Isaacs JD, Cohen SB, Emery P, Tak PP, Wang J, Lei G, et al. Effect of baseline rheumatoid factor and anticitrullinated peptide antibody serotype on rituximab clinical response: a meta-analysis. Ann Rheum Dis. 2013;72:329–36.
Wolbink GJ, Vis M, Lems W, Voskuyl AE, de Groot E, Nurmohamed MT, et al. Development of antiinfliximab antibodies and relationship to clinical response in patients with rheumatoid arthritis. Arthritis Rheum. 2006;54:711–5.
Bendtzen K, Geborek P, Svenson M, Larsson L, Kapetanovic MC, Saxne T. Individualized monitoring of drug bioavailability and immunogenicity in rheumatoid arthritis patients treated with the tumor necrosis factor alpha inhibitor infliximab. Arthritis Rheum. 2006;54:3782–9.
Bartelds GM, Wijbrandts CA, Nurmohamed MT, Stapel S, Lems WF, Aarden L, et al. Clinical response to adalimumab: relationship to anti-adalimumab antibodies and serum adalimumab concentrations in rheumatoid arthritis. Ann Rheum Dis. 2007;66:921–6.
Radstake TR, Svenson M, Eijsbouts AM, van den Hoogen FH, Enevold C, van Riel PL, et al. Formation of antibodies against infliximab and adalimumab strongly correlates with functional drug levels and clinical responses in rheumatoid arthritis. Ann Rheum Dis. 2009;68:1739–45.
Bartelds GM, Krieckaert CL, Nurmohamed MT, van Schouwenburg PA, Lems WF, Twisk JW, et al. Development of antidrug antibodies against adalimumab and association with disease activity and treatment failure during long-term follow-up. JAMA. 2011;305:1460–8.
Pascual-Salcedo D, Plasencia C, Ramiro S, Nuno L, Bonilla G, Nagore D, et al. Influence of immunogenicity on the efficacy of long-term treatment with infliximab in rheumatoid arthritis. Rheumatology. 2011;50:1445–52.
Maneiro JR, Salgado E, Gomez-Reino JJ. Immunogenicity of monoclonal antibodies against tumor necrosis factor used in chronic immune mediated inflammatory conditions: systematic review and metaanalysis. JAMA Intern Med. 2013;173:1416–28.
Krintel SB, Grunert VP, Hetland ML, Johansen JS, Rothfuss M, Palermo G, et al. The frequency of anti-infliximab antibodies in patients with rheumatoid arthritis treated in routine care and the associations with adverse drug reactions and treatment failure. Rheumatology. 2013;52:1245–53.
Plasencia C, Pascual-Salcedo D, Nuno L, Bonilla G, Villalba A, Peiteado D, et al. Influence of immunogenicity on the efficacy of long term treatment of spondyloarthritis with infliximab. Ann Rheum Dis. 2012;71:1955–60.
Jefferis R, Lefranc MP. Human immunoglobulin allotypes: possible implications for immunogenicity. MAbs. 2009;1:332–8.
Magdelaine-Beuzelin C, Vermeire S, Goodall M, Baert F, Noman M, Assche GV, et al. IgG1 heavy chain-coding gene polymorphism (G1m allotypes) and development of antibodies-to-infliximab. Pharmacogenet Genomics. 2009;19:383–7.
Stickler MM, Reddy A, Xiong JM, Hinton PR, DuBridge R, Harding FA. The human G1m1 allotype associates with CD4+ T-cell responsiveness to a highly conserved IgG1 constant region peptide and confers an asparaginyl endopeptidase cleavage site. Genes Immun. 2011;12:213–21.
Tatarewicz SM, Juan G, Swanson SJ, Moxness MS. Epitope characterization of pre-existing and developing antibodies to an aglycosylated monoclonal antibody therapeutic of G1m17,1 allotype. J Immunol Methods. 2012;382:93–100.
Bartelds GM, de Groot E, Nurmohamed MT, Hart MH, van Eede PH, Wijbrandts CA, et al. Surprising negative association between IgG1 allotype disparity and anti-adalimumab formation: a cohort study. Arthritis Res Ther. 2010;12:R221.
Lefranc MP, Lefranc G. The immunoglobulin facts book. London: Academic; 2001.
Isaacs JD, Watts RA, Hazleman BL, Hale G, Keogan MT, Cobbold SP, et al. Humanised monoclonal antibody therapy for rheumatoid arthritis. Lancet. 1992;340:748–52.
Pandey JP, Li Z. The forgotten tale of immunoglobulin allotypes in cancer risk and treatment. Exp Hematol Oncol. 2013;2:6.
Pandey JP. Immunoglobulin GM, and KM allotypes and vaccine immunity. Vaccine. 2000;19:613–7.
Pandey JP, Luo Y, Elston RC, Wu Y, Philp FH, Astemborski J, et al. Immunoglobulin allotypes influence IgG antibody responses to hepatitis C virus envelope proteins E1 and E2. Hum Immunol. 2008;69:158–64.
Kidd MJ, Chen Z, Wang Y, Jackson KJ, Zhang L, Boyd SD, et al. The inference of phased haplotypes for the immunoglobulin H chain V region gene loci by analysis of VDJ gene rearrangements. J Immunol. 2012;188:1333–40.
Watson CT, Steinberg KM, Huddleston J, Warren RL, Malig M, Schein J, et al. Complete haplotype sequence of the human immunoglobulin heavy-chain variable, diversity, and joining genes and characterization of allelic and copy-number variation. Am J Hum Genet. 2013;92:530–46.
Pandey JP. Genomewide association studies and assessment of risk of disease. N Engl J Med. 2010;363:2076–7. author reply 2077.
Cui J, Stahl EA, Saevarsdottir S, Miceli C, Diogo D, Trynka G, et al. Genome-wide association study and gene expression analysis identifies CD84 as a predictor of response to etanercept therapy in rheumatoid arthritis. PLoS Genet. 2013;9:e1003394.
Liu C, Batliwalla F, Li W, Lee A, Roubenoff R, Beckman E, et al. Genome-wide association scan identifies candidate polymorphisms associated with differential response to anti-TNF treatment in rheumatoid arthritis. Mol Med. 2008;14:575–81.
Plant D, Bowes J, Potter C, Hyrich KL, Morgan AW, Wilson AG, et al. Genome-wide association study of genetic predictors of anti-tumor necrosis factor treatment efficacy in rheumatoid arthritis identifies associations with polymorphisms at seven loci. Arthritis Rheum. 2011;63:645–53.
Umicevic Mirkov M, Cui J, Vermeulen SH, Stahl EA, Toonen EJ, Makkinje RR, et al. Genome-wide association analysis of anti-TNF drug response in patients with rheumatoid arthritis. Ann Rheum Dis. 2013;72:1375–81.
van Gestel AM, Haagsma CJ, van Riel PL. Validation of rheumatoid arthritis improvement criteria that include simplified joint counts. Arthritis Rheum. 1998;41:1845–50.
R: A Language and Environment for Statistical Computing. http://www.R-project.org.
Heath D. An introduction to experimental design and statistics for biology. Boca Raton, Florida: CRC Press; 1995.
Lander E, Kruglyak L. Genetic dissection of complex traits: guidelines for interpreting and reporting linkage results. Nat Genet. 1995;11:241–7.
Vieland VJ. The replication requirement. Nat Genet. 2001;29:244–5.
Guillemin F. Primer: the fallacy of subgroup analysis. Nat Clin Pract Rheumatol. 2007;3:407–13.
Rispens T, de Vrieze H, de Groot E, Wouters D, Stapel S, Wolbink GJ, et al. Antibodies to constant domains of therapeutic monoclonal antibodies: anti-hinge antibodies in immunogenicity testing. J Immunol Methods. 2012;375:93–9.
Namboodiri AM, Pandey JP. Differential inhibition of trastuzumab- and cetuximab-induced cytotoxicity of cancer cells by immunoglobulin G1 expressing different GM allotypes. Clin Exp Immunol. 2011;166:361–5.
Fransen J, van Riel PL. Outcome measures in inflammatory rheumatic diseases. Arthritis Res Ther. 2009;11:244.
American College of Rheumatology Committee to Reevaluate Improvement Criteria. A proposed revision to the ACR20: the hybrid measure of American College of Rheumatology response. Arthritis Rheum. 2007;57:193–202.
van Schouwenburg PA, Krieckaert CL, Rispens T, Aarden L, Wolbink GJ, Wouters D. Long-term measurement of anti-adalimumab using pH-shift-anti-idiotype antigen binding test shows predictive value and transient antibody formation. Ann Rheum Dis. 2013;72:1680–6.
Tracey D, Klareskog L, Sasso EH, Salfeld JG, Tak PP. Tumor necrosis factor antagonist mechanisms of action: a comprehensive review. Pharmacol Ther. 2008;117:244–79.
Drug In Focus. Tocilizumab. http://www.genericsweb.com/index.php?object_id=762.
Suzuki T, Ishii-Watabe A, Tada M, Kobayashi T, Kanayasu-Toyoda T, Kawanishi T, et al. Importance of neonatal FcR in regulating the serum half-life of therapeutic proteins containing the Fc domain of human IgG1: a comparative study of the affinity of monoclonal antibodies and Fc-fusion proteins to human neonatal FcR. J Immunol. 2010;184:1968–76.
Knox C, Law V, Jewison T, Liu P, Ly S, Frolkis A, et al. DrugBank 3.0: a comprehensive resource for ‘omics’ research on drugs. Nucleic Acids Res. 2011;39:D1035–41.
Shealy DJ, Cai A, Staquet K, Baker A, Lacy ER, Johns L, et al. Characterization of golimumab, a human monoclonal antibody specific for human tumor necrosis factor alpha. MAbs. 2010;2:428–39.
We thank Carmen Pena-Pena for technical assistance, Manuel Calaza for help with data management and members of Biologics in Rheumatoid Arthritis Genetics and Genomics Study Syndicate (BRAGGSS) for their patient data and sample contribution. Funding was provided by grants PI11/01048, PI12/01909 and by RETICS Program, RD08/0075/0019 and RD12/0009/0008 (RIER) of the Instituto de Salud Carlos III (Spain) that are partially financed by the European Regional Development Fund of the European Union.
The authors declare that they have no competing interests.
AM and AG designed the study; AM genotyped the samples, analyzed the data and wrote the first manuscript; EP-P, FN-S, VM, ARS, BM, YV, TS, AF-N, MCO, JN, JDC, AMa, DP-S, BJ, PC, MJM-R, RC, MAF, RG-P, FJB, CM, ER, LV, JJA-S, AB, JM, DP, JI, AWM, AB, AGW, JJG-R provided samples and clinical and demographic information of the subjects, contributed to the interpretation of the data and revised and approved the final version of the manuscript; JJG-R and AG obtained funding; AG supervised the study and has complete access to all data, revised the manuscript and approved its final version. All authors read and approved the final manuscript.
Primers and probes in 5′ → 3′ direction employed in this study. All primers were selected with Primer3 and Oligos software. They were checked to avoid formation of dimers between the oligonucleotides included in the reaction. Table S2. Clinical characteristics of patients from the discovery set of patients stratified by treatment and G1m1,17 allotype. Table S3. Association of G1m carrier status with treatment response according to the EULAR criteria at six months of treatment with anti-TNF in the discovery samples. Table S4. Association of the G1m carrier status with treatment response according to the EULAR criteria at six months of treatment with anti-TNF in the replication collections. Table S5. Comparison between the clinical characteristics of the three RA patient sets (discovery and two replications) with six months follow-up for infliximab treatment.
The following ethics committees approved sample collection and this study: Comité Ético de Investigación Clínica de Galicia (currently, Comité Autnómico de Ética de la Investigacion de Galicia), Comité Ético de Investigación Clínica del Hospital Universitario Virgen Macarena, Comité de Ética de la Investigación Clínica Fundación de Gestión Sanitaria del Hospital de la Santa Creu i Sant Pau, University of Thessally and Medical School Ethics Committee, Comité Ético de Investigación del Hospital Carlos Haya de Málaga, Comité Ético de Investigación Clínica de la Ciutat Sanitária i Universitária de Bellvitge, Comité Ético de Investigación Clínica del Hospital Clínic i Provincial de Barcelona, Comité de Bioética del Consejo Superior de Investigaciones Científicas, Cómite de Ética de Investigación Clínica del Hospital Universitario de la Paz, Comité Ético de Investigación Clínica del Hospital 12 de Octubre, Comité Ético de Investigación Clínica del Hospital Clínico Universitario Virgen de la Arrixaca, Comité Ético de Investigación Clínica del Hospital Universitario Virgen de las Nieves, Comité Ético de Investigación Cclínica del Hospital Virgen Victoria de Málaga, Cómite Ético de Investigación Clínica del Hospital San Cecilio, Cómite Ético de Investigacion Clínica del Instituto de Investigación Sanitaria Gregorio Marañón, Comité Ético de Investigación Clínica del Hospital Universitario Doctor Pesset, UK Central Office of Research Ethics Committee.