Marieke J. H. Coenen

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.

Genetics of rheumatoid arthritis contributes to biology and drug discovery

A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological data sets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA). Here we performed a genome-wide association study meta-analysis in a total of >100,000 subjects of European and Asian ancestries (29,880 RA cases and 73,758 controls), by evaluating ∼10 million single-nucleotide polymorphisms. We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 101 (refs 2, 3, 4). We devised an in silico pipeline using established bioinformatics methods based on functional annotation, cis-acting expression quantitative trait loci and pathway analyses—as well as novel methods based on genetic overlap with human primary immunodeficiency, haematological cancer somatic mutations and knockout mouse phenotypes—to identify 98 biological candidate genes at these 101 risk loci. We demonstrate that these genes are the targets of approved therapies for RA, and further suggest that drugs approved for other indications may be repurposed for the treatment of RA. Together, this comprehensive genetic study sheds light on fundamental genes, pathways and cell types that contribute to RA pathogenesis, and provides empirical evidence that the genetics of RA can provide important information for drug discovery.

Genome-wide association study of systemic sclerosis identifies CD247 as a new susceptibility locus

Systemic sclerosis (SSc) is an autoimmune disease characterized by fibrosis of the skin and internal organs that leads to profound disability and premature death. To identify new SSc susceptibility loci, we conducted the first genome-wide association study in a population of European ancestry including a total of 2,296 individuals with SSc and 5,171 controls. Analysis of 279,621 autosomal SNPs followed by replication testing in an independent case-control set of European ancestry (2,753 individuals with SSc (cases) and 4,569 controls) identified a new susceptibility locus for systemic sclerosis at CD247 (1q22–23, rs2056626, P = 2.09 × 10−7 in the discovery samples, P = 3.39 × 10−9 in the combined analysis). Additionally, we confirm and firmly establish the role of the MHC (P = 2.31 × 10−18), IRF5 (P = 1.86 × 10−13) and STAT4 (P = 3.37 × 10−9) gene regions as SSc genetic risk factors.

Risk HLA-DQA1 and PLA(2)R1 alleles in idiopathic membranous nephropathy.

BACKGROUND Idiopathic membranous nephropathy is a major cause of the nephrotic syndrome in adults, but its etiologic basis is not fully understood. We investigated the genetic basis of biopsy-proven cases of idiopathic membranous nephropathy in a white population. METHODS We performed independent genomewide association studies of single-nucleotide polymorphisms (SNPs) in patients with idiopathic membranous nephropathy from three populations of white ancestry (75 French, 146 Dutch, and 335 British patients). The patients were compared with racially matched control subjects; population stratification and quality controls were carried out according to standard criteria. Associations were calculated by means of a chi-square basic allele test; the threshold for significance was adjusted for multiple comparisons (with the Bonferroni method). RESULTS In a joint analysis of data from the 556 patients studied (398 men), we identified significant alleles at two genomic loci associated with idiopathic membranous nephropathy. Chromosome 2q24 contains the gene encoding M-type phospholipase A(2) receptor (PLA(2)R1) (SNP rs4664308, P=8.6×10(-29)), previously shown to be the target of an autoimmune response. Chromosome 6p21 contains the gene encoding HLA complex class II HLA-DQ alpha chain 1 (HLA-DQA1) (SNP rs2187668, P=8.0×10(-93)). The association with HLA-DQA1 was significant in all three populations (P=1.8×10(-9), P=5.6×10(-27), and P=5.2×10(-36) in the French, Dutch, and British groups, respectively). The odds ratio for idiopathic membranous nephropathy with homozygosity for both risk alleles was 78.5 (95% confidence interval, 34.6 to 178.2). CONCLUSIONS An HLA-DQA1 allele on chromosome 6p21 is most closely associated with idiopathic membranous nephropathy in persons of white ancestry. This allele may facilitate an autoimmune response against targets such as variants of PLA2R1. Our findings suggest a basis for understanding this disease and illuminate how adaptive immunity is regulated by HLA.

Meta-Analysis of Genome-Wide Association Studies in Celiac Disease and Rheumatoid Arthritis Identifies Fourteen Non-HLA Shared Loci

