Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Genetic Risk Score Predicting Risk of Rheumatoid Arthritis Phenotypes and Age of Symptom Onset

  • Lori B. Chibnik ,

    lchibnik@rics.bwh.harvard.edu

    Affiliation Program in Translational NeuroPsychiatric Genomics, Department of Neurology, Brigham and Women's Hospital, Boston, Massachusetts, United States of America

    ⨯
  • Brendan T. Keenan,

    Affiliation Program in Translational NeuroPsychiatric Genomics, Department of Neurology, Brigham and Women's Hospital, Boston, Massachusetts, United States of America

    ⨯
  • Jing Cui,

    Affiliation Section of Clinical Sciences, Division of Rheumatology, Immunology and Allergy, Department of Medicine, Brigham and Women's Hospital, Boston, Massachusetts, United States of America

    ⨯
  • Katherine P. Liao,

    Affiliation Section of Clinical Sciences, Division of Rheumatology, Immunology and Allergy, Department of Medicine, Brigham and Women's Hospital, Boston, Massachusetts, United States of America

    ⨯
  • Karen H. Costenbader,

    Affiliation Section of Clinical Sciences, Division of Rheumatology, Immunology and Allergy, Department of Medicine, Brigham and Women's Hospital, Boston, Massachusetts, United States of America

    ⨯
  • Robert M. Plenge,

    Affiliation Division of Rheumatology, Immunology and Allergy, Department of Medicine and Division of Genetics, Brigham and Women's Hospital, Boston, Massachusetts, United States of America

    ⨯
  • Elizabeth W. Karlson

    Affiliation Section of Clinical Sciences, Division of Rheumatology, Immunology and Allergy, Department of Medicine, Brigham and Women's Hospital, Boston, Massachusetts, United States of America

    ⨯

Abstract

Background

Cumulative genetic profiles can help identify individuals at high-risk for developing RA. We examined the impact of 39 validated genetic risk alleles on the risk of RA phenotypes characterized by serologic and erosive status.

Methods/Principal Findings

We evaluated single nucleotide polymorphisms at 31 validated RA risk loci and 8 Human Leukocyte Antigen alleles among 542 Caucasian RA cases and 551 Caucasian controls from Nurses' Health Study and Nurses' Health Study II. We created a weighted genetic risk score (GRS) and evaluated it as 7 ordinal groups using logistic regression (adjusting for age and smoking) to assess the relationship between GRS group and odds of developing seronegative (RFβˆ’ and CCPβˆ’), seropositive (RF+ or CCP+), erosive, and seropositive, erosive RA phenotypes. In separate case only analyses, we assessed the relationships between GRS and age of symptom onset.

In 542 RA cases, 317 (58%) were seropositive, 163 (30%) had erosions and 105 (19%) were seropositive with erosions. Comparing the highest GRS risk group to the median group, we found an OR of 1.2 (95% CIβ€Š=β€Š0.8–2.1) for seronegative RA, 3.0 (95% CIβ€Š=β€Š1.9–4.7) for seropositive RA, 3.2 (95% CIβ€Š=β€Š1.8–5.6) for erosive RA, and 7.6 (95% CIβ€Š=β€Š3.6–16.3) for seropositive, erosive RA. No significant relationship was seen between GRS and age of onset.

Conclusions/Significance

Results suggest that seronegative and seropositive/erosive RA have different genetic architecture and support the importance of considering RA phenotypes in RA genetic studies.

Citation: Chibnik LB, Keenan BT, Cui J, Liao KP, Costenbader KH, Plenge RM, et al. (2011) Genetic Risk Score Predicting Risk of Rheumatoid Arthritis Phenotypes and Age of Symptom Onset. PLoS ONE 6(9): e24380. https://doi.org/10.1371/journal.pone.0024380

Editor: Sunil K. Ahuja, South Texas Veterans Health Care System, United States of America

Received: April 5, 2011; Accepted: August 9, 2011; Published: September 12, 2011

Copyright: Β© 2011 Chibnik et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Funding: This work was supported by the National Institutes of Health (grants CA87969, CA49449, CA67262, CA50385, R049880, AR47782, AR0524). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Rheumatoid arthritis (RA) is a chronic systemic inflammatory disease that often leads to disability from joint damage and inflammation. Although RA is uncommon, with a worldwide prevalence of approximately 1%, it has a large economic and societal cost, primarily in terms of work-related disability [1]. Bony destruction, or erosions, are associated with work disability [2] and lower functional status [3], [4] and thus with a more debilitating disease phenotype. Before the more widespread use of biologics, the incidence of erosions in RA patients was ∼70% within the first 3 years after diagnosis [5]. However, with more advanced treatment the prevalence of erosions has decreased [6], [7].

Recently in genetic studies, RA has been divided into two phenotypes defined by presence or absence of serologic factors. Originally the subdivision was based on rheumatoid factor (RF), but more recently antibodies to cyclic citrullinated proteins (anti-CCP) have been used to define the two subtypes [8], [9], [10], [11]. Both RF and anti-CCP positivity have been linked to more severe disease features and outcomes [12], [13], [14], [15]. Presence of RF has also been shown to be a major predictor of both development and severity of joint erosions [16], [17], [18], [19]. Based on these previous results, we defined RA on a continuum of disease severity ranging from seronegative RA (least severe phenotype) to seropositive RA or erosive RA (more severe phenotypes) and finally to seropositive, erosive RA (most severe phenotype).

In addition to erosions and serologic status, age at onset of disease has been associated with RA outcomes, although the results have been varied, with older disease onset predicting worse outcomes in some [20], [21], [22], [23] and milder outcomes in other studies [22], [24]. Specifically, Bukhari and colleagues showed that those with older age of RA onset had higher odds of developing erosive disease and a worse severity of erosions as compared to those with onset <50 years old [25]. Moreover, Pease et.al. reported a slight increase in odds of erosions (although non-significant) for those with onset 65 years or older, but in contrast also found a 3 fold increase in odds of RA disease remission for the same age group [22]. Earlier studies have also shown a lower prevalence of HLA-DR4, the major genetic risk factor for RA, in patients with later onset of RA, although not always significant [26], [27], [28], [29]. Similar to serologic status, these and other studies have suggested that RA could potentially be divided into 2 subsets defined by earlier versus later age at onset [22], [26], [27], [28], [30].

Methods

Objectives

Karlson et al. showed that a weighted Genetic Risk Score (GRS) with 22 RA risk alleles showed a good discrimination between seropositive RA and controls. The addition of a weighted GRS score comprised of validated genetic risk factors showed improved discrimination when compared to a model with just clinical risk factors alone [31]. We extend this analysis in two ways, first, by adding the newly validated RA risk alleles to the GRS [32], [33], [34], [35] and second, by assessing the GRS in relation to the more specific phenotypes of RA along the severity continuum, including erosive status, seropositivity and age at first symptom onset of RA. We will show that the GRS is most applicable for the more extreme RA phenotypes defined by seropositive and erosive status, and consequently that these phenotypes have a different genetic architecture than the seronegative and non-erosive forms of the disease.

Participants

The Nurses' Health Study is a prospective cohort which enrolled 121,700 female nurses aged 30 to 55 years throughout the US in 1976. Of those, 32,826 (27%) participants provided blood samples for future studies and an additional 33,040 (27%) provided buccal cell samples for a total of 65,866 (54%) samples with available DNA. A similar prospective cohort, Nurses' Health Study II enrolled 116,609 female nurses aged 25 to 42 years in 1989, of which 29,611 (25%) provided blood samples for future studies. For these analyses the two cohorts will be combined and referred to as β€˜NHS’.

Ethics

All aspects of this study were approved by the Partners Human Research Committee, the Institutional Review Board of Partners Research Management. Three types of written informed consent were acquired for these studies. First, for questionnaires, NHS cohort participants were consented by paper at baseline in 1976 for repeated surveys and NHSII cohort participants were consented by paper at baseline in 1989 for repeated surveys. Second, for the sub-cohorts who contributed blood for the genetic analysis, the NHS participants were consented by paper in 1989 at blood draw and the NHS2 participants were consented in 1997 at blood draw. Finally, all RA self-reported cases signed informed consent to release medical records for review.

Phenotypic Rheumatoid Arthritis

A staged screening method was used to confirm cases of RA in the NHS cohort. A connective tissue disease (CTD) screening questionnaire was used to screen all self-reported cases for RA symptoms, followed by chart validation by two board-certified rheumatologists [36]. Four phenotypes of RA were defined using rheumatoid factor (RF) and/or CCP positivity and presence of radiographic changes/erosions. Rheumatoid factor was determined by chart review. Second generation CCP assays were performed among available pre-diagnosis or post-diagnosis blood samples for a subset of cases (nβ€Š=β€Š273 cases) as previously described [37], all others were obtained by chart reviews since the mid-2000's when the test became widely available. Erosions were determined by chart review [37]. The four RA phenotypes of primary interest included: 1) seronegative RA (both RF and CCP negative, nβ€Š=β€Š225, in supplemental); 2) seropositive RA (either RF+ or CCP+ or both, nβ€Š=β€Š317); 3) Erosive RA (presence of erosions, nβ€Š=β€Š163); and 4) seropositive, erosive RA (presence of erosions AND in group 2, nβ€Š=β€Š105). In addition, we also examined all RA (nβ€Š=β€Š542) and non-erosive RA (nβ€Š=β€Š379) as a supplemental analysis. The groups are not mutually exclusive, i.e. cases in the seropositive, erosive RA group are also included in the seropositive RA and the erosive RA groups. Age at RA onset was determined from chart reviews as either age at onset of RA symptoms, if available, or age at RA diagnosis.

Genotyping

Low resolution HLA-DRB1 genotyping was performed using PCR with sequence-specific primers (SSP) using OLERUP SSP kits (Qiagen, West Chester, Pennsylvania, USA), as previously described [38]. For samples with positive two-digit human leucocyte antigen (HLA) signals, SSP were used for high- resolution four-digit allele detection of DRB1*0401, *0404, *0405, *0408, *0101, *0102, *09 and *1001. All non-MHC risk alleles were genotyped using iPlex (Sequenom, San Diego, California, USA) at the Broad Institute, as previously described [33]. All SNP had call rates greater than 95% and Hardy–Weinberg equilibrium p values greater than 0.01.

Creation of Genetic Risk Score (GRS)

