Prospective meta-analysis of interleukin 1 gene complex polymorphisms confirms associations with ankylosing spondylitis

Objectives: The aim of the current study was to determine the contribution of interleukin (IL) 1 gene cluster polymorphisms previously implicated in susceptibility for ankylosing spondylitis (AS) to AS susceptibility in different populations worldwide. Methods: Nine polymorphisms in the IL1 gene cluster members IL1A (rs2856836, rs17561 and rs1894399), IL1B (rs16944), IL1F10 (rs3811058) and IL1RN (rs419598, the IL1RA VNTR, rs315952 and rs315951) were genotyped in 2675 AS cases and 2592 healthy controls recruited in 12 different centres in 10 countries. Association of variants with AS was tested by Mantel-Haenszel random effects analysis. Results: Strong association was observed with three single nucleotide polymorphisms (SNPs) in the IL1A gene (rs2856836, rs17561, rs1894399, p = 0.0036, 0.000019 and 0.0003, respectively). There was no evidence of significant heterogeneity of effects between centres, and no evidence of non-combinability of findings. The population attributable risk fraction of these variants in Caucasians is estimated at 4-6%. Conclusions: This study confirms that IL1A is associated with susceptibility to AS. Association of the other IL1 gene complex members could not be excluded in specific populations. Prospective meta-analysis is a useful tool in confirmation studies of genes associated with complex genetic disorders such as AS, providing sufficiently large sample sizes to produce robust findings often not achieved in smaller individual cohorts.

Genetic testing for haemochromatosis in patients with chondrocalcinosis

Hereditary haemochromatosis (HH) is the most common lethal monogenic human disease, affecting roughly 1 in 300 white northern Europeans. Homozygosity for the C282Y polymorphism within the HFE gene causes more than 80% of cases, with compound heterozygosity of the C282Y and H63D polymorphism also increasing susceptibility to disease. The aim of this study was to determine the frequency of the C282Y and H63D polymorphisms in the disease, and to assess the risk of HH in heterozygotes for the C282Y polymorphism. 128 patients were recruited because of either radiographic chondrocalcinosis (at least bicompartmental knee disease or joints other than the knee involved) or CPPD pseudogout. Genotyping of the HFE C282Y and H63D mutations was performed using PCR/SSP and genotypes for the C282Y polymorphism confirmed by PCR/RFLP. Historical white European control data were used for comparison. Two previously undiagnosed C282Y homozygotes (1.6%), and 16 C282Y heterozygotes (12.5%), including four (3.1%) C282Y/ H63D compound heterozygotes were identified. This represents a significant overrepresentation of C282Y homozygotes (relative risk 3.4, p-0.037), but the number of heterozygotes was not significantly increased. At a cost per test of £1 for each subject, screening all patients with chondrocalcinosis using the above ascertainment criteria costs only £64 for each case of haemochromatosis identified, clearly a highly cost effective test given the early mortality associated with untreated haemochromatosis. Routine screening for haemochromatosis in patients with appreciable chondrocatcinosis is recommended.

Author's reply to: Ankylosing spondylitis: evidence for a non-HLA-B*27 protective effect, Sjef Van Der Linden, Desieree Van Der Heijde

As mentioned in the letter by van der Linden and van der Heijde, Jurgen Braun’s excellent recent paper describing a survey of blood donors by questionnaire, clinical, and magnetic resonance imaging examinations revealed a prevalence of ankylosing spondylitis in B27 positive blood donors (6.4%)1-1 very similar to that reported by Gran et al(6.7%).1-2 It is probable that some of the differences in reported prevalence of ankylosing spondylitis by the various studies are because of methodological differences.

Identification of major loci controlling clinical manifestations of ankylosing spondylitis

