You are here

Article Text

Article menu
Electronic pages
Clinical and epidemiological research
Extended report
A genome-wide association study of rheumatoid arthritis without antibodies against citrullinated peptides
  1. L Bossini-Castillo1,
  2. C de Kovel2,
  3. H Kallberg3,4,
  4. R van ‘t Slot2,
  5. A Italiaander2,
  6. M Coenen5,
  7. P P Tak6,
  8. M D Posthumus7,
  9. C Wijmenga8,
  10. T Huizinga9,
  11. A H M van der Helm-van Mil9,
  12. G Stoeken-Rijsbergen9,
  13. Luis Rodriguez-Rodriguez10,
  14. Alejandro Balsa11,
  15. Isidoro González-Álvaro12,
  16. Miguel Ángel González-Gay13,
  17. Carmen Gómez-Vaquero14,
  18. B Franke5,
  19. LifeLines Cohort Study,
  20. S Vermeulen5,
  21. IE van der Horst-Bruinsma15,
  22. B A C Dijkmans15,
  23. G J Wolbink16,
  24. R A Ophoff2,
  25. M T Maehlen17,
  26. P van Riel18,
  27. M Merriman19,
  28. L Klareskog3,
  29. B A Lie17,
  30. T Merriman19,
  31. J B A Crusius20,
  32. E Brouwer7,
  33. J Martin1,
  34. N de Vries6,
  35. R Toes9,
  36. L Padyukov3,
  37. B P C Koeleman2
  1. 1Instituto de Parasitología y Biomedicina López-Neyra, Consejo Superior de Investigaciones Científicas (IPBLN-CSIC), Granada, Spain
  2. 2Department of Medical Genetics, UMCU, Utrecht, The Netherlands
  3. 3Rheumatology Unit, Department of Medicine, Karolinska University Hospital, Karolinska Institutet, Stockholm, Sweden
  4. 4Institute of Environmental Medicine, Karolinska Institutet, Sweden
  5. 5Department of Human Genetics, Radboud University Nijmegen Medical Centre, Nijmegen, The Netherlands
  6. 6Division of Clinical Immunology and Rheumatology, AMC, University of Amsterdam, Amsterdam, The Netherlands
  7. 7Department of Rheumatology, UMCG, Groningen, The Netherlands
  8. 8Department of Medical Genetics, UMCG, Groningen, The Netherlands
  9. 9Department of Rheumatology, LUMC, Leiden, The Netherlands
  10. 10Rheumatology Service, Hospital Clínico San Carlos, Madrid, Spain
  11. 11Rheumatology Service, Hospital Universitario La Paz, Madrid, Spain
  12. 12Rheumatology Service, Hospital Universitario La Princesa, Instituto de Investigación Sanitaria La Princesa, Madrid, Spain
  13. 13Rheumatology Service, Hospital Universitario Marqués de Valdecilla, IFIMAV, Santander, Spain
  14. 14Rheumatology Service, Hospital Universitari Bellvitge, Barcelona, Spain
  15. 15Department of Rheumatology, VUMC, Amsterdam, The Netherlands
  16. 16Jan van Breemen Research Institute, Amsterdam, The Netherlands
  17. 17Department of Medical Genetics; University of Oslo and Oslo University, hospital, Oslo, Norway; K. G. Jebsen Inflammation Research Centre, University of Oslo, Oslo, Norway
  18. 18Department of Rheumatology, Radboud University Nijmegen Medical Centre, Nijmegen, The Netherlands
  19. 19Department of Biochemistry, University of Otago, New Zealand
  20. 20Laboratory of Immunogenetics, Department of Medical Microbiology and Infection Control, VU University Medical Center, Amsterdam, The Netherlands

Abstract

Introduction Rheumatoid arthritis (RA) patients can be classified based on presence or absence of anticitrullinated peptide antibodies (ACPA) in their serum. This heterogeneity among patients may reflect important biological differences underlying the disease process. To date, the majority of genetic studies have focused on the ACPA-positive group. Therefore, our goal was to analyse the genetic risk factors that contribute to ACPA-negative RA.

Methods We performed a large-scale genome-wide association study (GWAS) in three Caucasian European cohorts comprising 1148 ACPA-negative RA patients and 6008 controls. All patients were screened using the Illumina Human Cyto-12 chip, and controls were genotyped using different genome-wide platforms. Population-independent analyses were carried out by means of logistic regression. Meta-analysis with previously published data was performed as follow-up for selected signals (reaching a total of 1922 ACPA-negative RA patients and 7087 controls). Imputation of classical HLA alleles, amino acid residues and single nucleotide polymorphisms was undertaken.

