ERAP1 association with ankylosing spondylitis is attributable to common genotypes rather than rare haplotype combinations.

We investigated the proposal that ankylosing spondylitis (AS) is associated with unusual ERAP1 genotypes. ERAP1 haplotypes were constructed for 213 AS cases and 46 rheumatoid arthritis controls using family data. Haplotypes were generated from five common ERAP1 single nucleotide polymorphisms (SNPs)-rs2287987 (M349V), rs30187 (K528R), rs10050860 (D575N), rs17482078 (R725Q), and rs27044 (Q730E). Haplotype frequencies were compared using Fisher's exact test. ERAP1 haplotypes imputed from the International Genetics of AS Consortium (IGAS) Immunochip study were also studied. In the family study, we identified only four common ERAP1 haplotypes ("VRNQE," "MKDRQ," "MRDRE," and "MKDRE") in both AS cases and controls apart from two rare (<0.5%) previously unreported haplotypes. There were no examples of the unusual ERAP1 haplotype combination ("*001/*005") previously reported by others in 53% of AS cases. As expected, K528-bearing haplotypes were increased in the AS family study (AS 43% vs. control 35%), due particularly to an increase in the MKDRQ haplotype (AS 35% vs. control 25%, P = 0.01). This trend was replicated in the imputed Immunochip data for the two K528-bearing haplotypes MKDRQ (AS 33% vs. controls 27%, P = 1.2 × 10-24) and MKDRE (AS 8% vs. controls 7%, P = 0.004). The ERAP1 association with AS is therefore predominantly attributable to common ERAP1 haplotypes and haplotype combinations.

Investigating the genetic association between ERAP1 and ankylosing spondylitis.

A strong association between ERAP1 and ankylosing spondylitis (AS) was recently identified by the Wellcome Trust Case Control Consortium and the Australo-Anglo-American Spondylitis Consortium (WTCCC-TASC) study. ERAP1 is highly polymorphic with strong linkage disequilibrium evident across the gene. We therefore conducted a series of experiments to try to identify the primary genetic association(s) with ERAP1. We replicated the original associations in an independent set of 730 patients and 1021 controls, resequenced ERAP1 to define the full extent of coding polymorphisms and tested all variants in additional association studies. The genetic association with ERAP1 was independently confirmed; the strongest association was with rs30187 in the replication set (P = 3.4 x 10(-3)). When the data were combined with the original WTCCC-TASC study the strongest association was with rs27044 (P = 1.1 x 10(-9)). We identified 33 sequence polymorphisms in ERAP1, including three novel and eight known non-synonymous polymorphisms. We report several new associations between AS and polymorphisms distributed across ERAP1 from the extended case-control study, the most significant of which was with rs27434 (P = 4.7 x 10(-7)). Regression analysis failed to identify a primary association clearly; we therefore used data from HapMap to impute genotypes for an additional 205 non-coding SNPs located within and adjacent to ERAP1. A number of highly significant associations (P < 5 x 10(-9)) were identified in regulatory sequences which are good candidates for causing susceptibility to AS, possibly by regulating ERAP1 expression.

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.

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.

Mutations in Known Monogenic High Bone Mass Loci Only Explain a Small Proportion of High Bone Mass Cases.

High bone mass (HBM) can be an incidental clinical finding; however, monogenic HBM disorders (eg, LRP5 or SOST mutations) are rare. We aimed to determine to what extent HBM is explained by mutations in known HBM genes. A total of 258 unrelated HBM cases were identified from a review of 335,115 DXA scans from 13 UK centers. Cases were assessed clinically and underwent sequencing of known anabolic HBM loci: LRP5 (exons 2, 3, 4), LRP4 (exons 25, 26), SOST (exons 1, 2, and the van Buchem's disease [VBD] 52-kb intronic deletion 3'). Family members were assessed for HBM segregation with identified variants. Three-dimensional protein models were constructed for identified variants. Two novel missense LRP5 HBM mutations ([c.518C>T; p.Thr173Met], [c.796C>T; p.Arg266Cys]) were identified, plus three previously reported missense LRP5 mutations ([c.593A>G; p.Asn198Ser], [c.724G>A; p.Ala242Thr], [c.266A>G; p.Gln89Arg]), associated with HBM in 11 adults from seven families. Individuals with LRP5 HBM (∼prevalence 5/100,000) displayed a variable phenotype of skeletal dysplasia with increased trabecular BMD and cortical thickness on HRpQCT, and gynoid fat mass accumulation on DXA, compared with both non-LRP5 HBM and controls. One mostly asymptomatic woman carried a novel heterozygous nonsense SOST mutation (c.530C>A; p.Ser177X) predicted to prematurely truncate sclerostin. Protein modeling suggests the severity of the LRP5-HBM phenotype corresponds to the degree of protein disruption and the consequent effect on SOST-LRP5 binding. We predict p.Asn198Ser and p.Ala242Thr directly disrupt SOST binding; both correspond to severe HBM phenotypes (BMD Z-scores +3.1 to +12.2, inability to float). Less disruptive structural alterations predicted from p.Arg266Cys, p.Thr173Met, and p.Gln89Arg were associated with less severe phenotypes (Z-scores +2.4 to +6.2, ability to float). In conclusion, although mutations in known HBM loci may be asymptomatic, they only account for a very small proportion (∼3%) of HBM individuals, suggesting the great majority are explained by either unknown monogenic causes or polygenic inheritance.

Genome-wide association study of ankylosing spondylitis identifies non-MHC susceptibility loci.

To identify susceptibility loci for ankylosing spondylitis, we undertook a genome-wide association study in 2,053 unrelated ankylosing spondylitis cases among people of European descent and 5,140 ethnically matched controls, with replication in an independent cohort of 898 ankylosing spondylitis cases and 1,518 controls. Cases were genotyped with Illumina HumHap370 genotyping chips. In addition to strong association with the major histocompatibility complex (MHC; P < 10(-800)), we found association with SNPs in two gene deserts at 2p15 (rs10865331; combined P = 1.9 x 10(-19)) and 21q22 (rs2242944; P = 8.3 x 10(-20)), as well as in the genes ANTXR2 (rs4333130; P = 9.3 x 10(-8)) and IL1R2 (rs2310173; P = 4.8 x 10(-7)). We also replicated previously reported associations at IL23R (rs11209026; P = 9.1 x 10(-14)) and ERAP1 (rs27434; P = 5.3 x 10(-12)). This study reports four genetic loci associated with ankylosing spondylitis risk and identifies a major role for the interleukin (IL)-23 and IL-1 cytokine pathways in disease susceptibility.