Objective. To identify genomic regions linked with determinants of age at symptom onset, disease activity, and functional impairment in ankylosing spondylitis (AS). Methods. A whole genome linkage scan was performed in 188 affected sibling pair families with 454 affected individuals. Traits assessed were age at symptom onset, disease activity assessed by the Bath Ankylosing Spondylitis Disease Activity Index (BASDAI), and functional impairment assessed by the Bath Ankylosing Spondylitis Functional Index (BASFI). Parametric and nonparametric quantitative linkage analysis was performed using parameters defined in a previous segregation study. Results. Heritabilities of the traits studied in this data set were as follows: BASDAI 0.49 (P = 0.0001, 95% confidence interval [95% CI] 0.23-0.75), BASFI 0.76 (P = 10-7, 95% CI 0.49-1.0), and age at symptom onset 0.33 (P = 0.005, 95% CI 0.04-0.62). No linkage was observed between the major histocompatibility complex (MHC) and any of the traits studied (logarithm of odds [LOD] score <1.0). "Significant" linkage (LOD score 4.0) was observed between a region on chromosome 18p and the BASDAI. Age at symptom onset showed "suggestive" linkage to chromosome 11p (LOD score 3.3). Maximum linkage with the BASFI was seen at chromosome 2q (LOD score 2.9). Conclusion. In contrast to the genetic determinants of susceptibility to AS, clinical manifestations of the disease measured by the BASDAI, BASFI, and age at symptom onset are largely determined by a small number of genes not encoded within the MHC.

The effect of HLA-DR genes on susceptibility to and severity of ankylosing spondylitis

Objective. To analyze the effect of HLA-DR genes on susceptibility to and severity of ankylosing spondylitis (AS). Methods. Three hundred sixty- three white British AS patients were studied; 149 were carefully assessed for a range of clinical manifestations, and disease severity was assessed using a structured questionnaire. Limited HLA class I typing and complete HLA-DR typing were performed using DNA-based methods. HLA data from 13,634 healthy white British bone marrow donors were used for comparison. Results. A significant association between DR1 and AS was found, independent of HLA-B27 (overall odds ratio [OR] 1.4, 95% confidence interval [95% CI] 1.1-1.8, P = 0.02; relative risk [RR] 2.7, 95% CI 1.5-4.8, P = 6 x 10-4 among homozygotes; RR 2.1, 95% CI 1.5-2.8, P = 5 x 10-6 among heterozygotes). A large but weakly significant association between DR8 and AS was noted, particularly among DR8 homozygotes (RR 6.8, 95% CI 1.6-29.2, P = 0.01 among homozygotes; RR 1.6, 95% CI 1.0-2.7, P = 0.07 among heterozygotes). A negative association with DR12 (OR 0.22, 95% CI 0.09-0.5, P = 0.001) was noted. HLA-DR7 was associated with younger age at onset of disease (mean age at onset 18 years for DR7-positive patients and 23 years for DR7-negative patients; Z score 3.21, P = 0.001). No other HLA class I or class H associations with disease severity or with different clinical manifestations of AS were found. Conclusion. The results of this study suggest that HLA-DR genes may have a weak effect on susceptibility to AS independent of HLA-B27, but do not support suggestions that they affect disease severity or different clinical manifestations.

Susceptibility to ankylosing spondylitis in twins: The role of genes, HLA, and the environment

Objective To determine the relative effects of genetic and environmental factors in susceptibility to ankylosing spondylitis (AS). Methods Twins with AS were identified from the Royal National Hospital for Rheumatic Diseases database. Clinical and radiographic examinations were performed to establish diagnoses, and disease severity was assessed using a combination of validated scoring systems. HLA typing for HLA-B27, HLA-B60, and HLA-DR1 was performed by polymerase chain reaction with sequence- specific primers, and zygosity was assessed using microsatellite markers. Genetic and environmental variance components were assessed with the program Mx, using data from this and previous studies of twins with AS. Results Six of 8 monozygotic (MZ) twin pairs were disease concordant, compared with 4 of 15 B27-positive dizygotic (DZ) twin pairs (27%) and 4 of 32 DZ twin pairs overall (12.5%). Nonsignificant increases in similarity with regard to age at disease onset and all of the disease severity scores assessed were noted in disease-concordant MZ twins compared with concordant DZ twins. HLA-B27 and B60 were associated with the disease in probands, and the rate of disease concordance was significantly increased among DZ twin pairs in which the co- twin was positive for both B27 and DR1. Additive genetic effects were estimated to contribute 97% of the population variance. Conclusion Susceptibility to AS is largely genetically determined, and the environmental trigger for the disease is probably ubiquitous. HLA-B27 accounts for a minority of the overall genetic susceptibility to AS.

A genome-wide screen for susceptibility loci in ankylosing spondylitis

