Investigation of type 1 diabetes and coeliac disease susceptibility loci for association with juvenile idiopathic arthritis
- Anne Hinks1,
- Paul Martin1,
- Edward Flynn1,
- Steve Eyre1,
- Jon Packham2,
- Anne Barton1,
- Childhood Arthritis Prospective Study (CAPS), UKRAG Consortium, BSPAR study group,
- Jane Worthington1,
- Wendy Thomson1
- 1Arthritis Research UK Epidemiology Unit, The University of Manchester, Manchester, UK
- 2Haywood Hospital, University Hospital of North Staffordshire, Stoke on Trent, Staffordshire, ST4 7LN
- Correspondence to Dr Anne Hinks, Arthritis Research UK Epidemiology Unit, Manchester Academic Health Science Centre, Stopford Building, The University of Manchester, Manchester M13 9PT, UK;
- Accepted 25 June 2010
- Published Online First 20 July 2010
Background There is strong evidence suggesting that juvenile idiopathic arthritis (JIA) shares many susceptibility loci with other autoimmune diseases.
Objective To investigate variants robustly associated with type 1 diabetes (T1D) or coeliac disease (CD) for association with JIA.
Methods Sixteen single-nucleotide polymorphisms (SNPs) already identified as susceptibility loci for T1D/CD were selected for genotyping in patients with JIA (n=1054) and healthy controls (n=3129). Genotype and allele frequencies were compared using the Cochrane–Armitage trend test implemented in PLINK.
Results One SNP in the LPP gene, rs1464510, showed significant association with JIA (ptrend=0.002, OR=1.18, 95% CI 1.06 to 1.30). A second SNP, rs653178 in ATXN2, also showed nominal evidence for association with JIA (ptrend=0.02, OR=1.13, 95% CI 1.02 to 1.25). The SNP, rs17810546, in IL12A showed subtype-specific association with enthesitis-related arthritis (ERA) subtype (ptrend=0.005, OR=1.88, 95% CI 1.2 to 2.94).
Conclusions Evidence for a novel JIA susceptibility locus, LPP, is presented. Association at the SH2B3/ATXN2 locus, previously reported to be associated with JIA in a US series, also supports this region as contributing to JIA susceptibility. In addition, a subtype-specific association of IL12A with ERA is identified. All findings will require validation in independent JIA cohorts.
In the past few years huge advances have been made in the understanding of the genes contributing to autoimmune disease susceptibility. It is also becoming clear that, despite the apparent clinical and phenotypic differences of the autoimmune diseases, they share a number of genetic risk factors.1 Of note, the genetic overlap between type 1 diabetes (T1D), coeliac disease (CD) and rheumatoid arthritis (RA)2 3 (Eyre et al, submitted).
Juvenile idiopathic arthritis (JIA) is the most common chronic rheumatic disease of childhood. It is defined as chronic inflammation of the synovial joints, with unknown aetiology, that starts before the age of 16 and persists for at least 6 weeks. Familial clustering of other autoimmune diseases in JIA families is well established, with an increased prevalence of autoimmunity in first- and second-degree relatives of patients.4
There is already strong evidence suggesting that JIA shares many susceptibility loci with other autoimmune diseases such as RA, T1D and multiple sclerosis.5 6 Therefore, the aim of this study was to test single-nucleotide polymorphisms (SNPs) robustly associated with T1D or CD in a large cohort of patients with JIA and controls to investigate the overlap between these diseases.
DNA was available for 1054 UK Caucasian patients with JIA (332 men, 715 women) as previously described.6 JIA cases were classified according to the International League of Associations for Rheumatology (ILAR) criteria7 (online supplementary table 1).
Three thousand one hundred and twenty-nine healthy UK Caucasian control DNA samples were available as described previously.8 All individuals were recruited with ethical approval and provided informed consent (North-West Multi-Centre Research Ethics Committee (MREC 99/8/84) and the University of Manchester Committee on the Ethics of Research on Human Beings).
SNPs that showed robust evidence for association (p<5×10−7) with T1D and/or CD2 and had not been investigated previously were identified. In total 16 SNPs were selected for genotyping.
All SNPs were genotyped, in UK JIA cases and controls, using the Sequenom iPlex MassARRAY platform according to manufacturers instructions (Sequenom, San Diego, California, USA, http://www.sequenom.com, accessed 9 July 2010). A 90% sample quality control rate and 90% SNP genotyping success rate was imposed on the analysis.
Power calculations were performed using QUANTO based on the effect sizes reported in previous studies of these SNPs in other autoimmune diseases. Calculations assumed a log-additive model and an α value corrected for the number of SNPs analysed. Genotype and allele frequencies were compared between cases with JIA and controls using PLINK.9 Comparisons were expressed as ORs and their 95% CIs as well as the Cochrane–Armitage test for trend. JIA is a phenotypically heterogeneous disease and can be classified into more clinically homogeneous diseases using the ILAR classification criteria.7 We examined whether there was evidence of a difference in allele frequencies between the seven ILAR subtypes. Differences between subtypes were assessed using χ2 tests on the 7×2 tables. Only if a significant difference was detected was stratification analysis by ILAR subtype performed and ORs and 95% CIs calculated for the subtypes.
Genotype counts for all SNPs were in Hardy–Weinberg equilibrium (p>0.05) in the control cohort. Three SNPs (rs12708716, rs45450798 and rs229541) failed SNP quality control thresholds and were excluded from further analysis, leaving 13 SNPs for analysis. A Bonferroni correction of 13 was applied to correct for the number of loci studied, resulting in a p value threshold of 0.004 for claims of significance (table 1).
One SNP in the LPP gene associated with CD, rs1464510, showed significant association with JIA at the corrected threshold (ptrend=0.002, OR=1.18, 95% CI 1.06 to 1.30) (table 1). A second SNP, rs653178, in the ATXN2 gene associated with CD, was nominally associated with JIA (ptrend=0.02, OR=1.13, 95% CI 1.02 to 1.25) but this was not statistically significant after correction for multiple testing (ptrend_corrected=0.26).
All SNPs were examined for differences in allele frequencies across the ILAR subtypes. The SNP, rs17810546, in IL12A associated with CD, showed a significant allele frequency difference (p=0.03) between the subtypes and therefore this SNP was re-analysed stratifying by ILAR subtype. The enthesitis-related arthritis (ERA) subtype showed strong association with IL12A SNP (ptrend=0.005, OR=1.88, 95% CI 1.2 to 2.94). There was no evidence for association with any other subtype.
The emerging overlap of autoimmune disease susceptibility loci makes many of the newly discovered T1D and CD loci appealing candidates for investigation as possible JIA susceptibility loci. There is already some evidence for the overlap of JIA and T1D susceptibility loci, with confirmed association of genes such as the major histocompatibility complex, PTPN22, IL2RA, IL2/21 region and CCR5 with both diseases.5 10,–,12
By investigating confirmed CD and T1D loci in JIA we have identified a novel locus conferring susceptibility to JIA. This SNP lies within the Lim domain containing preferred translocation partner in lipoma (LPP) gene. SNPs in LPP show validated association with CD.13 14 Interestingly, the rs1464510 SNP does not appear to be associated with T1D2 and studies in RA show conflicting results: in one large UK population study, no association was detected (ptrend=0.82, OR=1.01, 95% CI 0.94 to 1.08) (Eyre et al, submitted) while in a Dutch population, modest evidence for association was reported (p=0.012, OR=1.14, 95% CI 1.03 to 1.26).3
LPP has a number of known functions; it is important in cell migration, cell adhesion and acts as a versatile scaffolding and adaptor protein.15 Recently, LPP was identified as a substrate of the protein-tyrosine phosphatase 1B, a negative regulator of multiple signalling pathways downstream of receptor tyrosine kinases and functionally linked to Ras signalling.16 Further work will be required to determine which of the multiple potential pathways with which the protein interacts is affected by the variants that predispose to autoimmune disease.
We also found weak evidence for association of a SNP, rs653178, in the SH2B3/ATXN2 gene region with JIA, though it did not remain significant after correction for multiple testing. The SNP lies within intron 1 of the ATXN2 gene which lies adjacent to the SH2B3 gene. This SNP is in complete linkage disequilibrium (r2=1) with the SNP, rs3184504, in SH2B3, and both SNPs show validated association with CD.13 14 rs3184504 has been also been associated with both RA and CD in a Dutch population.3 The SH2B3 rs3184504 variant is a non-synonymous SNP, R262W, located in exon 3 of the gene. SH2B3 encodes the T-cell adapter protein Lnk, which regulates T-cell receptor-, growth factor- and cytokine receptor-mediated signalling. It has been shown to be a negative regulator of tumour necrosis factor α signalling in human endothelial cells.17 Another SNP, rs17696736, within the adjacent gene c12orf30, has also been reported to be associated with JIA in a US case–control cohort.18 There is only modest linkage disequilibrium between that SNP and the SNP tested in our study (r2=0.5). Hence, both may simply be markers of an, as yet, un-genotyped causal variant or there may be independent effects in the region.
We also found that the SNP, rs17810546, which lies 5′ to the IL12A gene, showed differences in allele frequencies across the ILAR subtypes and this was driven by a strong association with ERA. This is a subgroup comprising mainly male, HLA-B27-positive patients and it is characterised by the presence of enthesitis and arthritis. It may progressively involve the sacroiliac joints and spine, with symptoms similar to ankylosing spondylitis. It is important to note that the numbers in this subgroup are small (n=61) so this finding should be interpreted with caution and will require validation in an independent cohort. To date the IL12A gene has not been associated with ankylosing spondylitis.
Although none of the other SNPs studied showed significant association with JIA, we cannot unequivocally rule out association for many of these loci. This is owing to the limited power in this study to detect some of the loci that were identified in the studies for CD and T1D, which had much larger sample sizes. We had >80% power to detect an effect at RGS1, IL18RAP, IL12A, INS and ERBB3 and therefore feel more confident that the SNPs tested at those loci are not associated with JIA. OR plots for CD and T1D susceptibility genes across the different autoimmune diseases, JIA, RA, T1D and CD, are shown in figures 1 and 2, respectively.
In conclusion we have identified novel association of the LPP gene and show weak but confirmatory association of the ATXN2/SH2B3/c12orf30 region with JIA in a UK cohort. In addition, we report an ERA subtype-specific association with the IL12A gene. Investigation of all these loci in independent JIA case–control cohorts followed by meta-analysis will be required to confidently confirm or refute these as JIA susceptibility loci. Nonetheless, the approach of targeting variants associated with other autoimmune diseases is already yielding insights into the genetic complexity underlying susceptibility to this serious childhood disease.
Childhood arthritis prospective study (CAPS) group Eileen Baildam, Lynsey Brown, Joanne Buckley, Alice Chieng, Joyce Davidson, Michael Eltringham, Helen Foster, Mark Friswell, Janet Gardner-Medwin, Paul Gilbert, Kimme Hyrich, Julie Jones, Sham Lal, Mark Lay, Carol Lydon, Alexandra Meijer, Vicki Price, Jane Sim, Maureen Todd, Peter Ward, Lucy Wedderburn.
UKRAG Consortium Stephen Eyre, Anne Hinks, John Bowes, Gisela Orozco, Edward Flynn, Paul Martin, Wendy Thomson, Anne Barton, Jane Worthington: Arthritis Research UK Epidemiology Unit, Manchester Academic Health Science Centre, Stopford Building, The University of Manchester, Manchester, UK; Stephen Martin, James I Robinson, Ann W Morgan, Paul Emery: Leeds Institute of Molecular Medicine, Section of Musculoskeletal Disease, University of Leeds, Leeds LS9 7TF, UK; Anthony G Wilson: School of Medicine and Biomedical Sciences, The University of Sheffield, Sheffield S10 2JF, UK; Sophia Steer: Clinical and Academic Rheumatology, Kings College Hospital NHS Foundation Trust, Denmark Hill, London SE5 9RS, UK; Lynne Hocking, David M Reid: Bone Research Group, Department of Medicine and Therapeutics, University of Aberdeen AB25 2ZD, UK; Pille Harrison, Paul Wordsworth: University of Oxford Institute of Musculoskeletal Sciences, Botnar Research Centre, Oxford OX3 7LD, UK.
British Society of Paediatric and Adolescent Rheumatology (BSPAR) study group M Abinum, MD, M Becker, MD, A Bell, MD, A Craft, MD, E Crawley, MD, J David, MD, H Foster, MD, J Gardener-Medwin, MD, J Griffin, MD, A Hall, MD, M Hall, MD, A Herrick, MD, P Hollingworth, MD, L Holt, MD, S Jones, MD, G Pountain, MD, C Ryder, MD, T Southwood, MD, I Stewart, MD, H Venning, MD, L Wedderburn, MD, P Woo, MD, and S Wyatt, MD.
Funding This work was supported by Arthritis Research UK: AR UK grant reference No 17552. The authors would like to thank David Strachan for facilitating access to the 1958 Birth Cohort. The authors acknowledge use of genotype data from the British 1958 Birth Cohort DNA collection, funded by the Medical Research Council grant G0000934 and the Wellcome Trust grant 068545/Z/02. The authors acknowledge support from the NIHR Manchester Biomedical Research Council.
Ethics approval This study was conducted with the approval of the (North-West Multi-Centre Research Ethics Committee (MREC 99/8/84) and the University of Manchester Committee on the Ethics of Research on Human Beings).
Provenance and peer review Not commissioned; externally peer reviewed.
This paper is freely available online under the BMJ Journals unlocked scheme, see http://ard.bmj.com/info/unlocked.dtl