Article Text
Abstract
Objectives: To investigate whether HLA alleles, other than HLA-B27, influence predisposition to spondyloarthropathy (SpA) in multiplex families.
Methods: Seventy French families with at least two affected SpA members were recruited. Patients, and their first degree relatives were typed for HLA-A, B, C, and DR, and extended HLA haplotypes were determined. The distribution of HLA-A, C, and DR alleles carried on HLA-B27+ haplotypes in SpA families was compared with the distribution of these alleles among HLA-B27+ haplotypes in the French general population. Contribution to SpA susceptibility of HLA-A, B, C, and DR alleles, other than HLA-B27, was tested by transmission disequilibrium test. The contribution of HLA alleles to specific presentation features of SpA was examined.
Results: Frequencies of HLA-A, C, and DR alleles carried on HLA-B27+ haplotypes from SpA families were comparable with those seen in the French population, except for DR13 which was overrepresented among patients (pcorr<0.001). Most interestingly, the HLA-DR4 allele was transmitted in excess to patients with SpA, independently of linkage to HLA-B27 (pcorr=0.05), and in a direction opposite to that for HLA-B27+ unaffected siblings (pcorr=0.01). Finally, the distribution of HLA alleles was not related to the presentation feature of SpA.
Conclusion: HLA predisposition to familial SpA appears not to be limited to HLA-B27, but some HLA-DR alleles also have a significant influence. In particular, HLA-DR4 contributes significantly to a genetic predisposition to SpA, which may have implications in our understanding of SpA pathogenesis.
- HLA
- spondyloarthropathy
- ankylosing spondylitis
- transmission disequilibrium test
- AS, ankylosing spondylitis
- CI, confidence interval
- IBD, inflammatory bowel disease
- MHC, major histocompatibility complex
- PsA, psoriatic arthritis
- RA, rheumatoid arthritis
- RR, relative risk
- SpA, spondyloarthropathy
- TDT, transmission disequilibrium test
- uSpA, undifferentiated spondyloarthropathy
Statistics from Altmetric.com
- AS, ankylosing spondylitis
- CI, confidence interval
- IBD, inflammatory bowel disease
- MHC, major histocompatibility complex
- PsA, psoriatic arthritis
- RA, rheumatoid arthritis
- RR, relative risk
- SpA, spondyloarthropathy
- TDT, transmission disequilibrium test
- uSpA, undifferentiated spondyloarthropathy
Ankylosing spondylitis (AS), reactive arthritis, a subset of psoriatic arthritis (PsA), arthritis associated with idiopathic inflammatory bowel diseases (IBD), and undifferentiated spondyloarthropathy (uSpA) form a group of inflammatory rheumatic diseases identified as the spondyloarthropathies (SpA).1 It has long been recognised that different presentation features of SpA tend to combine within the same families.2,3 We previously conducted a systematic analysis of SpA phenotypic expression in a large number of multiplex families.4,5 Results of this study support the hypothesis that although familial SpA shows several distinct phenotypes, these all have predisposing genetic factors in common.4,5 This indicates that, at least in this particular context, genetic investigations should be conducted on the whole group of patients with SpA, rather than restricted to specific subsets.6
One genetic factor shared by all the varieties of SpA has already been identified—that is, HLA-B27.7 Studies on AS suggest that additional genetic factors, distinct from the HLA region, are involved.8,9 However, it has also been proposed that HLA-B27 itself may not explain all the predisposition to SpA borne within the major histocompatibility complex (MHC) region.10,11 In this study we took advantage of having recruited a large number of SpA multiplex families to investigate whether HLA genes, other than HLA-B27, might also contribute to SpA predisposition.
SUBJECTS AND METHODS
Subjects
Seventy multiplex families with SpA who have been recruited throughout France by the Groupe Français d'Etude Génétique des Spondylarthropathies, as previously described,4,5 were included in this study. Of the 188 patients identified, 187 met the criteria of Amor et al,12 and 182 met the European Spondylarthropathy Study Group (ESSG) criteria.1Table 1⇓ shows characteristics of the patients. One hundred and twenty (64%) of the patients had AS according to the modified New York criteria.13 Among the 67 (36%) patients with SpA without definite sacroiliitis, PsA was diagnosed in 12, arthritis associated with IBD in 4, Reiter's syndrome in 1,14 and uSpA in the remaining 51. An age matched group of 71 HLA-B27+ siblings belonging to the same families, for which a diagnosis of SpA could not be made, was used as a control group (table 1⇓).
Characteristics of the study group*
HLA typing
Standard serological methods were routinely used for HLA-A, B, and C typing.15 In several cases, HLA-A and B generic typing was performed on DNA extracted from peripheral venous blood leucocytes by a polymerase chain reaction-based sequence-specific primers method, using a commercial kit designed to recognise all the known broad specificities (Dynal France SA, Compiégne, France). HLA-DRB1 typing was performed by polymerase chain reaction-sequence-specific oligonucleotide probing in all cases, and frequently also by a standard serological method. HLA-A, B, and DRB1 typing was performed on all 188 patients with SpA, and their first degree healthy relatives. HLA-C typing was only performed on a first set of 38 families. MHC haplotypes were deduced by segregation analysis of HLA-A, B, C, and DRB1 genotypes within pedigrees.
Study design and statistical analysis
Influence of HLA alleles, other than HLA-B27, on predisposition to SpA was investigated differently for alleles encoded by an HLA-B27+ haplotype and for alleles encoded by an HLA-B27− haplotype.
The frequency of alleles encoded by an HLA-A, C, or DRB1 locus carried on an HLA-B27+ haplotype was compared with the distribution of those alleles in a reference panel of 198 HLA-B27+ haplotypes representing the French general population. This control panel was retrieved from a collection of 3424 haplotypes (Borot N, unpublished results), which have been previously established on a sample of 1356 French families collected equally from the different French regions (to constitute a representative sample), and fully typed across the MHC by serological methods.16 For this part of the analysis, each individual haplotype identified in an SpA family was counted only once, even if present in several patients. Allelic frequencies were compared between familial patients with SpA and the control French population, using Fisher's exact test. Odds ratio and its 95% confidence interval (CI) were also determined.
The contribution of individual alleles to SpA predisposition was also tested by transmission disequilibrium test (TDT), as previously described.17 Alleles were tested independently of each other, using binomial probability to determine whether transmission of a given allele by heterozygous parents to affected offsprings was significantly different than the expected frequency of 50%. In addition to patients, HLA-B27+ healthy siblings were also included in the analysis, to test for transmission distortion in the direction opposite to that for SpA siblings, using the unmatched 2×2 χ2 test.18 Alleles encoded by an HLA-B27− haplotype, could specifically be tested by this method to discriminate effects independent of linkage with HLA-B27.
Comparison of allele frequency between patients and unaffected HLA-B27+ siblings, and between different subgroups of patients with SpA, was performed, using Fisher's exact test. The relative risk (RR), and its 95% CI were also determined.
In all cases, pcorr refers to corrected p values which were obtained by multiplying crude p values by the number of comparisons. Definite statistical significance required a pcorr<0.05.
RESULTS
Frequency of HLA-A, C, and DR alleles carried on HLA-B27+ haplotypes
Typing of HLA-A, B, and DRB1 in 70 French families with multiple cases of SpA, and of HLA-C in a subset of 38 families resulted in the identification of 81 distinct HLA-B27+ extended haplotypes in patients. Noteworthy, is that in 10 families, two B27+ haplotypes were found in B27 homozygous patients. Of those haplotypes, 48 included HLA-C (A;C;B27;DR), whereas 33 did not (A;B27;DR). Allelic distribution at locus A (table 2⇓) and C (table 3⇓) in familial SpA did not differ from the distribution seen in a panel of HLA-B27+ haplotypes from the French general population. Distribution of HLA-DR alleles was also similar between familial SpA and the control group for most alleles, except for DR13, which was significantly overrepresented among SpA HLA-B27+ haplotypes (table 4⇓).
Comparison of the frequency of HLA-A alleles carried on HLA-B27+ haplotype between patients with familial spondyloarthropathy and the French control population*
Comparison of the frequency of HLA-C alleles carried on HLA-B27 positive haplotype between patients with familial spondyloarthropathy and the French control population*
Comparison of the frequency of HLA-DR alleles carried on HLA-B27 positive haplotype between patients with familial spondyloarthropathy and the French control population*
Transmission of HLA-A, B, and C alleles to patients
TDT was used to determine if class I MHC alleles, other than HLA-B27, were randomly transmitted to patients with SpA. For each HLA-A and C allele tested we performed a global TDT which included all haplotypes, and also a TDT restricted to HLA-B27− haplotypes, to search for transmission disequilibrium independent of linkage between A or C alleles and HLA-B27. For HLA-B alleles, only TDT applied to HLA-B27− haplotypes is shown.
Transmission of most HLA-A, B, and C alleles to patients was not significantly different from random. A limited number of alleles encoded by those class I MHC loci displayed some degree of unbalanced transmission, which failed to reach statistical significance after correction for the number of comparisons (table 5⇓). Such weak transmission disequilibrium, was rarely observed in TDT restricted to HLA-B27− haplotypes. Notably, a non-significant trend towards an excess of transmission to patients was seen for the HLA-B60 allele (table 5⇓).
Results of the transmission disequilibrium test performed for HLA-A, B, and C alleles in patients with familial spondyloarthropathy*
Transmission of HLA-DR alleles to patients and healthy HLA-B27+ siblings
TDT of HLA-DR alleles to patients with SpA was first performed with alleles carried on all haplotypes. Results of this TDT displayed unbalanced transmission of several alleles, reaching statistical significance after correction for the number of comparisons (table 6⇓). Hence, an excess of transmission to patients with SpA was observed for DR1 (pcorr=0.02) and DR4 (pcorr=0.002), whereas a reduced transmission was noticed for DR2 (pcorr=0.05). In addition, trends towards reduced transmission were found for DR3 and DR13, but failed to reach statistical significance after correction (table 6⇓). The transmission of HLA-DR alleles to unaffected HLA-B27+ siblings of patients with SpA was analysed in parallel using TDT, and compared with the pattern of transmission to patients with SpA, using an unmatched 2×2 χ2 test. We found a significant distortion of transmission of DR4 allele to HLA-B27+ unaffected siblings, in the direction opposite to that for SpA siblings (table 6⇓; pcorr= 0.006). No such distortion was seen with any other HLA-DR allele. These results suggest that an excess of transmission of DR4 to a patient with SpA might exist independently of linkage disequilibrium between DR4 and HLA-B27. A TDT restricted to HLA-DR alleles carried on HLA-B27− haplotypes was next performed to examine this possibility. This TDT showed an excess of transmission of the DR4 allele to patients with SpA (table 7⇓; pcorr=0.05). Furthermore, a distortion of transmission of DR4 allele to HLA-B27+ unaffected siblings in the direction opposite to that for SpA siblings was also confirmed with HLA-B27− haplotypes (table 7⇓; pcorr= 0.01).
Results of the transmission disequilibrium test performed for HLA-DR alleles carried on all haplotypes*
Results of the transmission disequilibrium test performed for HLA-DR alleles carried on HLA-B27 negative haplotypes*
Contribution of HLA alleles to specific presentation features of SpA
The foregoing results were obtained by examining all patients with SpA from mutiplex families. We also considered whether specific features of the SpA spectrum were preferentially associated with HLA alleles. No statistically significant association of HLA-A, B, and C alleles with specific manifestations was seen after correcting for the number of comparisons (data not shown). Likewise, most HLA-DR alleles were uniformly distributed among patients with SpA, whatever the manifestation examined (table 8⇓). Only the DR13 allele appeared to be less frequent among patients with sacroiliitis, and was conversely overrepresented among patients with IBD, but the differences failed to reach statistical significance (table 8⇓). Notably, the DR13 allele which was present in 50% of patients with IBD was carried on a B27− haplotype. Finally, the HLA-DR4 allele appeared to be uniquely overrepresented among HLA-B27+ patients, as compared with unaffected HLA-B27+ siblings (table 8⇓; pcorr=0.04).
Frequency of HLA-DR alleles in familial spondyloarthropathy (SpA) according to specific manifestations*
DISCUSSION
Genetic linkage between AS, the prototypical form of SpA, and the MHC region has been previously established.9,10,19 The strongest association between the MHC antigen and AS is unequivocally attributable to HLA-B27.11,20 However, there are some suggestions that other MHC alleles besides HLA-B27 may also contribute towards increasing susceptibility for AS.10,11 Numerous studies have attempted in the past to identify such non-B27 genes within the MHC which would directly contribute to AS susceptibility.11,21 However, such studies have commonly been hampered by the strong linkage disequilibrium between alleles encoded in separate loci within the MHC region, by the overwhelming presence of HLA-B27 in AS, and also by a lack of power with regard to highly polymorphic HLA loci. In this study we took advantage of having recruited a large number of mutiplex SpA families to examine separately the contribution of MHC alleles encoded by an HLA-B27+ haplotype and the contribution of alleles encoded by B27− haplotypes.
To analyse the genetic contribution of HLA alleles encoded by an HLA-B27+ haplotype, we compared the frequency of alleles encoded by SpA haplotypes to alleles encoded by HLA-B27+ haplotypes from the French general population. The distribution of alleles was remarkably similar between both populations, even though the control haplotypes were obtained from a historical study.16 This result confirms, on a larger scale, previous reports concerning HLA-A and C loci.22 It is also consistent with a priority role for HLA-B27 over putative secondary predisposing gene(s) located within the MHC region.20 The only remarkable difference concerned the HLA-DR13 allele, which was significantly represented among SpA B27+ haplotypes, whereas it was not detected among HLA-B27+ control haplotypes. This absence of DR13 in controls is unlikely to reflect a technical bias, because the DR13 allele was commonly detected among B27− haplotypes in the French population study (Borot N, unpublished results). Rather, it is consistent with a previously reported negative linkage disequilibrium between DR13 and HLA-B27 in the general population.23 Nevertheless, TDT failed to support linkage between DR13 and SpA within families. Hence, this increased prevalence of B27;DR13 haplotypes in SpA families may reflect population stratification (as a consequence of some putative degree of common origin, which would make patients more likely to share the DR13 allele with each other than with a control population), rather than a direct role for the DR13 allele in genetic predisposition to SpA. Alternatively, it might be secondary to the presence of another as yet unidentified predisposing gene carried on the B27;DR13 haplotype.
The most striking findings in this study were obtained with TDT, which allows one to determine both linkage and association between disease and given alleles. When all haplotypes were included in the TDT, unbalanced transmission of several HLA alleles to patients was found. In most cases this result was not confirmed by TDT restricted to HLA-B27− haplotypes, suggesting that it was accounted for by linkage disequilibrium between such alleles and HLA-B27. Notably, the excess transmission of HLA-DR1 to patients could entirely be attributed to linkage disequilibrium with B27 in this study, albeit DR1 has been suspected to predispose to AS independently of HLA-B27.23 In contrast, HLA-DR4 was transmitted in excess to patients independently of its positive linkage disequilibrium with HLA-B27. Furthermore, an opposite trend towards a lack of transmission of DR4 to unaffected HLA-B27+ siblings was seen. Together, these data show that the presence of DR4 contributed to SpA susceptibility, in addition to the HLA-B27 effect. For some of the other alleles, the small amount of meiosis available for TDT in this study might have been associated with a lack of power in demonstrating a statistically significant effect. This particularly might be the case for HLA-B60, which has been reported in several studies to increase the risk of AS.11
In this study we analysed the whole spectrum of SpA together, because we had previously shown that in familial SpA different subsets of SpA probably share major predisposition factors.4,5 Nevertheless, genes located within the MHC region may contribute to phenotypic diversity. In particular, HLA-DR4 might have preferentially contributed to peripheral arthritis, as previously suggested.24 Such an assumption was not confirmed by our data, which showed equivalent association between DR4 and SpA, whatever the presentation feature. Notably, this study had a statistical power >95% to detect a twofold increase in peripheral arthritis among patients with DR4. Furthermore, this study failed to repeat any of the previously reported preferential associations between specific features of SpA and HLA alleles.21 Several explanations might account for such discrepancies, among which are ethnic variations between different studies and the familial context, which was specific to this study. Another study in AS, which was largely based on familial cases, also failed to identify a significant influence of HLA-DR alleles on phenotypic presentation, including peripheral arthritis, uveitis, psoriasis, or IBD.23 Results of this study are suggestive of a negative association between HLA-DR13 and sacroiliitis and, conversely, of a positive association between DR13 and IBD, which was not statistically significant. However, this study had only a power of 12% to demonstrate an RR of IBD among DR13 patients, of a magnitude of two with 95% CI (that is, without correcting for the number of comparisons), because of the low proportion of IBD cases. Hence, we would have needed to include 4213 subjects to reach an 80% power of demonstrating an RR of such magnitude, or 1583 subjects to reach an 80% power of demonstrating an RR of three, after correcting for the number of comparisons. Nevertheless, these data highlight again HLA-DR13 as a candidate marker to be further investigated in a genetic study of SpA.
To our knowledge, this study is the first to show that HLA-DR4, another HLA allele within the MHC region besides HLA-B27, is associated with predisposition to SpA. The possibility that preferential association between SpA and a more important MHC gene in tight linkage disequilibrium with DR4 accounted for this result cannot entirely be excluded. Nevertheless, it is a further step towards improving our understanding of the role of the MHC in SpA. It has been proposed that CD4+ T cells might have a more important role in SpA pathogenesis than CD8+ T cells.25 In line with this hypothesis, HLA-DR4, which is a class II allele, might preferentially affect CD4+ T cell biology. Most interestingly, HLA-DR4 is also known as a major genetic factor predisposing to rheumatoid arthritis (RA).26 It has been reported that HLA-B27, in addition to HLA-DR4, may contribute to RA susceptibility.27,28 Such observations support the hypothesis that both alleles might contribute additively to chronic inflammatory joint disease, and that RA and SpA may share some pathogenic mechanism. Further work will be required to identify whether DR4 subtypes associated with RA are those preferentially transmitted to patients with SpA, which would strengthen this latter hypothesis.
Acknowledgments
We gratefully acknowledge the contribution of Dr Bernard Amor, Dr Francis Berenbaum, Dr Jean-Marie Berthelot, Dr Jean-Paul Blanquet, Dr Bernard Combe, Dr Emmanuelle Dernis-Labous, Dr Agnés Duché, Dr Jacques Fechtembaum, Dr Sandrine Guis, Dr Pierre Miossec, Dr Aleth Perdriger, Dr Anne-Marie Prieur, and Dr Alain Saraux, from Groupe Français d'Etude Génétique des Spondylarthropathies to the selection of patients, and the excellent technical assistance of Maryse Meziéres, Anne-Marie Poussy, and Gilles Michard.
This work was supported by grants from Programme Hospitalier de Recherche Clinique (PHRC; 1995), Société Française de Rhumatologie (SFR; 1997 and 1999), Association de Recherche Clinique en Rhumatologie (ARCR; 1996 and 1997), and Searle-Monsanto (1998 and 1999). Corinne Miceli-Richard was supported by a grant from Association de Recherche sur la Polyarthrite (ARP; 1999 and 2000).