Objective. To localize the regions containing genes that determine susceptibility to ankylosing spondylitis (AS). Methods. One hundred five white British families with 121 affected sibling pairs with AS were recruited, largely from the Royal National Hospital for Rheumatic Diseases AS database. A genome-wide linkage screen was undertaken using 254 highly polymorphic microsatellite markers from the Medical Research Council (UK) (MRC) set. The major histocompatibility complex (MHC) region was studied more intensively using 5 microsatellites lying within the HLA class III region and HLA-DRB1 typing. The Analyze package was used for 2-point analysis, and GeneHunter for multipoint analysis. Results. When only the MRC set was considered, 11 markers in 7 regions achieved a P value of ≤0.01. The maximum logarithm of odds score obtained was 3.8 (P = 1.4 x 10-5) using marker D6S273, which lies in the HLA class III region. A further marker used in mapping of the MHC class III region achieved a LOD score of 8.1 (P = 1 x 10-9). Nine of 118 affected sibling pairs (7.6%) did not share parental haplotypes identical by descent across the MHC, suggesting that only 31% of the susceptibility to AS is coded by genes linked to the MHC. The maximum non-MHC LOD score obtained was 2.6 (P = 0.0003) for marker D16S422. Conclusion. The results of this study confirm the strong linkage of the MHC with AS, and provide suggestive evidence regarding the presence and location of non-MHC genes influencing susceptibility to the disease.

Ectopic calcification among families in the Azores: Clinical and radiologic manifestations in families with diffuse idiopathic skeletal hyperostosis and chondrocalcinosis

Objective. Twelve families that were multiply affected with diffuse idiopathic skeletal hyperostosis (DISH) and/or chondrocalcinosis, were identified on the island of Terceira, The Azores, potentially supporting the hypothesis that the 2 disorders share common etiopathogenic factors. The present study was undertaken to investigate this hypothesis. Methods. One hundred three individuals from 12 unrelated families were assessed. Probands were identified from patients attending the Rheumatic Diseases Clinic, Hospital de Santo Espirito, in The Azores. Family members were assessed by rheumatologists and radiologists. Radiographs of all family members were obtained, including radiographs of the dorsolumbar spine, pelvis, knees, elbows, and wrists, and all cases were screened for known features of chondrocalcinosis. Results. Ectopic calcifications were identified in 70 patients. The most frequent symptoms or findings were as follows: axial pain, elbow, knee and metacarpophalangeal (MCP) joint pain, swelling, and/or deformity, and radiographic enthesopathic changes. Elbow and MCP joint periarticular calcifications were observed in 35 and 5 patients, respectively, and chondrocalcinosis was identified in 12 patients. Fifteen patients had sacroiliac disease (ankylosis or sclerosis) on computed tomography scans. Fifty-two patients could be classified as having definite (17%), probable (26%), or possible (31%) DISH. Concomitant DISH and chondrocalcinosis was diagnosed in 12 patients. Pyrophosphate crystals were identified from knee effusions in 13 patients. The pattern of disease transmission was compatible with an autosomal-dominant monogenic disease. The mean age at which symptoms developed was 38 years. Conclusion. These families may represent a familial type of pyrophosphate arthropathy with a phenotype that includes peripheral and axial enthesopathic calcifications. The concurrence of DISH and chondrocalcinosis suggests a shared pathogenic mechanism in the 2 conditions.

Role of NOD2 variants in spondylarthritis

