C Farrar

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.

Association between the interleukin 23 receptor and ankylosing spondylitis is confirmed by a new UK case–control study and meta-analysis of published series

OBJECTIVES It has been shown previously that IL-23R variants are associated with AS. We conducted an extended analysis in the UK population and a meta-analysis with the previously published studies, in order to refine these IL-23R associations with AS. METHODS The UK case-control study included 730 new cases and 1331 healthy controls. In the extended study, the 730 cases were combined with 1088 published cases. Allelic associations were analysed using contingency tables. In the meta-analysis, 3482 cases and 3150 controls from four different published studies and the new UK cases were combined. DerSimonian-Laird test was used to calculate random effects pooled odds ratios (ORs). RESULTS In the UK case-control study with new cases, four of the eight SNPs showed significant associations, whereas in the extended UK study, seven of the eight IL-23R SNPs showed significant associations (P < 0.05) with AS, maximal with rs11209032 (P < 10(-5), OR 1.3), when cases with IBD and/or psoriasis were excluded. The meta-analysis showed significant associations with all eight SNPs; the strongest associations were again seen not only with rs11209032 (P = 4.06 x 10(-9), OR approximately 1.2) but also with rs11209026 (P < 10(-10), OR approximately 0.6). CONCLUSIONS IL-23R polymorphisms are clearly associated with AS, but the primary causal association(s) is(are) still not established. These polymorphisms could contribute either increased or decreased susceptibility to AS; functional studies will be required for their full evaluation. Additionally, observed stronger associations with SNPs rs11209026 and rs11465804 upon exclusion of IBD and/or psoriasis cases may represent an independent association with AS.

Elucidating the chromosome 9 association with AS; CARD9 is a candidate gene

Ankylosing spondylitis (AS) is polygenic with contributions from the immunologically relevant genes HLA-B*27, ERAP1 and IL23R. A recent genome-wide association screen (GWAS) identified associations (P∼0.005) with the non-synonymous single-nucleotide polymorphisms (nsSNPs), rs4077515 and rs3812571, in caspase recruitment domain-containing protein 9 (CARD9) and small nuclear RNA-activating complex polypeptide 4 (SNAPC4) on chromosome 9q that had previously been linked to AS. We replicated these associations in a study of 730 AS patients compared with 2879 historic disease controls (rs4077515 P=0.0004, odds ratio (OR)=1.2, 95% confidence interval (CI)=1.1–1.4; rs3812571 P=0.0003, OR=1.2, 95% CI=1.1–1.4). Meta-analysis revealed strong associations of both SNPs with AS, rs4077515 P=0.000005, OR=1.2, 95% CI=1.1–1.3 and rs3812571 P=0.000006, OR=1.2, 95% CI=1.1–1.3. We then typed 1604 AS cases and 1020 controls for 13 tagging SNPs; 6 showed at least nominal association, 5 of which were in CARD9. We imputed genotypes for 13 additional SNPs but none was more strongly associated with AS than the tagging SNPs. Finally, interrogation of an mRNA expression database revealed that the SNPs most strongly associated with AS (or in strong linkage disequilibrium) were those most associated with CARD9 expression. CARD9 is a plausible candidate for AS given its central role in the innate immune response.

Analysis of killer immunoglobulin-like receptor genes in ankylosing spondylitis

Objectives: To assess the possible association of killer immunoglobulin-like receptor (KIR) genes, specifically KIR3DL1, KIR3DS1 and KIR3DL2, with ankylosing spondylitis (AS). Methods: 14 KIR genes were genotyped in 200 UK patients with AS and 405 healthy controls using multiplex polymerase chain reaction. Sequence-specific oligonucleotide probes were used to subtype 368 cases with AS and 366 controls for 12 KIR3DL2 alleles. Differences in KIR genotypes and KIR3DL2 allele frequencies were assessed using the χ2 test. Results: KIR3DL1 and KIR3DS1 gene frequencies were very similar in cases with AS and controls (odds ratio = 1.5, 95% confidence interval 0.8 to 3.0, and odds ratio = 1.02, 95% confidence interval 0.2 to 5.3, respectively). KIR3DL2 allele frequencies were not significantly different between cases with AS and controls. Conclusions: Neither the KIR gene content of particular KIR haplotypes nor KIR3DL2 polymorphisms contribute to AS.

The chromosome 16q region associated with ankylosing spondylitis includes the candidate gene tumour necrosis factor receptor type 1-associated death domain (TRADD)

