- Research article
- Open access
- Published:
Association of HLA-DRB1 amino acid residues with giant cell arteritis: genetic association study, meta-analysis and geo-epidemiological investigation
Arthritis Research & Therapy volume 17, Article number: 195 (2015)
Abstract
Introduction
Giant cell arteritis (GCA) is an autoimmune disease commonest in Northern Europe and Scandinavia. Previous studies report various associations with HLA-DRB1*04 and HLA-DRB1*01; HLA-DRB1 alleles show a gradient in population prevalence within Europe. Our aims were (1) to determine which amino acid residues within HLA-DRB1 best explained HLA-DRB1 allele susceptibility and protective effects in GCA, seen in UK data combined in meta-analysis with previously published data, and (2) to determine whether the incidence of GCA in different countries is associated with the population prevalence of the HLA-DRB1 alleles that we identified in our meta-analysis.
Methods
GCA patients from the UK GCA Consortium were genotyped by using single-strand oligonucleotide polymerization, allele-specific polymerase chain reaction, and direct sequencing. Meta-analysis was used to compare and combine our results with published data, and public databases were used to identify amino acid residues that may explain observed susceptibility/protective effects. Finally, we determined the relationship of HLA-DRB1*04 population carrier frequency and latitude to GCA incidence reported in different countries.
Results
In our UK data (225 cases and 1378 controls), HLA-DRB1*04 carriage was associated with GCA susceptibility (odds ratio (OR) = 2.69, P = 1.5×10−11), but HLA-DRB1*01 was protective (adjusted OR = 0.55, P = 0.0046). In meta-analysis combined with 14 published studies (an additional 691 cases and 4038 controls), protective effects were seen from HLA-DR2, which comprises HLA-DRB1*15 and HLA-DRB1*16 (OR = 0.65, P = 8.2×10−6) and possibly from HLA-DRB1*01 (OR = 0.73, P = 0.037). GCA incidence (n = 17 countries) was associated with population HLA-DRB1*04 allele frequency (P = 0.008; adjusted R2 = 0.51 on univariable analysis, adjusted R2 = 0.62 after also including latitude); latitude also made an independent contribution.
Conclusions
We confirm that HLA-DRB1*04 is a GCA susceptibility allele. The susceptibility data are best explained by amino acid risk residues V, H, and H at positions 11, 13, and 33, contrary to previous suggestions of amino acids in the second hypervariable region. Worldwide, GCA incidence was independently associated both with population frequency of HLA-DRB1*04 and with latitude itself. We conclude that variation in population HLA-DRB1*04 frequency may partly explain variations in GCA incidence and that HLA-DRB1*04 may warrant investigation as a potential prognostic or predictive biomarker.
Introduction
Genetic association studies are a well-established method for investigating genetic contributions to disease. In rheumatoid arthritis (RA) [1] and small-vessel vasculitis [2], genetically distinct subsets have been identified that have different associations with the major histocompatibility complex (MHC) region that encodes the HLA-DRB1 alleles. Comparison of HLA-DRB1 associations with RA in different ethnic groups helped to support the original “shared epitope” hypothesis of RA susceptibility [3] based on an amino acid risk motif at positions 67–74 in the third hypervariable region (HVR3) of the class II MHC molecule, encoded by HLA-DRB1. The group of RA “shared epitope” alleles now includes HLA-DRB1*01:01, HLA-DRB1*04:01, HLA-DRB1*04:04, HLA-DRB1*04:08, and HLA-DRB1*10:01; other alleles provide weaker protective effects, additional to the risk effects of the “shared epitope” [4]. Recently, it has been demonstrated that amino acid residues 11 and 13 in the first hypervariable region (HVR1) of class II MHC display the strongest associations with RA susceptibility [5].
Giant cell arteritis (GCA) incidence is highest in populations with Scandinavian ancestry [6–8] and this has led to suggestions that this might be due to genetic factors [9, 10]. Susceptibility to GCA has been reported to be associated with carriage of HLA-DRB1*04, but not all studies have shown an association and there are conflicting data as to whether there is an association with specific HLA-DRB1*04 alleles [11]. Under-representation of HLA-DRB1*01 in GCA patients from Rochester, Minnesota, led to a suggestion that the risk of GCA may be due to a DRYF motif at positions 28–31 in the second hypervariable region (HVR2) of MHC class II [12], but a Spanish study failed to replicate this finding [13]. To date, however, formal meta-analysis has not yet been performed to determine the major susceptibility and protective HLA-DRB1 alleles. The relative contribution of genetic and environmental factors as an explanation for geographical differences in GCA incidence also remains disputed [6, 7]. When genetic diversity within Europe is subjected to principal component analysis, the MHC is one of several genetic regions that are strongly associated with a component that runs along a north-south gradient from Norway/Sweden to Spain [14]. We therefore hypothesised that variations in the frequency of HLA-DRB1 GCA susceptibility alleles may partly explain geographical variations in GCA incidence.
Here, we report new GCA susceptibility data, combine these with the published data using meta-analysis, and propose a new hypothesis regarding a possible amino acid GCA 11-13-33 risk motif in HVR1 and HVR2 of class II MHC. This hypothesis fits the observed data better than previously proposed models.
Methods
Patients
The UK GCA Consortium was designed to support genetic association studies. Investigators, all experienced rheumatologists, recruited cases with a firm clinical diagnosis of GCA, based on all available information. Recruitment was retrospective. A positive temporal artery biopsy was not required as it was not always undertaken in classic presentations or could not be performed within an optimal time window or both. In some centres, the erythrocyte sedimentation rate was unavailable and so fulfilment of the 1990 American College of Rheumatology (ACR) criteria [15], which should not be used for clinical diagnosis of GCA [16], was not a requirement for inclusion. Clinical data on a subset of this cohort have already been published [17]. In this analysis, we included all patients who agreed to give a blood sample for genetic studies up to 2012 and where a sample was available. Written informed consent was provided by all patients, and the study was approved by the York Research Ethics Committee (reference 05/Q1108/28).
DNA extraction and genotyping
DNA was extracted from peripheral blood. HLA-DRB1 genotyping was performed by either single-stranded oligonucleotide polymerisation [18] or allele-specific polymerase chain reactions (standard primer sequences (HLA DRBplus Typing Kit, Amersham Biosciences, now part of GE Healthcare, Little Chalfont, UK), except for the forward primer of HLA-DRB1*10 which was redesigned as 5'-GCG GTT GCT GGA AAG ACG CG-3'). Direct sequencing was also performed to enable four-digit genotyping of HLA-DRB1*04 subtypes [18] because of previous reports of a HLA-DRB1*04 association of GCA at the two-digit level. The HaplotypeViewer program was developed to facilitate rapid four-digit genotyping from sequence electropherograms and is freely available [19].
Analysis of genotyping data
Control data from the UK Rheumatoid Arthritis Genetics (UKRAG) Consortium were used for this analysis. Initial logistic regression analyses were undertaken by assuming additive genetic models to estimate the effect of each potential susceptibility/protective allele. Adjustments for genetic effects already proposed in the literature (HLA-DRB1*04) were also performed.
Meta-analysis of giant cell arteritis susceptibility data
To identify case-control studies of HLA-DRB1 association with GCA susceptibility, a literature search was conducted in PubMed, without language restriction, by using the terms “HLA” and “(giant cell arteritis) OR (temporal arteritis)”. Reference lists of studies identified were also scanned. Publications were included if they provided sufficient detail on cases and controls to perform a meta-analysis. Where there were multiple publications with overlapping datasets, the report with the most complete dataset was chosen. Meta-analysis of the published summary carrier frequency data was performed (i.e., assuming a dominant mode of inheritance) because allele frequency and individual-level patient data were mostly unavailable from the authors of the studies. A random-effects model was used; the overall estimate was calculated by using as weights 1/(v i + τ), where v i is the variance of the estimated effect from the ith study and τ is the estimated between-study variance [20].
Worldwide giant cell arteritis incidence in relation to HLA-DRB1*04 population carrier frequencies
To identify reports of the incidence of GCA in different countries, a second literature search in PubMed was conducted with combinations of the medical subject heading terms “giant cell arteritis”, “temporal arteritis”, and “epidemiology”. Hand-searching was also performed in the reference lists of retrieved articles, review articles, and textbooks. Studies were included if they were available in full-text and included an estimate of the annual incidence, time period of the study, method of case definition, population studied, and geographical location of the study. Where necessary, the incidence figure was recalculated as number of new cases per 100,000 of the over-50 population per year. Studies that appeared to report duplicate or overlapping populations were excluded. Studies completing recruitment before 1980 were excluded in case of time trends in the incidence of GCA and because the quality of the reporting was generally lower for the older studies. Where more than one report existed for a single country (unless in ethnically distinct populations), the one with a later average period of recruitment was preferred. Where a single report included two separate sub-studies (regions or time periods), a weighted mean of the two sub-studies was used to arrive at an overall incidence figure.
We then sought data on ethnically matched HLA-DRB1 population allele frequencies at the two- and four-digit levels for each geographical region identified in the second literature search. We considered HLA-DRB1*04 alleles and those identified as being potential susceptibility/protective alleles in our own UK dataset. Methods have been reported elsewhere [21]; briefly, we first consulted the Allele Frequency Net Database [22] and then, if necessary, Ovid Medline and Embase. Carrier frequencies for control populations were converted to estimated allele frequencies by using the Hardy-Weinberg equation. Finally, the following predetermined rule was used to generate an estimate of population HLA-DRB1 allele frequencies: reports with over 500 (four-digit typing) or 1000 (two-digit typing) controls were identified and a weighted mean calculated. In the absence of large studies, studies with more than 100 (four-digit typing) or more than 200 (two-digit typing) were identified and a weighted mean calculated. Determination of geographical latitude and linear regression analysis were performed as previously described [21].
Development of amino acid risk motif model
After determination of susceptibility and protective HLA-DRB1 alleles, amino acid residues in HVR1, HVR2, and HVR3 were obtained from the IMGT (International ImMunoGeneTics Information System) database [23], accessed 31 January 2012). For samples with only two-digit typing, we estimated amino acid residues on the basis of geographically relevant population frequencies of the four-digit subtypes from the Allele Frequency Net Database [22], assigning a probability to residues when they varied within the four-digit subtypes (only necessary at positions 28, 32, 37, 67, 70, 71, and 74). For each of these, the expected misclassification rate when assigned by using population frequencies is less than 1 %, apart from positions 67 (2 %) and 71 (3 %). For each polymorphic position, samples were assigned a dosage (i.e., the expected number of copies) for each residue. Logistic regression was then used to test for association at each position separately, and degrees of freedom were equal to one less than the number of distinct residues. For the most significant positions, forward stepwise regression was used to identify the residues associated with disease risk at that position.
We used population HLA-DRB1 frequencies for inferring amino acid residues in both cases and controls. Under the null hypothesis of no association, the frequencies would be the same, and any bias introduced by using population frequencies for cases would be toward the null. Using HLA-DRB1 four-digit frequencies that have been observed in patients with GCA to infer the amino acid residues in the GCA cases could lead to a biased analysis with inflated false-positive rate.
Results that reach a nominal significance level of 0.05 are highlighted. For the exploratory hypotheses, these should be interpreted in the light of multiple testing. Analyses were performed in SPSS 15 (IBM Corporation, Armonk, NY, USA) and Stata SE (StataCorp LP, College Station, TX, USA).
Results
Patients
Two hundred twenty-five patients with GCA from 7 UK centres consented to analysis of genetic material for this study (125 from Leeds hospitals, including 38 from Otley; 33 from Harrogate; 23 from Southend; 17 from York; 16 from Dewsbury; 10 from Pontefract; and one from Ipswich). Their demographics and disease characteristics, including fulfilment of 1990 ACR criteria, are shown in Table 1. Of the 183 temporal artery biopsies performed, 140 (77 %) were positive. Patients were all European Caucasian.
Analysis of genotyping data
Allele frequencies in cases and 1378 UKRAG controls are shown in Table 2 with per-allele odds ratios with and without adjustment for HLA-DRB1*04, the previously proposed susceptibility allele. Initial analysis was performed at the two-digit level. Four-digit analysis was also performed for the common *04 subtypes, but statistical analysis was not performed on the rarer *04 subtypes.
In 225 patients with GCA and 1378 controls in the novel UK cohort, a susceptibility effect of HLA-DRB1*04 carriage was confirmed (odds ratio (OR) 2.69, 95 % confidence interval (CI) 2.02 to 3.58, P = 1.5×10−11) (Table 3). Possible protective effects from HLA-DRB1*01 and HLA-DRB1*15 were noted, but only HLA-DRB1*01 retained significance after adjusting for HLA-DRB1*04 (Table 2). The data were consistent with a dominant effect of HLA-DRB1*04 (OR for one copy 2.78, 95 % CI 2.07 to 3.72, OR for two copies 1.94, 95 % CI 0.95 to 3.96, compared with no copies).
The effect sizes for HLA-DRB1*04 carriage were similar when restricting analyses to biopsy-positive GCA cases (OR = 2.83, 1.99 to 4.03, P = 7.7×10−9); the number of biopsy-negative GCA cases was too small for a separate analysis.
Meta-analysis of giant cell arteritis susceptibility data
Meta-analysis of previously published data from 14 studies (691 cases, 4038 controls; Table 3) gave ORs of 2.45 (P = 9.2×10−24), 0.78 (P = 0.11), and 0.71 (P = 0.0019) for HLA-DRB1*04, *01, and *02, respectively (*02 is now reclassified as *15 and *16). There was more heterogeneity noted for some alleles than for others; the HLA-DRB1*04 meta-analysis had an I2 statistic of 0 % whereas the meta-analysis for HLA-DRB1*01 had an I2 statistic of 39.4 %. When our UK data was included, the meta-analysis still demonstrated protective effects for HLA-DRB1*01 and HLA-DRB1*02 (P = 0.037 and P = 8.2×10−6, respectively) (Table 3). In the six published articles with information on ethnicity, cases and controls were always stated to be either “white” or “Caucasian”. Therefore, subgroup meta-analysis by ethnic group was not possible.
Worldwide giant cell arteritis incidence in relation to HLA-DRB1*04 population carrier frequencies
Reliable population HLA-DRB1 data were not always available, notably for small, native tribes of Alaska and Saskatoon, where extrapolation from small and physically or genetically (or both) isolated communities was felt to be unwarranted. Table 4 summarises the GCA incidence articles included, together with the estimated population allele frequencies and the number of individuals on which these estimates are based. Substantial clinical heterogeneity was identified in the GCA incidence studies, including variations in the methods used to identify GCA cases and confirm GCA diagnosis.
At the two-digit level (Table 4 and Fig. 1a and b), 17 studies were included in the analysis of worldwide GCA incidence in relation to population HLA-DRB1 allele frequencies. In view of our meta-analysis findings, we extracted data for HLA-DRB1*04, HLA-DRB1*01, and HLA-DRB1*15 population frequencies (Table 4). The majority of these were from Europe and the Mediterranean area. Within this small dataset, HLA-DRB1*15 was more common in the general population at more northerly latitudes (r = 0.52, P = 0.038) whereas no significant association with latitude was seen for population HLA-DRB1*04 or HLA-DRB1*01 (r = 0.47, P = 0.057; r = 0.39, P = 0.133).
Predictors of GCA incidence in univariable analyses were HLA-DRB1*04 population allele frequency (P = 0.001, adjusted R2 = 0.51) and latitude (P = 0.004, adjusted R2 = 0.40), whereas HLA-DRB1*15 was non-significant. In multivariable analysis, both were significant predictors and each made independent contributions to the explanatory power of the model (HLA-DRB1*04, P = 0.008; latitude, P = 0.036; adjusted R2 of the model = 0.62). HLA-DRB1*15 made no additional contribution to the explanatory power of the model.
Development of model with 11-13-33 amino acid risk motif
Tests for association showed that the most significant position was at 13 (P = 1.2×10−9), followed by 11, 33, 37, and 9 (Fig. 1c). At position 13, the most significant residue was H (OR = 2.11, 95 % CI 1.61 to 2.77, P = 5.5×10−8, equivalent to H at 33 and also to the 04 allele). However, stepwise regression found additional contributions from residues S (OR = 1.38, 95 % CI 1.07 to 1.77, P = 0.014) and F (OR = 0.66, 95 % CI 0.44 to 0.99, P = 0.038) (Table 5). Similarly, multiple residues were found at positions 11 and 37. There is very strong linkage disequilibrium in this region, and so many of these residues at different positions almost always occur together, as illustrated in Additional file 1. For example, we did not have the power to distinguish between the risk effects of V, H, and H at positions 11, 13, and 33, respectively, since the *10 allele, which differs at residue 13, is very rare. The previously proposed DYF motif (positions 28, 30, and 31) in HVR2 (OR = 1.54, 95 % CI 1.21 to 1.96, P = 0.00038) did not explain the observed data as well as simple HLA-DRB1*04 carriage. Similarly, variation in amino acid residues within HVR3 is unlikely to explain the observed GCA susceptibility data (Fig. 1c), especially since the other alleles comprising the “RA shared epitope” were not over-represented in GCA.
Discussion
In this study, which includes both new UK data and the first formal meta-analysis of published data on HLA-DRB1 associations of GCA, we not only confirm a strong association of GCA with HLA-DRB1*04 allele carriage, including within our own UK data, but also identify possible protective effects of HLA-DRB1*01 and HLA-DRB1*15, supported by the meta-analysis of previous studies. We were able to impute amino acid residues quite reliably from published allele frequencies, enabling us to analyse amino acid residues even though four-digit typing was not available for every HLA-DRB1 allele. Based on this, it was the amino acid residues 11, 13, and 33 in the first and second hypervariable regions that best explained the observed HLA-DRB1 susceptibility and protective effects, rather than the previously proposed DRYF amino acid motif in the second hypervariable region [12]. We also observed that some non-HLA-DRB1*04 amino acid residues had additional effects (individual amino acid residues that were retained by a multivariable regression model for each separate amino acid position are shown in the last two columns in Table 5), suggesting additional genetic complexities that we did not have the power to investigate in depth. We then systematically extracted data on population prevalence of the identified susceptibility and protective HLA-DRB1 alleles and compared this with reports of GCA incidence in different countries. We found a significant and independent relationship of GCA incidence both with HLA-DRB1*04 and with latitude. Conversely, HLA-DRB1*15 was, if anything, protective and did not contribute to incidence of GCA in the geo-epidemiological study.
Strengths of this work include the presentation of the first UK HLA data in GCA, its presentation in the context of the international literature, the first meta-analysis of HLA-DRB1*04 GCA susceptibility studies, and the novel approach combining a traditional genetic association study with a geo-epidemiology approach. Using logistic regression for the UK, we could control for already-known HLA-DRB1 susceptibility effects in the per-allele analysis, which also has not been performed in other datasets, which mostly reported only carrier frequencies not allele frequencies. Based on this, we were able to suggest HLA-DRB1 amino acid residues that best fit the observed susceptibility/protective allele effects. This is the first synthesis of the literature on reported GCA incidence in relation to population HLA-DRB1 allele frequency and geographical latitude.
Our analysis is based on certain assumptions. Firstly, because many clinicians in the UK do not always request temporal artery biopsy except in cases of diagnostic doubt [24], we had prespecified in the analysis that GCA would be defined clinically rather than limiting inclusion to biopsy-positive cases only. The clinically diagnosed patients, however, had to be firmly diagnosed by an experienced consultant, and there had to be unequivocal clinical features and no alternative explanation for the symptoms after follow-up. Temporal artery biopsy is not 100 % sensitive for GCA; possible reasons for false-negative biopsies in our cohort included delays in obtaining biopsies resulting in resolution of inflammation, suboptimal biopsy length, and biopsy reporting based on the classic pathologic criteria, which may be overly stringent [25]; sometimes the temporal artery is spared in patients with GCA, particularly in those with predominant disease of the aorta and its proximal branches [26]. We conducted a sub-analysis of the biopsy-positive subset and found no difference in the observed effect size for HLA-DRB1*04 association compared with the whole group; with such a small number of biopsy-negative cases, no meaningful statement can be made about the effect size in that group. Our meta-analysis showed that the effect size in the cohort overall was also comparable to that observed in previous reports, some of which included only biopsy-positive cases. If not all the cases truly had GCA, this would have reduced the power of the study (“diluted out” the genetic association) but would have been highly unlikely to introduce artefactual genetic associations because the differential diagnosis of GCA is so wide. Similar pragmatic approaches to case definition for genetics studies, accepting a small, finite rate of misclassification in order to maximise recruitment, have been successfully used in other genetic association studies [27]. We also did not have the power to study whether there are differences in the effect size between regions of the UK, but regional variations in the incidence of diagnosed GCA have been described [28]; it remains unclear how far this is influenced by regional variations in population HLA-DRB1 frequency [29].
In regard to the HLA-DRB1 typing, it is recognised that HLA-DRB1 represents only a small part of the whole MHC and also that not having complete sequence-based four-digit typing may have resulted in some important information being missed. This study focused on HLA-DRB1 and we did not set out to analyse variation elsewhere in the MHC [30]. However, our finding that non-HLA-DRB1*04 residues also contributed significantly to GCA susceptibility/protection (Table 5) suggests that other alleles may also be involved. The MHC is a complex locus with extensive linkage disequilibrium, and an MHC-wide analysis requires larger datasets and specialised analysis methods. A concurrent international, collaborative large-scale genetic analysis of GCA (including samples from this study), using a different genotyping platform (Immunochip) with more extensive coverage of the HLA region [31], shows evidence of wider involvement of the MHC region while confirming the strong association with DRB1*04. Lastly, the literature reviews and meta-analyses are limited by the small number of studies in the literature, many of which were published some years ago, with corresponding variations in case ascertainment and in genotyping assays. Larger datasets using modern genotyping and statistical analysis methods will reveal further GCA susceptibility alleles within the whole HLA locus and allow their pattern of linkage disequilibrium to be analysed.
The P values reported here should be considered in the light of multiple testing, but owing to the a priori suggestion of HLA-DRB1*04 association and lack of consensus as to how to adjust for multiple testing at a multi-allelic locus where the different alleles are not independent of each other, we did not consider a Bonferroni correction to be appropriate here. Nevertheless, model over-fitting is a possibility, and it is essential that our findings be replicated in an independent dataset.
Conclusions
In summary, we report a novel approach to studying genetic influences of disease by combining traditional genetic association studies with geo-epidemiology methods that capitalise on publicly available data. Our new UK data and a synthesis of the published literature suggest that HLA-DRB1*04 might explain part of the observed geographical variation in GCA incidence. This is consistent with an autoimmune aetiology for GCA [32]. However, we found additional variation in susceptibility (Table 5) and incidence (Fig. 1a, b) that is not fully explained by HLA-DRB1*04 and is likely to relate to additional, unknown genetic and environmental factors. Previous studies of GCA have also demonstrated an association between HLA-DRB1*04 and visual loss [33] and also with glucocorticoid resistance [34]. Of interest, in Japan (where HLA-DRB1*04 population frequency is low), large-vessel vasculitis (Takayasu arteritis) is relatively more common than GCA. Takayasu arteritis was associated with alleles containing the 11-13-33 V-H-H motif (HLA-DRB1*0405) in a Turkish population but was not associated with another allele also containing V-H-H (HLA-DRB1*0401) in a European-American population; HLA-DRB1*1502, which was associated with Takayasu arteritis in both populations [35], does not contain the V-H-H motif. Very few patients in our dataset had large-vessel imaging, but genetic characterisation of the subset of GCA patients who have large-vessel involvement or temporal artery sparing or both [26] would be of interest in future studies. From a clinical perspective, further study of well-phenotyped cohorts is required to determine whether HLA-DRB1*04 may serve as a biomarker of pathophysiologically relevant phenotypic disease subsets in order to develop better risk stratification, prediction of response to glucocorticoids, and ultimately targeted therapies.
Abbreviations
- ACR:
-
American College of Rheumatology
- CI:
-
Confidence interval
- GCA:
-
Giant cell arteritis
- HLA:
-
Human leukocyte antigen
- HVR:
-
Hypervariable region
- MHC:
-
Major histocompatibility complex
- OR:
-
Odds ratio
- RA:
-
Rheumatoid arthritis
- UKRAG:
-
UK Rheumatoid Arthritis Genetics
References
Padyukov L, Seielstad M, Ong RT, Ding B, Ronnelid J, Seddighzadeh M, et al. A genome-wide association study suggests contrasting associations in ACPA-positive versus ACPA-negative rheumatoid arthritis. Ann Rheum Dis. 2011;70:259–65.
Lyons PA, Rayner TF, Trivedi S, Holle JU, Watts RA, Jayne DR, et al. Genetically distinct subsets within ANCA-associated vasculitis. N Engl J Med. 2012;367:214–23.
Gregersen PK, Silver J, Winchester RJ. The shared epitope hypothesis. An approach to understanding the molecular genetics of susceptibility to rheumatoid arthritis. Arthritis Rheum. 1987;30:1205–13.
Mackie SL, Taylor JC, Martin SG, Wordsworth P, Steer S, Wilson AG, et al. A spectrum of susceptibility to rheumatoid arthritis within HLA-DRB1: stratification by autoantibody status in a large UK population. Genes Immun. 2012;13:120–8.
Raychaudhuri S, Sandor C, Stahl EA, Freudenberg J, Lee HS, Jia X, et al. Five amino acids in three HLA proteins explain most of the association between MHC and seropositive rheumatoid arthritis. Nat Genet. 2012;44:291–6.
Gran JT. Some thoughts about the etiopathogenesis of temporal arteritis--a review. Scand J Rheumatol. 2002;31:1–5.
Gabriel SE, Michaud K. Epidemiological studies in incidence, prevalence, mortality, and comorbidity of the rheumatic diseases. Arthritis Res Ther. 2009;11:229.
Watts RA, Scott DG. Epidemiology of vasculitis. In: Ball GV, Fessler BJ, Bridges Jr SL, editors. Oxford Textbook of Vasculitis. 3rd ed. Oxford: Oxford University Press; 2014. p. 7–25.
Ball GV, Fessler BJ, Bridges Jr SL. Oxford Textbook of Vasculitis. 3rd ed. Oxford: Oxford University Press; 2014.
Rooney PJ, Rooney J, Balint G, Balint P. Polymyalgia rheumatica: 125 years of progress? Scott Med J. 2014;59:220–8.
Gonzalez-Gay MA, Amoli MM, Garcia-Porrua C, Ollier WE. Genetic markers of disease susceptibility and severity in giant cell arteritis and polymyalgia rheumatica. Semin Arthritis Rheum. 2003;33:38–48.
Weyand CM, Hicok KC, Hunder GG, Goronzy JJ. The HLA-DRB1 locus as a genetic component in giant cell arteritis. Mapping of a disease-linked sequence motif to the antigen binding site of the HLA-DR molecule. J Clin Invest. 1992;90:2355–61.
Martinez-Taboda VM, Bartolome MJ, Lopez-Hoyos M, Blanco R, Mata C, Calvo J, et al. HLA-DRB1 allele distribution in polymyalgia rheumatica and giant cell arteritis: influence on clinical subgroups and prognosis. Semin Arthritis Rheum. 2004;34:454–64.
Heath SC, Gut IG, Brennan P, McKay JD, Bencko V, Fabianova E, et al. Investigation of the fine structure of European populations with applications to disease association studies. Eur J Hum Genet. 2008;16:1413–29.
Hunder GG, Bloch DA, Michel BA, Stevens MB, Arend WP, Calabrese LH, et al. The American College of Rheumatology 1990 criteria for the classification of giant cell arteritis. Arthritis Rheum. 1990;33:1122–8.
Rao JK, Allen NB, Pincus T. Limitations of the 1990 American College of Rheumatology classification criteria in the diagnosis of vasculitis. Ann Intern Med. 1998;129:345–52.
Mackie SL, Dasgupta B, Hordon L, Gough A, Green M, Hollywood J, Dutta S, Bejarano V, Jarrett S, Morgan AW et al. Ischaemic manifestations in giant cell arteritis are associated with area level socio-economic deprivation, but not cardiovascular risk factors. Rheumatology (Oxford, England). 2011;50:2014–22.
Morgan AW, Thomson W, Martin SG, Carter AM, Erlich HA, Barton A, et al. Reevaluation of the interaction between HLA-DRB1 shared epitope alleles, PTPN22, and smoking in determining susceptibility to autoantibody-positive and autoantibody-negative rheumatoid arthritis in a large UK Caucasian population. Arthritis Rheum. 2009;60:2565–76.
Carr I. HaplotypeViewer. http://dna.leeds.ac.uk/HaplotypeViewer/. Accessed 17 Dec 2013.
DerSimonian R, Laird N. Meta-analysis in clinical trials. Control Clin Trials. 1986;7:177–88.
Watts RA, Mackie SL, MacGregor AJ. HLA allele variation as a potential explanation for the geographical distribution of granulomatosis with polyangiitis. Rheumatology (Oxford). 2015;54:359–62.
Gonzalez-Galarza FF, Christmas S, Middleton D, Jones AR. Allele frequency net: a database and online repository for immune gene frequencies in worldwide populations. Nucleic Acids Res. 2011;39:D913–9.
Robinson J, Mistry K, McWilliam H, Lopez R, Parham P, Marsh SGE. The IMGT/HLA Database. Nucleic Acids Res. 2011;39(Suppl 1):D1171–6. http://www.ebi.ac.uk/ipd/imgt/hla. Accessed 31 Jan 2012.
Davies C, Frost B, Eshan O, McLain AD, Shandall A. Temporal artery biopsy…who needs one? Postgrad Med J. 2006;82:476–8.
Zhou L, Luneau K, Weyand CM, Biousse V, Newman NJ, Grossniklaus HE. Clinicopathologic correlations in giant cell arteritis: a retrospective study of 107 cases. Ophthalmology. 2009;116:1574–80.
Brack A, Martinez-Taboada V, Stanson A, Goronzy JJ, Weyand CM. Disease pattern in cranial and large-vessel giant cell arteritis. Arthritis Rheum. 1999;42:311–7.
Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 2007;447:661–78.
Smeeth L, Cook C, Hall AJ. Incidence of diagnosed polymyalgia rheumatica and temporal arteritis in the United Kingdom, 1990–2001. Ann Rheum Dis. 2006;65:1093–8.
Winney B, Boumertit A, Day T, Davison D, Echeta C, Evseeva I, et al. People of the British Isles: preliminary analysis of genotypes and surnames in a UK-control population. Eur J Hum Genet. 2012;20:203–10.
Gonzalez-Gay MA, Rueda B, Vilchez JR, Lopez-Nevot MA, Robledo G, Ruiz MP, et al. Contribution of MHC class I region to genetic susceptibility for giant cell arteritis. Rheumatology (Oxford, England). 2007;46:431–34.
Carmona FD, Mackie SL, Martin J-E, Taylor JC, Vaglio A, Eyre S, et al. A large-scale genetic analysis reveals a strong contribution of the HLA Class II region to giant cell arteritis susceptibility. Am J Hum Genet. 2015;96:565–80.
Weyand CM, Schonberger J, Oppitz U, Hunder NN, Hicok KC, Goronzy JJ. Distinct vascular lesions in giant cell arteritis share identical T cell clonotypes. J Exp Med. 1994;179:951–60.
Gonzalez-Gay MA, Garcia-Porrua C, Llorca J, Hajeer AH, Branas F, Dababneh A, et al. Visual manifestations of giant cell arteritis. Trends and clinical spectrum in 161 patients. Medicine. 2000;79:283–92.
Rauzy O, Fort M, Nourhashemi F, Alric L, Juchet H, Ecoiffier M, et al. Relation between HLA DRB1 alleles and corticosteroid resistance in giant cell arteritis. Ann Rheum Dis. 1998;57:380–2.
Saruhan-Direskeneli G, Hughes T, Aksu K, Keser G, Coit P, Aydin SZ, et al. Identification of multiple genetic susceptibility loci in Takayasu arteritis. Am J Hum Genet. 2013;93:298–305.
Armstrong RD, Panayi GS, Welsh KI. Polymyalgia rheumatica and rheumatoid arthritis: similarity of HLA antigen frequencies. Arthritis Rheum. 1984;27:1438–9.
Bignon JD, Barrier J, Soulillou JP, Martin P, Grolleau JY. HLA DR4 and giant cell arteritis. Tissue Antigens. 1984;24:60–2.
Calamia KT, Moore SB, Elveback LR, Hunder GG. HLA-DR locus antigens in polymyalgia rheumatica and giant cell arteritis. J Rheumatol. 1981;8:993–6.
Cid MC, Ercilla G, Vilaseca J, Sanmarti R, Villalta J, Ingelmo M, et al. Polymyalgia rheumatica: a syndrome associated with HLA-DR4 antigen. Arthritis Rheum. 1988;31:678–82.
Combe B, Sany J, Le Quellec A, Clot J, Eliaou JF. Distribution of HLA-DRB1 alleles of patients with polymyalgia rheumatica and giant cell arteritis in a Mediterranean population. J Rheumatol. 1998;25:94–8.
Dababneh A, Gonzalez-Gay MA, Garcia-Porrua C, Hajeer A, Thomson W, Ollier W. Giant cell arteritis and polymyalgia rheumatica can be differentiated by distinct patterns of HLA class II association. J Rheumatol. 1998;25:2140–5.
Gouet D, Alcalay D, Azais I, Alcalay M, Becq-Giraudon B, Sudre Y, et al. HLA-DR antigens in polymyalgia rheumatica and giant cell arteritis. J Rheumatol. 1985;12:627–8.
Hansen JA, Healey LA, Wilske KR. Association between giant cell (temporal) arteritis and HLA-Cw3. Hum Immunol. 1985;13:193–8.
Jacobsen S, Baslund B, Madsen HO, Tvede N, Svejgaard A, Garred P. Mannose-binding lectin variant alleles and HLA-DR4 alleles are associated with giant cell arteritis. J Rheumatol. 2002;29:2148–53.
Ninet J, Gebuhrer L, Bonvoisin B, Mackiewicz R, Laurent H, Boussuge C, Rousset H, Le Petit JC, Betuel H. [Distribution of HLA-DR antigens in unrelated giant cell arteritis]. Presse medicale (Paris, France : 1983). 1987;16:1725–28.
Richardson JE, Gladman DD, Fam A, Keystone EC. HLA-DR4 in giant cell arteritis: association with polymyalgia rheumatica syndrome. Arthritis Rheum. 1987;30:1293–7.
Salvarani C, Boiardi L, Mantovani V, Ranzi A, Cantini F, Olivieri I, et al. HLA-DRB1, DQA1, and DQB1 alleles associated with giant cell arteritis in northern Italy. J Rheumatol. 1999;26:2395–9.
Weyand CM, Hunder NN, Hicok KC, Hunder GG, Goronzy JJ. HLA-DRB1 alleles in polymyalgia rheumatica, giant cell arteritis, and rheumatoid arthritis. Arthritis Rheum. 1994;37:514–20.
Franzen P, Sutinen S, von Knorring J. Giant cell arteritis and polymyalgia rheumatica in a region of Finland: an epidemiologic, clinical and pathologic study, 1984–1988. J Rheumatol. 1992;19:273–6.
Haimila K, Perasaari J, Linjama T, Koskela S, Saarinenl T, Lauronen J, et al. HLA antigen, allele and haplotype frequencies and their use in virtual panel reactive antigen calculations in the Finnish population. Tissue Antigens. 2013;81:35–43.
Nordborg C, Johansson H, Petursdottir V, Nordborg E. The epidemiology of biopsy-positive giant cell arteritis: special reference to changes in the age of the population. Rheumatology (Oxford, England). 2003;42:549–52.
Grunewald J, Brynedal B, Darlington P, Nisell M, Cederlund K, Hillert J, et al. Different HLA-DRB1 allele distributions in distinct clinical subgroups of sarcoidosis patients. Respir Res. 2010;11:25.
Hedstrom AK, Sundqvist E, Baarnhielm M, Nordin N, Hillert J, Kockum I, et al. Smoking and two human leukocyte antigen genes interact to increase the risk for multiple sclerosis. Brain. 2011;134:653–64.
Bjorkbacka H, Lavant EH, Fredrikson GN, Melander O, Berglund G, Carlson JA, et al. Weak associations between human leucocyte antigen genotype and acute myocardial infarction. J Intern Med. 2010;268:50–8.
van der Woude D, Lie BA, Lundstrom E, Balsa A, Feitsma AL, Houwing-Duistermaat JJ, et al. Protection against anti-citrullinated protein antibody-positive rheumatoid arthritis is predominantly associated with HLA-DRB1*1301: a meta-analysis of HLA-DRB1 associations with anti-citrullinated protein antibody-positive and anti-citrullinated protein antibody-negative rheumatoid arthritis in four European populations. Arthritis Rheum. 2010;62:1236–45.
Castro FA, Haimila K, Sareneva I, Schmitt M, Lorenzo J, Kunkel N, et al. Association of HLA-DRB1, interleukin-6 and cyclin D1 polymorphisms with cervical cancer in the Swedish population--a candidate gene approach. Int J Canc Suppl J Int Canc Suppl. 2009;125:1851–8.
Salvarani C, Macchioni P, Zizzi F, Mantovani W, Rossi F, Castri C, et al. Epidemiologic and immunogenetic aspects of polymyalgia rheumatica and giant cell arteritis in northern Italy. Arthritis Rheum. 1991;34:351–6.
Rendine S, Ferrero NM, Sacchi N, Costa C, Pollichieni S, Amoroso A. Estimation of human leukocyte antigen class I and class II high-resolution allele and haplotype frequencies in the Italian population and comparison with other European populations. Hum Immunol. 2012;73:399–404.
Invernizzi P, Selmi C, Poli F, Frison S, Floreani A, Alvaro D, Almasio P, Rosina F, Marzioni M, Fabris L et al. Human leukocyte antigen polymorphisms in Italian primary biliary cirrhosis: a multicenter study of 664 patients and 1992 healthy controls. Hepatology (Baltimore, Md). 2008:48:1906–12.
Gonzalez-Gay MA, Garcia-Porrua C, Rivas MJ, Rodriguez-Ledo P, Llorca J. Epidemiology of biopsy proven giant cell arteritis in northwestern Spain: trend over an 18 year period. Ann Rheum Dis. 2001;60:367–71.
Alcoceba M, Marin L, Balanzategui A, Sarasquete ME, Chillon MC, Martin-Jimenez P, et al. Frequency of HLA-A, −B and -DRB1 specificities and haplotypic associations in the population of Castilla y Leon (northwest-central Spain). Tissue Antigens. 2011;78:249–55.
Romero-Pinel L, Pujal JM, Martinez-Yelamos S, Gubieras L, Matas E, Bau L, et al. HLA-DRB1: genetic susceptibility and disability progression in a Spanish multiple sclerosis population. Eur J Neurol. 2011;18:337–42.
Elling P, Olsson AT, Elling H. Synchronous variations of the incidence of temporal arteritis and polymyalgia rheumatica in different regions of Denmark; association with epidemics of Mycoplasma pneumoniae infection. J Rheumatol. 1996;23:112–9.
Kruse C, Steffensen R, Varming K, Christiansen OB. A study of HLA-DR and -DQ alleles in 588 patients and 562 controls confirms that HLA-DRB1*03 is associated with recurrent miscarriage. Human reproduction (Oxford, England). 2004;19:1215–21.
Baldursson O, Steinsson K, Bjornsson J, Lie JT. Giant cell arteritis in Iceland. An epidemiologic and histopathologic analysis. Arthritis Rheum. 1994;37:1007–12.
Jonsson T, Gudmundsson KG, Bjarnadottir K, Hjalmarsdotti IB, Gudmundsson S, Arngrimsson R. Association of HLA-DRB1*01 with Dupuytren’s disease. Scand J Rheumatol. 2013;42:45–7.
Haugeberg G, Irgens KA, Thomsen RS. No major differences in incidence of temporal arteritis in northern and western Norway compared with reports from southern Norway. Scand J Rheumatol. 2003;32:318–9.
Munthe-Kaas MC, Carlsen KL, Carlsen KH, Egeland T, Haland G, Devulapalli CS, et al. HLA Dr-Dq haplotypes and the TNFA-308 polymorphism: associations with asthma and allergy. Allergy. 2007;62:991–8.
Karlsen TH, Boberg KM, Vatn M, Bergquist A, Hampe J, Schrumpf E, et al. Different HLA class II associations in ulcerative colitis patients with and without primary sclerosing cholangitis. Genes Immun. 2007;8:275–8.
Reinhold-Keller E, Herlyn K, Wagner-Bastmeyer R, Gross WL. Stable incidence of primary systemic vasculitides over five years: results from the German vasculitis register. Arthritis Rheum. 2005;53:93–9.
Reinshagen M, Loeliger C, Kuehnl P, Weiss U, Manfras BJ, Adler G, et al. HLA class II gene frequencies in Crohn’s disease: a population based analysis in Germany. Gut. 1996;38:538–42.
Dadoniene J, Kirdaite G, Mackiewicz Z, Rimkevicius A, Haugeberg G. Incidence of primary systemic vasculitides in Vilnius: a university hospital population based study. Ann Rheum Dis. 2005;64:335–6.
Ploski R, Butrimiene I, Kaminska E, Valiukiene K, Sliwinska P, Kubasiewicz E, et al. Rheumatoid arthritis in Poland and Lithuania: different clinical course and HLA associations despite similar genetic background. Ann Rheum Dis. 2005;64:165–6.
Bas-Lando M, Breuer GS, Berkun Y, Mates M, Sonnenblick M, Nesher G. The incidence of giant cell arteritis in Jerusalem over a 25-year period: annual and seasonal fluctuations. Clin Exp Rheumatol. 2007;25:S15–7.
Roitberg-Tambur A, Witt CS, Friedmann A, Safirman C, Sherman L, Battat S, et al. Comparative analysis of HLA polymorphism at the serologic and molecular level in Moroccan and Ashkenazi Jews. Tissue Antigens. 1995;46:104–10.
Sonnenblick M, Nesher G, Friedlander Y, Rubinow A. Giant cell arteritis in Jerusalem: a 12-year epidemiological study. Br J Rheumatol. 1994;33:938–41.
Kwon OJ, Brautbar C, Weintrob N, Sprecher E, Saphirman C, Bloch K, et al. Immunogenetics of HLA class II in Israeli Ashkenazi Jewish, Israeli non-Ashkenazi Jewish, and in Israeli Arab IDDM patients. Hum Immunol. 2001;62:85–91.
Zeitlin AA, Heward JM, Newby PR, Carr-Smith JD, Franklyn JA, Gough SC, et al. Analysis of HLA class II genes in Hashimoto’s thyroiditis reveals differences compared to Graves’ disease. Genes Immun. 2008;9:358–63.
Pamuk ON, Donmez S, Karahan B, Pamuk GE, Cakir N. Giant cell arteritis and polymyalgia rheumatica in northwestern Turkey: Clinical features and epidemiological data. Clin Exp Rheumatol. 2009;27:830–3.
Ucar F, Capkin E, Karkucak M, Yucel B, Sonmez M, Alver A, et al. Associations of HLA-DRB1 alleles with anti-citrullinated protein antibody-positive and anti-citrullinated protein antibody-negative rheumatoid arthritis in northern east part of Turkey. Int J Rheum Dis. 2012;15:538–45.
Saruhan-Direskeneli G, Uyar FA, Bakar S, Eraksoy M. Molecular analysis of HLA-DRB1, −DQA1 and -DQB1 polymorphism in Turkey. Tissue Antigens. 2000;55:171–4.
Kinikli G, Ates A, Turgay M, Akay G, Kinikli S, Tokgoz G. HLA-DRB1 genes and disease severity in rheumatoid arthritis in Turkey. Scand J Rheumatol. 2003;32:277–80.
Kermani TA, Schafer VS, Crowson CS, Hunder GG, Gabriel SE, Matteson EL, et al. Increase in age at onset of giant cell arteritis: a population-based study. Ann Rheum Dis. 2010;69:780–1.
Smith CA, Fidler WJ, Pinals RS. The epidemiology of giant cell arteritis. Report of a ten-year study in Shelby County, Tennessee. Arthritis Rheum. 1983;26:1214–9.
Leech NJ, Kitabchi AE, Gaur LK, Hagopian WA, Hansen J, Burghen GA, et al. Genetic and immunological markers of insulin dependent diabetes in Black Americans. Autoimmunity. 1995;22:27–32.
Ramstead CL, Patel AD. Giant cell arteritis in a neuro-ophthalmology clinic in Saskatoon, 1998–2003. Can J Ophthalmol. 2007;42:295–8.
Gladman DD, Kung TN, Siannis F, Pellett F, Farewell VT, Lee P. HLA markers for susceptibility and expression in scleroderma. J Rheumatol. 2005;32:1481–7.
Silverberg MS, Mirea L, Bull SB, Murphy JE, Steinhart AH, Greenberg GR, et al. A population- and family-based study of Canadian families reveals association of HLA DRB1*0103 with colonic involvement in inflammatory bowel disease. Inflamm Bowel Dis. 2003;9:1–9.
Acknowledgements
Earlier versions of this work, including the proposal about the contribution of the HVR1 region, were presented at the UKIVAS Lockwood Club in Manchester, May 2013, and at the Northern and Yorkshire Combined Regional Meeting in York, September 2013; many audience members provided constructive feedback and comments at that time. Membership of the UK GCA Consortium at the time these patients were recruited [17] and UKRAG Consortium [18] can be found in the references cited and in Additional file 2. We are grateful to all those who recruited patients at these sites, including Shouma Dutta and Jane Hollywood. We would like to thank Ian Carr for writing the HaplotypeViewer program, specifically for this study. Thanks to Clive Edelsten, ophthalmologist, Ipswich Hospital, for discussions about Scandinavian contributions to UK ancestry.
This work was supported by the Wellcome Trust/Academy of Medical Sciences Starter Grants for Clinical Lecturers scheme. SLM was funded by a Clinical Research Fellowship from Research into Ageing, then a Clinical Lectureship from the National Institute for Health Research (NIHR), and is currently funded by a Clinician Scientist Award from the NIHR. AWM, JCT, and JHB are supported in part by the NIHR-Leeds Musculoskeletal Biomedical Research Unit. The views expressed are those of the authors and not necessarily those of the National Health Service, the NIHR, or the Department of Health.
Grant funding
Grant funding was provided by the NIHR and the Wellcome Trust.
Author information
Authors and Affiliations
Consortia
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
SLM designed the study and wrote the case record form, recruited patients, and validated clinical data, including diagnoses, DNA extractions, and genotyping assays, and performed the searches and data extraction from the published literature and online databases for HLA-DRB1 data (GCA susceptibility, amino acid sequences, and population allele frequencies), participated in design of geo-epidemiology analyses and carried out the geo-epidemiology statistical analyses, and drafted the manuscript. JCT carried out all the other statistical analyses, including amino acid analyses, and helped to draft the manuscript. LH-R and SM participated in design and validation of the genotyping assays, including sequencing and analysis of electropherograms, including use of HaploViewer. BD, AG, MG, LH, SJ, and CTP recruited patients and validated clinical data, including diagnoses. UKRAG Consortium provided control HLA-DRB1 data. JHB supervised the statistical analyses, including amino acid analyses, and helped to draft the manuscript. RW supervised the design and conduct of the geo-epidemiology analyses, performed the literature search for incidence of GCA in different countries, and helped to draft the manuscript. AWM conceived of the study, supervised the design and conduct of the genetic studies, including amino acid analyses, and helped to draft the manuscript. All authors revised the draft for important intellectual content and read and approved the final manuscript.
Richard Watts and Ann W. Morgan contributed equally to this work.
Additional files
Additional file 1:
Appendices I and II (Membership of consortia). Description of contents: Names and institutional affiliations of members of the UK GCA Consortium (Appendix I) and the UKRAG Consortium (Appendix II).
Additional file 2:
Amino acids at hypervariable regions (HVR) in susceptibility, protective, and neutral GCA alleles. Table showing the amino acid residues in the hypervariable regions of HLA-DRB1, for each of the HLA-DRB1 alleles. Shaded columns denote the proposed 11-13-33 GCA risk motif.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.
The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.
To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Mackie, S.L., Taylor, J.C., Haroon-Rashid, L. et al. Association of HLA-DRB1 amino acid residues with giant cell arteritis: genetic association study, meta-analysis and geo-epidemiological investigation. Arthritis Res Ther 17, 195 (2015). https://doi.org/10.1186/s13075-015-0692-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13075-015-0692-4