Variants in CXCR4 associate with juvenile idiopathic arthritis susceptibility.

BACKGROUND: Juvenile idiopathic arthritis (JIA) is the most common chronic rheumatic disease among children, the etiology of which involves a strong genetic component, but much of the underlying genetic determinants still remain unknown. Our aim was to identify novel genetic variants that predispose to JIA. METHODS: We performed a genome-wide association study (GWAS) and replication in a total of 1166 JIA cases and 9500 unrelated controls of European ancestry. Correlation of SNP genotype and gene expression was investigated. Then we conducted targeted resequencing of a candidate locus, among a subset of 480 cases and 480 controls. SUM test was performed to evaluate the association of the identified rare functional variants. RESULTS: The CXCR4 locus on 2q22.1 was found to be significantly associated with JIA, peaking at SNP rs953387. However, this result is subjected to subpopulation stratification within the subjects of European ancestry. After adjusting for principal components, nominal significant association remained (p < 10(-4)). Because of its interesting known function in immune regulation, we carried out further analyses to assess its relationship with JIA. Expression of CXCR4 was correlated with CXCR4 rs953387 genotypes in lymphoblastoid cell lines (p = 0.014) and T-cells (p = 0.0054). In addition, rare non-synonymous and stop-gain sequence variants in CXCR4, putatively damaging for CXCR4 function, were significantly enriched in JIA cases (p = 0.015). CONCLUSION: Our results suggest the association of CXCR4 variants with JIA, implicating that this gene may be involved in the pathogenesis of autoimmune disease. However, because this locus is subjected to population stratification within the subjects of European ancestry, additional replication is still necessary for this locus to be considered a true risk locus for JIA. This cell-surface chemokine receptor has already been targeted in other diseases and may serve as a tractable therapeutic target for a specific subset of pediatric arthritis patients with additional replication and functional validation of the locus.

Dense genotyping of immune-related disease regions identifies 14 new susceptibility loci for juvenile idiopathic arthritis.

We used the Immunochip array to analyze 2,816 individuals with juvenile idiopathic arthritis (JIA), comprising the most common subtypes (oligoarticular and rheumatoid factor-negative polyarticular JIA), and 13,056 controls. We confirmed association of 3 known JIA risk loci (the human leukocyte antigen (HLA) region, PTPN22 and PTPN2) and identified 14 loci reaching genome-wide significance (P < 5 × 10(-8)) for the first time. Eleven additional new regions showed suggestive evidence of association with JIA (P < 1 × 10(-6)). Dense mapping of loci along with bioinformatics analysis refined the associations to one gene in each of eight regions, highlighting crucial pathways, including the interleukin (IL)-2 pathway, in JIA disease pathogenesis. The entire Immunochip content, the HLA region and the top 27 loci (P < 1 × 10(-6)) explain an estimated 18, 13 and 6% of the risk of JIA, respectively. In summary, this is the largest collection of JIA cases investigated so far and provides new insight into the genetic basis of this childhood autoimmune disease.

Susceptibility to childhood-onset rheumatoid arthritis: investigation of a weighted genetic risk score that integrates cumulative effects of variants at five genetic loci.

OBJECTIVE: Children with childhood-onset rheumatoid arthritis (RA) include those with rheumatoid factor or anti-citrullinated protein antibody-positive juvenile idiopathic arthritis. To test the hypothesis that adult-onset RA-associated variants are also associated with childhood-onset RA, we investigated RA-associated variants at 5 loci in a cohort of patients with childhood-onset RA. We also assessed the cumulative association of these variants in susceptibility to childhood-onset RA using a weighted genetic risk score (wGRS). METHODS: A total of 155 children with childhood-onset RA and 684 healthy controls were genotyped for 5 variants in the PTPN22, TRAF1/C5, STAT4, and TNFAIP3 loci. High-resolution HLA-DRB1 genotypes were available for 149 cases and 373 controls. We tested each locus for association with childhood-onset RA via logistic regression. We also computed a wGRS for each subject, with weights based on the natural log of the published odds ratios (ORs) for the alleles investigated, and used logistic regression to test the wGRS for association with childhood-onset RA. RESULTS: Childhood-onset RA was associated with TNFAIP3 rs10499194 (OR 0.60 [95% confidence interval 0.44-0.83]), PTPN22 rs2476601 (OR 1.61 [95% confidence interval 1.11-2.31]), and STAT4 rs7574865 (OR 1.41 [95% confidence interval 1.06-1.87]) variants. The wGRS was significantly different between cases and controls (P < 2 × 10(-16) ). Individuals in the third to fifth quintiles of wGRS had a significantly increased disease risk compared to baseline (individuals in the first quintile). Higher wGRS was associated with increased risk of childhood-onset RA, especially among males. CONCLUSION: The magnitude and direction of the association between TNFAIP3, STAT4, and PTPN22 variants and childhood-onset RA are similar to those observed in RA, suggesting that adult-onset RA and childhood-onset RA share common genetic risk factors. Using a wGRS, we have demonstrated the cumulative association of RA-associated variants with susceptibility to childhood-onset RA.

