Ann W. Morgan

Genome-wide association study meta-analysis identifies seven new rheumatoid arthritis risk loci

To identify new genetic risk factors for rheumatoid arthritis, we conducted a genome-wide association study meta-analysis of 5,539 autoantibody-positive individuals with rheumatoid arthritis (cases) and 20,169 controls of European descent, followed by replication in an independent set of 6,768 rheumatoid arthritis cases and 8,806 controls. Of 34 SNPs selected for replication, 7 new rheumatoid arthritis risk alleles were identified at genome-wide significance (P < 5 × 10−8) in an analysis of all 41,282 samples. The associated SNPs are near genes of known immune function, including IL6ST, SPRED2, RBPJ, CCR6, IRF5 and PXK. We also refined associations at two established rheumatoid arthritis risk loci (IL2RA and CCL21) and confirmed the association at AFF3. These new associations bring the total number of confirmed rheumatoid arthritis risk loci to 31 among individuals of European ancestry. An additional 11 SNPs replicated at P < 0.05, many of which are validated autoimmune risk alleles, suggesting that most represent genuine rheumatoid arthritis risk alleles.

Interaction between ERAP1 and HLA-B27 in ankylosing spondylitis implicates peptide handling in the mechanism for HLA-B27 in disease susceptibility

Ankylosing spondylitis is a common form of inflammatory arthritis predominantly affecting the spine and pelvis that occurs in approximately 5 out of 1,000 adults of European descent. Here we report the identification of three variants in the RUNX3, LTBR-TNFRSF1A and IL12B regions convincingly associated with ankylosing spondylitis (P < 5 × 10−8 in the combined discovery and replication datasets) and a further four loci at PTGER4, TBKBP1, ANTXR2 and CARD9 that show strong association across all our datasets (P < 5 × 10−6 overall, with support in each of the three datasets studied). We also show that polymorphisms of ERAP1, which encodes an endoplasmic reticulum aminopeptidase involved in peptide trimming before HLA class I presentation, only affect ankylosing spondylitis risk in HLA-B27–positive individuals. These findings provide strong evidence that HLA-B27 operates in ankylosing spondylitis through a mechanism involving aberrant processing of antigenic peptides.

Rheumatoid arthritis association at 6q23

The Wellcome Trust Case Control Consortium (WTCCC) identified nine single SNPs putatively associated with rheumatoid arthritis at P = 1 × 10−5 − 5 × 10−7 in a genome-wide association screen. One, rs6920220, was unequivocally replicated (trend P = 1.1 × 10−8) in a validation study, as described here. This SNP maps to 6q23, between the genes oligodendrocyte lineage transcription factor 3 (OLIG3) and tumor necrosis factor-α–induced protein 3 (TNFAIP3).

Genetic variants at CD28, PRDM1, and CD2/CD58 are associated with rheumatoid arthritis risk

To discover new rheumatoid arthritis (RA) risk loci, we systematically examined 370 SNPs from 179 independent loci with P < 0.001 in a published meta-analysis of RA genome-wide association studies (GWAS) of 3,393 cases and 12,462 controls. We used Gene Relationships Across Implicated Loci (GRAIL), a computational method that applies statistical text mining to PubMed abstracts, to score these 179 loci for functional relationships to genes in 16 established RA disease loci. We identified 22 loci with a significant degree of functional connectivity. We genotyped 22 representative SNPs in an independent set of 7,957 cases and 11,958 matched controls. Three were convincingly validated: CD2-CD58 (rs11586238, P = 1 × 10−6 replication, P = 1 × 10−9 overall), CD28 (rs1980422, P = 5 × 10−6 replication, P = 1 × 10−9 overall) and PRDM1 (rs548234, P = 1 × 10−5 replication, P = 2 × 10−8 overall). An additional four were replicated (P < 0.0023): TAGAP (rs394581, P = 0.0002 replication, P = 4 × 10−7 overall), PTPRC (rs10919563, P = 0.0003 replication, P = 7 × 10−7 overall), TRAF6-RAG1 (rs540386, P = 0.0008 replication, P = 4 × 10−6 overall) and FCGR2A (rs12746613, P = 0.0022 replication, P = 2 × 10−5 overall). Many of these loci are also associated to other immunologic diseases.