Thirty-nine validated risk alleles for RA are combined to form a continuous GRS. This GRS is a weighted combination of 8 HLA-DRB1 β€˜shared epitope’ (HLA-SE) alleles and 31 non-MHC risk alleles. Each allele is weighted by the natural log of the published OR and then summed over 39 alleles, as previously described [31], [39], [40]. Using a weighted risk score is important in this case since both PTPN22 and the HLA-SE have a stronger predictive relationship with RA as compared to the more recently discovered SNPs. The list of alleles and weights are presented in Supplemental Table 1. We evaluated the linkage disequilibrium (LD) structure of the risk alleles using HapMap release 22 and found little evidence of LD (largest R2β€Š=β€Š0.06) suggesting that LD has little to no affect on the variance of the GRS.

The ORs and weights for the HLA-SE alleles are from a meta-analysis of all published studies [41]. Odds ratios and weights for 5 out of the 31 non-MHC SNPs were taken from extensively replicated SNPs from published studies. These include PTPN22 (rs2476601) [42], TRAF1-C5 (rs3761847) [43], STAT4 (rs7574865) [44], TNFAIP3 (rs17066662, in LD with rs1099194, r2β€Š=β€Š1.0) [45] and TNFAIP3 (rs6920220) [45]. Odds ratios and weights for 9 of the 31 non-MHC alleles were taken from a meta-analysis of GWAS data for 3,393 cases and 12,462 controls with replication in 3,929 seropositive RA cases and 5,807 matched controls by Raychaudhuri et al, [33]. To avoid over-estimation of the true effect size we used the ORs from the replication phase of the study [32]. These SNPs include CD40 (rs4810485), CCL21 (rs2812378), CTLA4 (rs3087243), PADI4 (rs2240340), CDK6 (rs42041), TNFRSF14 (rs3890745), PRKCQ (rs4750316), KIF5A (rs1678542), and 4q27 (IL2/IL21) (rs6822844) [33]. The ORs and weights for 7 risk alleles were selected from the joint analysis from Raychaudhuri, et al [34] and were identified as functionally related to known RA risk loci by GRAIL, a bioinformatics analysis that identifies connections among genes in published abstracts. These include PTPRC (rs10919563), CD2 (rs11586238), CD28 (rs1980422), TAGAP (rs394581), RAG1 (rs540386), PRDM1 (rs548234), and FCGR2A (rs7552317) [34]. Finally, 10 SNPs were selected from the final combined analysis from a genome-wide association study (GWAS) meta-analysis of 5,539 autoantibody positive RA cases and 20,169 controls of European descent, followed by replication in an independent set of 6,768 RA cases and 8,806 controls [35]. These include SPRED2 (rs934734), ANKRD55/IL6ST (rs6859219), C5orf13/GIN1 (rs26232), PXK (rs13315591), RBPJ (rs874040), CCR6 (rs3093023), IRF5 (rs10488631), AAF3 (rs11676922), CCL21 (rs951005) and IL2RA (rs706778).

The continuous GRS score was then divided into 7 sub-groups. The thresholds for the groups were based on the Gaussian distribution in the controls. A more detailed description of methods is published elsewhere [31]. Briefly, dividing our score into 7 categories provided the most robust distribution, allowing us to parse out the highest and lowest risk groups while ensuring that there were sufficient numbers of cases and controls in these extreme categories of interest.

Study Sample Filtering

Each confirmed RA case was matched to one healthy control by cohort (NHS or NHSII), year of birth, menopausal status, and post-menopausal hormone use. Our initial nested case-control dataset consisted of 585 RA cases and 585 matched controls. To reduce the potential for population stratification we limited our analysis to self-reported Caucasian women, resulting in 564 RA cases and 571 controls. Since the HLA alleles have a large weight in the GRS we dropped any participant missing HLA data. Among the 564 RA cases, 22 (4%) were missing HLA, and among the 571 healthy controls, 20 (4%) were missing HLA. This left us with 542 RA cases and 551 healthy controls. For anyone missing other SNPs, we assigned them a value equal to the expected value (2*risk allele frequency defined in cases or controls separately).

Epidemiological Covariates

Smoking is the strongest environmental factor linked with RA, and its population attributable risk is 25% for all RA and 35% for seropositive RA [46], [47], [48]. Prospective, biennial questionnaires were used to collect covariate information from all NHS subjects. The questions include inquiries regarding diseases, lifestyle and health practices. Lifetime history of smoking was collected at the baseline questionnaire and data concerning current smoking status and number of cigarettes smoked per day were updated in each two year questionnaire cycle. Pack-years were calculated as number of packs per day smoked times number of years of smoking using the questionnaire cycle prior to the date of RA diagnosis or index date for matched controls.

Statistical Analysis

Demographic characteristics of the cohorts are described using means and standard deviations for continuous variables and frequency and proportions for categorical variables. Logistic regression analysis was used to calculate the odds of a phenotypic RA for each GRS risk group as compared to the median group (group 4). The odds of phenotypic RA for the most extreme risk group (group 7) as compared to the least extreme risk group (group 1) was calculated using an ordinal model that takes into account all the data in all the groups. A test for linear trend across all seven groups was performed using logistic regression, with each group equaling the median GRS level in that group. The discriminatory ability of the GRS to define case group vs. control group at different combinations of sensitivity and specificity was assessed using a Receiver Operating Characteristic (ROC) curve and computing the Area Under Curve (AUC). Finally, Pearson correlation coefficients were used to compare continuous GRS and age at RA symptom onset and Analysis of Variance (ANOVA) was used to calculate the mean age at RA symptom onset for each GRS risk group. Models were adjusted for year of birth and pack-years of smoking. All analyses were performed on SAS Version 9.1 (SAS Institute, Cary, NC).

