Autoimmune diseases are complex and have genetic and environmental susceptibility factors. The objective was to test the genetic association of systemic lupus erythematosus (SLE) and anti-neutrophil cytoplasmic antibody (ANCA) – associated systemic vasculitis (AAV) with SNPs in the IL2RA region and to correlate genotype with serum levels of IL-2RA.
Using a cohort of over 700 AAV patients, two SLE case-control studies and an SLE trio collection (totalling over 1000 SLE patients), and a TaqMan genotyping approach, we tested 3 SNPs in the IL2RA locus, rs11594656, rs2104286 & rs41295061, each with a prior association with autoimmune disease; rs11594656 and rs41295061 with type 1 diabetes (T1D) and rs2104286 with multiple sclerosis (MS) and T1D.
We show that SLE is associated with rs11594656 (P = 3.87 × 10-7) and there is some evidence of association of rs41295061 with AAV (P = 0.0122), which both have prior association with T1D. rs2104286, an MS and T1D – associated SNP in the IL2RA locus, is not associated with either SLE or AAV.
We have confirmed a previous suggestion that the IL2RA locus is associated with SLE and showed some evidence of association with AAV. Soluble IL-2RA concentrations correlate with rs11594656 genotype in quiescent disease in both AAV and SLE. Differential association of autoimmune diseases and SNPs within the IL2RA locus suggests that the IL2RA pathway may prove to play differing, as yet undefined, roles in each disease.
Most autoimmune diseases are polygenic, with genetic susceptibility being conferred by common variants at a number of loci. It is evident that some susceptibility loci are shared between autoimmune diseases. One locus that is associated with several diseases including T1D, MS and Graves' disease, is IL2RA [1-4].
The IL-2 pathway is critical for the maintenance of peripheral tolerance. Activated T cells and regulatory T cells express IL-2RA, the high affinity IL-2 receptor (or CD25). IL-2RA expression has also been reported on activated B cells, activated monocytes, and NK cells . IL-2 signalling supports the division and survival of effector T cells, plays a role in activation-induced cell death and is vital to regulatory T cell homeostasis . Its physiological role on other IL-2RA expressing cell types, however, is unclear. One human patient, homozygous for a 4 base pair exonic deletion in IL2RA, lacking a functional IL-2RA demonstrated increased susceptibility to viral, bacterial and fungal infections, but unlike Il2ra-deficient mice did not develop overt autoimmune disease . Thus, the IL-2 pathway has a profound effect on the T cell immune response, contributing to both activation and regulation depending on the immunological context.
In T1D, two disease-associated SNPs [rs11594656 and rs41295061] are required to explain association with this region of chromosome 10p15, indicating that there are multiple functional alleles of the IL2RA locus, presumably altering expression of the receptor [1,2]. In addition, associations have been established with multiple sclerosis (MS) [rs2104286] , and suggested with Graves' disease . The SNP, rs2104286, has also been associated recently with T1D, representing a third independent effect for the IL2RA region in this autoimmune disease . In systemic lupus erythematosus (SLE), SNPs without any prior autoimmune associations [rs7072793, rs4147539, rs7090530 and rs12251307] were examined in a candidate gene analysis of IL2RA, providing weak evidence for association . Soluble IL-2RA (sIL-2RA) is shed from human T cells following in vitro activation in response to a variety of polyclonal stimuli . Significantly elevated serum sIL-2RA has been described in infection, malignancy and in autoimmune disease, including T1D, rheumatoid arthritis, SLE and Wegener's granulomatosis [5,10], although the functional implications of this remain unclear. Genotype at rs11594656 and rs41295061 , and most recently with rs2104286 , has been shown to correlate with serum sIL-2RA concentrations. We therefore genotyped autoimmune disease-prone variants of the IL2RA region (rs11594656, rs41595061 and rs2104286) in two distinct autoimmune diseases associated with soluble IL-2RA elevation, using three independent SLE cohorts and one AAV cohort. SLE is a polygenic multi-system disease characterised by immune complex formation and deposition, frequently associated with anti-dsDNA antibodies and complement dysregulation. The clinical features vary between patients, both in severity and system involvement, and include joint, skin, renal and neurological manifestations. AAV can be divided into three clinical syndromes (Wegener's granulomatosis, microscopic polyangiitis and Churg-Strauss syndrome). All are characterised by inflammation of medium to small calibre blood vessels, and the presence of anti-neutrophil cytoplasmic antibody (ANCA). The clinical presentation of AAV is variable between patients. Severe manifestations include acute glomerulonephritis and granulomatous inflammation of the upper and lower respiratory tract. We also assessed serum titres of sIL-2RA in both diseases.
The AAV cohort (n = 744) comprises subjects from four sources, all meeting the Chapel Hill diagnostic criteria :
1. The MRC/Kidney Research UK National DNA Bank for Glomerulonephritis (KRUK-AAV). Individuals were between the ages of 18 and 70 years, were ANCA seropositive, and had biopsy-proven necrotizing glomerulonephritis.
2. The UK vasculitis cohort 2 was recruited from 9 centres in the UK and comprised patients seropositive for ANCA and/or with histological evidence of small vessel vasculitis.
3. Patients recruited from the University of Birmingham. All individuals were ANCA seropositive with firm clinical and/or histological evidence of vasculitis.
4. The Lupus and Vasculitis Service, Addenbrooke's Hospital, Cambridge. All individuals were ANCA seropositive with firm clinical and/or histological evidence of vasculitis.
The KRUK-SLE cohort (n = 220) was obtained from the MRC/Kidney Research UK National DNA Bank for Glomerulonephritis (KRUK) and the Lupus and Vasculitis Service, Addenbrooke's Hospital, Cambridge. All individuals were between the ages of 18 and 50 years with a definite diagnosis of lupus nephritis based on biopsy and on the clinical and serological features defined by the American College of Rheumatology (ACR) . The Imperial case-control cohort comprised 480 SLE patients and there were 471 trios in the Imperial family cohort (ie 471 affected sons/daughters, 942 parents) . Both Imperial cohorts used the ACR defined clinical and serological features. Since the KRUK-SLE cohort used the additional inclusion criterion of biopsy-proven renal involvement, it was analysed separately from the Imperial case-control cohort.
Genotyping was performed using TaqMan genotyping kits (Applied Biosystems) for each SNP, with fluorescence data captured using an ABI 7900HT (Applied Biosystems) after 40 cycles of PCR.
Assessment of disease activity
Birmingham Vasculitis Activity Scores (BVAS) were calculated for AAV patients seen by the Lupus and Vasculitis Service at Addenbrooke's hospital. BVAS assesses both constitutional and organ-specific symptoms, determining whether they are new or pre-existing . British Isles Lupus Assessment Group (BILAG) scores were used to assess disease activity in SLE patients seen by the Lupus and Vasculitis Service at Addenbrooke's hospital. BILAG scores assess 9 different organ systems based upon a physician's intention to treat . A flare was defined when three criteria were met: 1. Worsening disease activity scores (at least 1 major or 3 minor BVAS criteria for AAV; a new BILAG score A or B in any system for SLE), 2. clinical impression of active disease and 3. an increase in immunosuppressive treatment.
Control genotypes for 9115 individuals from the British 1958 Birth Cohort and UK Blood Service were obtained from the Juvenile Diabetes Research Foundation/Wellcome Trust Diabetes and Inflammation Laboratory . This cohort is an expansion of a dataset previously shown to be appropriate for us as UK-wide controls . The 1958 Birth Cohort DNA was collected as part of an ongoing study following all births in England, Scotland and Wales in one week in 1958 http://www.b58cgene.sgul.ac.uk webcite. For the SLE analysis, controls were split into two groups and one group used as controls for each case cohort: n = 4534 controls for the KRUK cohort and n = 4581 for the Imperial cohort. Some of these control individuals were used in the earlier T1D studies [1,2].
Serum soluble IL-2RA was measured in triplicate using an ELISA, according to the manufacturer's instructions (R&D Systems). Serum from both AAV and SLE patients was collected at presentation with a 'flare' (time 0, see Assessment of disease activity section) and three months later (3 months) by the Lupus and Vasculitis Service, Addenbrooke's Hospital, Cambridge.
Statistical analysis was performed using Prism (GraphPad), and R http://www.r-project.org webcite. Case control analyses were performed using χ2 tests for significance on 2 × 2 contingency tables. Odds ratios and 95% confidence intervals were calculated from the same 2 × 2 tables. The Imperial trio data was analysed using the transmission disequilibrium test (TDT). ELISA data were tested for significance using linear regression modelling or the Mann-Whitney test, as indicated in the figure legends.
This study was approved by the Cambridge Local Research Ethics Committee and by the Oversight Committee of the KRUK DNA Bank.
In all three of our SLE cohorts, rs11594656 was found to be significantly associated with disease (protective minor allele, as seen in T1D), combined P = 3.87 × 10-7 (KRUK case control P = 1 × 10-4, Imperial case control P = 0.0221 and Imperial family trio TDT P = 2.8 × 10-3) (Figure 1 and Table 1). Neither of the other two SNPs were associated with SLE (Figure 1 and Table 1). In AAV, in contrast, some evidence for an association was found with rs41295061 (P = 0.0122; minor allele odds ratio = 0.7739; 95% OR CI 0.6330 – 0.9461). The minor allele is protective as in T1D. However, neither rs2104286 or rs11594656 showed any association with AAV (Figure 1 and Table 1). At 80% power, we can detect effect sizes larger than 1.25 for rs11594656 in AAV and 1.4 for rs41295061 in SLE, across all the SLE cases and controls.
Figure 1. SLE and AAV are associated with different SNPs in the IL2RA region. The genomic context of rs2104286, rs41295061 and rs11594656 is shown, using the GBrowse implementation in T1DBase http://www.t1dbase.org webcite. All are non-coding SNPs. Negative decimal logarithms of P values (χ2 test; 1 degree of freedom) from each cohort from SLE and AAV are plotted for rs41295061 and rs11594656 but not for rs2104286, as an association of this MS-associated SNP was not detected with SLE or AAV (table 1). Open circles represent the KRUK case-control cohort; filled circles the Imperial case-control cohort and open squares the Imperial family trios (P value from TDT). The filled square represents the combined P value for rs11594656 in SLE. Crosses are used for the AAV cohort. The dashed lines correspond to P = 0.05.
Table 1. Genotype data at rs11594656, rs41295061 and rs2104286 for SLE and AAV
In both SLE and AAV, sIL-2RA concentrations were significantly higher than controls (Figure 2a). In SLE, disease activity as assessed by BILAG score did not show a significant correlation with sIL-2RA concentrations (Figure 2b), despite previous reports of a correlation between sIL-2RA concentrations and disease activity . However, in AAV, disease activity as assessed by BVAS showed a positive correlation with sIL-2RA concentration (Figure 2c, ). Three months after commencement of treatment, when both diseases were considerably less active (mean BILAG decreased from 16.1 to 4.7 and mean BVAS from 13.9 to 4.4), sIL-2RA concentrations correlated with the SLE associated rs11594656 genotype (Figures 2d &2e). In active AAV this correlation no longer held true (figure 2e).
Figure 2. Serum soluble IL-2RA in SLE and AAV and its relationship with disease activity and genotype at rs11594656. (a) A comparison of SLE and AAV soluble IL-2RA (sIL-2RA) at time 0 and 3 months compared with controls, as assessed by ELISA. Each dot represents a patient. (b, c) Correlation between disease activity and sIL-2RA. (BILAG for SLE; BVAS for AAV; see supplementary information). r2 and P values from linear regression modelling are shown. (d, e) sIL-2RA levels stratified by genotype at rs11594656 at both time 0 and 3 months for SLE (d) and AAV (e). P values for panels (a), (d) & (e) were generated using the Mann-Whitney test.
We have shown a novel IL2RA association with AAV and confirm the IL2RA locus association with SLE. This brings the number of autoimmune diseases associated with this gene to five, the others being T1D, MS and Graves' disease [1-4]. An association of SLE with IL2RA (rs7072793, P = 0.01) was generated by candidate gene re-analysis of genome-wide association data from 720 female patients . The present study uses 1130 SLE genotypes, of mixed gender, in three independent cohorts, and demonstrates a strong association with the IL2RA region (rs11594656, P = 3.87 × 10-7). Owing to strong co-segregation of these disorders in cases and families, susceptibility loci can be shared across different autoimmune diseases. Several examples at the gene-specific level are known, including CTLA-4  and PTPN22 . Our study underlines that this sharing can occur at the single gene level for IL2RA for T1D and SLE, two diseases in which co-segregation is less well recognised.
Our findings demonstrate that the association between IL2RA and autoimmune disease is not a simple one. SLE is associated with rs11594656, but not rs41295061 and in AAV there is some evidence for an association with rs41295061 but not with rs11594656. Moreover, in T1D, all three SNPs are associated with disease, but the strength of association varies with each (rs41295061, P = 6.43 × 10-25; rs2104286, P = 1.27 × 10-13; rs11594656, P = 3.37 × 10-6) [2,8]. In MS, whilst rs2104286 is associated with disease , rs11594656 or rs41295061 are not (rs11594656, P = 0.293; rs41295061, P = 0.625) .
There are differences at the functional level as well, with rs11594656 genotype correlated with sIL-2RA levels in both SLE and AAV despite only having genetic association with SLE. This correlation is present in both SLE and AAV three months following disease flare, but absent at the time of flare itself, perhaps because local inflammatory mechanisms provide a major contribution to sIL-2RA, overriding the influence of genotype. Thus the degree of autoimmune disease activity needs to be taken into account in studies attempting to correlate genotype with markers such as sIL-2RA.
There is precedent for such contrasting functional and genetic associations between SLE and AAV. FcγRIIIb, a neutrophil-specific low affinity activatory Fc receptor, has two allotypic variants, which have a significant effect on the affinity of the receptor for IgG. One of these, FCGR3B*NA1, is associated with SLE, and the other, FCGR3B*NA2, with AAV [20,21]. Moreover, copy number at FCGR3B is associated with expression and function of the receptor , and low copy number is associated with SLE  but not with AAV ( & data not shown). The association of a single gene with different diseases indicates the importance in disease of the biological pathway subserved by that gene, but it does not imply that defects in the gene lead to a similar functional perturbation in each case. Thus the same gene, IL2RA, is associated with both SLE and AAV, but may play quite distinct roles in these diseases.
We confirm the association of the IL2RA locus with SLE, by genotyping rs11594656 across three independent cohorts, using two different study designs. We report evidence for a novel association of the IL2RA locus with AAV. The MS-associated SNP, rs2104286, is not associated with either disease in our cohorts. Soluble IL-2RA concentrations were elevated in SLE and AAV. After three months of treatment, soluble IL-2RA concentrations in both diseases correlated with the SLE associated rs11594656 genotype. The contrasting genetic and functional associations found between SLE and AAV and the IL2RA locus presumably reflect differences in disease pathogenesis.
The authors declare that they have no competing interests.
KGCS and PAL conceived and designed the study. AW, TJV, JAT, KL & RAW provided access to DNA samples, and TJV & JAT assisted in data interpretation. EJC performed the SNP genotyping and its analysis. Additional genotyping and experimental oversight was provided by CEL. MRC contributed the ELISA data and its analysis. The manuscript was written by EJC & KGCS. Additionally, MRC, CEL, TJV, PAL & JAT actively contributed to discussion of the results of this manuscript. All authors read and approved the final manuscript.
We acknowledge the participation of all of the patients, control subjects and family members. We acknowledge use of DNA from the MRC/Kidney Research UK National DNA Bank for Glomerulonephritis and Neil Walker for his curation and supply of data from the British 1958 Birth Cohort. We thank Jane Hollis and Kirsty Townsend for serum sample collections, David Jayne and Afzal Chaudry for their clinical input, and Caroline Savage, Lorraine Harper and Dwomoa Adu for contributions to DNA collections.
This work was funded by the NIHR Cambridge Biomedical Research Centre, the Medical Research Council (MRC), the Wellcome Trust (Programme Grant Number 083650/Z/07/Z), and a combined Wellcome Trust – Juvenile Diabetes Research Foundation grant. The Cambridge Institute for Medical Research is in receipt of a Wellcome Trust Strategic Award (079895). EJC was funded by a MRC Doctoral Training Account studentship.
Vella A, Cooper JD, Lowe CE, Walker N, Nutland S, Widmer B, Jones R, Ring SM, McArdle W, Pembrey ME, Strachan D, Dunger DB, Twells RC, Clayton DG, Todd JA: Localization of a type 1 diabetes locus in the IL2RA/CD25 region by use of tag single-nucleotide polymorphisms.
Lowe CE, Cooper JD, Brusko T, Walker NM, Smyth DJ, Bailey R, Bourget K, Plagnol V, Field S, Atkinson M, Clayton DG, Wicker LS, Todd JA: Large-scale genetic fine mapping and genotype-phenotype associations implicate polymorphism in the IL2RA region in type 1 diabetes.
Brand OJ, Lowe CE, Heward JM, Franklyn JA, Cooper JD, Todd JA, Gough SC: Association of the interleukin-2 receptor alpha (IL-2Ralpha)/CD25 gene region with Graves' disease using a multilocus test and tag SNPs.
Maier LM, Lowe CE, Cooper J, Downes K, Anderson DE, Severson C, Clark PM, Healy B, Walker N, Aubin C, Oksenberg JR, Hauser SL, Compston A, Sawcer S, Consortium TIMSG, Jager PLD, Wicker LS, Todd JA, Hafler DA: IL2RA genetic heterogeneity in multiple sclerosis and type 1 diabetes susceptibility and soluble interleukin-2 receptor production.
The International Consortium for Systemic Lupus Erythematosus Genetics: Genome-wide association scan in women with systemic lupus erythematosus identifies susceptibility variants in ITGAM, PXK, KIAA1542 and other loci.
Jennette JC, Falk RJ, Andrassy K, Bacon PA, Churg J, Gross WL, Hagen EC, Hoffman GS, Hunder GG, Kallenberg CGM, Mccluskey RT, Sinico RA, Rees AJ, Van Es LA, Waldherr R, Wiik A: Nomenclature of Systemic Vasculitides.
Graham DS, Graham RR, Manku H, Wong AK, Whittaker JC, Gaffney PM, Moser KL, Rioux JD, Altshuler D, Behrens TW, Vyse TJ: Polymorphism at the TNF superfamily gene TNFSF4 confers susceptibility to systemic lupus erythematosus.
Stone JH, Hoffman GS, Merkel PA, Min YI, Uhlfelder ML, Hellmann DB, Specks U, Allen NB, Davis JC, Spiera RF, Calabrese LH, Wigley FM, Maiden N, Valente RM, Niles JL, Fye KH, McCune JW, St Clair EW, Luqmani RA: A disease-specific activity index for Wegener's granulomatosis: modification of the Birmingham Vasculitis Activity Score. International Network for the Study of the Systemic Vasculitides (INSSYS).
Isenberg DA, Rahman A, Allen E, Farewell V, Akil M, Bruce IN, D'Cruz D, Griffiths B, Khamashta M, Maddison P, McHugh N, Snaith M, Teh LS, Yee CS, Zoma A, Gordon C: BILAG 2004. Development and initial validation of an updated version of the British Isles Lupus Assessment Group's disease activity index for patients with systemic lupus erythematosus.
Todd JA, Walker NM, Cooper JD, Smyth DJ, Downes K, Plagnol V, Bailey R, Nejentsev S, Field SF, Payne F, Lowe CE, Szeszko JS, Hafler JP, Zeitels L, Yang JH, Vella A, Nutland S, Stevens HE, Schuilenburg H, Coleman G, Maisuria M, Meadows W, Smink LJ, Healy B, Burren OS, Lam AA, Ovington NR, Allen J, Adlem E, Leung HT, et al.: Robust associations of four new chromosome regions from genome-wide analyses of type 1 diabetes.
Ueda H, Howson JM, Esposito L, Heward J, Snook H, Chamberlain G, Rainbow DB, Hunter KM, Smith AN, Di Genova G, Herr MH, Dahlman I, Payne F, Smyth D, Lowe C, Twells RC, Howlett S, Healy B, Nutland S, Rance HE, Everett V, Smink LJ, Lam AC, Cordell HJ, Walker NM, Bordin C, Hulme J, Motzo C, Cucca F, Hess JF, et al.: Association of the T-cell regulatory gene CTLA4 with susceptibility to autoimmune disease.
Siriboonrit U, Tsuchiya N, Sirikong M, Kyogoku C, Bejrachandra S, Suthipinittharm P, Luangtrakool K, Srinak D, Thongpradit R, Fujiwara K, Chandanayingyong D, Tokunaga K: Association of Fcgamma receptor IIb and IIIb polymorphisms with susceptibility to systemic lupus erythematosus in Thais.
Willcocks LC, Lyons PA, Clatworthy MR, Robinson JI, Yang W, Newland SA, Plagnol V, McGovern NN, Condliffe AM, Chilvers ER, Adu D, Jolly EC, Watts R, Lau YL, Morgan AW, Nash G, Smith KGC: Copy number of FCGR3B, which is associated with systemic lupus erythematosus, correlates with protein expression and immune complex uptake.
Fanciulli M, Norsworthy PJ, Petretto E, Dong R, Harper L, Kamesh L, Heward JM, Gough SCL, de Smith A, Blakemore AIF, Froguel P, Owen CJ, Pearce SHS, Teixeira L, Guillevin L, Graham DSC, Pusey CD, Cook HT, Vyse TJ, Aitman TJ: FCGR3B copy number variation is associated with susceptibility to systemic, but not organ-specific, autoimmunity.
The pre-publication history for this paper can be accessed here: