Skip to main content
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

Genetic associations with radiological damage in rheumatoid arthritis: Meta-analysis of seven genome-wide association studies of 2,775 cases

  • Matthew Traylor,

    Roles Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

    Affiliations Department of Clinical Neurosciences, Stroke Research Group, University of Cambridge, Cambridge, United Kingdom, Department of Medical and Molecular Genetics, King’s College London, London, United Kingdom

  • Rachel Knevel,

    Roles Formal analysis, Investigation, Methodology, Writing – review & editing

    Affiliations Brigham and Women’s Hospital, Division of Genetics, Raychaudhuri Lab, Boston, MA, United States of America, Broad institute, Cambridge, MA, United States of America, Department of Rheumatology C1-R, Leiden University Medical Center, Albinusdreef, Leiden, the Netherlands

  • Jing Cui,

    Roles Formal analysis, Investigation, Methodology, Writing – review & editing

    Affiliation Division of Rheumatology Immunology and Allergy Brigham & Women's Hospital Harvard Medical School Boston, MA, United States of America

  • John Taylor,

    Roles Formal analysis, Investigation, Methodology, Writing – review & editing

    Affiliation Leeds Institute of Cancer & Pathology, Worsley Building Level 11 (LIDA), Clarendon Way, Leeds, United Kingdom

  • Westra Harm-Jan,

    Roles Formal analysis

    Affiliations Brigham and Women’s Hospital, Division of Genetics, Raychaudhuri Lab, Boston, MA, United States of America, Broad institute, Cambridge, MA, United States of America

  • Philip G. Conaghan,

    Roles Data curation, Writing – review & editing

    Affiliations Leeds Institute of Rheumatic and Musculoskeletal Medicine, University of Leeds, Leeds, United Kingdom, NIHR Leeds Biomedical Research Centre, Leeds Teaching Hospitals NHS Trust, Leeds, United Kingdom

  • Andrew P. Cope,

    Roles Conceptualization, Data curation, Funding acquisition, Writing – review & editing

    Affiliation Academic Department of Rheumatology, Centre for Molecular and Cellular Biology of Inflammation, King's College London, London, United Kingdom

  • Charles Curtis,

    Roles Data curation, Writing – review & editing

    Affiliations NIHR Biomedical Research Centre at South London and Maudsley NHS Foundation Trust and King’s College London, London, United Kingdom, SGDP Centre, Institute of Psychiatry, Psychology & Neuroscience, King’s College London, London, United Kingdom

  • Paul Emery,

    Roles Data curation, Writing – review & editing

    Affiliations Leeds Institute of Rheumatic and Musculoskeletal Medicine, University of Leeds, Leeds, United Kingdom, NIHR Leeds Biomedical Research Centre, Leeds Teaching Hospitals NHS Trust, Leeds, United Kingdom

  • Stephen Newhouse,

    Roles Data curation, Writing – review & editing

    Affiliations NIHR Biomedical Research Centre at South London and Maudsley NHS Foundation Trust and King’s College London, London, United Kingdom, Department of Biostatistics and Health Informatics, Institute of Psychiatry, Psychology and Neuroscience, King’s College London, London, United Kingdom, Farr Institute of Health Informatics Research, UCL Institute of Health Informatics, University College London, London, United Kingdom

  • Hamel Patel,

    Roles Data curation, Writing – review & editing

    Affiliation NIHR Biomedical Research Centre at South London and Maudsley NHS Foundation Trust and King’s College London, London, United Kingdom

  • Sophia Steer,

    Roles Conceptualization, Data curation, Writing – review & editing

    Affiliation Department of Rheumatology, King’s College Hospital, Denmark Hill, London, United Kingdom

  • Peter Gregersen,

    Roles Data curation, Writing – review & editing

    Affiliation The Feinstein Institute for Medical Research, Northwell Health, Manhasset, New York, United States of America

  • Nancy A. Shadick,

    Roles Data curation, Writing – review & editing

    Affiliation Division of Rheumatology Immunology and Allergy Brigham & Women's Hospital Harvard Medical School Boston, MA, United States of America

  • Michael E. Weinblatt,

    Roles Data curation, Writing – review & editing

    Affiliation Division of Rheumatology Immunology and Allergy Brigham & Women's Hospital Harvard Medical School Boston, MA, United States of America

  • Annette Van Der Helm-van Mil,

    Roles Data curation, Writing – review & editing

    Affiliation Department of Rheumatology C1-R, Leiden University Medical Center, Albinusdreef, Leiden, the Netherlands

  • Jennifer H. Barrett,

    Roles Conceptualization, Investigation, Methodology, Writing – original draft, Writing – review & editing

    Affiliations NIHR Leeds Biomedical Research Centre, Leeds Teaching Hospitals NHS Trust, Leeds, United Kingdom, School of Medicine, University of Leeds, Leeds, United Kingdom

  • Ann W. Morgan,

    Roles Conceptualization, Data curation, Investigation, Methodology, Writing – original draft, Writing – review & editing

    Affiliations Leeds Institute of Rheumatic and Musculoskeletal Medicine, University of Leeds, Leeds, United Kingdom, NIHR Leeds Biomedical Research Centre, Leeds Teaching Hospitals NHS Trust, Leeds, United Kingdom

  • Cathryn M. Lewis,

    Roles Conceptualization, Funding acquisition, Investigation, Methodology, Writing – original draft, Writing – review & editing

    Affiliations Department of Medical and Molecular Genetics, King’s College London, London, United Kingdom, SGDP Centre, Institute of Psychiatry, Psychology & Neuroscience, King’s College London, London, United Kingdom

  •  [ ... ],
  • Ian C. Scott

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Writing – original draft, Writing – review & editing

    i.scott@keele.ac.uk

    Affiliations Primary Care Centre Versus Arthritis, Research Institute for Primary Care and Health Sciences, Primary Care Sciences, Keele University, Keele, United Kingdom, Haywood Academic Rheumatology Centre, Haywood Hospital, Midlands Partnership NHS Foundation Trust, High Lane, Burslem, Staffordshire, United Kingdom

  • [ view all ]
  • [ view less ]

Abstract

Background

Previous studies of radiological damage in rheumatoid arthritis (RA) have used candidate-gene approaches, or evaluated single genome-wide association studies (GWAS). We undertook the first meta-analysis of GWAS of RA radiological damage to: (1) identify novel genetic loci for this trait; and (2) test previously validated variants.

Methods

Seven GWAS (2,775 RA cases, of a range of ancestries) were combined in a meta-analysis. Radiological damage was assessed using modified Larsen scores, Sharp van Der Heijde scores, and erosive status. Single nucleotide polymophsim (SNP) associations with radiological damage were tested at a single time-point using regression models. Primary analyses included age and disease duration as covariates. Secondary analyses also included rheumatoid factor (RF). Meta-analyses were undertaken in trans-ethnic and European-only cases.

Results

In the trans-ethnic primary meta-analysis, one SNP (rs112112734) in close proximity to HLA-DRB1, and strong linkage disequilibrium with the shared-epitope, attained genome-wide significance (P = 4.2x10-8). In the secondary analysis (adjusting for RF) the association was less significant (P = 1.7x10-6). In both trans-ethnic primary and secondary meta-analyses 14 regions contained SNPs with associations reaching P<5x10-6; in the European primary and secondary analyses 13 and 10 regions contained SNPs reaching P<5x10-6, respectively. Of the previously validated SNPs for radiological progression, only rs660895 (tagging HLA-DRB1*04:01) attained significance (P = 1.6x10-5) and had a consistent direction of effect across GWAS.

Conclusions

Our meta-analysis confirms the known association between the HLA-DRB1 shared epitope and RA radiological damage. The lack of replication of previously validated non-HLA markers highlights a requirement for further research to deliver clinically-useful prognostic genetic markers.

Background

Rheumatoid arthritis (RA) is a heterogeneous disease, exhibiting a variable course between patients. Prospectively identifying patients likely to develop severe phenotypes could allow early intensive therapy to be focussed on poor prognosis cases. This approach should optimise clinical and cost-effectiveness, but requires accurate prognostic markers.

Radiological damage is one measure of RA severity. It is moderately heritable (heritability from common variants estimated at 45%-58%) suggesting genetic markers could represent useful prognostic biomarkers [1]. The strongest genetic association with RA radiological progression is the shared epitope (SE), which represents consensus amino acid sequences (QRRAA, RRRAA and QKRAA) spanning positions 70–74 in the HLA-DRβ1 molecule, encoded by various HLA-DRB1 SE alleles [2]. The SE has been demonstrated to associate with greater radiological damage in a broad range of RA populations [3,4], although as it associates with the presence of rheumatoid factor (RF) and antibodies to citrullinated peptide antigens (ACPA) [5], both of which independently associate with radiological damage [6], its impact may be mediated by autoantibody status. Van Steenbergen et al previously identified a further 16 “validated” non-HLA variants (in 11 genes) for RA radiological progression [7]. In the Leiden Early Arthritis Cohort (EAC) these variants, in combination with the SE, explained up to 18% of the variance in radiological progression, although as this cohort was used to identify many of the variants, this finding requires replication.

The gold standard to identify genetic associations with a phenotype or trait is to perform a meta-analysis of all available genome-wide association studies (GWAS). To date, this has not been performed for RA radiological damage. Five individual GWAS have identified four non-HLA variants − rs7607479 in SPAG16, rs451066 in RAD51L1-ZFP36L1, rs11908352 in SLC12A5, and rs2833522 in HUNK/SCAF4 –associating with RA radiological damage that passed multiple testing correction thresholds, and replicated externally [811]. These studies were limited by modest sample sizes, with the largest containing 646 patients. Combining GWAS in a meta-analysis should increase statistical power to detect novel loci.

To this end, we have carried out the largest GWAS of RA radiological damage, by peforming and then combining seven independent GWAS (totalling 2,775 patients). We aimed to identify novel genetic loci for radiological damage. We also tested previously validated single nucleotide polymorphisms (SNPs) for their association with radiological damage.

Methods

Patients

GWAS were undertaken in: (1) Combination Anti-Rheumatic Drugs in Early RA (CARDERA) Genetics Cohort; (2) Yorkshire Early Arthritis Registry (YEAR); (3) Brigham and Women’s Hospital RA Sequential Study (BRASS); (4) Leiden Early Arthritis Clinic (EAC); (5) North American RA Consortium (NARAC); (6) GENetics of RA in individuals of African ancestry (GENRA) study; (7) South London RA Study (SLRAS). These cohorts have been described in detail previously [1217]. An overview is provided in Table 1 (with further details in S1 Table).

Radiological measures

Hand and feet X-rays were scored using either the Scott modification of the Larsen method (scoring 0–200 for joint space narrowing and erosive damage [18]), Sharp van Der Heidje Scores (SvHS; scoring 0–448 for joint space narrowing and erosive damage [19]) or erosions (present/absent). Larsen and SvHS are highly correlated; both evaluate erosions [20].

Genotyping, quality control and imputation

CARDERA was genotyped on the Immunochip, YEAR on the ImmunoChip and HumanOmniExpressExome Beadchip, BRASS on the Affymetrix 6.0, Leiden EAC on the Illumina iScan, GENRA on the Illumina MEGA array, NARAC on the Illumina Beadchip (HumanHap 550k), and SLRAS on the Illumina HumanOmniExpress. Prior to imputation all cohorts excluded SNPs with high levels of missingness, low minor allele frequency (MAF), and deviations from Hardy-weinberg equilibrium (HWE). Individuals were removed that did not segregate with reference populations, had high-levels of missingness, or whose phenotypic sex mismatched with sex inferred from genotype data. Genotype phasing and imputation were performed using SHAPEIT and IMPUTE2, respectively. Data were imputed to the 1,000 Genomes Phase 3 Panel, including all samples regardless of ancestry. Post-imputation markers with INFO scores <0.7 or MAF <0.05 were removed (the latter was undertaken as we had limited power to detect associations with low frequency variants). Genotypes were analysed as expected allele “dosages” in all analyses. Genotypes were converted to be reported on the forward strand. The number of variants available post-QC comprised: CARDERA: 3,181,676; YEAR: 6,209,766; Leiden EAC: 674,614; BRASS: 6,059,126; GENRA: 5,662,513; NARAC: 5,373,610; SLRAS: 5,715,821.”

Statistical analysis

X-ray time-point.

As RA X-ray progression can be non-linear [21], and several cohorts had cross-sectional X-ray scores, we tested SNP associations with radiological damage at a single time-point. In the repeated measures early RA cohorts (CARDERA, YEAR, and Leiden EAC) associations with end-point X-ray scores were tested (when radiological damage was greatest). In the repeated measures established RA cohort (BRASS) associations with baseline X-ray scores were tested (when sample size was greatest).