Epidemiology and candidate gene studies indicate a shared genetic basis for celiac disease (CD) and rheumatoid arthritis (RA), but the extent of this sharing has not been systematically explored. Previous studies demonstrate that 6 of the established non-HLA CD and RA risk loci (out of 26 loci for each disease) are shared between both diseases. We hypothesized that there are additional shared risk alleles and that combining genome-wide association study (GWAS) data from each disease would increase power to identify these shared risk alleles. We performed a meta-analysis of two published GWAS on CD (4,533 cases and 10,750 controls) and RA (5,539 cases and 17,231 controls). After genotyping the top associated SNPs in 2,169 CD cases and 2,255 controls, and 2,845 RA cases and 4,944 controls, 8 additional SNPs demonstrated P<5×10−8 in a combined analysis of all 50,266 samples, including four SNPs that have not been previously confirmed in either disease: rs10892279 near the DDX6 gene (Pcombined = 1.2×10−12), rs864537 near CD247 (Pcombined = 2.2×10−11), rs2298428 near UBE2L3 (Pcombined = 2.5×10−10), and rs11203203 near UBASH3A (Pcombined = 1.1×10−8). We also confirmed that 4 gene loci previously established in either CD or RA are associated with the other autoimmune disease at combined P<5×10−8 (SH2B3, 8q24, STAT4, and TRAF1-C5). From the 14 shared gene loci, 7 SNPs showed a genome-wide significant effect on expression of one or more transcripts in the linkage disequilibrium (LD) block around the SNP. These associations implicate antigen presentation and T-cell activation as a shared mechanism of disease pathogenesis and underscore the utility of cross-disease meta-analysis for identification of genetic risk factors with pleiotropic effects between two clinically distinct diseases.

Novel association in chromosome 4q27 region with rheumatoid arthritis and confirmation of type 1 diabetes point to a general risk locus for autoimmune diseases.

Recently, association of celiac disease with common single-nucleotide polymorphism (SNP) variants in an extensive linkage-disequilibrium block of 480 kb containing the KIAA1109, Tenr, IL2, and IL21 genes has been demonstrated in three independent populations (rs6822844P combined=1.3 x 10(-14)). The KIAA1109/Tenr/IL2/IL21 block corresponds to the Idd3 locus in the nonobese diabetic mouse model of type 1 diabetes (T1D). This block was recently found to be associated with T1D in a genomewide association study, although this finding lacks unequivocal confirmation. We therefore aimed to investigate whether the KIAA1109/Tenr/IL2/IL21 region is involved in susceptibility to multiple autoimmune diseases. We tested SNP rs6822844 for association with disease in 350 T1D-affected and 1,047 rheumatoid arthritis (RA)-affected Dutch patients and in 929 controls. We replicated the association with T1D (P=.0006; OR 0.64 [95% CI 0.50-0.83]), and revealed a similar novel association with RA (P=.0002; OR 0.72 [95% CI 0.61-0.86]). Our results replicate and extend the association found in the KIAA1109/Tenr/IL2/IL21 gene region with autoimmune diseases, implying that this locus is a general risk factor for multiple autoimmune diseases.

Genome-Wide Association Scan Identifies Candidate Polymorphisms Associated with Differential Response to Anti-TNF Treatment in Rheumatoid Arthritis

The prediction of response (or non-response) to anti-TNF treatment for rheumatoid arthritis (RA) is a pressing clinical problem. We conducted a genome-wide association study using the Illumina HapMap300 SNP chip on 89 RA patients prospectively followed after beginning anti-TNF therapy as part of Autoimmune Biomarkers Collaborative Network (ABCoN [Autoimmune Bio-markers Collaborative Network]) patient cohort. Response to therapy was determined by the change in Disease Activity Score (DAS28) observed after 14 wks. We used a two-part analysis that treated the change in DAS28 as a continuous trait and then incorporated it into a dichotomous trait of “good responder” and “nonresponder” by European League Against Rheumatism (EULAR) criteria.We corrected for multiple tests by permutation, and adjusted for potential population stratification using EIGENSTRAT. Multiple single nucleotide polymorphism (SNP) markers showed significant associations near or within loci including: the v-maf musculoaponeurotic fibrosarcoma oncogene homolog B (MAFB) gene on chromosome 20; the type I interferon gene IFNk on chromosome 9; and in a locus on chromosome 7 that includes the paraoxonase I (PON1) gene. An SNP in the IL10 promoter (rs1800896) that was previously reported as associated with anti-TNF response was weakly associated with response in this cohort. Replications of these results in independent and larger data sets clearly are required. We provide a reference list of candidate SNPs (P < 0.01) that can be investigated in future pharmacogenomic studies.

Identification of novel genetic markers associated with clinical phenotypes of systemic sclerosis through a genome-wide association strategy