Copy number of FCGR3B, which is associated with systemic lupus erythematosus, correlates with protein expression and immune complex uptake.

Copy number (CN) variation (CNV) has been shown to be common in regions of the genome coding for immune-related genes, and thus impacts upon polygenic autoimmunity. Low CN of FCGR3B has recently been associated with systemic lupus erythematosus (SLE). FcγRIIIb is a glycosylphosphatidylinositol-linked, low affinity receptor for IgG found predominantly on human neutrophils. We present novel data demonstrating that both in a family with FcγRIIIb-deficiency and in the normal population, FCGR3B CNV correlates with protein expression, with neutrophil uptake of and adherence to immune complexes, and with soluble serum FcγRIIIb. Reduced FcγRIIIb expression is thus likely to contribute to the impaired clearance of immune complexes, which is a feature of SLE, explaining the association between low FCGR3B CNV and SLE that we have confirmed in a Caucasian population. In contrast, antineutrophil cytoplasmic antibody–associated systemic vasculitis (AASV), a disease not associated with immune complex deposition, is associated with high FCGR3B CN. Thus, we define a role for FCGR3B CNV in immune complex clearance, a function that may explain why low FCGR3B CNV is associated with SLE, but not AASV. This is the first report of an association between disease-related gene CNV and variation in protein expression and function that may contribute to autoimmune disease susceptibility.

The Requirement for DNAM-1, NKG2D, and NKp46 in the Natural Killer Cell-Mediated Killing of Myeloma Cells

Recent evidence suggests a role for natural killer (NK) cells in the control of multiple myeloma. We show that expression of the NK cell receptor DNAM-1 (CD226) is reduced on CD56(dim) NK cells from myeloma patients with active disease compared with patients in remission and healthy controls. This suggested that this receptor might play a role in NK-myeloma interactions. The DNAM-1 ligands Nectin-2 (CD112) and the poliovirus receptor (PVR; CD155) were expressed by most patient myeloma samples analyzed. NK killing of patient-derived myelomas expressing PVR and/or Nectin-2 was DNAM-1 dependent, revealing a functional role for DNAM-1 in myeloma cell killing. In myeloma cell lines, cell surface expression of PVR was associated with low levels of NKG2D ligands, whereas cells expressing high levels of NKG2D ligands did not express PVR protein or mRNA. Furthermore, NK cell-mediated killing of myeloma cell lines was dependent on either DNAM-1 or NKG2D but not both molecules. In contrast, the natural cytotoxicity receptor NKp46 was required for the killing of all myeloma cell lines analyzed. Thus, DNAM-1 is important in the NK cell-mediated killing of myeloma cells expressing the cognate ligands. The importance of NKp46, NKG2D, and DNAM-1 in myeloma killing mirrors the differential expression of NK cell ligands by myeloma cells, reflecting immune selection during myeloma disease progression.

Association of rheumatoid factor and anti-cyclic citrullinated peptide positivity, but not carriage of shared epitope or PTPN22 susceptibility variants, with anti-tumour necrosis factor response in rheumatoid arthritis

Objective: To determine whether rheumatoid factor (RF), anti-cyclic citrullinated peptide (CCP) antibodies, or carriage of shared epitope (SE) and PTPN22 genetic susceptibility variants predict response to therapy in patients with rheumatoid arthritis (RA) treated with anti-tumour necrosis factor (TNF) agents. Methods: UK-wide multicentre collaborations were established to recruit a large cohort of patients treated with anti-TNF drugs for RA. Serum RF, anti-CCP antibody and SE status were determined using commercially available kits. PTPN22 R620W genotyping was performed by Sequenom MassArray. Linear regression analyses were performed to investigate the role of these four factors in predicting response to treatment by 6 months, defined as the absolute change in 28-joint Disease Activity Score (DAS28). Results: Of the 642 patients analysed, 46% received infliximab, 43% etanercept and 11% adalimumab. In all, 89% and 82% of patients were RF and anti-CCP positive, respectively. Patients that were RF negative had a 0.48 (95% CI 0.08 to 0.87) greater mean improvement in DAS28 compared to patients that were RF positive. A better response was also seen among patients that were anti-CCP negative. No association was demonstrated between drug response and SE or PTPN22 620W carriage. Conclusion: The presence of RF or anti-CCP antibodies was associated with a reduced response to anti-TNF drugs. However, these antibodies only account for a small proportion of the variance in treatment response. It is likely that genetic factors will contribute to treatment response, but these do not include the well established RA susceptibility loci, SE and PTPN22.