Results The combined analysis of the studied cohorts resulted in identification of a peak of association in the HLA-region and several suggestive non-HLA associations. Meta-analysis with previous reports confirmed the association of the HLA region with this subset and an observed association in the CLYBL locus remained suggestive. The imputation and deep interrogation of the HLA region led to identification of a two amino acid model (HLA-B at position 9 and HLA-DRB1 at position 11) that accounted for the observed genome-wide associations in this region.

Conclusions Our study shed light on the influence of the HLA region in ACPA-negative RA and identified a suggestive risk locus for this condition.

  • Ant-CCP
  • Gene Polymorphism
  • Rheumatoid Arthritis

http://dx.doi.org/10.1136/annrheumdis-2013-204591

Statistics from Altmetric.com

    Request Permissions

    If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

    Introduction

    Rheumatoid arthritis (RA) is a common autoimmune disease that is associated with a progressive loss of the joints induced by a chronic inflammation of the joint synovium.1 In this inflammatory environment, different products of cell apoptosis and necrosis accumulate, and citrullinated proteins can be detected.2 The production of anticitrullinated peptide/protein antibodies (ACPA) is a common but not essential characteristic of RA patients, which is thought to be influenced by the genetic background. Indeed, a strong correlation exists between ACPA and alleles of the HLA-DRB1 gene at chromosome 6 known as the shared epitope (SE).3 ,4 Hence, the presence or absence of ACPA divides RA patients into two serological and clinical subgroups, ACPA positive (ACPA+) and ACPA negative (ACPA−) RA.

    ACPA are highly specific for RA and appear years before the first clinical manifestation of RA.5 Moreover, a number of studies have shown that ACPA+ patients suffer more aggressive radiological joint damage.6–9 It is worth mentioning that, the most recent RA diagnostic criteria proposed by the American College of Rheumatology (ACR) includes ACPA among the classification factors.10

    Remarkably, different studies suggested the existence of a partial genetic overlap between the ACPA− and the ACPA+ phenotypes, but evidence supporting important genetic differences is increasing.11 In line with this, a study by Viatte et al12 reported that the SE, a HLA-DRB1*0401 tag-single nucleotide polymorphism (SNP) and variants located in TNFAIP3, GIN1/C5orf30, STAT4, ANKRD55/ IL6ST, BLK and PTPN22 showed significant associations with ACPA− RA patients. However, several RA susceptibility factors showed no association with the ACPA− subset, and those shared independently of the serotype revealed different strength of association and effect size. Two high-throughput genotyping efforts, a genome-wide association study (GWAS) and a SNP-dedicated Immunochip-based dense mapping, including ACPA− RA patients have been recently carried out. The GWAS focused on the genetic comparison of both serological RA subgroups.13 Although no genome-wide level association with the ACPA− subset was described, significant differences between ACPA+ and ACPA− patients in the HLA region were apparent13 The Immunochip approach revealed genome-wide level associations in the HLA region and the ANKRD55 locus with ACPA− RA, and supported differences with the ACPA+ subgroup in several loci.14 Additionally, a few loci have been associated with ACPA− RA but not with ACPA+ RA, such as CLEC16A and IRF5 that have been found following a candidate gene strategy.15 ,16 In spite of the increasing interest in the identification of genetic risk factors associated with ACPA− RA, the underlying genetic background in this subset of RA patients remains widely unknown.17

    The aim of this study was to identify novel susceptibility loci implicated in ACPA− RA susceptibility through a hypothesis-free GWA strategy. Therefore, we carried out a genome-wide combined analysis of three independent Caucasian European cohorts, and a follow-up phase using a previously studied GWAS cohort, comprising a total of 1922 ACPA− RA patients and 7087 unaffected controls.

    Patients and methods

    Populations

    All RA patients participating in this study met the 1987 ACR criteria for the classification of RA and were classified as ACPA− using standard methods as described elsewhere.13 ,18 Control populations comprised unaffected unrelated individuals from the same geographical and ethnical origin as cases. All individuals in this study were of European ancestry (self-reported and/or principal component analysis (PCA) derived, figure 1). This study was approved by the local ethics committees of the participating hospitals and all participating individuals gave written informed consent. A description of the analysed populations is provided in online supplementary note 1.

    Figure 1
    Figure 1

    Principal component analysis (PCA) plot of the analysed European populations. (A) The Netherlands I versus The Netherlands II populations PC1 versus PC3 plot; (B) The Netherlands I population PC1 versus PC2; (C) The Netherlands II population PC1 versus PC2; (D) Spanish population PC1 versus PC2 population.

    The power to detect an association with an OR of 1.40 considering 1922 cases and 7087 controls, at the p value <5×10−8 level, under an additive model and using a minor allele frequency (MAF) of 0.20 was of 96% (additional power calculations are provided in online supplementary table S1). For power estimation we used Power Calculator for Genome Wide Studies.19

    Genotyping and imputation

    The genotyping platforms for The Netherlands I, The Netherlands II and Spain cohorts are described in table 1. Quality control (QC) and imputation was conducted separately for each cohort and each chip type, using a common approach for all datasets (details in online supplementary note 2) using PLINK.20 The Swedish cohort genotyping and QC were performed as described in Padyukov et al.13

    View this table:
    Table 1

    Study design

    After QC, a total of 452 367 genotyped or imputed SNPs were shared between The Netherlands I, The Netherlands II and the Spanish case-control series; and 363 330 genotyped or imputed SNPs were shared between the previously mentioned cohorts and the Swedish population in Padyukov et al.13

    All SNP genotypes in the HLA region (chromosome 6 between the positions 29 000 000 and 34 000 000 reaching a total of 3882 SNPs) underwent a specific HLA imputation process as described in previous reports.21–23 This novel approach resulted in the imputation of classical HLA genotypes (HLA-A, HLA-B, HLA-C, HLA-DPA1, HLA-DPB1, HLA-DQA1, HLA-DQB1 and HLA-DRB1), their corresponding amino acid sequences and SNPs (see online supplementary note 3 and table S2).21 ,23

    Statistical methods

    The association of the imputed set of SNPs with ACPA− RA was tested by logistic regression, and this approach was followed separately in each cohort using PLINK (Genomic inflation factors are shown in online supplementary note 4 and quantile-quantile (QQ) plots in supplementary figure S1). We carried out individual population PCA using EIGENSTRAT software in order to detect population substructure.24 ,25 Then, the association analyses for The Netherlands I, The Netherlands II and Spain cohorts were adjusted for population substructure by including the first 10 PCs of each population as covariates in the logistic regression (λThe Netherlands I=1.07; λThe Netherlands II=1.02; λSpain=1.07; see online supplementary figure S1).

    A combined analysis using the inverse variance method under a fixed effects model was performed on the basis of the PCA-adjusted association results of The Netherlands I, The Netherlands II and Spain cohorts using PLINK. The inflation value for this analysis was λ=1.02. Heterogeneity across the datasets was evaluated using Cochran's Q test, and those loci showing a high heterogeneity (Q<0.05) were not considered for the validation step.

    Non-HLA SNPs and showing a significant p value at a tier 2 association level (p<5×10−5) in the combined analysis of The Netherlands I, The Netherlands II and Spain cohorts were selected for a validation step (table 2).

    View this table:
    Table 2

    Non-HLA SNPs associated with p value <5E-05 in the combined analysis of The Netherlands I, The Netherlands II and Spain populations

    In the validation phase, a combined analysis of the previously selected SNPs using the inverse variance method under a fixed effects model was performed using the association results of the three analysed populations and the ACPA− RA patients and controls of the Swedish cohort in Padyukov et al.13 A p value <5×10−6 was established as arbitrary threshold to consider a SNP as a suggestive association in the meta-analysis. Additionally, loci showing a high heterogeneity (Q<0.05) were discarded.

    For the analysis of the imputed data in the HLA region, we performed the association analyses by means of unconditional and conditional logistic regression analysis to account for dependency between the observed signals (details in online supplementary note 3). Finally, we searched recursively for models which better explained all association present in the HLA region (see online supplementary note 3).

    Regional plots were generated using LocusZoom, and the remaining figures were generated using R V.2.15.1.26 ,27

    Results

    Independent and combined analyses of three European ACPA− populations

    Independent analysis of The Netherlands I, The Netherlands II and Spain cohorts showed no associations at genome-wide significance level (p<5×10−8) (see online supplementary figures S2–S4). However, the combined analysis of the three cohorts identified a SNP in the CLYBL locus and two in SMIM21 showing significant risk associations at this level (figure 2). Moreover, 47 variants in 34 non-HLA loci showed a suggestive p value <5×10−5 and no heterogeneity, and were selected for follow-up (table 2 and figure 3). Additionally, the combined analysis of the three cohorts resulted in a peak of association in the HLA region in chromosome 6.

    Figure 2
    Figure 2

    Regional association plots showing (A) the CLYBL locus (B) the SMIM21 locus.

    Figure 3
    Figure 3

    Manhattan plot of the genome-wide association study (GWAS) of the three European cohorts comprising 1148 anti-citrullinated peptide antibodies-Rheumatoid arthritis (ACPA- RA) cases and 6008 healthy controls (The Netherlands I, The Netherlands II and Spanish populations). The −log10 of the inverse variance test p value of 452 367 SNPs is plotted against its physical chromosomal position. Chromosomes are shown in alternate colours. SNPs above the red line represent those with a p value <5×10−8 and variants in GWAS level loci are shown in orange.

    HLA region deep interrogation

    We further explored the association in the HLA region applying a new imputation method, which inferred the classical HLA alleles, polymorphic amino acid positions and SNPs in the studied populations.21 ,22 After imputation, the most significant association corresponded with amino acid position 67 of the HLA-DRβ1 molecule (p = 4.13×10−9, see online supplementary table S3). Three amino acids in this HLA-DRβ1 position were observed in our population (leucine, Leu; isoleucine, Ile; phenylalanine, Phe), but only the Leu67 variant showed genome-wide association (p=9.41×10−10. Furthermore, three additional amino acid residues showed a genome-wide level significant association in our study: two located in the HLA-DRB1 gene (the presence of threonine in the 181 amino acid position, Thr181, p value = 2.74×10−8, and the combination of serine, Ser, valine, Val and leucine, Leu, in the position 11, SerValLeu11, p value = 4.27×10−8) and one in the HLA-B locus (the presence of aspartic acid in the ninth position, Asp9, p value = 9.27×10−9). Moreover, these residues corresponded with the most associated amino acid positions in the omnibus test (see online supplementary note 3 and table S3). Then, we hypothesised that the association in the HLA region might be explained by polymorphic amino acid residues in the HLA molecules as observed in previous reports.14 ,21 ,22

    Therefore, we performed step-wise conditional regression analysis including each of the four amino acid positions that had an amino acid residue reaching genome-wide level association. These analyses identified two independent signals among the selected amino acid positions: HLA-B at position 9 and HLA-DRB1 at position 11(see online supplementary table S4). After controlling for the two previously mentioned associations, no signal showing a p<5×10−5 remained in the HLA region (figure 4).

    Figure 4
    Figure 4

    Manhattan plots showing the results of the unconditioned and conditioned analyses of the imputed HLA region. (A) Unconditioned analysis of the HLA region. (B) Conditional analysis controlling for the independent association of the 9th amino acid position at HLA-B. (C) Conditional analysis controlling for the previously mentioned association in the 9th amino acid position at HLA-B and the independent association of the 11th amino acid position at HLA-DRB1.

    Additionally, we observed that the HLA-B*0801 allele was the only classical HLA allele associated at genome-wide level in our study. However, this allele was indistinguishable from the HLA-B Asp9 variant using conditional logistic regression, due to their high linkage disequilibrium (r2=0.996). Regarding the HLA-DRB1 independent position, several classical HLA alleles share the genome-wide associated amino acid residues (see online supplementary note 5). Additionally, we confirmed that the hypothesis-free model was equivalent to the amino acid model, and addition of each amino acid variant to our model outperformed the model without it and achieved a better goodness-of-fit (table 3 and see online supplementary note 3). We also confirmed that the model including the two independent amino acid variants was the most parsimonious explanation for our data, and the addition of the remaining amino acids did not lead to a better goodness-of-fit (table 3).

    View this table:
    Table 3

    Comparison of goodness-of-fit of the different amino acid position models for the HLA region

    Follow-up meta-analysis

    After imputation, 363 330variants were shared among all three previously described European populations and the Swedish cohort described in Padyukov et al.13 The final set of 452 367 SNPs analysed in our cohorts tagged (with r2>0.8) a 52.37% of the SNPs with MAF>0.05 included in the HapMap phase 3 Caucasian of EUropean Ancestry (CEU) population and 51.63% of the SNPs in the CEU/Tuscans in Italy (TSI) populations (following the same parameters). As it can be observed in figure 5 and table 4, the variant in CLYB that previously showed genome-wide level association (rs9557321) did not reach the genome-wide significance threshold after meta-analysis but maintained a suggestive association (p value=5.82×10−8 OR=1.73). In the case of the genome-wide level associations in the SMIM21 locus (rs1943199, rs11663465), significant heterogeneity was found in the meta-analyses (Q<0.05) and the SNPs were discarded from analysis. Additionally, 30 out of 34 loci selected for the follow-up phase were shared between the genotyped and the Swedish populations, and two of them, CLYBL and SMIM21, were still significant considering a p value <5×10−6 in the meta-analysis (figure 3, regional plots are given in figure 2). The variants showing tier 2 associations were rs518167 (located in the 2nd intron of the GRM5 gene) and rs3790022 (located in the RNASEH2B-FAM124A intergenic region, ie, GUCY1B2 pseudogene) (table 4 and see online supplementary table S5)⇓.

    View this table:
    Table 4

    Previously selected non-HLA SNPs showing a p value <5×10−6 and HLA SNPs showing a p value <5×10−7 in the inverse variance meta-analysis of all the cohorts included the study

    Figure 5
    Figure 5

    Manhattan plot showing the meta-analysis comprising 1922 cases and 7087 controls from all the studied populations (The Netherlands I, The Netherlands II, Spain and Sweden populations). The −log10 of the inverse variance test p value of 363 330 SNPs is plotted against its physical chromosomal position. Chromosomes are shown in alternate colours. SNPs above the red line represent those with a p value <5×10−8 and variants in genome-wide association study level loci are shown in orange.

    Figure 6
    Figure 6

    Twin Manhattan plot showing the results of the combined analysis of the three European cohorts (left, dark blue), the Swedish cohort (right, green) and the meta-analysis of all the available cohorts (both sides, grey).

    Regarding the HLA region, no classical HLA allele information was available for the Swedish cohort, and unfortunately, our two-variant models could not be tested in this cohort. Nevertheless, the genome-wide associated variant in the HLA class I region, rs2596565, was included in the combined analysis, showed nominal association in the Swedish cohort (p=2.12×10−2 OR=1.24) and was the only SNP showing genome-wide significance in the HLA region (p=9.26×10−9 OR=1.4) (figure 6).

    Five out of 7 loci previously associated with ACPA− RA by Viatte et al12 were properly covered in our study (directly included or tagged by a variant showing r2>0.8). For these loci, we observed trends of association (p<0.05) in the TNFAIP3, GIN1/C5orf30 and PTPN22 loci (see online supplementary table S6). In the case of the HLA-DRB1* 04 : 01 tag-SNP, a trend association was observed, but the meta-analysis showed significant heterogeneity. Thus, we applied a random effects model and the initial association was lost (see online supplementary table S6).

    The previously reported variant of CLEC16A (rs6498169) was filtered out from the current analyses during QC, but the most strongly associated variant located in the gene showed a nominal association in the combined analysis (rs17803698 p=2.12×10−2 OR=1.14, r2=0.12).15 In the case of TNPO3-IRF5 region, we selected variants from Viatte et al12 and Sigurdsson et al16 that were included in the combined analysis. The rs12531711 showed a remarkably strong association with ACPA− RA in the combined analysis (p=4.35×10−5 OR=1.30).

    Discussion

    This report comprises the largest ACPA− RA cohort genotyped and analysed with a genome-wide platform to date (1922 ACPA− RA patients and 7087 non-affected controls). Our study identified a suggestive new risk factor for this condition, (CLYBL), confirmed the association of the HLA region with this subset and proposed a two-variant model including the HLA-DRB1 and HLA-B loci to explain the observed HLA association peak.

    The first GWAS that focused on the ACPA+ versus ACPA− genetic differences did not show evidence for an association of ACPA− RA with the HLA region in a cohort of Swedish RA patients.13 Nevertheless, recent reports have suggested a previously unidentified role for the HLA region.12 ,14 Taking into account that the Immunochip platform includes a dense mapping of the HLA region, and that the cohorts analysed in the present study are partially overlapping with the populations in Eyre et al, we consider that the resemblance of the results in genome-wide and Immunochip platforms confirm previous findings, and indicate that the coverage of the HLA region in our study was appropriate.14 Moreover, in Eyre et al, the five-amino acid model which was described by Raychaudhuri et al21 and accounted for the observed association in the HLA region was confirmed.14 Nevertheless, this novel HLA allele and amino acid residue imputation approach was not applied separately in the ACPA− group. In our ACPA− study, we used this approach to detect two-amino acid position (HLA-B at position 9 and HLA-DRB1 at position 11), as the most parsimonious model explaining the observed association in the HLA region in our data. The two of the identified amino acid position, HLA-DRB1 11th amino acid and HLA-B 9th amino acid, were located in the peptide binding groove of their corresponding HLA molecules.21 Moreover, both amino acid positions were shared between our model and the ACPA+ RA model proposed by Raychaudhuri et al21 and confirmed by Eyre et al14 However, no significant association was found for the SE variants or other ACPA+-related HLA-DRB1 amino acid positions (HLA-.DRB1 13, 71 and 74 positions), which can be due to real association divergences in these subgroups. Moreover, given the observed association of HLA-B, it can be hypothesised that HLA class I has a more relevant role in ACPA− RA than in ACPA+ RA, while the HLA class II association seems more predominant ACPA+ RA. However, the partial overlap and the differences between the proposed HLA models for both RA serological subgroups would require further studies to clarify if the HLA-DRB1 and HLA-B signals belong to both RA subgroups or are restricted to ACPA+ or ACPA−, respectively. Additionally, individuals showing borderline ACPA titres and seroconversion from negative to positive ACPA is a rare event (only 2%), but it may act as a confounding factor.28

    In relation to the novel candidate association with a SNP in the CLYBL gene, although the signal is slightly below the genome-wide significance level, our findings might support future studies to clarify the veracity of this association. CLYBL encodes a citrate lyase subunit β-like protein which has citrate (pro-3S)-lyase and ion-binding activities that are transported to the mitochondria. Interestingly, a SNP in the CLYBL locus has been recently associated with low serum levels of vitamin B12 (also known as cyanocobalamin).29 ,30 Although the relation between CLYBL and the vitamin B12 serum levels remains unclear, a link between vitamin B12 and RA has been described. In a report by Regal et al, 24% of RA patients had low vitamin B12 serum levels and 26% of them had true vitamin B12 deficiency.31 Besides, a previous study described that RA patients with vitamin B12 deficiency, folate deficiency, vitamin B6 deficiency and impaired renal function showed associated hyperhomocysteinemia that may have a role in promoting high cardiovascular morbidity in patients with RA.32 Moreover, significant associations with the familial risk of RA in offspring according to parental proband were reported for pernicious anaemia, which is usually the result of intrinsic factor insufficient secretion and consequent vitamin B12 deficiency.33

    Additionally, the design of our study led to some limitations mainly due to the variety of genotyping platforms used. Therefore, some associations might have been overlooked due to a limited coverage in the region despite the imputation step (ie, the ANKRD55/IL6ST locus). In spite of the large size of the ACPA− patient cohort, our study might have been underpowered to detect modest associations. Additionally, our conclusions may be applicable only in populations of European ancestry.

    The present report analysed the genetic component of a large cohort of ACPA− RA patients compared to non-affected controls following a genome-wide strategy, replicated previous findings in different non-HLA loci, such as TNFAIP3, GIN1/C5orf30, PTPN22, CLEC16A and TNPO3/IRF5, and revealed a novel suggestive susceptibility gene, CLYBL. Moreover, our study provided a deep insight into the influence of the HLA region in ACPA− RA and identified a two-amino acid residue model explaining this association. The present study together with previous evidence supported the existence of an ACPA-specific genetic component, and highlighted the importance of comprehensive genetic analysis of large ACPA− cohorts.

    References

    1. ↵
      1. McInnes IB,
      2. Schett G
      . The pathogenesis of rheumatoid arthritis. N Engl J Med 2011;365:2205–19.
      OpenUrlCrossRefPubMedWeb of Science
    2. ↵
      1. van Venrooij WJ,
      2. van Beers JJ,
      3. Pruijn GJ
      . Anti-CCP antibodies: the past, the present and the future. Nat Rev Rheumatol 2011;7:391–8.
      OpenUrlCrossRefPubMed
    3. ↵
      1. Huizinga TW,
      2. Amos CI,
      3. van der Helm-van Mil AH,
      4. et al
      . Refining the complex rheumatoid arthritis phenotype based on specificity of the HLA-DRB1 shared epitope for antibodies to citrullinated proteins. Arthritis Rheum 2005;52:3433–8.
      OpenUrlCrossRefPubMedWeb of Science
    4. ↵
      1. Lundstrom E,
      2. Kallberg H,
      3. Alfredsson L,
      4. et al
      . Gene-environment interaction between the DRB1 shared epitope and smoking in the risk of anti-citrullinated protein antibody-positive rheumatoid arthritis: all alleles are important. Arthritis Rheum 2009;60:1597–603.
      OpenUrlCrossRefPubMedWeb of Science
    5. ↵
      1. van Gaalen FA,
      2. Linn-Rasker SP,
      3. van Venrooij WJ,
      4. et al
      . Autoantibodies to cyclic citrullinated peptides predict progression to rheumatoid arthritis in patients with undifferentiated arthritis: a prospective cohort study. Arthritis Rheum 2004;50:709–15.
      OpenUrlCrossRefPubMedWeb of Science
    6. ↵
      1. Aho K,
      2. Heliovaara M,
      3. Maatela J,
      4. et al
      . Rheumatoid factors antedating clinical rheumatoid arthritis. J Rheumatol 1991;18:1282–4.
      OpenUrlPubMedWeb of Science
    7. ↵
      1. Nielen MM,
      2. van Schaardenburg D,
      3. Reesink HW,
      4. et al
      . Specific autoantibodies precede the symptoms of rheumatoid arthritis: a study of serial measurements in blood donors. Arthritis Rheum 2004;50:380–6.
      OpenUrlCrossRefPubMedWeb of Science
    8. ↵
      1. Rantapaa-Dahlqvist S,
      2. de Jong BA,
      3. Berglin E,
      4. et al
      . Antibodies against cyclic citrullinated peptide and IgA rheumatoid factor predict the development of rheumatoid arthritis. Arthritis Rheum 2003;48:2741–9.
      OpenUrlCrossRefPubMedWeb of Science
    9. ↵
      1. van der Helm-van Mil AH,
      2. Verpoort KN,
      3. Breedveld FC,
      4. et al
      . Antibodies to citrullinated proteins and differences in clinical progression of rheumatoid arthritis. Arthritis Res Ther 2005;7:R949–58.
      OpenUrlCrossRefPubMedWeb of Science
    10. ↵
      1. Aletaha D,
      2. Neogi T,
      3. Silman AJ,
      4. et al
      . 2010 rheumatoid arthritis classification criteria: an American College of Rheumatology/European League Against Rheumatism collaborative initiative. Ann Rheum Dis 2010;69:1580–8.
      OpenUrlAbstract/FREE Full Text
    11. ↵
      1. Kurreeman F,
      2. Liao K,
      3. Chibnik L,
      4. et al
      . Genetic basis of autoantibody positive and negative rheumatoid arthritis risk in a multi-ethnic cohort derived from electronic health records. Am J Hum Genet 2011;88:57–69.
      OpenUrlCrossRefPubMed
    12. ↵
      1. Viatte S,
      2. Plant D,
      3. Bowes J,
      4. et al
      . Genetic markers of rheumatoid arthritis susceptibility in anti-citrullinated peptide antibody negative patients. Ann Rheum Dis 2012;71:1984–90.
      OpenUrlAbstract/FREE Full Text
    13. ↵
      1. Padyukov L,
      2. Seielstad M,
      3. Ong RT,
      4. et al
      . A genome-wide association study suggests contrasting associations in ACPA-positive versus ACPA-negative rheumatoid arthritis. Ann Rheum Dis 2011;70:259–65.
      OpenUrlAbstract/FREE Full Text
    14. ↵
      1. Eyre S,
      2. Bowes J,
      3. Diogo D,
      4. et al
      . High-density genetic mapping identifies new susceptibility loci for rheumatoid arthritis. Nat Genet 2012;44:1336–40.
      OpenUrlCrossRefPubMed
    15. ↵
      1. Skinningsrud B,
      2. Lie BA,
      3. Husebye ES,
      4. et al
      . A CLEC16A variant confers risk for juvenile idiopathic arthritis and anti-cyclic citrullinated peptide antibody negative rheumatoid arthritis. Ann Rheum Dis 2010;69:1471–4.
      OpenUrlAbstract/FREE Full Text
    16. ↵
      1. Sigurdsson S,
      2. Padyukov L,
      3. Kurreeman FA,
      4. et al
      . Association of a haplotype in the promoter region of the interferon regulatory factor 5 gene with rheumatoid arthritis. Arthritis Rheum 2007;56:2202–10.
      OpenUrlCrossRefPubMedWeb of Science
    17. ↵
      1. Perricone C,
      2. Ceccarelli F,
      3. Valesini G
      . An overview on the genetic of rheumatoid arthritis: a never-ending story. Autoimmun Rev 2011;10:599–608.
      OpenUrlCrossRefPubMed
    18. ↵
      1. Arnett FC,
      2. Edworthy SM,
      3. Bloch DA,
      4. et al
      . The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum 1988;31:315–24.
      OpenUrlCrossRefPubMedWeb of Science
    19. ↵
      1. Skol AD,
      2. Scott LJ,
      3. Abecasis GR,
      4. et al
      . Joint analysis is more efficient than replication-based analysis for two-stage genome-wide association studies. Nat Genet 2006;38:209–13.
      OpenUrlCrossRefPubMedWeb of Science
    20. ↵
      1. Purcell S,
      2. Neale B,
      3. Todd-Brown K,
      4. et al
      . PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 2007;81:559–75.
      OpenUrlCrossRefPubMed
    21. ↵
      1. Raychaudhuri S,
      2. Sandor C,
      3. Stahl EA,
      4. et al
      . Five amino acids in three HLA proteins explain most of the association between MHC and seropositive rheumatoid arthritis. Nat Genet 2012;44:291–6.
      OpenUrlCrossRefPubMed
    22. ↵
      1. Achkar JP,
      2. Klei L,
      3. de Bakker PI,
      4. et al
      . Amino acid position 11 of HLA-DRbeta1 is a major determinant of chromosome 6p association with ulcerative colitis. Genes Immun 2012;13:245–52.
      OpenUrlCrossRefPubMedWeb of Science
    23. ↵
      1. Jia X,
      2. Han B,
      3. Onengut-Gumuscu S,
      4. et al
      . Imputing amino Acid polymorphisms in human leukocyte antigens. PLoS ONE 2013;8:e64683.
      OpenUrlCrossRefPubMed
    24. ↵
      1. Patterson N,
      2. Price AL,
      3. Reich D
      . Population structure and eigenanalysis. PLoS Genet 2006;2:e190.
      OpenUrlCrossRefPubMed
    25. ↵
      1. Price AL,
      2. Patterson NJ,
      3. Plenge RM,
      4. et al
      . Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet 2006;38:904–9.
      OpenUrlCrossRefPubMedWeb of Science
    26. ↵
      1. Pruim RJ,
      2. Welch RP,
      3. Sanna S,
      4. et al
      . LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics 2010;26:2336–7.
      OpenUrlAbstract/FREE Full Text
    27. ↵
      1. Team RC
      . R: A language and environment for statistical computing. R Foundation for Statistical Computing 2012.
    28. ↵
      1. Ursum J,
      2. Bos WH,
      3. van Dillen N,
      4. et al
      . Levels of anti-citrullinated protein antibodies and IgM rheumatoid factor are not associated with outcome in early arthritis patients: a cohort study. Arthritis Res Ther 2010;12:R8.
      OpenUrlCrossRefPubMed
    29. ↵
      1. Lin X,
      2. Lu D,
      3. Gao Y,
      4. et al
      . Genome-wide association study identifies novel loci associated with serum level of vitamin B12 in Chinese men. Hum Mol Genet 2012;21:2610–17.
      OpenUrlAbstract/FREE Full Text
    30. ↵
      1. Grarup N,
      2. Sulem P,
      3. Sandholt CH,
      4. et al
      . Genetic architecture of vitamin B12 and folate levels uncovered applying deeply sequenced large datasets. PLoS Genet 2013;9:e1003530.
      OpenUrlCrossRefPubMed
    31. ↵
      1. Segal R,
      2. Baumoehl Y,
      3. Elkayam O,
      4. et al
      . Anemia, serum vitamin B12, and folic acid in patients with rheumatoid arthritis, psoriatic arthritis, and systemic lupus erythematosus. Rheumatol Int 2004;24:14–19.
      OpenUrlCrossRefPubMed
    32. ↵
      1. Pettersson T,
      2. Friman C,
      3. Abrahamsson L,
      4. et al
      . Serum homocysteine and methylmalonic acid in patients with rheumatoid arthritis and cobalaminopenia. J Rheumatol 1998;25:859–63.
      OpenUrlPubMedWeb of Science
    33. ↵
      1. Hemminki K,
      2. Li X,
      3. Sundquist J,
      4. et al
      . Familial associations of rheumatoid arthritis with autoimmune diseases and related conditions. Arthritis Rheum 2009;60: 661–8.
      OpenUrlCrossRefPubMedWeb of Science
    View Abstract

    Supplementary materials

    • Supplementary Data

      This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

      Files in this Data Supplement:

    Footnotes

    • Handling editor Tore K Kvien

    • Correction notice This article has been corrected since it was published Online First. The initials and family names of some authors were in inverse order, such that family name was listed first followed by the initials; these have been corrected.

    • Acknowledgements We thank Sofía Vargas and Sonia García for their excellent technical assistance. We thank Banco Nacional de ADN (University of Salamanca, Spain), which supplied part of the control DNA samples, and we thank all patients and donors for their collaboration. Principal investigators of the Nijmegen Biomedical Study are LALM Kiemeney, M den Heijer, ALM Verbeek, DW Swinkels and B Franke.

    • Contributors MJHC, PPT, MDP, AZ, TH, AHMHM, GSR, LRR, AB, IGA, M AGG, CGV, SHV, BF, LLCS, IEHB, BACD, GJW, RAO, JF, PLCMR, LK, JBAC, EB, NV, RT, JM, and LP provided samples and clinical information, or control genotypes. BPCK, LP, LK, JM, NV, RT, designed the GWAS. LBC, CK, HK RS, AI performed genotyping quality control and association analyses of the GWAS data. The paper was written by LBC, CK and BPCK with contributions from all other authors. BPCK and CK coordinated the genotyping, analysis, and manuscript writing efforts for this study.

    • Funding This work was supported by the Dutch Arthritis Foundation (grant DAA 0901069) a Spanish grant from RETICS Program, RD12/0009/0004 (RIER) from Instituto de Salud Carlos III (ISCIII) and also by the grant BTCure for Rheumatoid Arthritis from the European IMI Program. Genotyping of the Dutch control samples was sponsored by US National Institutes of Mental Health funding, R01 MH078075 (ROA)

    • Competing interests None.

    • Patient consent Obtained.

    • Ethics approval Ethics Committee of participating academic centers.

    • Provenance and peer review Not commissioned; externally peer reviewed.