Statistical models.

Log-transformed Larsen/SvHS, or binary erosive status were used as response variables. Outcomes were regressed on imputed genotype dosages assuming an additive genetic model using linear (in CARDERA, Leiden EAC, BRASS, SLRAS, and NARAC), negative binomial (in YEAR, due to over-dispersion) or logistic (in GENRA) regression models.

The same modelling covariates were included to provide consistency. Sex, age, disease duration, rheumatoid factor (RF) and (where available) treatment were tested for associations with X-ray scores in each cohort to determine the best predictors across studies. Only age and disease duration were associated with radiological scores in every cohort, and these were included as covariates. Ancestry-informative principal components (PCs) were included to account for population stratification. These were calculated using EIGENSTRAT. Different numbers were used in each study, with a minimum of 2 used [22].

Two analyses were undertaken: a primary analysis using these covariates, and a secondary analysis, including RF as an additional covariate to account for genetic variants exerting their effects through RF formation (which associates with both the SE and X-ray damage).

Meta-analysis.

SNP association results from GWAS were combined using a Z-score weighted approach in METAL [23]. METAL performs a meta-analysis using information on P-values, effect directions, and the number of cases evaluated at each SNP. This method was used because other meta-analysis methods, such as the inverse variance-weighted average method [24], require effect size data at each marker, and as our various GWAS utlised different statistical models and X-ray scoring systems, their beta-values were non-comparable.

We assessed QQ-plots and genomic inflation factors to confirm that genome-wide statistics had the expected distributions. Genomic control correction was used to adjust for any residual inflation of test statistics. In the overall meta-analysis of all SNPs, we only analysed SNPs that had been genotyped or imputed in at least half of all individuals in the total sample.

Two meta-analyses were performed: (1) a trans-ethnic meta-analysis (maximising sample size); (2) a European ancestry meta-analysis (reducing heterogeneity across populations).

Associations with previously validated SNPs.

Van Steenbergen et al reported 17 validated variants for radiological progression. These comprised variants in the following 12 loci: SPAG16, IL-15, C5orf30, HLA-DRB1, OPG, DKK-1, IL2RA, RAD51L1/ZFP36L1, GRZB, IL-4R, MMP-9, and CD40. They defined “validated” as being a SNP studied in several cohorts, with the association independently replicated or found significant in a meta-analysis of all published data [7]. We tested their association with radiological damage in our meta-analysis. We included all available SNPs, regardless of if they were present in less than half of all individuals. We used the HLA-DRB1*04:01 tagging SNP (rs660895) identified in the Eyre et al RA susceptibility meta-analysis [25] to represent the SE, as this is the commonest SE-encoding allele in RA cases [2].

Significance thresholds.

For the testing of previously validated SNPs we used a Bonferroni corrected P-value of 0.004167 (12 loci tested). For the genome-wide analyses we used the standard P-value threshold of <5x10-8 for genome-wide significance, and <5x10-6 for “suggestive” significance. In the genome-wide trans-ethnic analysis, based on sampling from a Chi-squared distribution with , where q is trait variance and N is sample size, we had 80% power to detect a variant explaining 1.4% of the variance of X-ray damage at P<5x10-8 [26]. We made assumptions to enable us to perform a power calculation for this heterogeneous dataset. As most datasets used the quantitative Larsen scale, we simulated chi-squared statistics with an NCP given by N*qtl_variance/(1-qtl_variance) and calculated power based on the number of observations reaching genome-wide significance, using the principles reported by Yang et al [26]. We assumed we had 80% study power for qtl_variance values which produced genome-wide significant values for 80% of simulations.

Ethics approval and consent to participate

Ethical approval was granted for each of the genetic studies as follows: CARDERA genetics cohort was approved by the National Research Ethics Service Committee East of England—Essex (reference: 11/EE/0544); YEAR by the Multi-Centre Research Ethics Committee (MREC) (reference 99/3/48); Leiden EAC by the Medical Ethics Committee Leiden University Medical Center; BRASS by the Institutional Review Board of Brigham and Women’s Hospital, Partners Health Care, Boston, MA; GENRA by the National Research Ethics Service Committee London—Dulwich (reference: 11/LO/1244); NARAC by the Via Christi Institutional Review Board and North Shore-LIJ Health System Institutional Review Board; SLRAS by the Guy's Hospital Local Research Committee (reference 99/11/06) and Lewisham Hospital Local Research Committee (reference 01/05/02). Informed consent was provided by all participating patients.