Rheumatoid arthritis susceptibility loci at chromosomes 10p15, 12q13 and 22q13.

The WTCCC study identified 49 SNPs putatively associated with rheumatoid arthritis at P = 1 × 10−4 − 1 × 10−5 (tier 3). Here we show that three of these SNPs, mapping to chromosome 10p15 (rs4750316), 12q13 (rs1678542) and 22q13 (rs3218253), are also associated (trend P = 4 × 10−5, P = 4 × 10−4 and P = 4 × 10−4, respectively) in a validation study of 4,106 individuals with rheumatoid arthritis and an expanded reference group of 11,238 subjects, confirming them as true susceptibility loci in individuals of European ancestry.

Re-evaluation of putative rheumatoid arthritis susceptibility genes in the post-genome wide association study era and hypothesis of a key pathway underlying susceptibility

Rheumatoid arthritis (RA) is an archetypal, common, complex autoimmune disease with both genetic and environmental contributions to disease aetiology. Two novel RA susceptibility loci have been reported from recent genome-wide and candidate gene association studies. We, therefore, investigated the evidence for association of the STAT4 and TRAF1/C5 loci with RA using imputed data from the Wellcome Trust Case Control Consortium (WTCCC). No evidence for association of variants mapping to the TRAF1/C5 gene was detected in the 1860 RA cases and 2930 control samples tested in that study. Variants mapping to the STAT4 gene did show evidence for association (rs7574865, P = 0.04). Given the association of the TRAF1/C5 locus in two previous large case–control series from populations of European descent and the evidence for association of the STAT4 locus in the WTCCC study, single nucleotide polymorphisms mapping to these loci were tested for association with RA in an independent UK series comprising DNA from >3000 cases with disease and >3000 controls and a combined analysis including the WTCCC data was undertaken. We confirm association of the STAT4 and the TRAF1/C5 loci with RA bringing to 5 the number of confirmed susceptibility loci. The effect sizes are less than those reported previously but are likely to be a more accurate reflection of the true effect size given the larger size of the cohort investigated in the current study.

Genome-Wide Association Study of Genetic Predictors of Anti–Tumor Necrosis Factor Treatment Efficacy in Rheumatoid Arthritis Identifies Associations with Polymorphisms at Seven Loci

Objective Anti–tumor necrosis factor (anti-TNF) agents are successful therapies in rheumatoid arthritis (RA); however, inadequate response occurs in 30–40% of patients treated. Knowledge of the genetic factors that influence response may facilitate personalized therapy. The purpose of this study was to identify genetic predictors of response to anti-TNF therapy in RA and to validate our findings in independent cohorts. Methods Data from genome-wide association (GWA) studies were available from the Wellcome Trust Case Control Consortium for 566 anti-TNF–treated RA patients. Multivariate linear regression analysis of changes in the Disease Activity Score in 28 joints at 6 months was conducted at each single-nucleotide polymorphism (SNP) using an additive model. Associated markers (P < 10−3) were genotyped in 2 independent replication cohorts (n = 379 and n = 341), and a combined analysis was performed. Results Of 171 successfully genotyped markers demonstrating association with treatment response in the GWA data, 7 were corroborated in the combined analysis. The strongest effect was at rs17301249, mapping to the EYA4 gene on chromosome 6: the minor allele conferred improved response to treatment (coefficient −0.27, P = 5.67−05). The minor allele of rs1532269, mapping to the PDZD2 gene, was associated with a reduced treatment response (coefficient 0.20, P = 7.37−04). The remaining associated SNPs mapped to intergenic regions on chromosomes 1, 4, 11, and 12. Conclusion Using a genome-wide strategy, we have identified and validated the association of 7 genetic loci with response to anti-TNF treatment in RA. Additional confirmation of these findings in further cohorts will be required.