Objective. To investigate the role of the gene NOD2 in susceptibility to, and clinical manifestations of, ankylosing spondylitis (AS). Methods. A case-control study of NOD2 polymorphisms known to be associated with Crohn's disease (CD) (Pro268 Ser, Arg702 Trp, GlY908 Arg, and Len1007fsinsC) was performed in 229 cases of primary AS with no diagnosed inflammatory bowel disease (IBD), 197 cases of AS associated with IBD (referred to as colitic spondylarthritis; comprising 78 with CD and 119 with ulcerative colitis [UC]), and 229 ethnically matched, healthy controls. Associations between NOD2 polymorphisms and several clinical features of AS, including disease severity assessed by questionnaire and age at spondylarthritis onset, were also investigated. Exclusion linkage mapping of chromosome 16 was performed in a separate group of 185 multicase families with AS. Results. An association was identified between Gly908 Arg and UC spondylarthritis (P = 0.016, odds ratio [OR] 4.6, 95% confidence interval [95% CI] 1.316), and a nonsignificant trend with a similar magnitude was observed in association with CD spondylarthritis (P = 0.08, OR 3.9, 95% CI 0.8-18). The Pro268Ser variant was inversely associated with UC spondylarthritis (P = 0.003, OR 0.55, 95% CI 0.37-0.82), but not with CD spondylarthritis. No association was demonstrated between NOD2 variants and primary AS, or between other variants of NOD2 and either UC or CD spondylarthritis. Carriage of the Pro268 Ser polymorphism was associated with greater disease activity as measured by the Bath Ankylosing Spondylitis Disease Activity Index (P = 0.002). Although patients with CD had a younger age at spondylarthritis onset than did those with UC (22.4 years versus 26.4 years; P = 0.01), no association was noted between the NOD2 variants linked with CD and age at spondylarthritis onset. In primary AS, the presence of a gene with a magnitude of association >2.0 was excluded (exclusion logarithm of odds score less than -2.0), and no association was observed with the microsatellite D16S3136. Conclusion. NOD2 variants do not significantly affect the risk of developing primary AS, but may influence susceptibility to, and clinical manifestations of, colitic spondylarthritis.

High-throughput sequencing of IL23R reveals a low-frequency, nonsynonymous single-nucleotide polymorphism that is associated with ankylosing spondylitis in a Han Chinese population

Objective Ankylosing spondylitis (AS) is a highly heritable common inflammatory arthritis that targets the spine and sacroiliac joints of the pelvis, causing pain and stiffness and leading eventually to joint fusion. Although previous studies have shown a strong association of IL23R with AS in white Europeans, similar studies in East Asian populations have shown no association with common variants of IL23R, suggesting either that IL23R variants have no role or that rare genetic variants contribute. The present study was undertaken to screen IL23R to identify rare variants associated with AS in Han Chinese. Methods A 170-kb region containing IL23R and its flanking regions was sequenced in 50 patients with AS and 50 ethnically matched healthy control subjects from a Han Chinese population. In addition, the 30-kb region of peak association in white Europeans was sequenced in 650 patients with AS and 1,300 healthy controls. Validation genotyping was undertaken in 846 patients with AS and 1,308 healthy controls. Results We identified 1,047 variants, of which 729 were not found in the dbSNP genomic build 130. Several potentially functional rare variants in IL23R were identified, including one nonsynonomous single-nucleotide polymorphism (nsSNP), Gly149Arg (position 67421184 GA on chromosome 1). Validation genotyping showed that the Gly149Arg variant was associated with AS (odds ratio 0.61, P = 0.0054). Conclusion This is the first study to implicate rare IL23R variants in the pathogenesis of AS. The results identified a low-frequency nsSNP with predicted loss-of-function effects that was protectively associated with AS in Han Chinese, suggesting that decreased function of the interleukin-23 (IL-23) receptor protects against AS. These findings further support the notion that IL-23 signaling has an important role in the pathogenesis of AS.

Omission of author name in the article by Davidson et al (Arthritis Rheum, July 2013)

This article corrects: Brief Report: High-Throughput Sequencing of IL23R Reveals a Low-Frequency, Nonsynonymous Single-Nucleotide Polymorphism That Is Associated With Ankylosing Spondylitis in a Han Chinese Population Vol. 65, Issue 7, 1747–1752, Article first published online: 2 JUL 2013

Association of ERAP1, but not IL23R, with ankylosing spondylitis in a Han Chinese population