Availability of data and materials

Individual-level genotype and phenotype data from each of the separate GWAS are not publically available, as the centres from which patients were recruited are reported in the primary publications outlining the individual genetic datasets, leading to the potential for patient confidentiality to be affected. Summary statistics will be made available through the NHGRI-EBI GWAS Catalog (https://www.ebi.ac.uk/gwas/downloads/summary-statistics). Further data requests can be made to the following non-author contacts: Ms Isabel Sinha (isabel.sinha@kcl.ac.uk) for CARDERA, GENRA, and SLRAS; Dr James Robinson (J.I.Robinson@leeds.ac.uk) for YEAR; Vivi Feathers (vfeathers@bwh.harvard.edu) for BRASS; Dr Annette Lee (Anlee@northwell.edu) for NARAC; and Dr E Niemantsverdriet (E.Niemantsverdriet@lumc.nl) for the Leiden EAC.

Results

Cohort characteristics

Cohort characteristics are presented in Table 1. The majority of cases were female (67.2–80.8%), and seropositive (57.4–100% RF-positive; 52.8–100% ACPA-positive). The mean age in all cohorts was between 54.5 and 60.9 years, except NARAC, which had a mean age of 40.8 years. As expected, the cohorts containing established RA patients had higher rates of radiological damage than those evaluating early RA patients. The total number of patients included in the trans-ethnic and European meta-analyses comprised 2,775 and 2,527 cases, respectively. The total number of SNPs assessed in the trans-ethnic meta-analysis was 4,802,696, and in the European meta-analysis was 2,723,488.

Trans-ethnic meta-analysis

In the primary analysis (age and disease duration as covariates), one SNP (rs112112734, on chromosome 6, in close proximity to HLA-DRB5, HLA-DRB6 and HLA-DRB1) attained genome-wide significance (P = 4.2x10-8, regional association plot in Fig 1). The same direction of effect was observed in the 6 GWAS in which this marker was available (it was absent in the Leiden EAC). This SNP is in strong linkage disequilibrium (LD) with the HLA-DRB1 SE (R2 = 0.93 with the lead HLA-DRB1 SNP, rs9268839, from the Okada et al RA susceptibility meta-analysis [27]). Fourteen regions contained SNPs with associations reaching P<5x10-6 (S2 Table).

thumbnail
Fig 1. Regional association plot for region surrounding rs112112734 showing evidence of significance at P<5x10-8 in the primary trans-ethnic meta-analysis.

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

In the secondary analysis (age, disease duration and RF as covariates), no SNPs attained genome-wide significance. Fourteen regions contained SNPs with associations reaching P<5x10-6, one of which was rs112112734 (P = 1.7x10-6) (S3 Table). Plots of–log10(P-value) by genomic position for the trans-ethnic and European analyses are provided in Fig 2, and QQ plots in S1 Fig.

thumbnail
Fig 2. Manhattan plots for radiological damage.

Panel A = trans-ethnic primary analysis; Panel B = European primary analysis; Panel C = trans-ethnic secondary analysis; Panel D = European secondary analysis; primary analysis includes age, disease duration and genotypes as explanatory variables; secondary analysis also includes RF as an explanatory variable.

https://doi.org/10.1371/journal.pone.0223246.g002

European meta-analysis

In the primary analysis, thirteen regions contained SNPs reaching P<5x10-6 (S2 Table), with the most significant SNP being rs112112734 (P = 8.6x10-8).

In the secondary analysis, ten regions contained SNPs with associations reaching P<5x10-6 (S3 Table). rs112112734 was not amongst these.

Previously validated SNPs

Two SNPs associated with radiological damage in our European meta-analysis (Table 2). These comprised: (1) rs660895 (tagging HLA-DRB1*04:01; P = 1.6x10-5) which had the same direction of effect across all 6 GWAS (the G allele, indicating an *04:01 copy, associating with increased damage), and (2) rs7607479 (in SPAG16), which was available in four GWAS, two of which had increased and two reduced radiological damage associated with the T allele. In the trans-ethnic meta-analysis, only rs660895 was significant (P = 7.0x10-5) with the same direction of effect observed in all but one GWAS. Only 2 SNPs in the trans-ethnic meta-analysis, and 4 SNPs in the European meta-analysis had a shared direction of effect across all available GWAS.

thumbnail
Table 2. Previously validated genetic variants for radiological progression reported by Van Steenbergen et al [7], and their association with radiological damage.

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

Discussion

Our study has three key findings. Firstly, it confirms the known association between this trait and the HLA-region, with one SNP in this region (rs112112734, in high LD with the HLA-DRB1 SE) attaining genome-wide significance. Secondly, it demonstrates the challenges in replicating genetic variants for traits across GWAS. Of the 17 previously validated SNPs, only the HLA-DRB1*04:01 tagging SNP had a significant association and consistent direction of effect across GWAS. Thirdly, it highlights the difficulties in performing meta-analysis of GWAS of continuous disease outcomes like RA radiological damage, which involves combining summary statistics from highly heterogeneous patient cohorts.

The association between the SE and radiological damage in RA is well established [3]. As the SE associates with RF and ACPA, which themselves are linked with radiological damage, we undertook a secondary analysis including RF as a covariate (ACPA was not used, as ACPA status was unavailable in SLRAS, and two cohorts only included ACPA-positive cases). The association observed at chromosome 6 was less significant in this secondary analysis, and rs112112734 no longer reached genome-wide significance. This suggests the observed association between the HLA-region and X-ray damage in our meta-analysis may be, at least in part, mediated by autoantibody production. More recently, Viatte et al reported that amino acid polymorphisms at position 11 in HLA-DRβ1 also associate with erosions, independently of the shared epitope [28]; owing to an absence of imputed amino acid polymorphism data, and the fact that this amino acid polymorphism is not “tagged” by a single SNP, we were unable to evaluate this association within our meta-analysis. Future research should focus on fine mapping of the HLA locus to provide a more comprehensive understanding of the association between HLA alleles and X-ray damage.

Our failure to replicate associations for previously validated non-HLA markers is disappointing. There are several potential explanations for this. Firstly, the validated variants reported by Van Steenbergen et al were for radiological progression, and we tested their association with radiological damage (although only radiological “progressors” will accrue damage). Secondly, several of these variants were identified in ACPA-positive cases only, and our cohorts included cases with a range of frequencies of ACPA. Thirdly, not all SNPs were available in all cohorts, although only one SNP (rs8192916 in GRZB) had a sample size of <1,000 cases. Irrespective of the explanation, current data on seven GWAS strengthen the notion that these variants lack clinical utility at identifying severe phenotypes [7].

The strengths of our study are that it represents the largest analysis of genetic predictors of RA radiological damage, includes individuals from a range of ethnicities (increasing the generalisability of its results), and is the first meta-analysis of an RA prognostic trait. Its weaknesses are the marked heterogeneity of the included GWAS (which varied in their disease durations, X-ray scoring systems used, and statistical models), and its modest sample size, limiting our ability to detect variants of a small effect size. As our meta-analysis replicated the genetic locus with the largest effect size for radiological damage (the HLA region), we consider that our modest sample size is the most likely reason we did not identify non-HLA associations attaining genome-wide significance. An alternative explanation is that the trait of radiological damage has low heritability. Due to a modest sample size, and clinical heterogeneity, we could not determine heritability in our GWAS meta-analysis with confidence, although previous research has reported RA radiological damage to be moderately heritable (estimated heritability rate of 45–58%) [1].

Conclusions

In conclusion, our meta-analysis replicates the association between the HLA-region and radiological damage in RA. It demonstrates the complexities of undertaking meta-analysis of GWAS of RA characteristics. For these to be more precisely defined, a collaborative international approach is required, standardising the collection of phenotypic data in genotyped patients, creating large homogeneous cohorts that are suitable for meta-analysis.

Supporting information

S1 Table. Genome-Wide association study details.

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

(PDF)

S2 Table. Lead SNPs from regions showing evidence of significance at P<5x10-6 in the primary trans-ethnic and European meta-analysis.

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

(PDF)

S3 Table. Lead SNPs from regions showing evidence of significance at P<5x10-6 in the secondary trans-ethnic and European meta-analysis.

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

(PDF)

S1 Fig. QQ plots for meta-analyses.

A = primary trans-ethnic meta-analysis (inflation factor 1.00); B = primary European meta-analysis (inflation factor 1.01); C = secondary trans-ethnic meta-analysis (inflation factor 1.00); D = secondary European meta-analysis (inflation factor 1.00).

https://doi.org/10.1371/journal.pone.0223246.s004

(TIF)

References

  1. 1. Knevel R, Grondal G, Huizinga TWJ, Visser AW, Jonsson H, Vikingsson A, et al. Genetic predisposition of the severity of joint destruction in rheumatoid arthritis: a population-based study. Ann Rheum Dis. 2012;71: 707–709. pmid:22219137
  2. 2. Raychaudhuri S, Sandor C, Stahl EA, Freudenberg J, Lee H-S, Jia X, 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–296. pmid:22286218
  3. 3. Gorman J, MD MPH, Lum R, Chen J, Suarez-Almazor M, MD P, et al. Impact of Shared Epitope Genotype and Ethnicity on Erosive Disease: A Meta-Analysis of 3,240 Rheumatoid Arthritis Patients. Arthritis Rheum. 2004;50: 400–412. pmid:14872482
  4. 4. Gonzalez-Gay MA, Garcia-Porrua C, Hajeer AH. Influence of human leukocyte antigen-DRB1 on the susceptibility and severity of rheumatoid arthritis. Semin Arthritis Rheum. 2002;31: 355–360. pmid:12077707
  5. 5. Balsa A, Cabezon A, Orozco G, Cobo T, Miranda-Carus E, Lopez-Nevot MA, et al. Influence of HLA DRB1 alleles in the susceptibility of rheumatoid arthritis and the regulation of antibodies against citrullinated proteins and rheumatoid factor. Arthritis Res Ther. 2010;12: R62. pmid:20370905
  6. 6. Mewar D, Coote A, Moore DJ, Marinou I, Keyworth J, Dickson MC, et al. Independent associations of anti-cyclic citrullinated peptide antibodies and rheumatoid factor with radiographic severity of rheumatoid arthritis. Arthritis Res Ther. 2006;8: R128. pmid:16859535
  7. 7. van Steenbergen HW, Tsonaka R, Huizinga TW, le Cessie S, van der Helm-van Mil AH. Predicting the severity of joint damage in rheumatoid arthritis; the contribution of genetic factors. Ann Rheum Dis. 2015;74: 876–882. pmid:24431388
  8. 8. Arya R, Del Rincon I, Farook VS, Restrepo JF, Winnier DA, Fourcaudot MJ, et al. Genetic Variants Influencing Joint Damage in Mexican Americans and European Americans With Rheumatoid Arthritis. Genet Epidemiol. 2015;39: 678–688. pmid:26498133
  9. 9. de Rooy DPC, Zhernakova A, Tsonaka R, Willemze A, Kurreeman BAS, Trynka G, et al. A genetic variant in the region of MMP-9 is associated with serum levels and progression of joint damage in rheumatoid arthritis. Ann Rheum Dis. 2014;73: 1163–1169. pmid:23696630
  10. 10. Knevel R, Klein K, Somers K, Ospelt C, Houwing-Duistermaat JJ, van Nies JAB, et al. Identification of a genetic variant for joint damage progression in autoantibody-positive rheumatoid arthritis. Ann Rheum Dis. 2014;73: 2038–2046. pmid:23956247
  11. 11. de Rooy DPC, Tsonaka R, Andersson MLE, Forslind K, Zhernakova A, Frank-Bertoncelj M, et al. Genetic Factors for the Severity of ACPA-negative Rheumatoid Arthritis in 2 Cohorts of Early Disease: A Genome-wide Study. J Rheumatol. 2015;42: 1383–1391. pmid:26077402
  12. 12. Scott IC, Rijsdijk F, Walker J, Quist J, Spain SL, Tan R, et al. Do Genetic Susceptibility Variants Associate with Disease Severity in Early Active Rheumatoid Arthritis?. J Rheumatol. 2015;42: 1131–1140. pmid:25979711
  13. 13. Teare MD, Knevel R, Morgan MD, Kleszcz A, Emery P, Moore DJ, et al. Allele-dose association of the C5orf30 rs26232 variant with joint damage in rheumatoid arthritis. Arthritis Rheum. 2013;65: 2555–2561. pmid:23817893
  14. 14. Iannaccone CK, Lee YC, Cui J, Frits ML, Glass RJ, Plenge RM, et al. Using genetic and clinical data to understand response to disease-modifying anti-rheumatic drug therapy: data from the Brigham and Women’s Hospital Rheumatoid Arthritis Sequential Study. Rheumatology (Oxford). 2011;50: 40–46. pmid:20847201
  15. 15. van der Linden MPM, Feitsma AL, le Cessie S, Kern M, Olsson LM, Raychaudhuri S, et al. Association of a single-nucleotide polymorphism in CD40 with the rate of joint destruction in rheumatoid arthritis. Arthritis Rheum. 2009;60: 2242–2247. pmid:19644859
  16. 16. Traylor M, Curtis C, Patel H, Breen G, Hyuck Lee S, Xu X, et al. Genetic and environmental risk factors for rheumatoid arthritis in a UK African ancestry population: the GENRA case-control study. Rheumatology (Oxford). 2017;56: 1282–1292. pmid:28407095
  17. 17. Steer S, Lad B, Grumley JA, Kingsley GH, Fisher SA. Association of R602W in a protein tyrosine phosphatase gene with a high risk of rheumatoid arthritis in a British population: evidence for an early onset/disease severity effect. Arthritis Rheum. 2005;52: 358–360. pmid:15641088
  18. 18. Scott DL, Houssien DA, Laasonen L. Proposed modification to Larsen’s scoring methods for hand and wrist radiographs. Br J Rheumatol. 1995;34: 56. pmid:7881840
  19. 19. van der Heijde D, Dankert T, Nieman F, Rau R, Boers M. Reliability and sensitivity to change of a simplification of the Sharp/van der Heijde radiological assessment in rheumatoid arthritis. Rheumatology (Oxford). 1999;38: 941–947. pmid:10534543
  20. 20. Pincus T, Larsen A, Brooks RH, Kaye J, Nance EP, Callahan LF. Comparison of 3 quantitative measures of hand radiographs in patients with rheumatoid arthritis: Steinbrocker stage, Kaye modified Sharp score, and Larsen score. J Rheumatol. 1997;24: 2106–2112. pmid:9375867
  21. 21. Plant MJ, Jones PW, Saklatvala J, Ollier WE, Dawes PT. Patterns of radiological progression in early rheumatoid arthritis: results of an 8 year prospective study. J Rheumatol. 1998;25: 417–426. pmid:9517757
  22. 22. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006;38: 904–909. pmid:16862161
  23. 23. Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010;26: 2190–2191. pmid:20616382
  24. 24. Lee CH, Cook S, Lee JS, Han B. Comparison of Two Meta-Analysis Methods: Inverse-Variance-Weighted Average and Weighted Sum of Z-Scores. Genomics Inform. 2016;14: 173–180. pmid:28154508
  25. 25. Eyre S, Bowes J, Diogo D, Lee A, Barton A, Martin P, et al. High-density genetic mapping identifies new susceptibility loci for rheumatoid arthritis. Barton A Morgan A, Wilson G, Hyrich K, Moitra RK, Prouse PJ, Shawe DJ, Crisp AJ, Gaston JS, Hall FC, Hazleman BL, Jenner JR, Lillicrap MS, Ostor A, Silverman B, Speed C, Bruce IN, Hyrich K, Ho P, Gorodkin R, Armstrong D, Chuck AJ, Hailwood S, Kumar N, Ba IJ, editor. Nat Genet. 2012;44: 1336–1340. pmid:23143596
  26. 26. Yang J, Wray NR, Visscher PM. Comparing apples and oranges: equating the power of case-control and quantitative trait association studies. Genet Epidemiol. 2010;34: 254–257. pmid:19918758
  27. 27. Okada Y, Wu D, Trynka G, Raj T, Terao C, Ikari K, et al. Genetics of rheumatoid arthritis contributes to biology and drug discovery. Eyre S 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, Andrew BJ, editor. Nature. 2014;506: 376–381. pmid:24390342
  28. 28. Viatte S, Plant D, Han B, Fu B, Yarwood A, Thomson W, et al. Association of HLA-DRB1 haplotypes with rheumatoid arthritis severity, mortality, and treatment response. JAMA. 2015;313: 1645–1656. pmid:25919528