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.

The X-chromosome and susceptibility to ankylosing spondylitis

Objective. Ankylosing spondylitis (AS) affects 0.25-1.0% of the population, and its etiology is incompletely understood. Susceptibility to this highly familial disease (λ(s) = 58) is primarily genetically determined. There is a significant sex bias in AS, and there are differences in recurrence risk to the offspring of affected mothers and fathers, suggesting that there may be an X-linked recessive effect. We undertook an X- chromosome linkage study to determine any contribution of the X-chromosome to AS susceptibility. Methods. A linkage study of the X-chromosome using 234 affected sibling pairs was performed to investigate this hypothesis. Results. No linkage of the X-chromosome with susceptibility to AS was found. Model- free multipoint linkage analysis strongly excluded any significant genetic contribution (λ ≥1.5) to AS susceptibility encoded on the X-chromosome (logarithm of odds [LOD] <-2.0). Smaller genetic effects (A ≥1.3) were also found to be unlikely (LOD <-1.0). Conclusion. The sex bias in AS is not explained by X-chromosome-encoded genetic effects. The disease model best explaining the sex bias in occurrence and transmission of AS is a polygenic model with a higher susceptibility threshold in females.

Enrichment of circulating interleukin-17-secreting interleukin-23 receptor-positive γ/δ T cells in patients with active ankylosing spondylitis

Objective Ankylosing spondylitis (AS) is a common inflammatory arthritis affecting primarily the axial skeleton. IL23R is genetically associated with AS. This study was undertaken to investigate and characterize the role of interleukin-23 (IL-23) signaling in AS pathogenesis. Methods The study population consisted of patients with active AS (n = 17), patients with psoriatic arthritis (n = 8), patients with rheumatoid arthritis, (n = 9), and healthy subjects (n = 20). IL-23 receptor (IL-23R) expression in T cells was determined in each subject group, and expression levels were compared. Results The proportion of IL-23R-expressing T cells in the periphery was 2-fold higher in AS patients than in healthy controls, specifically driven by a 3-fold increase in IL-23R-positive γ/δ T cells in AS patients. The proportions of CD4+ and CD8+ cells that were positive for IL-17 were unchanged. This increased IL-23R expression on γ/δ T cells was also associated with enhanced IL-17 secretion, with no observable IL-17 production from IL-23R-negative γ/δ T cells in AS patients. Furthermore, γ/δ T cells from AS patients were heavily skewed toward IL-17 production in response to stimulation with IL-23 and/or anti-CD3/CD28. Conclusion Recently, mouse models have shown IL-17-secreting γ/δ T cells to be pathogenic in infection and autoimmunity. Our data provide the first description of a potentially pathogenic role of these cells in a human autoimmune disease. Since IL-23 is a maturation and growth factor for IL-17-producing cells, increased IL-23R expression may regulate the function of this putative pathogenic γ/δ T cell population.

Whole-genome linkage analysis of rheumatoid arthritis susceptibility loci in 252 affected sibling pairs in the United Kingdom

Objective. To undertake a systematic wholegenome screen to identify regions exhibiting genetic linkage to rheumatoid arthritis (RA). Methods. Two hundred fifty-two RA-affected sibling pairs from 182 UK families were genotyped using 365 highly informative microsatellite markers. Microsatellite genotyping was performed using fluorescent polymerase chain reaction primers and semiautomated DNA sequencing technology. Linkage analysis was undertaken using MAPMAKER/SIBS for single-point and multipoint analysis. Results. Significant linkage (maximum logarithm of odds score 4.7 [P = 0.000003] at marker D6S276, 1 cM from HLA-DRB1) was identified around the major histocompatibility complex (MHC) region on chromosome 6. Suggestive linkage (P < 7.4 × 10-4) was identified on chromosome 6q by single- and multipoint analysis. Ten other sites of nominal linkage (P < 0.05) were identified on chromosomes 3p, 4q, 7p, 2 regions of 10q, 2 regions of 14q, 16p, 21q, and Xq by single-point analysis and on 3 sites (1q, 14q, and 14q) by multipoint analysis. Conclusion. Linkage to the MHC region was confirmed. Eleven non-HLA regions demonstrated evidence of suggestive or nominal linkage, but none reached the genome-wide threshold for significant linkage (P = 2.2 × 10-5). Results of previous genome screens have suggested that 6 of these regions may be involved in RA susceptibility.

The effect of HLA-DR on susceptibility to rheumatoid arthritis is influenced by the associated lymphotoxin α-tumor necrosis factor haplotype

Objective. HLA-DRB1, a major genetic determinant of susceptibility to rheumatoid arthritis (RA), is located within 1,000 kb of the gene encoding tumor necrosis factor (TNF). Because certain HLA-DRB1*04 subtypes increase susceptibility to RA, investigation of the role of the TNF gene is complicated by linkage disequilibrium (LD) between TNF and DRB1 alleles. By adequately controlling for this LD, we aimed to investigate the presence of additional major histocompatibility complex (MHC) susceptibility genes. Methods. We identified 274 HLA-DRB1*04-positive cases of RA and 271 HLA-DRB1*04-positive population controls. Each subject was typed for 6 single-nucleotide polymorphisms within a 4.5-kb region encompassing TNF and lymphotoxin a (LTA). LTA-TNF haplotypes in these unrelated individuals were determined using a combination of family data and the PHASE software program. Results. Significant differences in LTA-TNF haplotype frequencies were observed between different subtypes of HLA-DRB1*04. The LTA-TNF haplotypes observed were very restricted, with only 4 haplotypes constituting 81% of all haplotypes present. Among individuals carrying DRB1*0401, the LTA-TNF 2 haplotype was significantly underrepresented in cases compared with controls (odds ratio 0.5 [95% confidence interval 0.3-0.8], P = 0.007), while in those with DRB1*0404, the opposite effect was observed (P = 0.007). Conclusion. These findings suggest that the MHC contains genetic elements outside the LTA-TNF region that modify the effect of HLA-DRB1 on susceptibility to RA.