The aim of this study was to determine, through a genome-wide association study (GWAS), the genetic components contributing to different clinical sub-phenotypes of systemic sclerosis (SSc). We considered limited (lcSSc) and diffuse (dcSSc) cutaneous involvement, and the relationships with presence of the SSc-specific auto-antibodies, anti-centromere (ACA), and anti-topoisomerase I (ATA). Four GWAS cohorts, comprising 2,296 SSc patients and 5,171 healthy controls, were meta-analyzed looking for associations in the selected subgroups. Eighteen polymorphisms were further tested in nine independent cohorts comprising an additional 3,175 SSc patients and 4,971 controls. Conditional analysis for associated SNPs in the HLA region was performed to explore their independent association in antibody subgroups. Overall analysis showed that non-HLA polymorphism rs11642873 in IRF8 gene to be associated at GWAS level with lcSSc (P = 2.32×10−12, OR = 0.75). Also, rs12540874 in GRB10 gene (P = 1.27 × 10−6, OR = 1.15) and rs11047102 in SOX5 gene (P = 1.39×10−7, OR = 1.36) showed a suggestive association with lcSSc and ACA subgroups respectively. In the HLA region, we observed highly associated allelic combinations in the HLA-DQB1 locus with ACA (P = 1.79×10−61, OR = 2.48), in the HLA-DPA1/B1 loci with ATA (P = 4.57×10−76, OR = 8.84), and in NOTCH4 with ACA P = 8.84×10−21, OR = 0.55) and ATA (P = 1.14×10−8, OR = 0.54). We have identified three new non-HLA genes (IRF8, GRB10, and SOX5) associated with SSc clinical and auto-antibody subgroups. Within the HLA region, HLA-DQB1, HLA-DPA1/B1, and NOTCH4 associations with SSc are likely confined to specific auto-antibodies. These data emphasize the differential genetic components of subphenotypes of SSc.

Common and different genetic background for rheumatoid arthritis and coeliac disease

Recent genome-wide association studies (GWAS) have revealed genetic risk factors in autoimmune and inflammatory disorders. Several of the associated genes and underlying pathways are shared by various autoimmune diseases. Rheumatoid arthritis (RA) and coeliac disease (CD) are two autoimmune disorders which have commonalities in their pathogenesis. We aimed to replicate known RA loci in a Dutch RA population, and to investigate whether the effect of known RA and CD risk factors generalize across the two diseases. We selected all loci associated to either RA or CD in a GWAS and confirmed in an independent cohort, with a combined P-value cut-off P < 5 x 10(-6). We genotyped 11 RA and 11 CD loci in 1368 RA patients, 795 CD patients and 1683 Dutch controls. We combined our results in a meta-analysis with UK GWAS on RA (1860 cases; 2938 controls) and CD (767 cases; 1422 controls). In the Dutch RA cohort, the PTPN22 and IL2/IL21 variants showed convincing association (P = 3.4 x 10(-12) and P = 2.8 x 10(-4), respectively). Association of RA with the known CD risk variant in the SH2B3 was also observed, predominantly in the subgroup of rheumatoid factor-positive RA patients (P = 0.0055). In a meta-analysis of Dutch and UK data sets, shared association with six loci (TNFAIP3, IL2/IL21, SH2B3, LPP, MMEL1/TNFRSF14 and PFKFB3/PRKCQ) was observed in both RA and CD cohorts. We confirmed two known loci and identified four novel ones for shared CD-RA genetic risk. Most of the shared loci further emphasize a role for adaptive and innate immunity in these diseases.

Genome-wide association study and gene expression analysis identifies CD84 as a predictor of response to etanercept therapy in rheumatoid arthritis.

Anti-tumor necrosis factor alpha (anti-TNF) biologic therapy is a widely used treatment for rheumatoid arthritis (RA). It is unknown why some RA patients fail to respond adequately to anti-TNF therapy, which limits the development of clinical biomarkers to predict response or new drugs to target refractory cases. To understand the biological basis of response to anti-TNF therapy, we conducted a genome-wide association study (GWAS) meta-analysis of more than 2 million common variants in 2,706 RA patients from 13 different collections. Patients were treated with one of three anti-TNF medications: etanercept (n = 733), infliximab (n = 894), or adalimumab (n = 1,071). We identified a SNP (rs6427528) at the 1q23 locus that was associated with change in disease activity score (ΔDAS) in the etanercept subset of patients (P = 8×10−8), but not in the infliximab or adalimumab subsets (P>0.05). The SNP is predicted to disrupt transcription factor binding site motifs in the 3′ UTR of an immune-related gene, CD84, and the allele associated with better response to etanercept was associated with higher CD84 gene expression in peripheral blood mononuclear cells (P = 1×10−11 in 228 non-RA patients and P = 0.004 in 132 RA patients). Consistent with the genetic findings, higher CD84 gene expression correlated with lower cross-sectional DAS (P = 0.02, n = 210) and showed a non-significant trend for better ΔDAS in a subset of RA patients with gene expression data (n = 31, etanercept-treated). A small, multi-ethnic replication showed a non-significant trend towards an association among etanercept-treated RA patients of Portuguese ancestry (n = 139, P = 0.4), but no association among patients of Japanese ancestry (n = 151, P = 0.8). Our study demonstrates that an allele associated with response to etanercept therapy is also associated with CD84 gene expression, and further that CD84 expression correlates with disease activity. These findings support a model in which CD84 genotypes and/or expression may serve as a useful biomarker for response to etanercept treatment in RA patients of European ancestry.