Results

Subjects

Five hundred and forty two RA cases were identified with a mean age at RA symptom onset of 56 (SD, 11). Of these, 317 (58%) were seropositive, 163 (30%) had evidence of erosions and 105 (19%) had seropositive, erosive RA. Five hundred and fifty-one controls were selected among NHS participants who gave a blood or buccal cell sample. The mean age at time of blood sample was 55 (SD, 8) years for cases and 56 (SD, 8) years for controls. Demographic information for cases and controls are presented in Table 1.

thumbnail
Download:
Table 1. Characteristics of RA cases and controls in the Nurses' Health Study.

https://doi.org/10.1371/journal.pone.0024380.t001

Relationship between Risk Alleles and Seropositive RA

Bivariate associations between single risk alleles and odds of seropositive RA in the NHS data are presented in Supplemental Table S1. In most cases the direction of the association in the NHS data and the published ORs is the same, although most confidence intervals cross the null value of 1.0. This is to be expected since the individual effect sizes are small, and thus we do not have the power to see significant bivariate associations.

GRS and Odds of Phenotypic RA

The results of the association analysis of three of the outcomes, seropositive RA, erosive RA and seropositive, erosive RA are presented in Table 2. The additional outcomes, all RA, seronegative RA and non-erosive RA are presented in Supplemental Table S1. Group 4 (the median level of risk) was used as the referent group in this analysis. Those with a GRS in group 7 had a significantly increased odds of seropositive RA (ORβ€Š=β€Š3.0; 95%CI 1.9–4.7), erosive RA (ORβ€Š=β€Š3.2; 95%CI 1.8–5.6) and seropositive/erosive RA (ORβ€Š=β€Š7.6; 95%CI 3.6–16.3), with the highest increased odds being for seropositive/erosive RA. Those in top GRS group (group 7) had no significant increase in odds of seronegative RA (ORβ€Š=β€Š1.2; 95%CI 0.8–2.1) (Supplementary Table S2). The discrimination ability of the model, measured as AUC, for predicting seronegative RA was 0.563, only slightly above the null value of 0.50 (a null model would have an AUC of 0.500, whereas a perfect model would have an AUC of 1.0). The other 3 phenotypes of RA showed better discrimination with AUCs of 0.654, 0.644 and 0.712 for seropositive RA, erosive RA and seropositive, erosive RA respectively. The ROC curves for 4 outcomes, seronegative, seropositive, erosive and seropositive, erosive RA are represented in Figure 1. Since the outcomes vary across the models we cannot directly compare the AUCs using the known methods [49].

thumbnail
Download:
Figure 1. Receiver Operating Characteristic (ROC) curve for 4 phenotypes of RA.

https://doi.org/10.1371/journal.pone.0024380.g001

thumbnail
Download:
Table 2. Weighted GRS groups and odd ratios of seronegative, Seropositive, Erosive and Seropositive/Erosive RA in NHS.

https://doi.org/10.1371/journal.pone.0024380.t002

In the ordinal model, which takes into account all the data in all the groups, we see a significant increase odds of seropositive RA (ORβ€Š=β€Š5.7; 95%CI 3.5–9.3), erosive RA (ORβ€Š=β€Š4.9; 95%CI 3.2–10.8) and seropositive/erosive RA (ORβ€Š=β€Š14.2; 95%CI 6.5–30.9) for group 7 as compared to group 1 (Table 2). In addition, we see an increased odds of seronegative RA (ORβ€Š=β€Š2.0; 95%CI 1.2–3.5) for those with a GRS in the top group (group 7) compared to the lowest risk group (group 1) (Supplementary Table S2). Strongly significant linear trends were seen in the seropositive, erosive and seropositive/erosive RA case groups, with all p-values<0.0001 (Table 2). For seronegative RA the p for trend was 0.007 (Supplementary Table S2).

Association between GRS and Age at RA Symptom Onset

The results for the association between the GRS and age at RA symptom onset are presented in Table 3 and Supplementary Table S3. The adjusted mean ages at RA symptom onset were not significantly different among the 7 GRS groups for any of the phenotype subgroups (p>0.05 for all). The correlations between continuous GRS for seropositive RA, erosive RA and seropositive, erosive RA all were negative (βˆ’0.09, βˆ’0.08, βˆ’0.11 respectively) indicating that the larger the GRS, the younger the age at first symptom; however none of these correlations were significant.

thumbnail
Download:
Table 3. Relationship between weighted GRS as groups and as continuous and age at RA symptom start.

https://doi.org/10.1371/journal.pone.0024380.t003

Discussion