Objective The results of a recent genome-wide association study have shown that ERAP1 and IL23R are associated with ankylosing spondylitis (AS) in Caucasian populations from North America and the UK. Based on these findings, we undertook the current study to investigate whether single-nucleotide polymorphisms (SNPs) covering the genes ERAP1 and IL23R are associated with AS in a Han Chinese population. Methods A case-control study was performed in Han Chinese patients with AS (n = 527) and controls (n = 945) from Shanghai and Nanjing. All patients met the modified New York criteria for AS. The Sequenom iPlex platform was used to genotype cases and controls for 21 tag SNPs covering IL23R and 38 tag SNPs covering ERAP1. Statistical analysis was performed using the Cochran-Armitage test for trend. Results Multiple SNPs in ERAP1 were significantly associated with AS (for rs27980, P = 0.0048; for rs7711564, P = 0.0081). However, no association was observed between IL23R and AS (for all SNPs, P > 0.1). The nonsynonymous SNP in IL23R, rs11209026, widely thought to be the primary AS-associated SNP in IL23R in Europeans, was found not to be polymorphic in Chinese. Conclusion Our results demonstrate that genetic polymorphisms in ERAP1 are associated with AS in Han Chinese, suggesting a common pathogenic mechanism for the disease in Chinese and Caucasian populations, and that IL23R is not associated with AS in Chinese, indicating a difference in the mechanism of disease pathogenesis between Chinese and Caucasian populations. This may result from the fact that rs11209026, the nonsynonymous SNP in IL23R, is not polymorphic in Chinese patients, providing further evidence that rs11209026 is the key polymorphism associated with AS (and likely inflammatory bowel disease and psoriasis) in this gene.

A linkage study across the T cell receptor A and T cell receptor B loci in families with rheumatoid arthritis

Objective. To examine whether the T cell receptor (TCR) A or TCRB loci exhibit linkage with disease in multiplex rheumatoid arthritis (RA) families. Methods. A linkage study was performed in 184 RA families from the UK Arthritis and Rheumatism Council Repository, each containing at least 1 affected sibpair. The microsatellites D14S50, TCRA, and D14S64 spanning the TCRA locus and D7S509, Vβ6.7, and D7S688 spanning the TCRB locus were used as DNA markers. The subjects were genotyped using a semiautomated polymerase chain reaction-based method. Two-point and multipoint linkage analyses were performed. Results. Nonparametric single-marker likelihood odds (LOD) scores were 0.49 (P = 0.07) for D14S50, 0.65 (P = 0.04) for TCRA, 0.07 (P = 0.29) for D14S64, 0.01 (P = 0.43) for D7S509, 0.0 (P = 0.50) for Vβ6.7, and 0.0 (P = 0.50) for D7S688. By multipoint analysis, there was no evidence of linkage at TCRB (LOD score 0), and the maximum LOD score at the TCRA locus was 0.37 (at D14S50). The presence of a susceptibility locus (LOD score < -2.0) was excluded, with lambda ≤ 1.8 at TCRA and ≤1.4 at TCRB. Conclusion. These linkage studies provide no significant evidence of a major germline-encoded TCRA or TCRB component of susceptibility to RA.

Is disease severity in ankylosing spondylitis genetically determined?

Objective. To assess the role of genes and the environment in determining the severity of ankylosing spondylitis. Methods: One hundred seventy-three families with >1 case of ankylosing spondylitis were recruited (120 affected sibling pairs, 26 affected parent-child pairs, 20 families with both first- and second-degree relatives affected, and 7 families with only second-degree relatives affected), comprising a total of 384 affected individuals. Disease severity was assessed by the Bath Ankylosing Spondylitis Disease Activity Index (BASDAI) and functional impairment was determined using the Bath Ankylosing Spondylitis Functional Index (BASFI). Disease duration and age at onset were also studied. Variance-components modeling was used to determine the genetic and environmental components Contributing to familiality of the traits examined, and complex segregation analysis was performed to assess different disease models. Results. Both the disease activity and functional capacity as assessed by the BASDAI and the BASFI, respectively, were found to be highly familial (BASDAI familiality 0.51 [P = 10-4], BASFI familiality 0,68 [P = 3 × 10-7]). No significant shared environmental component was demonstrated to be associated with either the BASDAI or the BASFI. Including age at disease onset and duration of disease as covariates made no difference in the heritability assessments. A strong correlation was noted between the BASDAI and the BASFI (genetic correlation 0.9), suggesting the presence of shared determinants of these 2 measures. However, there was significant residual heritability for each measure independent of the other (BASFI residual heritability 0.48, BASDAI 0,36), perhaps indicating that not all genes influencing disease activity influence chronicity. No significant heritability of age at disease onset was found (heritability 0.18; P = 0.2). Segregation studies suggested the presence of a single major gene influencing the BASDAI and the BASFI. Conclusion. This study demonstrates a major genetic contribution to disease severity in ankylosing spondylitis. As with susceptibility to ankylosing spondylitis, shared environmental factors play little role in determining the disease severity.