Multiple juvenile idiopathic arthritis subtypes demonstrate proinflammatory IgG glycosylation.

OBJECTIVE: Rheumatoid arthritis is associated with an excess of agalactosylated (G0) IgG that is considered relatively proinflammatory. Assessment of this association in juvenile idiopathic arthritis (JIA) is complicated by age-dependent IgG glycan variation. The aim of this study was to conduct the first large-scale survey of IgG glycans in healthy children and patients with JIA, with a focus on early childhood, the time of peak JIA incidence. METHODS: IgG glycans from healthy children and disease-modifying antirheumatic drug-naive patients with JIA were characterized using high-performance liquid chromatography. Agalactosylated glycans were quantitated with reference to monogalactosylated (G1) species. Associations were sought between the G0:G1 ratio and disease characteristics. RESULTS: Among healthy children ages 9 months to 16 years (n = 165), the G0:G1 ratio was highly age dependent, with the ratio peaking to 1.19 in children younger than age 3 years and declining to a nadir of 0.83 after age 10 years (Spearman's ρ = 0.60, P < 0.0001). In patients with JIA (n = 141), the G0:G1 ratio was elevated compared with that in control subjects (1.32 versus 1.02; P < 0.0001). The G0:G1 ratio corrected for age was abnormally high in all JIA subtypes (enthesitis-related arthritis was not assessed), most strikingly in systemic JIA. Glycosylation aberrancy was comparable in patients with and those without antinuclear antibodies and in both early- and late-onset disease and exhibited at most a weak correlation with markers of inflammation. CONCLUSION: IgG glycosylation is skewed toward proinflammatory G0 variants in healthy children, in particular during the first few years of life. This deviation is exaggerated in patients with JIA. The role for IgG glycan variation in immune function in children, including the predilection of JIA for early childhood, remains to be defined.

Investigation of rheumatoid arthritis susceptibility loci in juvenile idiopathic arthritis confirms high degree of overlap.

OBJECTIVES: Rheumatoid arthritis (RA) shares some similar clinical and pathological features with juvenile idiopathic arthritis (JIA); indeed, the strategy of investigating whether RA susceptibility loci also confer susceptibility to JIA has already proved highly successful in identifying novel JIA loci. A plethora of newly validated RA loci has been reported in the past year. Therefore, the aim of this study was to investigate these single nucleotide polymorphisms (SNP) to determine if they were also associated with JIA. METHODS: Thirty-four SNP that showed validated association with RA and had not been investigated previously in the UK JIA cohort were genotyped in JIA cases (n=1242), healthy controls (n=4281), and data were extracted for approximately 5380 UK Caucasian controls from the Wellcome Trust Case-Control Consortium 2. Genotype and allele frequencies were compared between cases with JIA and controls using PLINK. A replication cohort of 813 JIA cases and 3058 controls from the USA was available for validation of any significant findings. RESULTS: Thirteen SNP showed significant association (p<0.05) with JIA and for all but one the direction of association was the same as in RA. Of the eight loci that were tested, three showed significant association in the US cohort. CONCLUSIONS: A novel JIA susceptibility locus was identified, CD247, which represents another JIA susceptibility gene whose protein product is important in T-cell activation and signalling. The authors have also confirmed association of the PTPN2 and IL2RA genes with JIA, both reaching genome-wide significance in the combined analysis.

Genome-wide association analysis of juvenile idiopathic arthritis identifies a new susceptibility locus at chromosomal region 3q13.

OBJECTIVE: In a genome-wide association study of Caucasian patients with juvenile idiopathic arthritis (JIA), we have previously described findings limited to autoimmunity loci shared by JIA and other diseases. The present study was undertaken to identify novel JIA-predisposing loci using genome-wide approaches. METHODS: The discovery cohort consisted of Caucasian JIA cases (n = 814) and local controls (n = 658) genotyped on the Affymetrix Genome-Wide SNP 6.0 Array, along with 2,400 out-of-study controls. In a replication study, we genotyped 10 single-nucleotide polymorphisms (SNPs) in 1,744 cases and 7,010 controls from the US and Europe. RESULTS: Analysis within the discovery cohort provided evidence of associations at 3q13 within C3orf1 and near CD80 (rs4688011) (odds ratio [OR] 1.37, P = 1.88 × 10(-6) ) and at 10q21 near JMJD1C (rs647989 [OR 1.59, P = 6.1 × 10(-8) ], rs12411988 [OR 1.57, P = 1.16 × 10(-7) ], and rs10995450 [OR 1.31, P = 6.74 × 10(-5) ]). Meta-analysis provided further evidence of association for these 4 SNPs (P = 3.6 × 10(-7) for rs4688011, P = 4.33 × 10(-5) for rs6479891, P = 2.71 × 10(-5) for rs12411988, and P = 5.39 × 10(-5) for rs10995450). Gene expression data on 68 JIA cases and 23 local controls showed cis expression quantitative trait locus associations for C3orf1 SNP rs4688011 (P = 0.024 or P = 0.034, depending on the probe set) and JMJD1C SNPs rs6479891 and rs12411988 (P = 0.01 or P = 0.04, depending on the probe set and P = 0.008, respectively). Using a variance component liability model, it was estimated that common SNP variation accounts for approximately one-third of JIA susceptibility. CONCLUSION: Genetic association results and correlated gene expression findings provide evidence of JIA association at 3q13 and suggest novel genes as plausible candidates in disease pathology.

Hierarchy of risk of childhood-onset rheumatoid arthritis conferred by HLA-DRB1 alleles encoding the shared epitope.

OBJECTIVE: Associations between shared epitope (SE)-encoding HLA-DRB1 alleles and rheumatoid arthritis (RA) are well established. However, only a limited number of studies have investigated these alleles in patients with childhood-onset RA, which is defined as rheumatoid factor- and/or anti-citrullinated protein antibody-positive juvenile idiopathic arthritis. The aims of this study were to investigate the largest cohort of patients with childhood-onset RA for association with SE alleles and to determine whether there is a hierarchy of risk based on the amino acid sequence of the SE. METHODS: High-resolution HLA-DRB1 genotypes were obtained for 204 patients with childhood-onset RA and 373 healthy control subjects. Odds ratios (ORs) and 95% confidence intervals (95% CIs) were calculated for different SE-encoding HLA-DRB1 alleles. In addition, genotype ORs were calculated for combinations of SE alleles classified into S(2) , S(3P) , or L alleles, based on amino acid sequences in position 70-74 of the DRß1 chain, as proposed by Tezenas du Montcel et al. RESULTS: We confirmed associations between HLA-DRB1 SE alleles and childhood-onset RA (76% of patients carried 1 or 2 SE alleles compared with 46% of control subjects; OR 3.81, 95% CI 2.4-6.0, P < 1 × 10(-7) ). We also observed associations between individual SE alleles (HLA-DRB1*0101, *0401, *0404, *0405, *0408, and *1001) and childhood-onset RA. Genotype-specific risk estimates suggested a hierarchy of risk, with the highest risk among individuals heterozygous for S(2) /S(3P) (OR 22.3, 95% CI 9.9-50.5, P < 0.0001). CONCLUSION: We confirm the association between SE-encoding HLA-DRB1 alleles and susceptibility to childhood-onset RA. The excess risk conferred by carriage of the combination of S(2) and S(3P) risk alleles suggests that children with DRß1 chains containing the KRAA and QRRAA or RRRAA sequences are especially susceptible to RA.

Association of the 5-aminoimidazole-4-carboxamide ribonucleotide transformylase gene with response to methotrexate in juvenile idiopathic arthritis.

OBJECTIVES: Methotrexate (MTX) is the mainstay treatment for juvenile idiopathic arthritis (JIA), however approximately 30% of children will fail to respond to the drug. Identification of genetic predictors of response to MTX would be invaluable in developing optimal treatment strategies for JIA. Using a candidate gene approach, single nucleotide polymorphisms (SNPs) within genes in the metabolic pathway of MTX, were investigated for association with response to treatment in JIA cases. METHODS: Tagging SNPs were selected across 13 MTX metabolic pathway genes and were genotyped using Sequenom genotyping technology in subjects recruited from the Sparks Childhood Arthritis Response to Medication Study. Response to MTX was defined using the American College of Rheumatology (ACR) paediatric response criteria and SNP genotype frequencies were compared between the worst and best responders (ACR-Ped70) to MTX. An independent cohort of US JIA cases was available for validation of initial findings. RESULTS: One SNP within the inosine triphosphate pyrophosphatase gene (ITPA) and two SNPs within 5-aminoimidazole-4-carboxamide ribonucleotide transformylase gene (ATIC) were significantly associated with a poor response to MTX. One of the ATIC SNPs showed a trend towards association with MTX response in an independent cohort of US JIA cases. Meta-analysis of the two studies strengthened this association (combined p value=0.002). CONCLUSIONS: This study presents association of a SNP in the ATIC gene with response to MTX in JIA. There is now growing evidence to support a role of the ATIC gene with response to MTX treatment. These results could contribute towards a better understanding of and ability to predict MTX response in JIA.

Generation of novel pharmacogenomic candidates in response to methotrexate in juvenile idiopathic arthritis: correlation between gene expression and genotype.

OBJECTIVES: Little is known about the mechanisms of efficacy of methotrexate (MTX) in childhood arthritis, or genetic influences upon response to MTX. The aims of this study were to use gene expression profiling to identify novel pathways/genes altered by MTX and then investigate these genes for genotype associations with response to MTX treatment. METHODS: Gene expression profiling before and after MTX treatment was performed on 11 children with juvenile idiopathic arthritis (JIA) treated with MTX, in whom response at 6 months of treatment was defined. Genes showing the most differential gene expression after the treatment were selected for single nucleotide polymorphism (SNP) genotyping. Genotype frequencies were compared between nonresponders and responders (ACR-Ped70). An independent cohort was available for validation. RESULTS: Gene expression profiling before and after MTX treatment revealed 1222 differentially expressed probes sets (fold change >1.7, P<0.05) and 1065 when restricted to full responder cases only. Six highly differentially expressed genes were analyzed for genetic association in response to MTX. Three SNPs in the SLC16A7 gene showed significant association with MTX response. One SNP showed validated association in an independent cohort. CONCLUSION: This study is the first, to our knowledge, to evaluate gene expression profiles in children with JIA before and after MTX, and to analyze genetic variation in differentially expressed genes. We have identified a gene, which may contribute to genetic variability in MTX response in JIA, and established as proof of principle that genes that are differentially expressed at mRNA level after drug administration may also be good candidates for genetic analysis.

The susceptibility loci juvenile idiopathic arthritis shares with other autoimmune diseases extend to PTPN2, COG6, and ANGPT1.

OBJECTIVE: To test for associations between non-major histocompatibility complex susceptibility loci previously reported in autoimmune diseases and juvenile idiopathic arthritis (JIA). METHODS: Published autoimmune disease genome-wide association studies were reviewed, and 519 single-nucleotide polymorphisms (SNPs) were selected for association testing. The initial cohort included 809 JIA cases and 3,535 controls of non-Hispanic, European ancestry. Of the SNPs, 257 were successfully genotyped, while 168 were imputed with quality. Based on findings in the initial cohort, replication was sought for 21 SNPs in a second cohort of 1,015 JIA cases and 1,569 controls collected in the US and Germany. For the initial cohort, tests for association were adjusted for potential confounding effects of population structure by including principal components derived from a genome-wide association study as covariates in logistic regression models. Odds ratios (ORs) and 95% confidence intervals were calculated. RESULTS: Testing for association of previously reported autoimmune disease genetic associations in the initial cohort suggested associations with JIA in 13 distinct loci. Of these, 7 were validated in the replication cohort. Meta-analysis results for the replicating loci included PTPN22 (rs6679677 [OR 1.58, P = 1.98 × 10(-12) ], rs2476601 [OR 1.64, P = 1.90 × 10(-13) ], and rs2488457 [OR 1.32, P = 6.74 × 10(-8) ]), PTPN2 (rs1893217 [OR = 1.33, P = 1.60 × 10(-9) ] and rs7234029 [OR 1.35, P = 1.86 × 10(-10) ]), ADAD1-IL2-IL21 (rs17388568 [OR 1.24, P = 1.13 × 10(-6) ] and rs13143866 [OR 0.83, P = 1.95 × 10(-4) ]), STAT4 (rs3821236 [OR = 1.27, P = 2.36 × 10(-6) ] and rs7574865 [OR = 1.31, P = 2.21 × 10(-6) ]), C12orf30 (rs17696736 [OR = 1.19, P = 2.59 × 10(-5) ]), COG6 (rs7993214 [OR = 0.76, P = 1.10 × 10(-5) ]), and ANGPT1 (rs1010824 [OR = 0.79, P = 2.91 × 10(-4) ]). These polymorphisms have been reported in diseases such as rheumatoid arthritis, type 1 diabetes mellitus, Crohn's disease, and multiple sclerosis. CONCLUSION: General susceptibility loci for autoimmunity are shared across diseases, including JIA, suggesting the potential for common therapeutic targets and mechanisms.

Biologic similarities based on age at onset in oligoarticular and polyarticular subtypes of juvenile idiopathic arthritis.

OBJECTIVE: To explore biologic correlates to age at onset in patients with juvenile idiopathic arthritis (JIA) using peripheral blood mononuclear cell (PBMC) gene expression analysis. METHODS: PBMCs were isolated from 56 healthy controls and 104 patients with recent-onset JIA (39 with persistent oligoarticular JIA, 45 with rheumatoid factor-negative polyarticular JIA, and 20 with systemic JIA). RNA was amplified and labeled using NuGEN Ovation, and gene expression was assessed with Affymetrix HG-U133 Plus 2.0 GeneChips. RESULTS: A total of 832 probe sets revealed gene expression differences (false discovery rate 5%) in PBMCs from children with oligoarticular JIA whose disease began before age 6 years (early-onset disease) compared with those whose disease began at or after age 6 years (late-onset disease). In patients with early-onset disease, there was greater expression of genes related to B cells and less expression of genes related to cells of the myeloid lineage. Support vector machine analyses identified samples from patients with early- or late-onset oligoarticular JIA (with 97% accuracy) or from patients with early- or late-onset polyarticular JIA (with 89% accuracy), but not from patients with systemic JIA or healthy controls. Principal components analysis showed that age at onset was the major classifier of samples from patients with oligoarticular JIA and patients with polyarticular JIA. CONCLUSION: PBMC gene expression analysis reveals biologic differences between patients with early-and late-onset JIA, independent of classification based on the number of joints involved. These data suggest that age at onset may be an important parameter to consider in JIA classification. Furthermore, pathologic mechanisms may vary with age at onset, and understanding these processes may lead to improved treatment of JIA.

A functional RANKL polymorphism associated with younger age at onset of rheumatoid arthritis.

OBJECTIVE: We previously observed the association of the co-occurrence of the HLA-DRB1 shared epitope (SE) and RANKL single-nucleotide polymorphisms (SNPs) with younger age at the onset of rheumatoid arthritis (RA) in 182 rheumatoid factor (RF)-positive European American patients with early-onset RA. The aim of this study was to fine-map the 48-kb RANKL region in the extended cohort of 210 European American RF-positive patients with early RA, to seek replication of RA-associated SNPs in additional RA cohorts of 501 European Americans and 298 African Americans, and to explore the functional consequences of RA-associated SNPs. METHODS: SNP genotyping was conducted using pyrosequencing or TaqMan polymerase chain reaction (PCR) assays. Associations of rs7984870 with RANKL expression in plasma, peripheral blood mononuclear cells, and isolated T cells were quantified using enzyme-linked immunosorbent assay and reverse transcription-PCR. Site-directed mutagenesis of rs7984870 within the 2-kb RANKL promoter was performed to drive the luciferase reporter gene in osteoblast and stromal cell lines. Interaction of DNA and protein was determined by electrophoretic mobility shift assay. RESULTS: A single promoter SNP, rs7984870, was consistently significantly associated with earlier age at the onset of RA in 3 independent seropositive (RF or anti-cyclic citrullinated peptide antibody) RA cohorts but not in seronegative RA patients. The C risk allele of rs7984870 conferred 2-fold higher plasma RANKL levels in RF-positive patients with RA, significantly elevated RANKL messenger RNA expression in activated normal T cells, and increased promoter activity after stimulation in vitro via differential binding to the transcription factor SOX5. CONCLUSION: The RANKL promoter allele that increased transcription levels upon stimulation might promote interaction between activated T cells and dendritic cells, predisposing to a younger age at the onset of RA in seropositive European American and African American patients.

Immature cell populations and an erythropoiesis gene-expression signature in systemic juvenile idiopathic arthritis: implications for pathogenesis.

INTRODUCTION: Previous observations suggest that active systemic juvenile idiopathic arthritis (sJIA) is associated with a prominent erythropoiesis gene-expression signature. The aim of this study was to determine the association of this signature with peripheral blood mononuclear cell (PBMC) subpopulations and its specificity for sJIA as compared with related conditions. METHODS: The 199 patients with JIA (23 sJIA and 176 non-sJIA) and 38 controls were studied. PBMCs were isolated and analyzed for multiple surface antigens with flow cytometry and for gene-expression profiles. The proportions of different PBMC subpopulations were compared among sJIA, non-sJIA patients, and controls and subsequently correlated with the strength of the erythropoiesis signature. Additional gene-expression data from patients with familial hemophagocytic lymphohistiocytosis (FHLH) and from a published sJIA cohort were analyzed to determine whether the erythropoiesis signature was present. RESULTS: Patients with sJIA had significantly increased proportions of immature cell populations, including CD34+ cells, correlating highly with the strength of the erythropoiesis signature. The erythropoiesis signature strongly overlapped with the gene-expression pattern in purified immature erythroid precursors. The expansion of immature cells was most prominently seen in patients with sJIA and anemia, even in the absence of reticulocytosis. Patients with non-sJIA and anemia did not exhibit the erythropoiesis signature. The erythropoiesis signature was found to be prominent in patients with FHLH and in a published cohort of patients with active sJIA, but not in patients with inactive sJIA. CONCLUSIONS: An erythropoiesis signature in active sJIA is associated with the expansion of CD34+ cells, also is seen in some patients with FHLH and infection, and may be an indicator of ineffective erythropoiesis and hemophagocytosis due to hypercytokinemia.

Juvenile idiopathic arthritis and HLA class I and class II interactions and age-at-onset effects.

OBJECTIVE: The aim of this study was to quantitate risk and to examine heterogeneity for HLA at high resolution in patients with the most common subtypes of juvenile idiopathic arthritis (JIA), IgM rheumatoid factor-negative polyarticular JIA and oligoarticular JIA. Use of 4-digit comprehensive HLA typing enabled great precision, and a large cohort allowed for consideration of both age at disease onset and disease subtype. METHODS: Polymerase chain reaction-based high-resolution HLA typing for class I and class II loci was accomplished for 820 patients with JIA and 273 control subjects. Specific HLA epitopes, potential interactions of alleles at specific loci and between loci (accounting for linkage disequilibrium and haplotypic associations), and an assessment of the current International League of Associations for Rheumatology classification criteria were considered. RESULTS: An HLA-DRB1/DQB1 effect was shown to be exclusively attributable to DRB1 and was similar between patients with oligoarticular JIA and a younger subgroup of patients with polyarticular JIA. Furthermore, patients with polyarticular JIA showed age-specific related effects, with disease susceptibility in the group older than age 6 years limited to an effect of the HLA-DRB1*08 haplotype, which is markedly different from the additional susceptibility haplotypes, HLA-DRB1*1103/1104, found in the group with oligoarticular JIA and the group of younger patients with polyarticular JIA. Also in contrast to findings for oligoarticular JIA, patients with polyarticular arthritis had no evidence of an HLA class I effect. Markers associated with a reduced risk of disease included DRB1*1501, DRB1*0401, and DRB1*0701. DRB1*1501 was shown to reduce risk across the whole cohort, whereas DRB1*0401 and DRB1*0701 were protective for selected JIA subtypes. Surprisingly, the disease predisposition mediated by DPB1*0201 in individuals without any disease-predisposing DRB1 alleles was great enough to overcome even the very strong protective effect observed for DRB1*1501. CONCLUSION: Inherited HLA factors in JIA show similarities overall as well as differences between JIA subtypes.

Sequence feature variant type (SFVT) analysis of the HLA genetic association in juvenile idiopathic arthritis.

The immune response HLA class II DRB1 gene provides the major genetic contribution to Juvenile Idiopathic Arthritis (JIA), with a hierarchy of predisposing through intermediate to protective effects. With JIA, and the many other HLA associated diseases, it is difficult to identify the combinations of biologically relevant amino acid (AA) residues directly involved in disease due to the high level of HLA polymorphism, the pattern of AA variability, including varying degrees of linkage disequilibrium (LD), and the fact that most HLA variation occurs at functionally important sites. In a subset of JIA patients with the clinical phenotype oligoarticular-persistent (OP), we have applied a recently developed novel approach to genetic association analyses with genes/proteins sub-divided into biologically relevant smaller sequence features (SFs), and their "alleles" which are called variant types (VTs). With SFVT analysis, association tests are performed on variation at biologically relevant SFs based on structural (e.g., beta-strand 1) and functional (e.g., peptide binding site) features of the protein. We have extended the SFVT analysis pipeline to additionally include pairwise comparisons of DRB1 alleles within serogroup classes, our extension of the Salamon Unique Combinations algorithm, and LD patterns of AA variability to evaluate the SFVT results; all of which contributed additional complementary information. With JIA-OP, we identified a set of single AA SFs, and SFs in which they occur, particularly pockets of the peptide binding site, that account for the major disease risk attributable to HLA DRB1. These are (in numeric order): AAs 13 (pockets 4 and 6), 37 and 57 (both pocket 9), 67 (pocket 7), 74 (pocket 4), and 86 (pocket 1), and to a lesser extent 30 (pockets 6 and 7) and 71 (pockets 4, 5, and 7).

Gene expression signatures in polyarticular juvenile idiopathic arthritis demonstrate disease heterogeneity and offer a molecular classification of disease subsets.

OBJECTIVE: To determine whether peripheral blood mononuclear cells (PBMCs) from children with recent-onset polyarticular juvenile idiopathic arthritis (JIA) exhibit biologically or clinically informative gene expression signatures. METHODS: Peripheral blood samples were obtained from 59 healthy children and 61 children with polyarticular JIA prior to treatment with second-line medications, such as methotrexate or biologic agents. RNA was extracted from isolated mononuclear cells, fluorescence labeled, and hybridized to commercial gene expression microarrays (Affymetrix HG-U133 Plus 2.0). Data were analyzed using analysis of variance at a 5% false discovery rate threshold after robust multichip analysis preprocessing and distance-weighted discrimination normalization. RESULTS: Initial analysis revealed 873 probe sets for genes that were differentially expressed between polyarticular JIA patients and healthy controls. Hierarchical clustering of these probe sets distinguished 3 subgroups within the polyarticular JIA group. Prototypical patients within each subgroup were identified and used to define subgroup-specific gene expression signatures. One of these signatures was associated with monocyte markers, another with transforming growth factor beta-inducible genes, and a third with immediate early genes. Correlation of gene expression signatures with clinical and biologic features of JIA subgroups suggested relevance to aspects of disease activity and supported the division of polyarticular JIA into distinct subsets. CONCLUSION: Gene expression signatures in PBMCs from patients with recent-onset polyarticular JIA reflect discrete disease processes and offer a molecular classification of disease.

Subtype-specific peripheral blood gene expression profiles in recent-onset juvenile idiopathic arthritis.

OBJECTIVE: To identify differences in peripheral blood gene expression between patients with different subclasses of juvenile idiopathic arthritis (JIA) and healthy controls in a multicenter study of patients with recent-onset JIA prior to treatment with disease-modifying antirheumatic drugs (DMARDs) or biologic agents. METHODS: Peripheral blood mononuclear cells (PBMCs) from 59 healthy children and 136 patients with JIA (28 with enthesitis-related arthritis [ERA], 42 with persistent oligoarthritis, 45 with rheumatoid factor [RF]-negative polyarthritis, and 21 with systemic disease) were isolated from whole blood. Poly(A) RNA was labeled using a commercial RNA amplification and labeling system (NuGEN Ovation), and gene expression profiles were obtained using commercial expression microarrays (Affymetrix HG-U133 Plus 2.0). RESULTS: A total of 9,501 differentially expressed probe sets were identified among the JIA subtypes and controls (by analysis of variance; false discovery rate 5%). Specifically, 193, 1,036, 873, and 7,595 probe sets were different in PBMCs from the controls compared with those from the ERA, persistent oligoarthritis, RF-negative polyarthritis, and systemic JIA patients, respectively. In patients with persistent oligoarthritis, RF-negative polyarthritis, and systemic JIA subtypes, up-regulation of genes associated with interleukin-10 (IL-10) signaling was prominent. A hemoglobin cluster was identified that was underexpressed in ERA patients but overexpressed in systemic JIA patients. The influence of JAK/STAT, ERK/MAPK, IL-2, and B cell receptor signaling pathways was evident in patients with persistent oligoarthritis. In systemic JIA, up-regulation of innate immune pathways, including IL-6, Toll-like receptor/IL-1 receptor, and peroxisome proliferator-activated receptor signaling, were noted, along with down-regulation of gene networks related to natural killer cells and T cells. Complement and coagulation pathways were up-regulated in systemic JIA, with a subset of these genes being differentially expressed in other subtypes as well. CONCLUSION: Expression analysis identified differentially expressed genes in PBMCs obtained early in the disease from patients with different subtypes of JIA and in healthy controls, providing evidence of immunobiologic differences between these forms of childhood arthritis.

Association of the IL2RA/CD25 gene with juvenile idiopathic arthritis.

OBJECTIVE: IL2RA/CD25, the gene for interleukin-2 receptor alpha, is emerging as a general susceptibility gene for autoimmune diseases because of its role in the development and function of regulatory T cells and the association of single-nucleotide polymorphisms (SNPs) within this gene with type 1 diabetes mellitus (DM), Graves' disease, rheumatoid arthritis (RA), and multiple sclerosis (MS). The aim of this study was to determine whether SNPs within the IL2RA/CD25 gene are associated with juvenile idiopathic arthritis (JIA). METHODS: Three SNPs within the IL2RA/CD25 gene, that previously showed evidence of an association with either RA, MS, or type 1 DM, were selected for genotyping in UK JIA cases (n=654) and controls (n=3,849). Data for 1 SNP (rs2104286) were also available from North American JIA cases (n=747) and controls (n=1,161). Association analyses were performed using Plink software. Odds ratios (ORs) and 95% confidence intervals (95% CIs) were calculated. RESULTS: SNP rs2104286 within the IL2RA/CD25 gene was significantly associated with UK JIA cases (OR for the allele 0.76 [95% CI 0.66-0.88], P for trend=0.0002). A second SNP (rs41295061) also showed modest evidence for association with JIA (OR 0.80 [95% CI 0.63-1.0], P=0.05). Association with rs2104286 was convincingly replicated in the North American JIA cohort (OR 0.84 [95% CI 0.65-0.99], P for trend=0.05). Meta-analysis of the 2 cohorts yielded highly significant evidence of association with JIA (OR 0.76 [95% CI 0.62-0.88], P=4.9x10(-5)). CONCLUSION: These results provide strong evidence that the IL2RA/CD25 gene represents a JIA susceptibility locus. Further investigation of the gene using both genetic and functional approaches is now required.