We found that a weighted genetic risk score was associated with development of seropositive RA, erosive RA and seropositive, erosive RA phenotypes. Although there was a significant linear trend with a continuous GRS39 measure predicting seronegative RA, with the exception of group 7 compared to group 1, there was no significant relationship when the score was divided into groups. In contrast, we found a strong and significant association between both continuous and grouped GRS39 and the erosive and/or seropositive phenotypes. Subjects with the highest GRS score (group 7) had a 3.2 times increase of odds of erosive RA as compared to the median group. This odds ratio increased to 7.6 when limiting the phenotype to those with seropositive, erosive RA. We observed similar results when comparing extreme GRS scores (group 7 vs. group 1), where we found a 5 times increased odds of erosive RA and a 14 times increased odds for seropositive, erosive RA. This suggests that the GRS has a stronger association with the more severe phenotype; however narrowing the phenotype definition resulted in a widened confidence interval. Thus, although we detected a stronger effect size (i.e. larger OR), there was also greater variability in the association, most likely due to the small sample size in this group.

One interesting result is the association between the GRS with 39 risk alleles and seropositive RA. We found that group 7 had a 3.0 times increased odds of seropositive RA as compared to group 4. This is similar to the 2.9 times increased odds found by Karlson et al [31] with the GRS based on 22 risk alleles. In addition, we observed a similar increase in ORs in the ordinal model when comparing group 7 to group 1, where the OR was 6.3 (from Karlson et al, with 22 risk alleles) and 5.7 in our analysis that included 17 additional risk alleles. Similarly, the combination of risk alleles also displayed a good ability to discriminate between an RA case and control when the case is defined as seropositive RA, erosive RA or seropositive, erosive RA. However, the GRS showed very little, if any, ability to discriminate between seronegative RA and controls with an AUC of 0.563. When we compare the seropositive RA model with the 39 alleles to the one from Karlson et al. with 22 alleles we see no improvement from 0.660 (GRS22) to 0.654 (GRS39). This suggests that the addition of these 17 newly discovered RA alleles, whose individual ORs range from 1.10 to 1.23, does not improve the predictive ability of the GRS. As genetic discoveries progress with next generation sequencing, it is likely that cumulative GRS will improve in its predictive ability.

Our results for seronegative RA should be viewed in the context of prior research. The loci used in the GRS were discovered and the weights determined using only studies that include seropositive RA cases. Although there have been a few genetic markers that have associated with an increase risk of seronegative RA, such as HLA-DR1*03 [50], HLA-DR3 [51], and allelic forms of DCIR [52] and IRF5 [53], there may be as yet undiscovered loci that predict the seronegative RA phenotype. In a dataset containing 1500 cases/1500 controls, Kurreeman et al [54] demonstrated that a GRS based on 28 non-HLA risk alleles was associated with seronegative RA with an AUC of 0.55 and a p-value for a linear association of 0.0008 also suggesting only a very modest association for these risk alleles with the seronegative phenotype.

It has been shown that the HLA-SE is strongly associated with both RF status and presence of anti-CCP antibodies [9], [55], [56]. More specifically, anti-CCP antibodies play a vital role in the causal pathway between HLA-SE and erosions [29]. This is one explanation of the results demonstrating that GRS39 performs similarly when using erosive status to define severe disease, rather than seropositivity. In addition, the observation of an AUC of 0.712 for GRS39 identifying seropositive, erosive RA cases suggests that a narrower definition of RA leads to better discriminative ability. This lends support to the argument that RA falls along a severity continuum starting with seronegative as least severe and leading to seropositive, erosive RA as most severe.

We found that earlier age at onset of RA may potentially be associated with increased GRS. While the correlations were weak and not statistically significant, this does suggest that perhaps those with earlier age at RA onset have a higher β€œload” of genetic risk factors than those with later onset. Previous studies have shown an earlier age of diagnosis of RA both for those having any HLA-SE compared to none [26], [27], [28], and for any PTPN22 T allele compared to CCP [29]. Since both HLA and PTPN22 have a strong influence on the GRS score, this may be one explanation for the inverse relationship between the GRS and age at onset. The strongest effects that we detected for GRS and age at onset were with the seronegative and seropositive phenotypes. With this number of subjects, we had 37% and 35% statistical power to detect a significant ρ of 0.11 in seronegative and a ρ of βˆ’0.09 in seropositive RA. It is possible with more subjects in all phenotype groups we might have been able to detect significant relationships.

One limitation of our study is that we only have anti-CCP status tested at one time point, either up to 12 years prior to time of RA diagnosis or after diagnosis for the subset of cases without blood sample collected. The lack of information for anti-CCP results in the medical records due to the recent development of this test limits our ability to study anti-CCP results after diagnosis in all cases. We have not systematically collected outcome data after diagnosis of RA in this cohort, thus we do not know if some of the subjects defined as seronegative at diagnosis will later go on to convert to seropositive. This could lead to misclassification bias, with some truly seropositive RA subjects being misclassified as seronegative, which would bias us away from the null in the analysis. However, as we have found only modest associations within the seronegative group we do not believe that this has affected our analysis. This is also the case with erosive disease status, which based on chart data included notes ranging from the date of diagnosis where subjects have not had time to develop erosions to many years of follow-up. Another possible limitation to this study is the lack of data to test for population stratification. However, a subset of this sample (437 RA cases and 437 controls) [38] was genotyped for the lactase gene (rs4988235), known to exhibit substantial variation in allele frequency from Northern to Southern Europe [57], [58]. No significant differences were found between cases and controls, arguing strongly against any significant population stratification in this dataset.

