Short Report | Open | Published:
Independent replication analysis of genetic loci with previous evidence of association with juvenile idiopathic arthritis
Pediatric Rheumatologyvolume 11, Article number: 12 (2013)
Over the last five years, there have been numerous reports of association of juvenile idiopathic arthritis with single nucleotide polymorphisms (SNPs) at various loci outside the major histocompatibility complex (MHC) region. However, the majority of these association findings have been generated using a limited number of international cohorts, and thus there is benefit in further independent replication. To address this, we examined a total of 56 SNPs in 42 non-MHC gene regions previously reported to be associated with JIA, in the ChiLdhood Arthritis Risk factor Identification sTudY (CLARITY), a new Australian collection of cases and healthy child controls.
Genotyping was performed on a total of 324 JIA cases (mean age 9.7 years, 67.3% female) and 568 controls (mean age 7.8 years, 40.7% female). We demonstrated clear evidence for replication of association of JIA with SNPs in or around c12orf30, c3orf1, PTPN22, STAT4, and TRAF1-C5, confirming the involvement of these loci in disease risk. Further, we generated evidence supportive of replication of association of JIA with loci containing AFF3, CD226, MBL2, PSTPIP1, and RANTES (CCL5). These results were robust to sensitivity analyses for ethnicity.
We have provided valuable independent data as to the underlying genetic architecture of this understudied pediatric autoimmune disease, further confirming five loci outside the MHC, and supporting a role for a further five loci in determining disease risk.
Our understanding of the genetic basis of juvenile idiopathic arthritis (JIA) has recently increased, but still lags behind many other autoimmune diseases. This is largely due to the paucity of DNA collections internationally. While association of JIA with variation in the major histocompatibility complex (MHC) is well-established, over the last five years, there have been a number of reports of new JIA susceptibility loci that lie outside this region. These findings have resulted from candidate gene approaches (for example, examining genes known to be associated with rheumatoid arthritis), and more recently, from a limited number of genome-wide association study (GWAS) approaches. The vast majority of reports describe discovery and replication findings generated from two large sample collections, from the UK and the US. Thus, although solid evidence for association with JIA has usually been described, for many of the identified loci, further replication in an entirely independent JIA sample would be beneficial in confirming their contribution to disease risk.
To address this, we performed a single nucleotide polymorphism (SNP) replication study within a new sample of JIA cases and healthy hospital-based child controls, collected so far as part of the ongoing Australian CLARITY study (ChiLdhood Arthritis Risk factor Identification sTudY). We selected 56 SNPs from 42 gene regions. Candidacy was based on published evidence of association, or a trend towards association, with total JIA, or one or more of its subtypes. These genes included ADAD1-IL2-IL21[4–6], AFF3, ANGPT1, ATXN2, BACH2, c12orf30[6, 8], c3orf1, CCR5, CD14, CD226, CD247, CLEC16A, COG6, CTLA4, ERAP1, IL12A, IL15, IL23R[6, 13], IL2RA[2, 6, 14], IL7R, JMJD1C-REEP3, KIF5A, LPP, MBL2, MEFV, NLRP3, NOD2, NRBF2-EGR2, PRKCQ, PSTPIP1, PTPN2[2, 6], PTPN22[6, 18], RANTES (CCL5), STAT4[6, 8, 15], TNFA, TNFAIP3[8, 15], TRAF1-C5[15, 21], and VTCN1. Additional SNPs from four other loci not attributed to any gene in the original publication but lying closest to the genes DCN1, FHIT, HUNK, and SLITRK5 were also selected. Genotyping was performed on a total of 324 JIA cases (mean age 9.7 years, 67.3% female) and 568 controls (mean age 7.8 years, 40.7% female) using the Sequenom MassARRAY system (assay design details available from the authors). After quality control pruning (removal of SNPs with Hardy Weinberg Equilibrium p < 0.01, and SNPs or samples with < 90% genotyping success rate) genotypes were successfully generated for 51 SNPs at 38 gene loci in 318 cases and 556 controls. SNPs excluded from analysis are listed in Additional file1. Allelic, genotypic, additive (Cochrane-Armitage test for trend), dominant and recessive association analyses were performed using PLINK v1.07. Given the a priori evidence for association of each SNP, no adjustments for multiple testing were made.
Outcomes of the association analyses in CLARITY cases and controls, including p values for the most highly significant test and allelic test for all SNPs analysed is shown in Table 1. Figure 1 shows a forest plot of the allelic odds ratios generated in the CLARITY sample compared to previously published odds ratios (ORs) for the majority of SNPs tested. Some SNPs were not included in this figure, including SNPs for which an OR was not presented in the original report, or where SNPs were previously associated specifically with rarer JIA subtypes. Our data demonstrates clear evidence of replication (defined as p < 0.05 for any test, ORs or case-control allele frequency differences in a direction consistent with previous reports) for SNPs at loci containing the genes ATXN2, c12orf30, c3orf1, PTPN22, STAT4, and TRAF1-C5. Possible evidence of replication (defined as p < 0.2 for any test, ORs or case-control allele frequency differences in a direction consistent with previous reports) was also generated for SNPs near AFF3, CD226, MBL2, PSTPIP1, and RANTES (CCL5). For IL15, we found a significant association, but in the opposite direction to that previously reported. The IL15 SNP is an A/T transversion with frequencies of both alleles close to 50%, and thus there is the possibility of allele reversal. Although our minor T allele frequency in controls was entirely consistent with that found in the prior study, it is difficult to be sure that the IL15 assays used in both studies were based on the same DNA strand. We therefore cannot be sure if our data supports replication. Additionally for CTLA4, considered a general autoimmunity susceptibility gene but with conflicting association results for JIA[5, 24], we did not generate any evidence of replication.
We then performed a sensitivity analysis in which we included only cases (n = 200) and controls (n = 341) of European ancestry using the stringent definition of self-reported European ancestry of all four of the child’s grandparents. Non-European participants, along with participants for whom full grandparent data was not provided, were excluded. The outcome of this re-analysis is shown in Additional file1: Table S1 and Additional file1: Figure S1. In general, the results of this European-only analysis were not materially different to the full sample analysis (taking into account reduced statistical power resulting from a reduction of sample size), with one exception. ATXN2 appeared influenced by ethnicity, with an opposite, non-significant, direction of effect in the European subgroup.
In conclusion, we have provided independent replication data for JIA susceptibility loci that have previously been identified using a generally limited number of international sample collections. We have confirmed association of JIA with SNPs close-by to c12orf30, c3orf1, PTPN22, STAT4, and TRAF1-C5; and we have provided further support for the association of SNPs close-by to AFF3, CD226, MBL2, PSTPIP1, and CCL5. A limitation of our study was our relatively small sample size; our full dataset had 80% power to detect an OR of 1.4 for an allele at 20% frequency in the population at an alpha of 0.05. Given that many of the published ORs were less than 1.4, and that a number of the SNPs analysed had minor allele frequencies less than 20%, we cannot, from our current data, exclude association of any of the SNPs examined. Our current sample size also precluded detailed subtype-specific association analyses. Large collaborative GWAS efforts would be beneficial in confirming the outstanding genes, and providing further novel insights into the breadth of genetic loci involved in JIA susceptibility.
Hollenbach JA, Thompson SD, Bugawan TL, Ryan M, Sudman M, Marion M, Langefeld CD, Thomson G, Erlich HA, Glass DN: Juvenile idiopathic arthritis and HLA class I and class II interactions and age-at-onset effects. Arthritis Rheum. 2010, 62: 1781-1791. 10.1002/art.27424.
Hinks A, Cobb J, Sudman M, Eyre S, Martin P, Flynn E, Packham J, Barton A, Worthington J, Langefeld CD: Investigation of rheumatoid arthritis susceptibility loci in juvenile idiopathic arthritis confirms high degree of overlap. Ann Rheum Dis. 2012, 71: 1117-1121. 10.1136/annrheumdis-2011-200814.
Ellis JA, Ponsonby AL, Pezic A, Chavez RA, Allen RC, Akikusa JD, Munro JE: CLARITY: ChiLdhood Arthritis Risk factor Identification sTudY. Pediatr Rheumatol. 2012, 10: 37-10.1186/1546-0096-10-37.
Albers HM, Kurreeman FA, Stoeken-Rijsbergen G, Brinkman DM, Kamphuis SS, van Rossum MA, Girschick HJ, Wouters C, Saurenmann RK, Hoppenreijs E: Association of the autoimmunity locus 4q27 with juvenile idiopathic arthritis. Arthritis Rheum. 2009, 60: 901-904. 10.1002/art.24296.
Hinks A, Eyre S, Ke X, Barton A, Martin P, Flynn E, Packham J, Worthington J, Thomson W: Association of the AFF3 gene and IL2/IL21 gene region with juvenile idiopathic arthritis. Genes Immun. 2010, 11: 194-198. 10.1038/gene.2009.105.
Thompson SD, Sudman M, Ramos PS, Marion MC, Ryan M, Tsoras M, Weiler T, Wagner M, Keddache M, Haas JP: The susceptibility loci juvenile idiopathic arthritis shares with other autoimmune diseases extend to PTPN2, COG6, and ANGPT1. Arthritis Rheum. 2010, 62: 3265-3276. 10.1002/art.27688.
Hinks A, Martin P, Flynn E, Eyre S, Packham J, Barton A, Worthington J, Thomson W: Investigation of type 1 diabetes and coeliac disease susceptibility loci for association with juvenile idiopathic arthritis. Ann Rheum Dis. 2010, 69: 2169-2172. 10.1136/ard.2009.126938.
Prahalad S, Hansen S, Whiting A, Guthery SL, Clifford B, McNally B, Zeft AS, Bohnsack JF, Jorde LB: Variants in TNFAIP3, STAT4, and C12orf30 loci associated with multiple autoimmune diseases are also associated with juvenile idiopathic arthritis. Arthritis Rheum. 2009, 60: 2124-2130. 10.1002/art.24618.
Thompson SD, Marion MC, Sudman M, Ryan M, Tsoras M, Howard TD, Barnes MG, Ramos PS, Thomson W, Hinks A: Genome-wide association analysis of juvenile idiopathic arthritis identifies a new susceptibility locus at chromosomal region 3q13. Arthritis Rheum. 2012, 64: 2781-2791. 10.1002/art.34429.
Hinks A, Martin P, Flynn E, Eyre S, Packham J, Barton A, Worthington J, Thomson W: Association of the CCR5 gene with juvenile idiopathic arthritis. Genes Immun. 2010, 11: 584-589. 10.1038/gene.2010.25.
Zeng HS, Chen XY, Luo XP: The association with the -159C/T polymorphism in the promoter region of the CD14 gene and juvenile idiopathic arthritis in a Chinese Han population. J Rheumatol. 2009, 36: 2025-2028. 10.3899/jrheum.081093.
Skinningsrud B, Lie BA, Husebye ES, Kvien TK, Forre O, Flato B, Stormyr A, Joner G, Njolstad PR, Egeland T, Undlien DE: A CLEC16A variant confers risk for juvenile idiopathic arthritis and anti-cyclic citrullinated peptide antibody negative rheumatoid arthritis. Ann Rheum Dis. 2010, 69: 1471-1474. 10.1136/ard.2009.114934.
Hinks A, Martin P, Flynn E, Eyre S, Packham J, Barton A, Worthington J, Thomson W: Subtype specific genetic associations for juvenile idiopathic arthritis: ERAP1 with the enthesitis related arthritis subtype and IL23R with juvenile psoriatic arthritis. Arthritis Res Ther. 2011, 13: R12-10.1186/ar3235.
Hinks A, Ke X, Barton A, Eyre S, Bowes J, Worthington J, Thompson SD, Langefeld CD, Glass DN, Thomson W: Association of the IL2RA/CD25 gene with juvenile idiopathic arthritis. Arthritis Rheum. 2009, 60: 251-257. 10.1002/art.24187.
Hinks A, Eyre S, Ke X, Barton A, Martin P, Flynn E, Packham J, Worthington J, Thomson W: Overlap of disease susceptibility loci for rheumatoid arthritis and juvenile idiopathic arthritis. Ann Rheum Dis. 2010, 69: 1049-1053. 10.1136/ard.2009.110650.
Gergely P, Pazar B, Nagy ZB, Gombos T, Rajczy K, Balogh Z, Orban I, Sevcic K, Poor G: Structural polymorphisms in the mannose-binding lectin gene are associated with juvenile idiopathic arthritis. J Rheumatol. 2009, 36: 843-847. 10.3899/jrheum.080681.
Day TG, Ramanan AV, Hinks A, Lamb R, Packham J, Wise C, Punaro M, Donn RP: Autoinflammatory genes and susceptibility to psoriatic juvenile idiopathic arthritis. Arthritis Rheum. 2008, 58: 2142-2146. 10.1002/art.23604.
Hinks A, Barton A, John S, Bruce I, Hawkins C, Griffiths CE, Donn R, Thomson W, Silman A, Worthington J: Association between the PTPN22 gene and rheumatoid arthritis and juvenile idiopathic arthritis in a UK population: further support that PTPN22 is an autoimmunity gene. Arthritis Rheum. 2005, 52: 1694-1699. 10.1002/art.21049.
Yao TC, Tsai YC, Huang JL: Association of RANTES promoter polymorphism with juvenile rheumatoid arthritis. Arthritis Rheum. 2009, 60: 1173-1178. 10.1002/art.24422.
Jimenez-Morales S, Velazquez-Cruz R, Ramirez-Bello J, Bonilla-Gonzalez E, Romero-Hidalgo S, Escamilla-Guerrero G, Cuevas F, Espinosa-Rosales F, Martinez-Aguilar NE, Gomez-Vera J: Tumor necrosis factor-alpha is a common genetic risk factor for asthma, juvenile rheumatoid arthritis, and systemic lupus erythematosus in a Mexican pediatric population. Hum Immunol. 2009, 70: 251-256. 10.1016/j.humimm.2009.01.027.
Behrens EM, Finkel TH, Bradfield JP, Kim CE, Linton L, Casalunovo T, Frackelton EC, Santa E, Otieno FG, Glessner JT: Association of the TRAF1-C5 locus on chromosome 9 with juvenile idiopathic arthritis. Arthritis Rheum. 2008, 58: 2206-2207. 10.1002/art.23603.
Hinks A, Barton A, Shephard N, Eyre S, Bowes J, Cargill M, Wang E, Ke X, Kennedy GC, John S: Identification of a novel susceptibility locus for juvenile idiopathic arthritis by genome-wide association analysis. Arthritis Rheum. 2009, 60: 258-263. 10.1002/art.24179.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81: 559-575. 10.1086/519795.
Prahalad S, Bohnsack JF, Whiting A, Clifford B, Jorde LB, Guthery SL, Thompson SD, Glass DN, Bamshad MJ: Lack of association of functional CTLA4 polymorphisms with juvenile idiopathic arthritis. Arthritis Rheum. 2008, 58: 2147-2152. 10.1002/art.23602.
We gratefully thank the families who kindly participated in this study. We thank Betty Lim and Sarah Macnee for their involvement in the initiation of the study and recruitment of participants, Amanda Hawker and William Siero for assistance with control recruitment and data entry, and Susan Matthyz-Rosa for assistance with data entry. We acknowledge support from the National Health and Medical Research Council of Australia, Arthritis Australia, Rebecca L Cooper Foundation, LEW Carty Charitable Fund, ANZ Medical Research and Technologies in Victoria Fund, Equity Trustees Lynne Quayle Charitable Trust Fund, The Australian Academy of Science, and The Victorian State Government Operational Infrastructure Support Program. JE was supported by a National Health and Medical Research Council of Australia Capacity Building Grant in Population Health, and an Australian Research Council Future Fellowship.
The authors declare that they have no competing interests.
JE and JM conceived, designed and led the study. RC performed the SNP genotyping. AP managed the study data. ALP advised and assisted with study design and collection of epidemiological data. RA and JA assisted with case study design, and JM, RA and JA assisted with recruitment through their paediatric rheumatology clinics. JE wrote the manuscript, and all authors participated in drafting the manuscript to the final version. All authors read and approved the final manuscript.