Objective To replicate and refine the reported association of ankylosing spondylitis (AS) with two non-synonymous single nucleotide polymorphisms (nsSNPs) on chromosome 16q22.1. Methods Firstly, 730 independent UK patients with AS were genotyped for rs9939768 and rs6979 and allele frequencies were compared with 2879 previously typed historic disease controls. Secondly, the two data sets were combined in meta-analyses. Finally, 5 tagging SNPs, located between rs9939768 and rs6979, were analysed in 1604 cases and 1020 controls. Results The association of rs6979 with AS was replicated, p=0.03, OR=1.14 (95% CI 1.01 to 1.28), and a trend for association with rs9939768 detected, p=0.06, OR=1.25 (95% CI 0.99 to 1.57). Meta-analyses revealed association of both SNPs with AS, p=0.0008, OR=1.31 (95% CI 1.12 to 1.54) and p=0.0009, OR=1.15 (95% CI 1.06 to 1.23) for rs9939768 and rs6979, respectively. New associations with rs9033 and rs868213 (p=0.00002, OR=1.23 (95% CI 1.12 to 1.36) and p=0.00002 OR=1.45 (95% CI 1.22 to 1.72), respectively, were identified. Conclusions The region on chromosome 16 that has been replicated in the present work is interesting as the highly plausible candidate gene, tumour necrosis factor receptor type 1 (TNFR1)-associated death domain (TRADD), is located between rs9033 and rs868213. It will require additional work to identify the primary genetic association(s) with AS.

The histone demethylase JARID1A is associated with susceptibility to ankylosing spondylitis.

Associations with disease identified by genome-wide association studies (GWAS) must be replicated and refined to validate causative variants. In the Wellcome Trust Case Control Consortium (WTCCC) GWAS using 14 500 non-synonymous single nucleotide polymorphisms (nsSNP), rs11062385 (a nsSNP in JARID1A) showed nominal association with ankylosing spondylitis (AS) (P=0.0006, odds ratio (OR)=1.26, 95% confidence interval (95% CI)=1.1–1.4). To replicate and refine the association of JARID1A, rs11062385 was genotyped in 730 further cases and compared with allele frequencies in non-AS disease cohorts typed by WTCCC. We replicated the initial association (P=0.04, OR=1.16, 95% CI=1.01–1.34) and identified a strengthened association with AS in a meta-analysis of this new study combined with the original WTCCC study (P=0.0001, OR=1.21, 95% CI=1.10–1.33). We also genotyped nine further intronic tagging SNPs in JARID1A in 1604 AS cases and 1020 new control samples, but none was associated with AS. JARID1A or a locus in strong linkage disequilibrium with it is a positional candidate for susceptibility to AS.

A common functional variant of endoplasmic reticulum aminopeptidase 2 (ERAP2) that reduces major histocompatibility complex class I expression is not associated with ankylosing spondylitis

Sir, The strong genetic association between AS and HLA-B27 has defied explanation for nearly 40 years. However, the additional discovery of a strong association between AS and ERAP1 (endoplasmic reticulum aminopeptidase 1), a gene that almost certainly operates in the trimming of peptides for optimal binding to MHC class I molecules, has rekindled hopes of rapid advances in this field [1]. It has been suggested that another aminopeptidase, endoplasmic reticulum aminopeptidase 2 (ERAP2), may act in concert with ERAP1, trimming residues inefficiently removed by ERAP1 [2]. The association described recently between ERAP2 and Crohn’s disease, which shares many clinical and genetic overlaps with AS [3], suggests that ERAP2 is worthy of further study in AS. An experiment of nature allows us to do this relatively simply. ERAP2 has evolved under balancing selection, similar to the MHC, and includes a high-frequency variant (∼50%) that influences antigen presentation [4]. A single nucleotide polymorphism (SNP), rs2248374 (A to G), located within the 5′ canonical splice site of exon 10, results in an alternatively spliced ERAP2 mRNA that is degraded by nonsense-mediated decay (NMD). Homozygosity for the minor G allele (carried by ∼25% of the population) results in failure to express ERAP2 protein; in turn this genotype is also associated with reduced surface MHC Class I expression on human B cells [4]. Such a dramatic phenotypic effect provides a great opportunity for studying the potential role of ERAP2 in AS, analogous to using a knockout model. We therefore specifically tested for differences in frequency of rs2248374 in 470 sporadic AS cases and 420 healthy, ethnically matched blood donors to determine whether this functionally important variant is implicated in AS. All patients were Caucasians, of UK origin and fulfilled the modified New York criteria for AS. Eighty-four per cent of the AS patients were HLA-B27 positive. All patients gave informed consent and ethical approval was obtained (Multicentre Research Ethics Committee 98/5/23). Genotyping involved restriction fragment length polymorphism (RFLP) analysis of PCR products performed under standard conditions. Primers (forward 5′-GCATCCATGGCTAATGTGCR and reverse 5′-GTTGTGGGAAAGCCGAACTA) amplified a 370-bp product that, in the presence of the G allele, was digested to 214 and 156 bp products by HphI. Genotype and allele frequencies in AS cases and controls were compared using the Cochrane–Armitage test of trend and the chi-squared test, respectively. This study had 80% power under a log-additive model to detect an odds ratio (OR) of 1.3 with a population risk of AS of 0.04% and a risk allele frequency of 46% (frequency in our controls) at a significance level of 0.05. There were no significant differences between allele or genotype frequencies in AS cases and controls at rs2248374 (Table 1). Additionally, due to the potential functional significance of homozygosity at the minor G allele (i.e. no functional ERAP2 protein), we also compared the frequency of GG homozygotes in AS cases and controls, but these proved almost identical again (Table 1). Table 1 Genotypes, minor allele frequency (MAF), OR, 95% CI and P-value for the ERAP2 SNP rs2248374 The association between ERAP1 and AS is the second only to HLA-B27, but formal testing for additional association with the neighbouring gene, ERAP2, had not been previously undertaken. A study of familial AS that typed one marker within ERAP2 reported that an ERAP1/ERAP2 haplotype was over-represented in affected family members [5]. However, such studies are frequently confounded by linkage disequilibrium effects. If genetic variation at ERAP2 were to play a role in susceptibility to AS, then a dramatic loss-of-function variant, such as rs2248374, would appear to be a prime candidate for the ‘causal variant’ at this locus. However, we have now excluded any significant effect arising from this important functional ERAP2 variant. Nonetheless, the possibility of other associations between ERAP2 and AS has not been excluded by this study. This could be done only by systematic mapping of the gene with tagging SNPs in a very large sample after controlling for the significant linkage disequilibrium between ERAP1 and ERAP2 would this be achievable. Our study was insufficiently powered to do this.