In summary, many arguments have been made in the last few years for subdividing RA into different phenotypes [8], [9], [10], [11]. The analyses here add credence to these arguments. We demonstrate different genetic associations for the different RA sub-types, with only a modest relationship seen in the least severe phenotype, seronegative and the strongest relationship seen with the most severe phenotype, seropositive, erosive RA. This suggests that seropositive RA has a different underlying genetic basis than seronegative RA and thus, in future research, studying the two phenotypes separately would lead to greater understanding of the genetic and functional make-up of the disease.

Supporting Information

Table S1.

Genotype frequencies and association with seropositive RA in for 39 RA risk alleles.

https://doi.org/10.1371/journal.pone.0024380.s001

(DOCX)

Table S2.

Weighted GRS groups and odd ratios of All, seronegative and non-Erosive RA in NHS.

https://doi.org/10.1371/journal.pone.0024380.s002

(DOCX)

Table S3.

Relationship between weighted GRS as groups and as continuous and age at RA symptom start.

https://doi.org/10.1371/journal.pone.0024380.s003

(DOCX)

Acknowledgments

The authors would like to thank all the participants and staff of the Nurses' Health Study and Nurses' Health Study II.

Author Contributions

Conceived and designed the experiments: LBC KHC EWK. Performed the experiments: RMP. Analyzed the data: LBC BTK JC. Contributed reagents/materials/analysis tools: RMP EWK. Wrote the paper: LBC BTK KPL KHC EWK. Chart reviews: KPL KHC EWK.

