Genetic susceptibility to complex diseases has been intensively studied during the last decade, yet only signals with small effect have been found leaving open the possibility that subgroups within complex traits show stronger association signals. In rheumatoid arthritis (RA), autoantibody production serves as a helpful discriminator in genetic studies and today anti-citrullinated cyclic peptide (anti-CCP) antibody positivity is employed for diagnosis of disease. The HLA-DRB1 locus is known as the most important genetic contributor for the risk of RA, but is not sufficient to drive autoimmunity and additional genetic and environmental factors are involved. Hence, we addressed the association of previously discovered RA loci with disease-specific autoantibody responses in RA patients stratified by HLA-DRB1*04.
We investigated 2178 patients from three RA cohorts from Sweden and Spain for 41 genetic variants and four autoantibodies, including the generic anti-CCP as well as specific responses towards citrullinated peptides from vimentin, alpha-enolase and type II collagen.
Our data demonstrated different genetic associations of autoantibody-positive disease subgroups in relation to the presence of DRB1*04. Two specific subgroups of autoantibody-positive RA were identified. The SNP in PTPN22 was associated with presence of anti-citrullinated enolase peptide antibodies in carriers of HLA-DRB1*04 (Cochran-Mantel-Haenszel test P = 0.0001, Pcorrected <0.05), whereas SNPs in CDK6 and PADI4 were associated with anti-CCP status in DRB1*04 negative patients (Cochran-Mantel-Haenszel test P = 0.0004, Pcorrected <0.05 for both markers). Additionally we see allelic correlation with autoantibody titers for PTPN22 SNP rs2476601 and anti-citrullinated enolase peptide antibodies in carriers of HLA-DRB1*04 (Mann Whitney test P = 0.02) and between CDK6 SNP rs42041 and anti-CCP in non-carriers of HLA-DRB1*04 (Mann Whitney test P = 0.02).
These data point to alternative pathways for disease development in clinically similar RA subgroups and suggest an approach for study of genetic complexity of disease with strong contribution of HLA.
The study of complex diseases has revealed complicated patterns of inheritance  in which involvement of multiple variants and environmental conditions, as well as gene-gene and gene-environmental interactions, have made the discovery of genetic causes for complex diseases a challenging task. One main methodological caveat is the definition of phenotypes. As clinical features of complex diseases are heterogeneous, one cannot expect simple genetic patterns to explain such complex heterogeneity. By stratifying for more homogeneous disease subgroups/phenotypes one may try to achieve stronger association and better understanding of disease mechanisms ,.
Rheumatoid arthritis (RA) is a common complex disease with incidence around 0.5 to 1.0% in different populations , and a steadily growing prevalence in aging Western and Eastern societies. Due to implementation of several targeted treatments in combination with immunosuppressive therapy the overall quality of life for RA patients has dramatically improved and many patients may stay in remission for long periods of time . However, current treatments mostly decelerate rather than cure disease, while disease-driving mechanisms remain poorly understood ,. Despite the discovery of multiple genetic associations with RA, variants within the HLA-DRB1 locus, (that is, HLA-DRB1*04 alleles), remain the most significant contributors to the risk of autoantibody-positive RA –.
Antibody responses against citrullinated epitopes of vimentin, fibrinogen, type-II collagen, alpha-enolase represent specific features of RA ,,,. Patients may display antibodies to one or several of these modified self-proteins. Most of these autoantibody specificities are confined within, and could thus be recognized on testing as anti-citrullinated cyclic peptide (CCP) antibodies ,, which is a generic test recognizing most anti-citrulinated peptide antibody (ACPAs). A strong association between the shared epitope (SE) HLA-DRB1 alleles, specifically for SE HLA-DRB1*04 alleles, and development of both anti-CCP and anti-citrullinated alpha-enolase peptide-1 (CEP-1) antibodies has been reported . Subsets of RA defined by other combinations of antibodies to citrullinated autoantigens have been shown to display very different degrees of association with the common HLA-DRB1 risk alleles that constitute the group of shared epitope alleles ,. These different profiles of ACPAs may reflect distinct biological and immunological courses that are determined by genetic niches of susceptibility. Although it is clinically challenging to differentiate these subgroups due to rather similar symptoms, it is essential to dissect these in view of diagnosis and treatment, and ultimately for the understanding of disease mechanisms and possible prevention.
To find additional links between genetics and serology of RA, which will allow for studies of the relationship between genotypes and phenotypes, we employed a large population-based study from Sweden, the Epidemiological Investigation of Rheumatoid Arthritis (EIRA) with incident cases of RA and two smaller cohorts of cases with well-established/chronic RA. We also considered known relative stability of anti-CCP levels through the RA development ,. Our hypothesis is that a polymorphism outside the HLA-DRB1 locus may contribute and shape the development of certain serological subtypes of RA, which are otherwise clinically indistinguishable. We found that contribution of different non-HLA single nucleotide polymorphisms (SNPs) associated with RA in the development of distinct ACPAs does not overlap, and may define distinct subgroups of disease with an SE-positive or an SE-negative background.
Patients and healthy subjects
DNA and serum samples were collected from three independent cohorts from Sweden and Spain (Additional file 1: Table S13); all experiments were performed in accordance with the Declaration of Helsinki and were approved by Stockholm Ethical Review board or Comite Etico de Investigacion Clinica de Galicia and all subjects gave informed consent. Serum and DNA samples were stored at −80°C until use. This was a case-case study to test the hypothesis of the contribution of different genetic factors in the development of serologically determined subgroups of RA.
Additional file 1: Table S1.. Association between anti-citrullinated cyclic peptide (anti-CCP) and single nucleotide polymorphisms (SNPs) in cohort 1. Table S2. Association between anti-Cit-Vimentin and SNP in cohort 1. Table S3. Association between anti-citrullinated alpha-enolase peptide-1 (CEP-1) and SNP in cohort 1. Table S4. Association between anti-CitC1III and SNP in cohort 1. Table S5. Association between anti-CCP and SNP in cohort 2. Table S6. Association between anti-Cit-Vimentin and SNP in cohort 2. Table S7. Association between anti-CEP1 and SNP in cohort 2. Table S8. Association between anti-CitC1 and SNP in cohort 2. Table S9. Association between anti-CCP and SNP in cohort 3. Table S10. Association between anti-Cit-Vimentin and SNP in cohort 3. Table S11. Association between anti-CEP-1 and SNP in cohort 3. Table S12. Complete meta-analysis results. Table S13. Baseline characteristics of three analyzed RA cohorts. Table S14. Genetic variants used in the study.
Format: XLSX Size: 514KB Download file
The cohorts included were: 1) For initial study we analyzed a cohort of 1,362 patients with incident RA (cohort 1) from a population-based case–control study (EIRA) ,. The details of the EIRA study have been described previously . Briefly, a case was defined as a person in the study base who received a new diagnosis of RA from a rheumatologist (within 1 year after the onset of symptoms in 85% of the cases) and fulfilled the American College of Rheumatology 1987 criteria for the classification of RA . Cases were recruited from all public and a majority of private rheumatology units in the study area; 2) cohort 2 comprised 379 patients with established RA, who all fulfilled the American College of Rheumatology criteria  and were attending the Rheumatology Clinic at the Karolinska University Hospital, Stockholm, Sweden; 3) cohort 3 comprised 437 patients with established RA classified according to the 1987 American College of Rheumatology criteria and of Spanish ancestry; DNA and serum samples from these patients were obtained from a single hospital. Their clinical characteristics have been already described in full .
Detection of IgG anti-CCP antibodies
ELISA Anti-CCP2 test (Immunoscan RA, Mark 2, Euro-Diagnostica, Malmö, Sweden) was used to determine the levels of anti-CCP IgG antibodies. Quantification of the results and use of a cut-off value of 25 U/ml were according to the manufacturer’s instructions.
Detection of IgG antibodies against citrullinated alpha-enolase peptide-1 (CEP-1, citrullinated type-II collagen (citC1III) and citrullinated vimentin (cit-vim)
DNA was extracted from ethylenediaminetetraacetic acid (EDTA)-treated blood by the salting-out method . Genotyping for HLA-DRB1 allotypes was conducted using the sequence-specific primer-PCR method (DR low-resolution kit; Olerup SSP, Saltsjöbaden, Sweden) as previously described . DRB1*04- and DRB1*01-positive patients were further subtyped by Olerup SSP DRB1*04 (Olerup SSP) subtyping kits, respectively ,. In replication cohort 3, HLA-DBR1 alleles were genotyped by a sequencing-based typing method using the AlleleSEQR HLA-DRB1 Typing kit (Abbott Diagnostics, Abbott Park, Germany), which includes bidirectional sequencing of the second exon of DRB1. Ambiguous samples were additionally sequenced with group-specific primers (AlleleSEQR HLA-DRB1 GSSP, Abbott Diagnostics, Abbott Park, Germany).
Genotyping for RA-predisposing allelic variants
The 41 SNPs in the first phase of this study were selected from previous studies of genetic risk factors for RA; references for each SNP are presented in Additional file 1: Table S14. These 41 SNPs were studied in the two first cohorts of patients. Genotypes corresponding to the EIRA cohort were retrieved from a genome-wide association study (GWAS) (26 SNPs) , from dense mapping of the HLA region (4 SNPs)  or from additional genotyping (11 SNPs) that was performed by TaqMan allelic discrimination assay (Applied Biosystems, Foster City, CA, USA). Genotyping of cohort-2 was performed using a 64-OpenArray platform (Applied Biosystems) with the chip-based TaqMan genotyping technology. Genotyping was performed according to manufacturer’s instructions, and genotype calls were made using AutoCaller (Applied Biosystems). All SNPs that were analyzed in this study are listed in Additional file 1: Table S14.
With the results of the two first cohorts we selected 12 SNPs for further analysis. The criteria of selection were either a P-value <0.01 for association of the SNP with a specific subgroup of autoantibody-positive patients in any of the two cohorts, or to show the same direction of association and a P-value <0.2 in the two cohorts. Genotyping for the 12 selected SNPs in cohort 3 (Spanish patients) was performed by PCR amplification followed by single-base extension with the SNaPshot Multiplex Kit (Applied Biosystems). Details of the protocol are available by the authors upon request. SNP rs2476601 was genotyped with a Taqman SNP genotyping assay (Applied Biosystems, Carlsbad, CA, USA).
We performed analyses either for the whole patient cohort, or for separate HLA-DRB1*04-positive patients and HLA-DRB1*04-negative patients and further calculated the odds ratios (OR) and 95% CI for the association with anti-CCP and for three ACPA specificities. A dominant genetic model or carrier model that assess the effect of the presence of the risk allele, irrespective of the presence of one or two copies, was used for all analyses to avoid too few counts of some patient subgroups. The OR were calculated by median-unbiased estimation and exact CI using the mid-P method  with patients without the specific type of antibody as the reference group.
The three cohorts were jointly analyzed (citC1III data were only available for two cohorts) using the Cochran-Mantel-Haenszel (CMH) test. An adjusted P-value of 0.05 after Bonferroni correction for multiple testing was considered statistically significant in the meta-analysis. Homogeneity between cohorts was tested by Tarone's test. As we detected no significant heterogeneity between cohorts a CMH fixed-effect model was fitted using the restricted maximum-likelihood estimator. The metafor R package was used for these analyses . Antibody levels and disease activity score in 28 joints (DAS28) were compared by the Mann–Whitney test.
We determined the associations between genotypes and phenotypes by performing a three-step analysis that is represented in the analytical workflow chart in Figure 1. At the first step, by analysis of cohort 1 and replication in cohort 2, we selected 12 SNPs for further replication. In a second step we validated the results in an independent cohort (cohort 3) and in the third step we estimated OR for RA with certain serologic profiles, based on data from all three cohorts, by meta-analysis. Our main goal was to detect possible associations of non-DRB1 SNPs with autoantibody-defined subgroups in relation to presence/absence of HLA-DRB1*04 alleles. However, we also performed analysis for the whole group of patients without this stratification.
Figure 1. Visual description of the analysis. SNP, single nucleotide polymorphism.
Association study in cohort 1 and cohort 2
Cohort 1 comprised 1,362 RA patients from the Swedish EIRA study: 70% of the patients were positive for anti-CCP antibodies, whereas defined ACPA fine specificities were less frequent; 50% of the patients were positive for anti-CEP-1 and 42% were positive for anti-cit-Vim and for anti-citC1III antibodies (Additional file 1: Table S13). We performed a case-case association study to test the hypothesis of different genetic backgrounds between RA subgroups identified by autoantibody specificity. Multiple SNPs were associated in the whole cohort or after stratification by the major SE allele, HLA-DRB1*04, but only few passed conservative corrections for multiple testing (P ≤0.0004, Additional file 1: Tables S1-S4), namely: SNP rs6457617 (HLA-DQB1 gene) for several specificities, for the whole cohort and also after stratification for HLA-DRB1*04, and rs2476601 (PTPN22 gene) for the HLA-DRB1*04-positive RA group with anti-CEP-1 autoantibody.
With the goal to gain better statistical power and to validate the associations identified, we performed analysis in cohort 2, an independent cohort of 379 Swedish RA patients with established disease (Additional file 1: Table S13). Most patients in this cohort were positive for anti-CCP (73%) and for anti-citVim (58%). The other two ACPAs were less frequent (46% for anti-CEP-1 and 38% for anti-citC1III). Again, although multiple SNPs demonstrated a moderate trend towards association, very few passed correction for multiple comparisons (Additional file 1: Tables S5-S8). Therefore, we choose to select the best candidate SNPs from association studies in cohorts 1 and 2 for further validation. Two criteria were used for selection, either a P-value <0.01 in any of the two cohorts, or a change in the same direction with P <0.2 in the two cohorts for the same comparison. This procedure generated a list of 12 SNPs that were used at the next step. Extracted data for association with these SNPs from analyses of cohorts 1 and 2 are presented in Table 1.
Table 1. Selected associations found in cohorts 1 and 2
Association study in cohort 3
The 12 SNPs were genotyped in an independent collection of patients of Spanish ancestry, who had established RA (n = 437). They showed a similar fraction of anti-CCP positivity, but a lower frequency of the ACPAs (Additional file 1: Table S13). The analysis revealed moderate levels of association for SNPs from PTPN22, PADI4, CDK6 and MICA in some disease subgroups (Additional file 1: Tables S9-S11). Importantly, the direction of associations followed the same pattern as in cohorts 1 and 2.
Meta-analysis of the three cohorts
Following examination for heterogeneity, a fixed-effect meta-analysis was conducted for the association of these SNPs with anti-citrulline immunity in RA using data available from the three cohorts; Table 2 and Figure 2 provide summarized results, which were Bonferroni-corrected for the number of tests done. Additional file 1: Table S12 provides the full disclosure of the analysis.
Table 2. Association of genetic markers with serological subgroups of rheumatoid arthritis: meta-analysis summary
Figure 2. Overview of the most significant associations in the study. The bars represent the odds ratio on the horizontal scale with the 99.999% CI for comparison between group of persons with and without autoantibody (distinguished by color of the bar) and carriers and non-carriers of the major allele with separation to HLA-DRB1*04-positive and HLA-DRB1*04-negative. Results are presented for the final selection of SNPs and meta-analysis of three cohorts (two cohorts for anti-CitC1). Indicated P-values are the values after Bonferroni correction. Full details are presented in Additional file 1: Table S12.
The results following stratification into carriers of HLA-DRB1*04 and non-carriers of HLA-DRB1*04 showed associations of the SNP in PTPN22 with antibodies to CEP-1 in the HLA-DRB1*04-positive group (adjusted P-value 0.007), and of two SNPs, in PADI4 and in CDK6, with anti-CCP response in the HLA-DRB1*04-negative group (adjusted p-value 0.038 and 0.034, respectively). In addition, two SNPs showed a trend to association with P-value <0.1 in the non-carriers of HLA-DRB1*04: CDK6 and PRKCQ with anti-citVim and anti-CEP-1 respectively (Table 2, Additional file 1: Table S12, and Figure 2). These data additionally endorse the association of RA subgroups with the CDK6 polymorphism and suggest PRKCQ as a relevant candidate for follow-up studies.
In the analyses without stratification according to HLA-DRB1*04, the SNP assigned to MICA was associated with response to cit-Vim (adjusted P-value 0.008), the PTPN22 SNP was associated with antibodies to CEP-1 (adjusted P-value 0.012), and the HLA-DQB1 SNP was associated with the anti-citC1III antibodies (adjusted P = 9 × 10−5) (Additional file 1: Table S12). These findings suggest a relative independence of the association with the MICA SNP from HLA-DRB1*04 as the OR are very uniform in the two groups of patients stratified by this allele (OR = 1.4 in the carriers and in non-carriers). This independence was not present in the association with the PTPN22 SNP that showed a stronger significant association in the carriers of DRB1*04 than in the non-carriers (OR = 1.7 and 1.1, respectively).
Analysis of serological and clinical data
As our data showed that SNPs in PADI4, PTPN22 and CDK6 together with HLA-DRB1*04 carrier status were associated with the presence of some autoantibodies, we tested if the genotypes were correlated with the titers of these antibodies. We saw a coherent significant difference for anti-CCP antibody levels between homozygous CC patients and those carrying the G allele of rs42041 from CDK6 in cohort 1 and 2 for patients without DRB1*04 alleles (P = 0.0193 and 0.0290 for cohorts 1 and 2, Additional file 2: Figure S1a). Similarly, patients having the AA/AG variant of rs2476601 from PTPN22 had higher levels of antibodies towards cit-Vim in cohort 1 and 2 but only in patients with DRB1*04 alleles (P = 0.0238, Additional file 2: Figure S1c).
Additional file 2: Figure S1.. Levels of autoantibodies in relation to CDK6, PADI4 and PTPN22 genotypes. Subjects from cohort 1 (upper panel) and cohort 2 (lower panel) were classified negative or positive to HLA-DRB1*04 and according to their CDK6, PADI4 or PTPN22 genotypes. Antibody levels against CCP (CDK6, PADI4) and cit-Vim (PTPN22) were compared between the different genotype groups. Bars indicate the median levels of antibody response. *P <0.05.
Format: JPEG Size: 77KB Download file
To compare genetic subgroups defined by certain autoantibody specificity, HLA-DRB1*04 and non-HLA alleles (PTPN22, PADI4 and CDK6), we analyzed data on disease activity at baseline in cohort 1 measured by the DAS28, but no differences were found (data not shown).
In this study we have demonstrated that specificity of autoantibody responses in seropositive RA is regulated not only by HLA-DRB1*04, but also by at least three additional non-HLA genetic variants: PADI4, PTPN22 and CDK6. Our data suggest that PTPN22 involvement is restricted to HLA-DRB1*04 carriers, while associations of PADI4 and CDK6 are restricted to the non-carriers, and they also differ in autoantibody specificity. While concordance of HLA-DRB1*04 and PTPN22 SNP in anti-CEP-1 positive RA was demonstrated previously in a single cohort ,, the involvement of genetic variants of PADI4 and CDK6 in risk of development of RA subgroups only in the absence of HLA-DRB1*04 is a new finding. We also have observed a difference in the autoantibody titers corresponding with the previous associations, (Additional file 2: Figure S1) while no difference was found in relation to clinical activity of disease at baseline measured by the DAS28.
Genetic and environmental factors are involved in the etiology of complex diseases. In RA, such susceptibility factors have been primarily found for ACPA-positive disease ,,. Among RA-predisposing genes, the involvement of the HLA-DRB1 polymorphism has been extensively studied, as it was the first identified and it was also the strongest to date ,. Association of the shared epitope alleles of HLA-DRB1 with ACPA (anti-CCP)-positive RA indicates probable mechanisms for disease development. In this work, we further examined the significance of other, mainly non-HLA variants, for the development of different autoantibodies in a case-case study. We investigated the effects of these genetic variants alone or in combination with HLA-DRB1*04, the most frequent SE alleles in European Caucasians, in three independent RA cohorts. Although the selection of SNPs was based on RA susceptibility studies our results are not related to susceptibility by itself, because no comparison with healthy individuals was done.
Despite the fact that HLA-DRB1*04 alleles are strongly associated with RA and anti-citrulline immunity, few studies have so far addressed the contribution of these allelic variants to RA development via antigen presentation of citrullinated peptides –. It is not clear why and how HLA-DRB1*04-positive RA patients develop ACPAs with different specificities towards non-overlapping epitopes. It is therefore possible that the HLA-DR*04 alleles (mainly *04:01) together with additional genes outside this region work collectively to shape ACPA responses in RA. In this study, we found three allelic variants outside HLA-DRB1 that associate with a particular type of autoantibody in RA. These genetic differences may provide a plausible explanation for the phenotypic heterogeneity of RA, for example, differences in the autoantibody profile of RA patients. However, the precise functional contribution of these variants to RA development remains to be determined.
The PTPN22 gene (MIM ID*600716) encodes an 807 amino acid lymphoid tyrosine phosphatase of non-receptor type (LYP). An SNP in this gene was found to be associated with several autoimmune diseases . This SNP is a non-synonymous change of arginine for tryptophan at position 620 of the protein. The biochemical effect and the resulting disease risk of this SNP have been very difficult to study and there is not yet consensus about its functional consequences. It has been suggested that this SNP determines a gain of function of the LYP phosphatase, but there other studies showing a loss of function of LYP when analyzed together with the c-Src tyrosine kinase . Therefore, it is still unclear how the R620W change might contribute to autoimmunity. In addition, most studies of PTPN22 function refer to T cell function either at the level of thymic selection or of mature T cell activation, but B cells also express LYP and the removal of autoreactive B cells is deficient in subjects bearing the risk allele of PTPN22, making for alternative mechanisms of disease risk .
The CDK6 gene (MIM ID* 603368) encodes a 326 amino acid serine/threonine-protein kinase, which is ubiquitously expressed and involved in cell cycle control and differentiation by promoting G1/S transition. A genetic polymorphism at CDK6 was recently found to be associated with RA  and with joint destruction during RA . Interestingly, polymorphisms in this gene are also associated with adult height  and with white blood cells counts in some populations ,.
The PADI4 gene (MIM ID *605347) codes for peptidyl arginine deiminase, type 4, which mediates citrullination of different proteins, including histones. In histones, citrullination affects methylation and this way PADI4 influences gene expression . PADI4 variants have previously been observed to be associated with RA in the Japanese population  and recently in European Caucasians . Although the enzyme is directly involved in citrullination of proteins, until now genetic association has not been restricted to RA patients positive for autoantibodies to citrullinated antigens .
There were three loci showing significant association with particular autoantibodies in the analysis non-stratified by HLA-DRB1*04. The associations between the anti-citC1III antibodies and the HLA-DQB1 SNP and between the anti-citVim antibodies and the SNP in MICA could reflect linkage disequilibrium with the shared epitope alleles of HLA-DRB1 because these SNPs are in the HLA region. However, we cannot exclude an effect of other variants in this region as has been shown for RA susceptibility . The third association, between anti-CEP-1 antibodies and the PTPN22 SNP and anti-CEP-1 antibodies, although observed in the non-stratified analysis was not independent of the presence of HLA-DRB1*04. It was much stronger in the carriers of HLA-DRB1*04 than in the non-carriers, as already discussed.
What are the possible mechanisms underlying the identified associations and how will our findings affect the understanding of, or the management of RA? It is clear that HLA-DRB1*04, as well as other SE alleles, do not represent a single determinant for development of autoimmune responses. We previously demonstrated the possibility of gene-gene and gene-environment interactions in the development of RA which suggested that smoking and several non-HLA genes could serve as possible second determinants for a break in tolerance to self antigens and the development of RA. Two of the genes in our study, CDK6 and PTPN22, are directly involved in signaling pathways that are essential for immune cell function. However, due to redundancy of these pathways we cannot point to the specific type of cells that are involved in the development of disease: it could be downstream mechanisms with T helper cells as primer players, or upstream mechanisms related to antigen presenting cells, B cells, macrophages/monocytes or dendritic cells. As the combination of HLA-DRB1*04 and variations from these genes are unlikely to be sufficient for RA development, one would expect a third determinant to be involved in this interaction. Our findings warrant future studies of the development of autoimmune responses in individuals at risk of RA, with certain subphenotypes, based on serological testing. It would guide both understanding of specific mechanisms and probable treatment towards precision medicine in contrast to current treatment based on relatively non-specific immunosuppressive therapy (that is, RA symptoms). Such studies may guide both understanding of specific mechanisms and future development of more specific therapies for patients from these different subsets. One such option that may be discussed based on our current observations in ACPA-negative RA patients is to test the selective cyclin-dependent kinase 4/6 (CDK4/6) inhibitor (Momilactone B) which is currently being investigated for treatment of colon cancer .
There are several strengths in our study. First, we performed analysis in three independent cohorts of RA patients. Second, we used cohorts from two different countries, with different frequencies of autoantibodies, and different types of RA. This variability between sample collections is important for generalization of our findings. Third, our final conclusions were based on the Bonferroni correction and this is a very conservative approach when analyzing correlated variables, such as the four autoantibodies in our study.
There are also some weaknesses in the current study: 1) as the phenotypes we have studied involved a combination of traits, they comprised relatively small numbers of observations even with several thousands of patients included in the study; 2) inclusion of incident cases (EIRA) and established RA (cohorts 2 and 3) may have decreased the power of our study to discover associations that are specific for a particular disease phase; 3) we did not address the contribution of other SE alleles beyond DRB1*04 in our study due to a relatively low frequency of HLA-DRB1*01 and *10 in our cohorts; 4) it also would be beneficial to address the influence of smoking, but we did not have enough statistical power for inclusion of an additional covariate in our analyses.
Autoimmune diseases display strong association with HLA alleles, but a substantial number of non-HLA variants are also associated at genome-wide level. Specific subphenotype analyses of these associated loci have been rare, and have usually not been stratified by HLA alleles. Our study provides a prototype that could be useful for other studies aiming to disentangle RA complexity. The ongoing replications and extensions of such studies should provide a basis for a molecular re-classification of RA based on genetic and serological features, rather than mainly on clinical data . Specifically, our study confirms the hypothesis that polymorphisms outside the HLA-DRB1 locus contribute to the development of specific serological subphenotypes of RA, which may have different mechanisms of disease development.
ACPA: anti-citrullinated peptide antibody
anti-CCP: anti-citrullinated cyclic peptide
CEP-1: citrullinated alpha-enolase peptide-1
citC1III: citrullinated type-II collagen trimer
cit-Vim: citrullinated vimentin
CMH test: Cochran-Mantel-Haenszel test
DAS28: disease activity score based on 28 joints
EIRA: epidemiological investigation of rheumatoid arthritis
ELISA: enzyme-linked immunosorbent assay
GWAS: genome-wide association study
HLA: human leucocyte antigen (locus)
LYP: lymphocyte phosphatase
OR: odds ratio
PCR: polymerase chain reaction
RA: rheumatoid arthritis
SE: shared epitope
SNP: single nucleotide polymorphism. Gene names are not listed and follow current genetic nomenclature
The authors declare that they have no competing interests.
OS, DGC, AG, VM and LP made substantial contributions to conception and design, acquisition of data, analysis and interpretation of data; AM, EPP, JJGR, MS, KUK, LI and AIC made substantial contributions in acquisition of data; BD made substantial contributions to analysis and interpretation of data; RH, LA, LK and JT made substantial contributions to conception and design of the study; OS, DGC, AG, VM and LP have been involved in drafting the manuscript and revising it critically for important intellectual content; AM, EPP, JJGR, MS, KUK, LI, AIC, RH, LA, LK and JT have been involved in revising the manuscript critically for important intellectual content. All authors gave final approval of the version to be published. LP agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
Current affiliation of Omri Snir is Department of Immunology, Centre for Immune Regulation, Oslo University Hospital-Rikshospitalet, University of Oslo, Oslo, Norway. Current affiliation of Katharina Ute Klich is Department of Cell and Developmental Biology, Max Planck Institute for Molecular Biomedicine, Münster, Germany. Current affiliation of Bo Ding is Observational Research Center, AstraZeneca, Mölndal, Sweden.
We wish to acknowledge financial support from EU-funded projects BTCure and Masterswitch, and support from the Swedish innovation agency Vinnova (COMBINE project). LP was also supported by The Swedish Research Council (Vetenskapsrådet) and King Gustaf V 80-year fund. DGC participation was supported by the Bioinformatic Infrastructure for Life Sciences (BILS). Work in the Spanish group was supported by grant PI11/01048 of the Instituto de Salud Carlos III (Spain) that is partially financed by the European Regional Development Fund of the EU. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Klareskog L, Stolt P, Lundberg K, Kallberg H, Bengtsson C, Grunewald J, Ronnelid J, Harris HE, Ulfgren AK, Rantapaa-Dahlqvist S, Eklund A, Padyukov L, Alfredsson L: A new model for an etiology of rheumatoid arthritis: smoking may trigger HLA-DR (shared epitope)-restricted immune reactions to autoantigens modified by citrullination.
Huizinga TW, Amos CI, van der Helm-van Mil AH, Chen W, van Gaalen FA, Jawaheer D, Schreuder GM, Wener M, Breedveld FC, Ahmad N, Lum RF, de Vries RR, Gregersen PK, Toes RE, Criswell LA: Refining the complex rheumatoid arthritis phenotype based on specificity of the HLA-DRB1 shared epitope for antibodies to citrullinated proteins.
Snir O, Widhe M, Hermansson M, von Spee C, Lindberg J, Hensen S, Lundberg K, Engstrom A, Venables PJ, Toes RE, Holmdahl R, Klareskog L, Malmstrom V: Antibodies to several citrullinated antigens are enriched in the joints of rheumatoid arthritis patients.
Snir O, Widhe M, von Spee C, Lindberg J, Padyukov L, Lundberg K, Engstrom A, Venables PJ, Lundeberg J, Holmdahl R, Klareskog L, Malmstrom V: Multiple antibody reactivities to citrullinated antigens in sera from patients with rheumatoid arthritis: association with HLA-DRB1 alleles.
Raychaudhuri S, Sandor C, Stahl EA, Freudenberg J, Lee HS, Jia X, Alfredsson L, Padyukov L, Klareskog L, Worthington J, Siminovitch KA, Bae SC, Plenge RM, Gregersen PK, de Bakker PI: Five amino acids in three HLA proteins explain most of the association between MHC and seropositive rheumatoid arthritis.
van der Helm-van Mil AH, Verpoort KN, Breedveld FC, Huizinga TW, Toes RE, de Vries RR: The HLA-DRB1 shared epitope alleles are primarily a risk factor for anti-cyclic citrullinated peptide antibodies and are not an independent risk factor for development of rheumatoid arthritis.
Arthritis and Rheum 2006, 54:1117-1121. Publisher Full Text
Verpoort KN, Cheung K, Ioan-Facsinay A, van der Helm-van Mil AH, de Vries-Bouwstra JK, Allaart CF, Drijfhout JW, de Vries RR, Breedveld FC, Huizinga TW, Pruijn GJ, Toes RE: Fine specificity of the anti-citrullinated protein antibody response is influenced by the shared epitope alleles.
Mahdi H, Fisher BA, Kallberg H, Plant D, Malmstrom V, Ronnelid J, Charles P, Ding B, Alfredsson L, Padyukov L, Symmons DP, Venables PJ, Klareskog L, Lundberg K: Specific interaction between genotype, smoking and autoimmunity to citrullinated alpha-enolase in the etiology of rheumatoid arthritis.
van der Woude D, Alemayehu WG, Verduijn W, de Vries RR, Houwing-Duistermaat JJ, Huizinga TW, Toes RE: Gene-environment interaction influences the reactivity of autoantibodies to citrullinated antigens in rheumatoid arthritis.
Nat Genet 2010, 42:814-816.
Author reply: 816PubMed Abstract | Publisher Full Text
Lundberg K, Bengtsson C, Kharlamova N, Reed E, Jiang X, Kallberg H, Pollak-Dorocic I, Israelsson L, Kessel C, Padyukov L, Holmdahl R, Alfredsson L, Klareskog L: Genetic and environmental determinants for disease risk in subsets of rheumatoid arthritis defined by the anticitrullinated protein/peptide antibody fine specificity profile.
Ronnelid J, Wick MC, Lampa J, Lindblad S, Nordmark B, Klareskog L, van Vollenhoven RF: Longitudinal analysis of citrullinated protein/peptide antibodies (anti-CP) during 5 year follow up in early rheumatoid arthritis: anti-CP status predicts worse disease activity and greater radiological progression.
Stolt P, Bengtsson C, Nordmark B, Lindblad S, Lundberg I, Klareskog L, Alfredsson L: Quantification of the influence of cigarette smoking on rheumatoid arthritis: results from a population based case–control study, using incident cases.
Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, Cooper NS, Healey LA, Kaplan SR, Liang MH, Luthra HS, Medsger TA, Mitchell DM, Neustadt DH, Pinals RS, Schaller JG, Sharp JT, Wilder RL, Hunder GG: The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis.
Montes A, Dieguez-Gonzalez R, Perez-Pampin E, Calaza M, Mera-Varela A, Gomez-Reino JJ, Gonzalez A: Particular association of clinical and genetic features with autoimmunity to citrullinated alpha-enolase in rheumatoid arthritis.
Olerup O, Zetterquist H: HLA-DR typing by PCR amplification with sequence-specific primers (PCR-SSP) in 2 hours: an alternative to serological DR typing in clinical practice including donor-recipient matching in cadaveric transplantation.
Lundstrom E, Kallberg H, Alfredsson L, Klareskog L, Padyukov L: 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.
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.
Padyukov L, Seielstad M, Ong RT, Ding B, Ronnelid J, Seddighzadeh M, Alfredsson L, Klareskog L: A genome-wide association study suggests contrasting associations in ACPA-positive versus ACPA-negative rheumatoid arthritis.
Ding B, Padyukov L, Lundstrom E, Seielstad M, Plenge RM, Oksenberg JR, Gregersen PK, Alfredsson L, Klareskog L: Different patterns of associations with anti-citrullinated protein antibody-positive and anti-citrullinated protein antibody-negative rheumatoid arthritis in the extended major histocompatibility complex region.
van der Woude D, Lie BA, Lundstrom E, Balsa A, Feitsma AL, Houwing-Duistermaat JJ, Verduijn W, Nordang GB, Alfredsson L, Klareskog L, Pascual-Salcedo D, Gonzalez-Gay MA, Lopez-Nevot MA, Valero F, Roep BO, Huizinga TW, Kvien TK, Martin J, Padyukov L, de Vries RR, Toes RE: Protection against anti-citrullinated protein antibody-positive rheumatoid arthritis is predominantly associated with HLA-DRB1*1301: a meta-analysis of HLA-DRB1 associations with anti-citrullinated protein antibody-positive and anti-citrullinated protein antibody-negative rheumatoid arthritis in four European populations.
Feitsma AL, van der Voort EI, Franken KL, el Bannoudi H, Elferink BG, Drijfhout JW, Huizinga TW, de Vries RR, Toes RE, Ioan-Facsinay A: Identification of citrullinated vimentin peptides as T cell epitopes in HLA-DR4-positive patients with rheumatoid arthritis.
Snir O, Rieck M, Gebe JA, Yue BB, Rawlings CA, Nepom G, Malmström V, JH B: Identification and functional characterization of T cells reactive to citrullinated-vimentin in HLA-DRB1*0401 humanized mice and RA patients.
James E, Rieck M, Pieper J, Gebe JA, Yue BB, Tatum M, Peda M, Sandin C, Klareskog L, Malmstrom V, Buckner JH: Citrulline specific Th1 cells are increased in rheumatoid arthritis and their frequency is influenced by disease duration and therapy.
Arthritis Rheum 2014, 66:1712-1722. Publisher Full Text
Scally SW, Petersen J, Law SC, Dudek NL, Nel HJ, Loh KL, Wijeyewickrema LC, Eckle SB, van Heemst J, Pike RN, McCluskey J, Toes RE, La Gruta NL, Purcell AW, Reid HH, Thomas R, Rossjohn J: A molecular basis for the association of the HLA-DRB1 locus, citrullination, and rheumatoid arthritis.
Begovich AB, Carlton VE, Honigberg LA, Schrodi SJ, Chokkalingam AP, Alexander HC, Ardlie KG, Huang Q, Smith AM, Spoerke JM, Conn MT, Chang M, Chang SY, Saiki RK, Catanese JJ, Leong DU, Garcia VE, McAllister LB, Jeffery DA, Lee AT, Batliwalla F, Remmers E, Criswell LA, Seldin MF, Kastner DL, Amos CI, Sninsky JJ, Gregersen PK: A missense single-nucleotide polymorphism in a gene encoding a protein tyrosine phosphatase (PTPN22) is associated with rheumatoid arthritis.
Menard L, Saadoun D, Isnardi I, Ng YS, Meyers G, Massad C, Price C, Abraham C, Motaghedi R, Buckner JH, Gregersen PK, Meffre E: The PTPN22 allele encoding an R620W variant interferes with the removal of developing autoreactive B cells in humans.
Raychaudhuri S, Remmers EF, Lee AT, Hackett R, Guiducci C, Burtt NP, Gianniny L, Korman BD, Padyukov L, Kurreeman FA, Chang M, Catanese JJ, Ding B, Wong S, van der Helm-van Mil AH, Neale BM, Coblyn J, Cui J, Tak PP, Wolbink GJ, Crusius JB, van der Horst-Bruinsma IE, Criswell LA, Amos CI, Seldin MF, Kastner DL, Ardlie KG, Alfredsson L, Costenbader KH, Altshuler D, et al.: Common variants at CD40 and other loci confer risk of rheumatoid arthritis.
van der Linden MP, Feitsma AL, le Cessie S, Kern M, Olsson LM, Raychaudhuri S, Begovich AB, Chang M, Catanese JJ, Kurreeman FA, van Nies J, van der Heijde DM, Gregersen PK, Huizinga TW, Toes RE, van der Helm-Van Mil AH: Association of a single-nucleotide polymorphism in CD40 with the rate of joint destruction in rheumatoid arthritis.
Weedon MN, Lango H, Lindgren CM, Wallace C, Evans DM, Mangino M, Freathy RM, Perry JR, Stevens S, Hall AS, Samani NJ, Shields B, Prokopenko I, Farrall M, Dominiczak A, Diabetes Genetics I, Wellcome Trust Case Control C, Johnson T, Bergmann S, Beckmann JS, Vollenweider P, Waterworth DM, Mooser V, Palmer CN, Morris AD, Ouwehand WH, Cambridge GEMC, Zhao JH, Li S, Loos RJ, et al.: Genome-wide association analysis identifies 20 loci that influence adult height.
Okada Y, Hirota T, Kamatani Y, Takahashi A, Ohmiya H, Kumasaka N, Higasa K, Yamaguchi-Kabata Y, Hosono N, Nalls MA, Chen MH, van Rooij FJ, Smith AV, Tanaka T, Couper DJ, Zakai NA, Ferrucci L, Longo DL, Hernandez DG, Witteman JC, Harris TB, O'Donnell CJ, Ganesh SK, Matsuda K, Tsunoda T, Tanaka T, Kubo M, Nakamura Y, Tamari M, Yamamoto K, et al.: Identification of nine novel loci associated with white blood cell subtypes in a Japanese population.
Reiner AP, Lettre G, Nalls MA, Ganesh SK, Mathias R, Austin MA, Dean E, Arepalli S, Britton A, Chen Z, Couper D, Curb JD, Eaton CB, Fornage M, Grant SF, Harris TB, Hernandez D, Kamatini N, Keating BJ, Kubo M, LaCroix A, Lange LA, Liu S, Lohman K, Meng Y, Mohler ER 3rd, Musani S, Nakamura Y, O'Donnell CJ, Okada Y, et al.: Genome-wide association study of white blood cell count in 16,388 African Americans: the continental origins and genetic epidemiology network (COGENT).
Wang Y, Wysocka J, Sayegh J, Lee YH, Perlin JR, Leonelli L, Sonbuchner LS, McDonald CH, Cook RG, Dou Y, Roeder RG, Clarke S, Stallcup MR, Allis CD, Coonrod SA: Human PAD4 regulates histone arginine methylation levels via demethylimination.
Suzuki A, Yamada R, Chang X, Tokuhiro S, Sawada T, Suzuki M, Nagasaki M, Nakayama-Hamada M, Kawaida R, Ono M, Ohtsuki M, Furukawa H, Yoshino S, Yukioka M, Tohma S, Matsubara T, Wakitani S, Teshima R, Nishioka Y, Sekine A, Iida A, Takahashi A, Tsunoda T, Nakamura Y, Yamamoto K: Functional haplotypes of PADI4, encoding citrullinating enzyme peptidylarginine deiminase 4, are associated with rheumatoid arthritis.
Eyre S, Bowes J, Diogo D, Lee A, Barton A, Martin P, Zhernakova A, Stahl E, Viatte S, McAllister K, Amos CI, Padyukov L, Toes RE, Huizinga TW, Wijmenga C, Trynka G, Franke L, Westra HJ, Alfredsson L, Hu X, Sandor C, de Bakker PI, Davila S, Khor CC, Heng KK, Andrews R, Edkins S, Hunt SE, Langford C, Symmons D, et al.: High-density genetic mapping identifies new susceptibility loci for rheumatoid arthritis.
Too CL, Murad S, Dhaliwal JS, Larsson P, Jiang X, Ding B, Alfredsson L, Klareskog L, Padyukov L: Polymorphisms in peptidylarginine deiminase associate with rheumatoid arthritis in diverse Asian populations: evidence from MyEIRA study and meta-analysis.
Park C, Jeong NY, Kim GY, Han MH, Chung IM, Kim WJ, Yoo YH, Choi YH: Momilactone B induces apoptosis and G1 arrest of the cell cycle in human monocytic leukemia U937 cells through downregulation of pRB phosphorylation and induction of the cyclin-dependent kinase inhibitor p21Waf1/Cip1.