References

  1. 1. Mikuls TR, Saag KG (2001) Comorbidity in rheumatoid arthritis. Rheum Dis Clin North Am 27: 283–303.
  2. 2. Kirwan JR (2001) Links between radiological change, disability, and pathology in rheumatoid arthritis. J Rheumatol 28: 881–886.
  3. 3. Clarke AE, St-Pierre Y, Joseph L, Penrod J, Sibley JT, et al. (2001) Radiographic damage in rheumatoid arthritis correlates with functional disability but not direct medical costs. J Rheumatol 28: 2416–2424.
  4. 4. Sokka T, Pincus T (2001) Markers for work disability in rheumatoid arthritis. J Rheumatol 28: 1718–1722.
  5. 5. Lee DM, Weinblatt ME (2001) Rheumatoid arthritis. Lancet 358: 903–911.
  6. 6. Finckh A, Choi HK, Wolfe F (2006) Progression of radiographic joint damage in different eras: trends towards milder disease in rheumatoid arthritis are attributable to improved treatment. Ann Rheum Dis 65: 1192–1197.
  7. 7. Finckh A, Liang MH, van Herckenrode CM, de Pablo P (2006) Long-term impact of early treatment on radiographic progression in rheumatoid arthritis: A meta-analysis. Arthritis Rheum 55: 864–872.
  8. 8. Ding B, Padyukov L, Lundstrom E, Seielstad M, Plenge RM, et al. (2009) 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. Arthritis Rheum 60: 30–38.
  9. 9. Huizinga TW, Amos CI, van der Helm-van Mil AH, Chen W, van Gaalen FA, et al. (2005) Refining the complex rheumatoid arthritis phenotype based on specificity of the HLA-DRB1 shared epitope for antibodies to citrullinated proteins. Arthritis Rheum 52: 3433–3438.
  10. 10. Klareskog L, Catrina AI, Paget S (2009) Rheumatoid arthritis. Lancet 373: 659–672.
  11. 11. Klareskog L, Wedren S, Alfredsson L (2009) On the origins of complex immune-mediated disease: the example of rheumatoid arthritis. J Mol Med 87: 357–362.
  12. 12. De Rycke L, Peene I, Hoffman IE, Kruithof E, Union A, et al. (2004) Rheumatoid factor and anticitrullinated protein antibodies in rheumatoid arthritis: diagnostic value, associations with radiological progression rate, and extra-articular manifestations. Ann Rheum Dis 63: 1587–1593.
  13. 13. Rantapaa-Dahlqvist S, de Jong BA, Berglin E, Hallmans G, Wadell G, et al. (2003) Antibodies against cyclic citrullinated peptide and IgA rheumatoid factor predict the development of rheumatoid arthritis. Arthritis Rheum 48: 2741–2749.
  14. 14. Ronnelid J, Wick MC, Lampa J, Lindblad S, Nordmark B, et al. (2005) 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. Ann Rheum Dis 64: 1744–1749.
  15. 15. van Leeuwen MA, Westra J, van Riel PL, Limburg PC, van Rijswijk MH (1995) IgM, IgA, and IgG rheumatoid factors in early rheumatoid arthritis predictive of radiological progression? Scand J Rheumatol 24: 146–153.
  16. 16. Bukhari M, Lunt M, Harrison BJ, Scott DG, Symmons DP, et al. (2002) Rheumatoid factor is the major predictor of increasing severity of radiographic erosions in rheumatoid arthritis: results from the Norfolk Arthritis Register Study, a large inception cohort. Arthritis Rheum 46: 906–912.
  17. 17. Jansen LM, van der Horst-Bruinsma IE, van Schaardenburg D, Bezemer PD, Dijkmans BA (2001) Predictors of radiographic joint damage in patients with early rheumatoid arthritis. Ann Rheum Dis 60: 924–927.
  18. 18. Reneses S, Gonzalez-Escribano MF, Fernandez-Suarez A, Pestana L, Davila B, et al. (2009) The value of HLA-DRB1 shared epitope, βˆ’308 tumor necrosis factor-alpha gene promoter polymorphism, rheumatoid factor, anti-citrullinated peptide antibodies, and early erosions for predicting radiological outcome in recent-onset rheumatoid arthritis. J Rheumatol 36: 1143–1149.
  19. 19. Vittecoq O, Pouplin S, Krzanowska K, Jouen-Beades F, Menard JF, et al. (2003) Rheumatoid factor is the strongest predictor of radiological progression of rheumatoid arthritis in a three-year prospective study in community-recruited patients. Rheumatology (Oxford) 42: 939–946.
  20. 20. Kaarela K (1985) Prognostic factors and diagnostic criteria in early rheumatoid arthritis. Scand J Rheumatol Suppl 57: 1–54.
  21. 21. Luukkainen R, Kaarela K, Isomaki H, Martio J, Kiviniemi P, et al. (1983) The prediction of radiological destruction during the early stage of rheumatoid arthritis. Clin Exp Rheumatol 1: 295–298.
  22. 22. Pease CT, Bhakta BB, Devlin J, Emery P (1999) Does the age of onset of rheumatoid arthritis influence phenotype?: a prospective study of outcome and prognostic factors. Rheumatology (Oxford) 38: 228–234.
  23. 23. Terkeltaub R, Esdaile J, Decary F, Tannenbaum H (1983) A clinical study of older age rheumatoid arthritis with comparison to a younger onset group. J Rheumatol 10: 418–424.
  24. 24. van der Heijde DM, van Riel PL, van Leeuwen MA, van 't Hof MA, van Rijswijk MH, et al. (1991) Older versus younger onset rheumatoid arthritis: results at onset and after 2 years of a prospective followup study of early rheumatoid arthritis. J Rheumatol 18: 1285–1289.
  25. 25. Bukhari M, Lunt M, Barton A, Bunn D, Silman A, et al. (2007) Increasing age at symptom onset is associated with worse radiological damage at presentation in patients with early inflammatory polyarthritis. Ann Rheum Dis 66: 389–393.
  26. 26. Hazes JM, Dijkmans BA, Hoevers JM, Janson JJ, de Vries RR, et al. (1989) DR4 prevalence related to the age at disease onset in female patients with rheumatoid arthritis. Ann Rheum Dis 48: 406–408.
  27. 27. Jaraquemada D, Ollier W, Awad J, Young A, Silman A, et al. (1986) HLA and rheumatoid arthritis: a combined analysis of 440 British patients. Ann Rheum Dis 45: 627–636.
  28. 28. Yukioka M, Wakitani S, Murata N, Toda Y, Ogawa R, et al. (1998) Elderly-onset rheumatoid arthritis and its association with HLA-DRB1 alleles in Japanese. Br J Rheumatol 37: 98–101.
  29. 29. Karlson EW, Chibnik LB, Cui J, Plenge RM, Glass RJ, et al. (2008) Associations between human leukocyte antigen, PTPN22, CTLA4 genotypes and rheumatoid arthritis phenotypes of autoantibody status, age at diagnosis and erosions in a large cohort study. Ann Rheum Dis 67: 358–363.
  30. 30. Bajocchi G, La Corte R, Locaputo A, Govoni M, Trotta F (2000) Elderly onset rheumatoid arthritis: clinical aspects. Clin Exp Rheumatol 18: S49–50.
  31. 31. Karlson EW, Chibnik LB, Kraft P, Cui J, Keenan BT, et al. (2010) Cumulative association of 22 genetic variants with seropositive rheumatoid arthritis risk. Ann Rheum Dis 69: 1077–1085.
  32. 32. Lohmueller KE, Pearce CL, Pike M, Lander ES, Hirschhorn JN (2003) Meta-analysis of genetic association studies supports a contribution of common variants to susceptibility to common disease. Nat Genet 33: 177–182.
  33. 33. Raychaudhuri S, Remmers EF, Lee AT, Hackett R, Guiducci C, et al. (2008) Common variants at CD40 and other loci confer risk of rheumatoid arthritis. Nat Genet 40: 1216–1223.
  34. 34. Raychaudhuri S, Thomson BP, Remmers EF, Eyre S, Hinks A, et al. (2009) Genetic variants at CD28, PRDM1 and CD2/CD58 are associated with rheumatoid arthritis risk. Nat Genet 41: 1313–1318.
  35. 35. Stahl EA, Raychaudhuri S, Remmers EF, Xie G, Eyre S, et al. (2010) Genome-wide association study meta-analysis identifies seven new rheumatoid arthritis risk loci. Nature Genetics. in press.
  36. 36. Liang MH, Meenan RF, Cathcart ES, Schur PH (1980) A screening strategy for population studies in systemic lupus erythematosus. Series design. Arthritis Rheum 23: 153–157.
  37. 37. Karlson EW, Chibnik LB, Tworoger SS, Lee IM, Buring JE, et al. (2009) Biomarkers of inflammation and development of rheumatoid arthritis in women from two prospective cohort studies. Arthritis Rheum 60: 641–652.
  38. 38. Costenbader KH, Chang SC, De Vivo I, Plenge R, Karlson EW (2008) Genetic polymorphisms in PTPN22, PADI-4, and CTLA-4 and risk for rheumatoid arthritis in two longitudinal cohort studies: evidence of gene-environment interactions with heavy cigarette smoking. Arthritis Res Ther 10: R52.
  39. 39. De Jager PL, Chibnik LB, Cui J, Reischl J, Lehr S, et al. (2009) Integration of genetic risk factors into a clinical algorithm for multiple sclerosis susceptibility: a weighted genetic risk score. Lancet Neurol 8: 1111–1119.
  40. 40. Karlson EW, Chibnik LB, Kraft P, Cui J, Keenan BT, et al. (2010) Cumulative association of 22 genetic variants with seropositive rheumatoid arthritis risk. Ann Rheum Dis. in press.
  41. 41. Fernando MM, Stevens CR, Walsh EC, De Jager PL, Goyette P, et al. (2008) Defining the role of the MHC in autoimmunity: a review and pooled analysis. PLoS Genet 4: e1000024.
  42. 42. Begovich AB, Carlton VE, Honigberg LA, Schrodi SJ, Chokkalingam AP, et al. (2004) A missense single-nucleotide polymorphism in a gene encoding a protein tyrosine phosphatase (PTPN22) is associated with rheumatoid arthritis. Am J Hum Genet 75: 330–337.
  43. 43. Plenge RM, Seielstad M, Padyukov L, Lee AT, Remmers EF, et al. (2007) TRAF1-C5 as a risk locus for rheumatoid arthritis–a genomewide study. N Engl J Med 357: 1199–1209.
  44. 44. Remmers EF, Plenge RM, Lee AT, Graham RR, Hom G, et al. (2007) STAT4 and the risk of rheumatoid arthritis and systemic lupus erythematosus. N Engl J Med 357: 977–986.
  45. 45. Plenge RM, Cotsapas C, Davies L, Price AL, de Bakker PI, et al. (2007) Two independent alleles at 6q23 associated with risk of rheumatoid arthritis. Nat Genet 39: 1477–1482.
  46. 46. Costenbader KH, Feskanich D, Mandl LA, Karlson EW (2006) Smoking intensity, duration, and cessation, and the risk of rheumatoid arthritis in women. Am J Med 119: 503.e501–509
  47. 47. Kallberg H, Ding B, Padyukov L, Bengtsson C, Ronnelid J, et al. (2010) Smoking is a major preventable risk factor for rheumatoid arthritis: estimations of risks after various exposures to cigarette smoke. Ann Rheum Dis.
  48. 48. Stolt P, Bengtsson C, Nordmark B, Lindblad S, Lundberg I, et al. (2003) Quantification of the influence of cigarette smoking on rheumatoid arthritis: results from a population based case-control study, using incident cases. Ann Rheum Dis 62: 835–841.
  49. 49. DeLong ER, DeLong DM, Clarke-Pearson DL (1988) Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics 44: 837–845.
  50. 50. Lundstrom E, Kallberg H, Smolnikova M, Ding B, Ronnelid J, et al. (2009) Opposing effects of HLA-DRB1*13 alleles on the risk of developing anti-citrullinated protein antibody-positive and anti-citrullinated protein antibody-negative rheumatoid arthritis. Arthritis Rheum 60: 924–930.
  51. 51. Verpoort KN, van Gaalen FA, van der Helm-van Mil AH, Schreuder GM, Breedveld FC, et al. (2005) Association of HLA-DR3 with anti-cyclic citrullinated peptide antibody-negative rheumatoid arthritis. Arthritis Rheum 52: 3058–3062.
  52. 52. Lorentzen JC, Flornes L, Eklow C, Backdahl L, Ribbhammar U, et al. (2007) Association of arthritis with a gene complex encoding C-type lectin-like receptors. Arthritis Rheum 56: 2620–2632.
  53. 53. Sigurdsson S, Padyukov L, Kurreeman FA, Liljedahl U, Wiman AC, et al. (2007) Association of a haplotype in the promoter region of the interferon regulatory factor 5 gene with rheumatoid arthritis. Arthritis Rheum 56: 2202–2210.
  54. 54. Kurreeman F, Liao K, Chibnik L, Hickey B, Stahl E, et al. (2010) Genetic basis of autoantibody positive and negative rheumatoid arthritis risk in a multi-ethnic cohort derived from electronic health records. Am J Hum Genet 88: 57–69.
  55. 55. Gorman JD, David-Vaudey E, Pai M, Lum RF, Criswell LA (2004) Particular HLA-DRB1 shared epitope genotypes are strongly associated with rheumatoid vasculitis. Arthritis Rheum 50: 3476–3484.
  56. 56. van der Helm-van Mil AH, Verpoort KN, Breedveld FC, Huizinga TW, Toes RE, et al. (2006) 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 Rheum 54: 1117–1121.
  57. 57. Bersaglieri T, Sabeti PC, Patterson N, Vanderploeg T, Schaffner SF, et al. (2004) Genetic signatures of strong recent positive selection at the lactase gene. Am J Hum Genet 74: 1111–1120.
  58. 58. Campbell CD, Ogburn EL, Lunetta KL, Lyon HN, Freedman ML, et al. (2005) Demonstrating stratification in a European American population. Nat Genet 37: 868–872.