Search

Objective: We previously reported an analysis of single-nucleotide polymorphisms (SNPs) in 3 validated European rheumatoid arthritis (RA) susceptibility loci, TAGAP, TNFAIP3, and CCR6, in African American patients with RA. Unexpectedly, the disease-associated alleles were different in African Americans from those in Europeans. In an effort to better define their contribution, we performed additional SNP genotyping in these genes., Methods: Seven SNPs were genotyped in 446 African American patients with RA and in 733 African American control subjects. Differences in minor allele frequency between the RA cases and controls were analyzed after controlling for the global proportion of European admixture, and pairwise linkage disequilibrium (LD) was estimated among the SNPs., Results: Three SNPs were significantly associated with RA: the TNFAIP3 rs719149 A allele (OR 1.22 [95% confidence interval (95% CI) 1.03-1.44], P = 0.02), the TAGAP rs1738074 G allele (OR 0.75 [95% CI 0.63-0.89, P = 0.0012), and the TAGAP rs4709267 G allele (OR 0.74 [95% CI 0.60-0.91], P = 0.004). Pairwise LD between the TAGAP SNPs was low (r(2) = 0.034). The haplotype containing minor alleles for both TAGAP SNPs was uncommon (4.5%). After conditional analysis of each TAGAP SNP, its counterpart remained significantly associated with RA (rs1738074 for rs4709267 P = 0.00001 and rs4709267 for rs1738074 P = 0.00005), suggesting independent effects., Conclusion: SNPs in regulatory regions of TAGAP and an intronic SNP (TNFAIP3) are potential susceptibility loci in African Americans. Pairwise LD, haplotype analysis, and SNP conditioning analysis suggest that these 2 SNPs in TAGAP are independent susceptibility alleles. Additional fine-mapping of this gene and functional genomic studies of these SNPs should provide further insight into the role of these genes in RA., (Copyright © 2012 by the American College of Rheumatology.)

A new study reports the mapping of gene expression in primary immune cell subsets, showing the presence of cell type-specific cis and trans expression quantitative trait loci (eQTLs). The identification of cell type-specific trans-regulated networks can inform functional studies of susceptibility loci identified from genome-wide association studies for human complex diseases.

Mutations in GBA1 gene result in defective acid β-glucosidase and the complex phenotype of Gaucher disease (GD) related to the accumulation of glucosylceramide-laden macrophages. The phenotype is highly variable even among patients harboring identical GBA1 mutations. We hypothesize that modifier gene(s) underlie phenotypic diversity in GD and performed a GWAS study in Ashkenazi Jewish patients with type 1 GD (GD1), homozygous for N370S mutation. Patients were assigned to mild, moderate, or severe disease categories using composite disease severity scoring systems. Whole-genome genotyping for >500,000 SNPs was performed to search for association signals using OQLS algorithm in 139 eligible patients. Several SNPs in linkage disequilibrium within the CLN8 gene locus were associated with the GD1 severity: SNP rs11986414 was associated with GD1 severity at P value 1.26 × 10(-6) . Compared to mild disease, risk allele A at rs11986414 conferred an odds ratio of 3.72 for moderate/severe disease. Loss of function mutations in CLN8 causes neuronal ceroid-lipofuscinosis, but our results indicate that its increased expression may protect against severe GD1. In cultured skin fibroblasts, the relative expression of CLN8 was higher in mild GD compared to severely affected patients, in whom CLN8 risk alleles were overrepresented. In an in vitro cell model of GD, CLN8 expression was increased, which was further enhanced in the presence of bioactive substrate, glucosylsphingosine. Taken together, CLN8 is a candidate modifier gene for GD1 that may function as a protective sphingolipid sensor and/or in glycosphingolipid trafficking. Future studies should explore the role of CLN8 in pathophysiology of GD., (Copyright © 2012 Wiley Periodicals, Inc.)

Rheumatoid arthritis is a common autoimmune disease characterized by chronic inflammation. We report a meta-analysis of genome-wide association studies (GWAS) in a Japanese population including 4,074 individuals with rheumatoid arthritis (cases) and 16,891 controls, followed by a replication in 5,277 rheumatoid arthritis cases and 21,684 controls. Our study identified nine loci newly associated with rheumatoid arthritis at a threshold of P < 5.0 × 10(-8), including B3GNT2, ANXA3, CSF2, CD83, NFKBIE, ARID5B, PDE2A-ARAP1, PLD4 and PTPN2. ANXA3 was also associated with susceptibility to systemic lupus erythematosus (P = 0.0040), and B3GNT2 and ARID5B were associated with Graves' disease (P = 3.5 × 10(-4) and 2.9 × 10(-4), respectively). We conducted a multi-ancestry comparative analysis with a previous meta-analysis in individuals of European descent (5,539 rheumatoid arthritis cases and 20,169 controls). This provided evidence of shared genetic risks of rheumatoid arthritis between the populations.

We have previously shown that rheumatoid arthritis (RA) risk alleles overlap between different ethnic groups. Here, we utilize a multiethnic approach to show that we can effectively discover RA risk alleles. Thirteen putatively associated SNPs that had not yet exceeded genome-wide significance (p < 5 × 10(-8)) in our previous RA genome-wide association study (GWAS) were analyzed in independent sample sets consisting of 4,366 cases and 17,765 controls of European, African American, and East Asian ancestry. Additionally, we conducted an overall association test across all 65,833 samples (a GWAS meta-analysis plus the replication samples). Of the 13 SNPs investigated, four were significantly below the study-wide Bonferroni corrected p value threshold (p < 0.0038) in the replication samples. Two SNPs (rs3890745 at the 1p36 locus [p = 2.3 × 10(-12)] and rs2872507 at the 17q12 locus [p = 1.7 × 10(-9)]) surpassed genome-wide significance in all 16,659 RA cases and 49,174 controls combined. We used available GWAS data to fine map these two loci in Europeans and East Asians, and we found that the same allele conferred risk in both ethnic groups. A series of bioinformatic analyses identified TNFRSF14-MMEL1 at the 1p36 locus and IKZF3-ORMDL3-GSDMB at the 17q12 locus as the genes most likely associated with RA. These findings demonstrate empirically that a multiethnic approach is an effective strategy for discovering RA risk loci, and they suggest that combining GWASs across ethnic groups represents an efficient strategy for gaining statistical power., (Copyright © 2012 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

The genetic association of the major histocompatibility complex (MHC) to rheumatoid arthritis risk has commonly been attributed to alleles in HLA-DRB1. However, debate persists about the identity of the causal variants in HLA-DRB1 and the presence of independent effects elsewhere in the MHC. Using existing genome-wide SNP data in 5,018 individuals with seropositive rheumatoid arthritis (cases) and 14,974 unaffected controls, we imputed and tested classical alleles and amino acid polymorphisms in HLA-A, HLA-B, HLA-C, HLA-DPA1, HLA-DPB1, HLA-DQA1, HLA-DQB1 and HLA-DRB1, as well as 3,117 SNPs across the MHC. Conditional and haplotype analyses identified that three amino acid positions (11, 71 and 74) in HLA-DRβ1 and single-amino-acid polymorphisms in HLA-B (at position 9) and HLA-DPβ1 (at position 9), which are all located in peptide-binding grooves, almost completely explain the MHC association to rheumatoid arthritis risk. This study shows how imputation of functional variation from large reference panels can help fine map association signals in the MHC.

Background: Relatively small, reproductively isolated populations with reduced genetic diversity may have advantages for genomewide association mapping in disease genetics. The Ashkenazi Jewish population represents a unique population for study based on its recent (< 1,000 year) history of a limited number of founders, population bottlenecks and tradition of marriage within the community. We genotyped more than 1,300 Ashkenazi Jewish healthy volunteers from the Hebrew University Genetic Resource with the Illumina HumanOmni1-Quad platform. Comparison of the genotyping data with that of neighboring European and Asian populations enabled the Ashkenazi Jewish-specific component of the variance to be characterized with respect to disease-relevant alleles and pathways., Results: Using clustering, principal components, and pairwise genetic distance as converging approaches, we identified an Ashkenazi Jewish-specific genetic signature that differentiated these subjects from both European and Middle Eastern samples. Most notably, gene ontology analysis of the Ashkenazi Jewish genetic signature revealed an enrichment of genes functioning in transepithelial chloride transport, such as CFTR, and in equilibrioception, potentially shedding light on cystic fibrosis, Usher syndrome and other diseases over-represented in the Ashkenazi Jewish population. Results also impact risk profiles for autoimmune and metabolic disorders in this population. Finally, residual intra-Ashkenazi population structure was minimal, primarily determined by class 1 MHC alleles, and not related to host country of origin., Conclusions: The Ashkenazi Jewish population is of potential utility in disease-mapping studies due to its relative homogeneity and distinct genomic signature. Results suggest that Ashkenazi-associated disease genes may be components of population-specific genomic differences in key functional pathways.

Selective IgA deficiency (IgAD; serum IgA<0.07 g/l) is the most common form of human primary immune deficiency, affecting approximately 1∶600 individuals in populations of Northern European ancestry. The polygenic nature of IgAD is underscored by the recent identification of several new risk genes in a genome-wide association study. Among the characterized susceptibility loci, the association with specific HLA haplotypes represents the major genetic risk factor for IgAD. Despite the robust association, the nature and location of the causal variants in the HLA region remains unknown. To better characterize the association signal in this region, we performed a high-density SNP mapping of the HLA locus and imputed the genotypes of common HLA-B, -DRB1, and -DQB1 alleles in a combined sample of 772 IgAD patients and 1,976 matched controls from 3 independent European populations. We confirmed the complex nature of the association with the HLA locus, which is the result of multiple effects spanning the entire HLA region. The primary association signal mapped to the HLA-DQB1*02 allele in the HLA Class II region (combined P = 7.69×10(-57); OR = 2.80) resulting from the combined independent effects of the HLA-B*0801-DRB1*0301-DQB1*02 and -DRB1*0701-DQB1*02 haplotypes, while additional secondary signals were associated with the DRB1*0102 (combined P = 5.86×10(-17); OR = 4.28) and the DRB1*1501 (combined P = 2.24×10(-35); OR = 0.13) alleles. Despite the strong population-specific frequencies of HLA alleles, we found a remarkable conservation of these effects regardless of the ethnic background, which supports the use of large multi-ethnic populations to characterize shared genetic association signals in the HLA region. We also provide evidence for the location of association signals within the specific extended haplotypes, which will guide future sequencing studies aimed at characterizing the precise functional variants contributing to disease pathogenesis., Competing Interests: RRG, WO, and TWB are full-time employees of Genentech.

To measure the strength of natural selection that acts upon single nucleotide variants (SNVs) in a set of human genes, we calculate the ratio between nonsynonymous SNVs (nsSNVs) per nonsynonymous site and synonymous SNVs (sSNVs) per synonymous site. We transform this ratio with a respective factor f that corrects for the bias of synonymous sites towards transitions in the genetic code and different mutation rates for transitions and transversions. This method approximates the relative density of nsSNVs (rdnsv) in comparison with the neutral expectation as inferred from the density of sSNVs. Using SNVs from a diploid genome and 200 exomes, we apply our method to immune system genes (ISGs), nervous system genes (NSGs), randomly sampled genes (RSGs), and gene ontology annotated genes. The estimate of rdnsv in an individual exome is around 20% for NSGs and 30-40% for ISGs and RSGs. This smaller rdnsv of NSGs indicates overall stronger purifying selection. To quantify the relative shift of nsSNVs towards rare variants, we next fit a linear regression model to the estimates of rdnsv over different SNV allele frequency bins. The obtained regression models show a negative slope for NSGs, ISGs and RSGs, supporting an influence of purifying selection on the frequency spectrum of segregating nsSNVs. The y-intercept of the model predicts rdnsv for an allele frequency close to 0. This parameter can be interpreted as the proportion of nonsynonymous sites where mutations are tolerated to segregate with an allele frequency notably greater than 0 in the population, given the performed normalization of the observed nsSNV to sSNV ratio. A smaller y-intercept is displayed by NSGs, indicating more nonsynonymous sites under strong negative selection. This predicts more monogenically inherited or de-novo mutation diseases that affect the nervous system.

Crohn's disease (CD) is a complex disorder resulting from the interaction of intestinal microbiota with the host immune system in genetically susceptible individuals. The largest meta-analysis of genome-wide association to date identified 71 CD-susceptibility loci in individuals of European ancestry. An important epidemiological feature of CD is that it is 2-4 times more prevalent among individuals of Ashkenazi Jewish (AJ) descent compared to non-Jewish Europeans (NJ). To explore genetic variation associated with CD in AJs, we conducted a genome-wide association study (GWAS) by combining raw genotype data across 10 AJ cohorts consisting of 907 cases and 2,345 controls in the discovery stage, followed up by a replication study in 971 cases and 2,124 controls. We confirmed genome-wide significant associations of 9 known CD loci in AJs and replicated 3 additional loci with strong signal (p<5×10⁻⁶). Novel signals detected among AJs were mapped to chromosomes 5q21.1 (rs7705924, combined p = 2×10⁻⁸; combined odds ratio OR = 1.48), 2p15 (rs6545946, p = 7×10⁻⁹; OR = 1.16), 8q21.11 (rs12677663, p = 2×10⁻⁸; OR = 1.15), 10q26.3 (rs10734105, p = 3×10⁻⁸; OR = 1.27), and 11q12.1 (rs11229030, p = 8×10⁻⁹; OR = 1.15), implicating biologically plausible candidate genes, including RPL7, CPAMD8, PRG2, and PRG3. In all, the 16 replicated and newly discovered loci, in addition to the three coding NOD2 variants, accounted for 11.2% of the total genetic variance for CD risk in the AJ population. This study demonstrates the complementary value of genetic studies in the Ashkenazim., Competing Interests: The authors have declared that no competing interests exist.

Complement receptor 1 (CR1) levels have been associated with malarial susceptibility and/or severity of the disease in different population groups, and CR1 is a receptor for Plasmodium falciparum. In this study, multiple CR1 single-nucleotide polymorphisms (SNPs) showed strong evidence of population differentiation between Sardinian and other European ethnic groups. Cross population algorithms comparing haplotype structure and differences in haplotype and allele frequency distribution provided additional support for natural selection of CR1 in Sardinia. The predominant Sardinian CR1 haplotype included SNPs that are associated with decreased CR1 levels in Europeans and other population groups. Previous studies have shown that the SNPs within the dominant Sardinian haplotype have a significantly higher frequency in a malaria endemic compared with non-endemic regions in India. Together with the historical evidence of the prevalence of malaria in Sardinia, these data support the role of malaria leading to positive selection of this CR1 haplotype in Sardinia.

The immune and coagulation systems are both implicated in the pathogenesis of rheumatoid arthritis (RA). Plasma carboxypeptidase B (CPB), which is activated by the thrombin/thrombomodulin complex, plays a procoagulant role during fibrin clot formation. However, an antiinflammatory role for CPB is suggested by the recent observation that CPB can cleave proinflammatory mediators, such as C5a, bradykinin, and osteopontin. Here, we show that CPB plays a central role in downregulating C5a-mediated inflammatory responses in autoimmune arthritis. CPB deficiency exacerbated inflammatory arthritis in a mouse model of RA, and cleavage of C5a by CPB suppressed the ability of C5a to recruit immune cells in vivo. In human patients with RA, genotyping of nonsynonymous SNPs in the CPB-encoding gene revealed that the allele encoding a CPB variant with longer half-life was associated with a lower risk of developing radiographically severe RA. Functionally, this CPB variant was more effective at abrogating the proinflammatory properties of C5a. Additionally, expression of both CPB and C5a in synovial fluid was higher in patients with RA than in those with osteoarthritis. These findings suggest that CPB plays a critical role in dampening local, C5a-mediated inflammation and represents a molecular link between inflammation and coagulation in autoimmune arthritis.

Protein tyrosine phosphatase nonreceptor type 22 (PTPN22) gene polymorphisms are associated with many autoimmune diseases. The major risk allele encodes an R620W amino acid change that alters B cell receptor (BCR) signaling involved in the regulation of central B cell tolerance. To assess whether this PTPN22 risk allele affects the removal of developing autoreactive B cells, we tested by ELISA the reactivity of recombinant antibodies isolated from single B cells from asymptomatic healthy individuals carrying one or two PTPN22 risk allele(s) encoding the PTPN22 R620W variant. We found that new emigrant/transitional and mature naive B cells from carriers of this PTPN22 risk allele contained high frequencies of autoreactive clones compared with those from non-carriers, revealing defective central and peripheral B cell tolerance checkpoints. Hence, a single PTPN22 risk allele has a dominant effect on altering autoreactive B cell counterselection before any onset of autoimmunity. In addition, gene array experiments analyzing mature naive B cells displaying PTPN22 risk allele(s) revealed that the association strength of PTPN22 for autoimmunity may be due not only to the impaired removal of autoreactive B cells but also to the upregulation of genes such as CD40, TRAF1, and IRF5, which encode proteins that promote B cell activation and have been identified as susceptibility genes associated with autoimmune diseases. These data demonstrate that early B cell tolerance defects in autoimmunity can result from specific polymorphisms and precede the onset of disease.

A marked increase in concomitant autoimmune diseases has previously been noted in patients with myasthenia gravis (MG). We show that these diseases occur both before and after the onset of MG and that the process is not influenced by thymectomy. IgA deficiency (IgAD), which is strongly associated with the same HLA haplotype as early onset MG, has recently been suggested to be an autoimmune disease. However, there was no increase in the prevalence of IgAD in a large cohort of Swedish MG patients., (Copyright © 2011 Elsevier B.V. All rights reserved.)

The aim of this study was to determine, through a genome-wide association study (GWAS), the genetic components contributing to different clinical sub-phenotypes of systemic sclerosis (SSc). We considered limited (lcSSc) and diffuse (dcSSc) cutaneous involvement, and the relationships with presence of the SSc-specific auto-antibodies, anti-centromere (ACA), and anti-topoisomerase I (ATA). Four GWAS cohorts, comprising 2,296 SSc patients and 5,171 healthy controls, were meta-analyzed looking for associations in the selected subgroups. Eighteen polymorphisms were further tested in nine independent cohorts comprising an additional 3,175 SSc patients and 4,971 controls. Conditional analysis for associated SNPs in the HLA region was performed to explore their independent association in antibody subgroups. Overall analysis showed that non-HLA polymorphism rs11642873 in IRF8 gene to be associated at GWAS level with lcSSc (P = 2.32×10(-12), OR = 0.75). Also, rs12540874 in GRB10 gene (P = 1.27 × 10(-6), OR = 1.15) and rs11047102 in SOX5 gene (P = 1.39×10(-7), OR = 1.36) showed a suggestive association with lcSSc and ACA subgroups respectively. In the HLA region, we observed highly associated allelic combinations in the HLA-DQB1 locus with ACA (P = 1.79×10(-61), OR = 2.48), in the HLA-DPA1/B1 loci with ATA (P = 4.57×10(-76), OR = 8.84), and in NOTCH4 with ACA P = 8.84×10(-21), OR = 0.55) and ATA (P = 1.14×10(-8), OR = 0.54). We have identified three new non-HLA genes (IRF8, GRB10, and SOX5) associated with SSc clinical and auto-antibody subgroups. Within the HLA region, HLA-DQB1, HLA-DPA1/B1, and NOTCH4 associations with SSc are likely confined to specific auto-antibodies. These data emphasize the differential genetic components of subphenotypes of SSc., Competing Interests: The authors have declared that no competing interests exist.

A common allele at the TAGAP gene locus demonstrates a suggestive, but not conclusive association with risk of rheumatoid arthritis (RA). To fine map the locus, we conducted comprehensive imputation of CEU HapMap single-nucleotide polymorphisms (SNPs) in a genome-wide association study (GWAS) of 5,500 RA cases and 22,621 controls (all of European ancestry). After controlling for population stratification with principal components analysis, the strongest signal of association was to an imputed SNP, rs212389 (P=3.9 × 10(-8), odds ratio=0.87). This SNP remained highly significant upon conditioning on the previous RA risk variant (rs394581, P=2.2 × 10(-5)) or on a SNP previously associated with celiac disease and type I diabetes (rs1738074, P=1.7 × 10(-4)). Our study has refined the TAGAP signal of association to a single haplotype in RA, and in doing so provides conclusive statistical evidence that the TAGAP locus is associated with RA risk. Our study also underscores the utility of comprehensive imputation in large GWAS data sets to fine map disease risk alleles.

Chronic lymphocytic leukemia (CLL) is an incurable adult disease of unknown etiology. Understanding the biology of CLL cells, particularly cell maturation and growth in vivo, has been impeded by lack of a reproducible adoptive transfer model. We report a simple, reproducible system in which primary CLL cells proliferate in nonobese diabetes/severe combined immunodeficiency/γc(null) mice under the influence of activated CLL-derived T lymphocytes. By co-transferring autologous T lymphocytes, activated in vivo by alloantigens, the survival and growth of primary CFSE-labeled CLL cells in vivo is achieved and quantified. Using this approach, we have identified key roles for CD4(+) T cells in CLL expansion, a direct link between CD38 expression by leukemic B cells and their activation, and support for CLL cells preferentially proliferating in secondary lymphoid tissues. The model should simplify analyzing kinetics of CLL cells in vivo, deciphering involvement of nonleukemic elements and nongenetic factors promoting CLL cell growth, identifying and characterizing potential leukemic stem cells, and permitting preclinical studies of novel therapeutics. Because autologous activated T lymphocytes are 2-edged swords, generating unwanted graph-versus-host and possibly autologous antitumor reactions, the model may also facilitate analyses of T-cell populations involved in immune surveillance relevant to hematopoietic transplantation and tumor cytoxicity.

Objective: To perform a genome-wide association study (GWAS) in Koreans in order to identify susceptibility loci for rheumatoid arthritis (RA)., Methods: We generated high-quality genotypes for 441,398 single-nucleotide polymorphisms (SNPs) in 801 RA cases and 757 controls. We then tested 79 markers from 46 loci for replication in an independent sample of 718 RA cases and 719 controls., Results: Genome-wide significance (P < 5 × 10(-08) ) was attained by markers from the major histocompatibility complex region and from the PADI4 gene. The replication data showed nominal association signals (P < 5 × 10(-02) ) for markers from 11 of the 46 replicated loci, greatly exceeding random expectation. Genes that were most significant in the replication stage and in the combined analysis include the known European RA loci BLK, AFF3, and CCL21. Thus, in addition to the previously associated STAT4 alleles, variants at these three loci may contribute to RA not only among Europeans, but also among Asians. In addition, we observed replication signals near the genes PTPN2, FLI1, ARHGEF3, LCP2, GPR137B, TRHDE, and CGA1. Based on the excess of small P values in the replication stage study, we estimate that more than half of these loci are genuine RA susceptibility genes. Finally, we systematically analyzed the presence of association signals in Koreans at established European RA loci, which showed a significant enrichment of European RA loci among the Korean RA loci., Conclusion: Genetic risk for RA involves both population-specific loci as well as many shared genetic susceptibility loci in comparisons of Asian and European populations., (Copyright © 2011 by the American College of Rheumatology.)

The major histocompatibility complex (MHC) class II transactivator gene (CIITA) encodes an important transcription factor regulating genes required for human leukocyte antigen (HLA) class II MHC-restricted antigen presentation. MHC genes, particularly HLA class II, are strongly associated with risk of developing rheumatoid arthritis (RA). Given the strong biological relationship between CIITA and HLA class II genes, a comprehensive investigation of CIITA variation in RA was conducted. This study tested 31 CIITA single-nucleotide polymorphisms in 2542 RA cases and 3690 controls (N=6232). All individuals were of European ancestry, as determined by ancestry informative genetic markers. No evidence for association between CIITA variation and RA was observed after a correction for multiple testing was applied. This is the largest study to fully characterize common genetic variation in CIITA, including an assessment of haplotypes. Results exclude even a modest role for common CIITA polymorphisms in susceptibility to RA.

Systemic lupus erythematosus (SLE) is a clinically heterogeneous, systemic autoimmune disease characterized by autoantibody formation. Previously published genome-wide association studies (GWAS) have investigated SLE as a single phenotype. Therefore, we conducted a GWAS to identify genetic factors associated with anti-dsDNA autoantibody production, a SLE-related autoantibody with diagnostic and clinical importance. Using two independent datasets, over 400,000 single nucleotide polymorphisms (SNPs) were studied in a total of 1,717 SLE cases and 4,813 healthy controls. Anti-dsDNA autoantibody positive (anti-dsDNA +, n = 811) and anti-dsDNA autoantibody negative (anti-dsDNA -, n = 906) SLE cases were compared to healthy controls and to each other to identify SNPs associated specifically with these SLE subtypes. SNPs in the previously identified SLE susceptibility loci STAT4, IRF5, ITGAM, and the major histocompatibility complex were strongly associated with anti-dsDNA + SLE. Far fewer and weaker associations were observed for anti-dsDNA - SLE. For example, rs7574865 in STAT4 had an OR for anti-dsDNA + SLE of 1.77 (95% CI 1.57-1.99, p = 2.0E-20) compared to an OR for anti-dsDNA - SLE of 1.26 (95% CI 1.12-1.41, p = 2.4E-04), with p(heterogeneity)<0.0005. SNPs in the SLE susceptibility loci BANK1, KIAA1542, and UBE2L3 showed evidence of association with anti-dsDNA + SLE and were not associated with anti-dsDNA - SLE. In conclusion, we identified differential genetic associations with SLE based on anti-dsDNA autoantibody production. Many previously identified SLE susceptibility loci may confer disease risk through their role in autoantibody production and be more accurately described as autoantibody propensity loci. Lack of strong SNP associations may suggest that other types of genetic variation or non-genetic factors such as environmental exposures have a greater impact on susceptibility to anti-dsDNA - SLE., Competing Interests: The authors have declared that no competing interests exist.

Systemic lupus erythematosus (SLE, MIM152700) is an autoimmune disease characterized by self-reactive antibodies resulting in systemic inflammation and organ failure. TNFAIP3, encoding the ubiquitin-modifying enzyme A20, is an established susceptibility locus for SLE. By fine mapping and genomic re-sequencing in ethnically diverse populations, we fully characterized the TNFAIP3 risk haplotype and identified a TT>A polymorphic dinucleotide (deletion T followed by a T to A transversion) associated with SLE in subjects of European (P = 1.58 × 10(-8), odds ratio = 1.70) and Korean (P = 8.33 × 10(-10), odds ratio = 2.54) ancestry. This variant, located in a region of high conservation and regulatory potential, bound a nuclear protein complex composed of NF-κB subunits with reduced avidity. Further, compared with the non-risk haplotype, the haplotype carrying this variant resulted in reduced TNFAIP3 mRNA and A20 protein expression. These results establish this TT>A variant as the most likely functional polymorphism responsible for the association between TNFAIP3 and SLE.

Systemic lupus erythematosus (SLE) is a genetically complex disease with heterogeneous clinical manifestations. Recent studies have greatly expanded the number of established SLE risk alleles, but the distribution of multiple risk alleles in cases versus controls and their relationship to subphenotypes have not been studied. We studied 22 SLE susceptibility polymorphisms with previous genome-wide evidence of association (p < 5 x 10⁻¹²⁸) in 1919 SLE cases from 9 independent Caucasian SLE case series and 4813 independent controls. The mean number of risk alleles in cases was 15.1 (SD 3.1) while the mean in controls was 13.1 (SD 2.8), with trend p = 4 x 10⁻⁸. We defined a genetic risk score (GRS) for SLE as the number of risk alleles with each weighted by the SLE risk odds ratio (OR). The OR for high-low GRS tertiles, adjusted for intra-European ancestry, sex, and parent study, was 4.4 (95% CI 3.8-5.1). We studied associations of individual SNPs and the GRS with clinical manifestations for the cases: age at diagnosis, the 11 American College of Rheumatology classification criteria, and double-stranded DNA antibody (anti-dsDNA) production. Six subphenotypes were significantly associated with the GRS, most notably anti-dsDNA (OR(high-low) = 2.36, p = 9e-9), the immunologic criterion (OR(high-low) = 2.23, p = 3e-7), and age at diagnosis (OR(high-low) = 1.45, p = 0.0060). Finally, we developed a subphenotype-specific GRS (sub-GRS) for each phenotype with more power to detect cumulative genetic associations. The sub-GRS was more strongly associated than any single SNP effect for 5 subphenotypes (the above plus hematologic disorder and oral ulcers), while single loci are more significantly associated with renal disease (HLA-DRB1, OR = 1.37, 95% CI 1.14-1.64) and arthritis (ITGAM, OR = 0.72, 95% CI 0.59-0.88). We did not observe significant associations for other subphenotypes, for individual loci or the sub-GRS. Thus our analysis categorizes SLE subphenotypes into three groups: those having cumulative, single, and no known genetic association with respect to the currently established SLE risk loci., Competing Interests: RRG and TWB are full-time employees of Genentech.

Affected relatives are essential for pedigree linkage analysis, however, they cause a violation of the independent sample assumption in case-control association studies. To avoid the correlation between samples, a common practice is to take only one affected sample per pedigree in association analysis. Although several methods exist in handling correlated samples, they are still not widely used in part because these are not easily implemented, or because they are not widely known. We advocate the effective sample size method as a simple and accessible approach for case-control association analysis with correlated samples. This method modifies the chi-square test statistic, p-value, and 95% confidence interval of the odds-ratio by replacing the apparent number of allele or genotype counts with the effective ones in the standard formula, without the need for specialized computer programs. We present a simple formula for calculating effective sample size for many types of relative pairs and relative sets. For allele frequency estimation, the effective sample size method captures the variance inflation exactly. For genotype frequency, simulations showed that effective sample size provides a satisfactory approximation. A gene which is previously identified as a type 1 diabetes susceptibility locus, the interferon-induced helicase gene (IFIH1), is shown to be significantly associated with rheumatoid arthritis when the effective sample size method is applied. This significant association is not established if only one affected sib per pedigree were used in the association analysis. Relationship between the effective sample size method and other methods - the generalized estimation equation, variance of eigenvalues for correlation matrices, and genomic controls - are discussed., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

Epidemiology and candidate gene studies indicate a shared genetic basis for celiac disease (CD) and rheumatoid arthritis (RA), but the extent of this sharing has not been systematically explored. Previous studies demonstrate that 6 of the established non-HLA CD and RA risk loci (out of 26 loci for each disease) are shared between both diseases. We hypothesized that there are additional shared risk alleles and that combining genome-wide association study (GWAS) data from each disease would increase power to identify these shared risk alleles. We performed a meta-analysis of two published GWAS on CD (4,533 cases and 10,750 controls) and RA (5,539 cases and 17,231 controls). After genotyping the top associated SNPs in 2,169 CD cases and 2,255 controls, and 2,845 RA cases and 4,944 controls, 8 additional SNPs demonstrated P<5 × 10(-8) in a combined analysis of all 50,266 samples, including four SNPs that have not been previously confirmed in either disease: rs10892279 near the DDX6 gene (P(combined) =  1.2 × 10(-12)), rs864537 near CD247 (P(combined) =  2.2 × 10(-11)), rs2298428 near UBE2L3 (P(combined) =  2.5 × 10(-10)), and rs11203203 near UBASH3A (P(combined) =  1.1 × 10(-8)). We also confirmed that 4 gene loci previously established in either CD or RA are associated with the other autoimmune disease at combined P<5 × 10(-8) (SH2B3, 8q24, STAT4, and TRAF1-C5). From the 14 shared gene loci, 7 SNPs showed a genome-wide significant effect on expression of one or more transcripts in the linkage disequilibrium (LD) block around the SNP. These associations implicate antigen presentation and T-cell activation as a shared mechanism of disease pathogenesis and underscore the utility of cross-disease meta-analysis for identification of genetic risk factors with pleiotropic effects between two clinically distinct diseases., Competing Interests: YL has employment and stock interest in Celera.

Systemic lupus erythematosus (SLE) is considered to be the prototypic autoimmune disease, with a complex genetic architecture influenced by environmental factors. We sought to replicate a putative association at 11p13 not yet exceeding genome-wide significance (p < 5 × 10(-8)) identified in a genome-wide association study (GWAS). Our GWA scan identified two intergenic SNPs located between PDHX and CD44 showing suggestive evidence of association with SLE in cases of European descent (rs2732552, p = 0.004, odds ratio [OR] = 0.78; rs387619, p = 0.003, OR = 0.78). The replication cohort consisted of >15,000 subjects, including 3562 SLE cases and 3491 controls of European ancestry, 1527 cases and 1811 controls of African American (AA) descent, and 1265 cases and 1260 controls of Asian origin. We observed robust association at both rs2732552 (p = 9.03 × 10(-8), OR = 0.83) and rs387619 (p = 7.7 × 10(-7), OR = 0.83) in the European samples with p(meta) = 1.82 × 10(-9) for rs2732552. The AA and Asian SLE cases also demonstrated association at rs2732552 (p = 5 × 10(-3), OR = 0.81 and p = 4.3 × 10(-4), OR = 0.80, respectively). A meta-analysis of rs2732552 for all racial and ethnic groups studied produced p(meta) = 2.36 × 10(-13). This locus contains multiple regulatory sites that could potentially affect expression and functions of CD44, a cell-surface glycoprotein influencing immunologic, inflammatory, and oncologic phenotypes, or PDHX, a subunit of the pyruvate dehydrogenase complex.

Selective immunoglobulin A deficiency (IgAD) is the most common primary immunodeficiency in Caucasians. It has previously been suggested to be associated with a variety of concomitant autoimmune diseases. In this review, we present data on the prevalence of IgAD in patients with Graves disease (GD), systemic lupus erythematosus (SLE), type 1 diabetes (T1D), celiac disease (CD), myasthenia gravis (MG) and rheumatoid arthritis (RA) on the basis of both our own recent large-scale screening results and literature data. Genetic factors are important for the development of both IgAD and various autoimmune disorders, including GD, SLE, T1D, CD, MG and RA, and a strong association with the major histocompatibility complex (MHC) region has been reported. In addition, non-MHC genes, such as interferon-induced helicase 1 (IFIH1) and c-type lectin domain family 16, member A (CLEC16A), are also associated with the development of IgAD and some of the above diseases. This indicates a possible common genetic background. In this review, we present suggestive evidence for a shared genetic predisposition between these disorders.

Prevention of disease can in principle be accomplished by identification of environmental and/or lifestyle risk and protective factors followed by public health measures (such as for smoking and lung cancer), or by modification of the individual's reactions to disease-inducing factors (such as in vaccinations against microbes). This review discusses both options based on emerging understanding of aetiologies in inflammatory rheumatic diseases such as rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE). The major current opportunity for public health-based prevention lies in avoiding smoking. In RA, recent studies have calculated that, in Sweden (a country characterised by a low frequency of smoking), 20% of all RA cases and 33% of all cases of ACPA-positive RA would not have occurred in a smoke-free society. Smoking is also a major risk factor for SLE but no population attribution is yet available. New avenues for individualised and biology-based prevention are provided by the demonstration that several autoimmune rheumatic diseases are preceded by emergence of subclinical autoimmunity followed by laboratory-based signs of inflammation and finally overt disease. Examples of this process are provided from studies of autoimmunity to citrullinated proteins (in RA), to dsDNA (in SLE in general) and to Ro52 epitopes (in the case of neonatal heart block). The recognition of this sequence of events provides opportunities to intervene specifically and potentially curatively before onset of full-blown disease. Such prevention can be accomplished by modification of inciting antigens (environment), by modification of immunity (more or less specific immunomodulation) or by modification of specific gene functions. In all cases, prevention will be different in different subsets of disease and differ at different time points of disease development. Thus, the road map towards prevention of autoimmune rheumatic diseases includes increased understanding of how genes, environment and immunity interact.

To pinpoint true positive single-nucleotide polymorphism (SNP) associations in a genome-wide association study (GWAS) of rheumatoid arthritis (RA), we categorize genetic loci by external knowledge. We test both the 'enrichment of associated loci' in a locus category and the 'combined association' of a locus category. The former is quantified by the odds ratio for the presence of SNP associations at the loci of a category, whereas the latter is quantified by the number of loci in a category that have SNP associations. These measures are compared with their expected values as obtained from the permutation of the affection status. To account for linkage disequilibrium (LD) among SNPs, we view each LD block as a genetic locus. Positional candidates were defined as loci implicated by earlier GWAS results, whereas functional candidates were defined by annotations regarding the molecular roles of genes, such as gene ontology categories. As expected, immune-related categories show the largest enrichment signal, although it is not very strong. The intersection of positional and functional candidate information predicts novel RA loci near the genes TEC/TXK, MBL2 and PIK3R1/CD180. Notably, a combined association signal is not only produced by immune-related categories, but also by most other categories and even randomly defined categories. The unspecific quality of these signals limits the possible conclusions from combined association tests. It also reduces the magnitude of enrichment test results. These unspecific signals might result from common variants of small effect and hardly concentrated in candidate categories, or an inflated size of associated regions from weak LD with infrequent mutations.

To understand the genetic predisposition to selective immunoglobulin A deficiency (IgAD), we performed a genome-wide association study in 430 affected individuals (cases) from Sweden and Iceland and 1,090 ethnically matched controls, and we performed replication studies in two independent European cohorts. In addition to the known association of HLA with IgAD, we identified association with a nonsynonymous variant in IFIH1 (rs1990760G>A, P = 7.3 x 10(-10)) which was previously associated with type 1 diabetes and systemic lupus erythematosus. Variants in CLEC16A, another known autoimmunity locus, showed suggestive evidence for association (rs6498142C>G, P = 1.8 x 10(-7)), and 29 additional loci were identified with P < 5 x 10(-5). A survey in IgAD of 118 validated non-HLA autoimmunity loci indicated a significant enrichment for association with autoimmunity loci as compared to non-autoimmunity loci (P = 9.0 x 10(-4)) or random SNPs across the genome (P < 0.0001). These findings support the hypothesis that autoimmune mechanisms may contribute to the pathogenesis of IgAD.

CLEC16A, a putative immunoreceptor, was recently established as a susceptibility locus for type I diabetes and multiple sclerosis. Subsequently, associations between CLEC16A and rheumatoid arthritis (RA), Addison's disease and Crohn's disease have been reported. A large comprehensive and independent investigation of CLEC16A variation in RA was pursued. This study tested 251 CLEC16A single-nucleotide polymorphisms in 2542 RA cases (85% anti-cyclic citrullinated peptide (anti-CCP) positive) and 2210 controls (N=4752). All individuals were of European ancestry, as determined by ancestry informative genetic markers. No evidence for significant association between CLEC16A variation and RA was observed. This is the first study to fully characterize common genetic variation in CLEC16A including assessment of haplotypes and gender-specific effects. The previously reported association between RA and rs6498169 was not replicated. Results show that CLEC16A does not have a prominent function in susceptibility to anti-CCP-positive RA.

Previous work has demonstrated that Northern and Southern European ancestries are associated with specific systemic lupus erythematosus (SLE) manifestations. In this study, 1855 SLE cases of European descent were genotyped for 4965 single-nucleotide polymorphisms and principal components analysis of genotype information was used to define population substructure. The first principal component (PC1) distinguished Northern from Southern European ancestry, PC2 differentiated Eastern from Western European ancestry and PC3 delineated Ashkenazi Jewish ancestry. Compared with Northern European ancestry, Southern European ancestry was associated with autoantibody production (odds ratio (OR)=1.40, 95% confidence interval (CI) 1.07-1.83) and renal involvement (OR 1.41, 95% CI 1.06-1.87), and was protective for discoid rash (OR=0.51, 95% CI 0.32-0.82) and photosensitivity (OR=0.74, 95% CI 0.56-0.97). Both serositis (OR=1.46, 95% CI 1.12-1.89) and autoantibody production (OR=1.38, 95% CI 1.06-1.80) were associated with Western compared to Eastern European ancestry. Ashkenazi Jewish ancestry was protective against neurologic manifestations of SLE (OR=0.62, 95% CI 0.40-0.94). Homogeneous clusters of cases defined by multiple PCs demonstrated stronger phenotypic associations. Genetic ancestry may contribute to the development of SLE endophenotypes and should be accounted for in genetic studies of disease characteristics.

A genome-wide association screen for primary biliary cirrhosis risk alleles was performed in an Italian cohort. The results from the Italian cohort replicated IL12A and IL12RB associations, and a combined meta-analysis using a Canadian dataset identified newly associated loci at SPIB (P = 7.9 x 10(-11), odds ratio (OR) = 1.46), IRF5-TNPO3 (P = 2.8 x 10(-10), OR = 1.63) and 17q12-21 (P = 1.7 x 10(-10), OR = 1.38).

Multiple sclerosis (MS) is a common demyelinating disease of the central nervous system mediated by autoimmune and neurodegenerative pathogenic mechanisms. Multiple genes account for its moderate heritability, but the only genetic region shown to have a large replicable effect on MS susceptibility is the major histocompatibility complex (MHC). Strong linkage disequilibrium (LD) across the MHC has made it difficult to fully characterize individual genetic contributions of this region to MS risk in previous studies. African Americans are at a lower risk for MS when compared with northern Europeans and Americans of European descent, but greater haplotypic diversity and distinct patterns of LD suggest that this population may be particularly informative for fine-mapping efforts. To examine the role of the MHC in African American MS, a case-control association study was performed with 499 African American MS patients and 750 African American controls that were genotyped for 6040 MHC region single nucleotide polymorphisms (SNPs). A replication data set consisting of 451 African American patients and 718 African American controls was genotyped for selected SNPs. Two MHC class II SNPs, rs2647040 and rs3135021, were significant in the replication cohort and partially tagged DRB1*15 alleles. Surprisingly, in comparison to similar studies of individuals of European descent, the MHC seems to play a smaller role in MS susceptibility in African Americans, consistent with pervasive genetic heterogeneity across ancestral groups, and may explain the difference in MS susceptibility between African Americans and individuals of European descent.

We genotyped individuals with primary biliary cirrhosis and unaffected controls for suggestive risk loci (genome-wide association P < 1 x 10(-4)) identified in a previous genome-wide association study. Combined analysis of the genome-wide association and replication datasets identified IRF5-TNPO3 (combined P = 8.66 x 10(-13)), 17q12-21 (combined P = 3.50 x 10(-13)) and MMEL1 (combined P = 3.15 x 10(-8)) as new primary biliary cirrhosis susceptibility loci. Fine-mapping studies showed that a single variant accounts for the IRF5-TNPO3 association. As these loci are implicated in other autoimmune conditions, these findings confirm genetic overlap among such diseases.

Sialic acid acetylesterase (SIAE) is an enzyme that negatively regulates B lymphocyte antigen receptor signalling and is required for the maintenance of immunological tolerance in mice. Heterozygous loss-of-function germline rare variants and a homozygous defective polymorphic variant of SIAE were identified in 24/923 subjects of European origin with relatively common autoimmune disorders and in 2/648 controls of European origin. All heterozygous loss-of-function SIAE mutations tested were capable of functioning in a dominant negative manner. A homozygous secretion-defective polymorphic variant of SIAE was catalytically active, lacked the ability to function in a dominant negative manner, and was seen in eight autoimmune subjects but in no control subjects. The odds ratio for inheriting defective SIAE alleles was 8.6 in all autoimmune subjects, 8.3 in subjects with rheumatoid arthritis, and 7.9 in subjects with type I diabetes. Functionally defective SIAE rare and polymorphic variants represent a strong genetic link to susceptibility in relatively common human autoimmune disorders.

Alopecia areata (AA) is among the most highly prevalent human autoimmune diseases, leading to disfiguring hair loss due to the collapse of immune privilege of the hair follicle and subsequent autoimmune attack. The genetic basis of AA is largely unknown. We undertook a genome-wide association study (GWAS) in a sample of 1,054 cases and 3,278 controls and identified 139 single nucleotide polymorphisms that are significantly associated with AA (P

Objective: Anti-tumor necrosis factor alpha (anti-TNF) therapy is a mainstay of treatment in rheumatoid arthritis (RA). The aim of the present study was to test established RA genetic risk factors to determine whether the same alleles also influence the response to anti-TNF therapy., Methods: A total of 1,283 RA patients receiving etanercept, infliximab, or adalimumab therapy were studied from among an international collaborative consortium of 9 different RA cohorts. The primary end point compared RA patients with a good treatment response according to the European League Against Rheumatism (EULAR) response criteria (n = 505) with RA patients considered to be nonresponders (n = 316). The secondary end point was the change from baseline in the level of disease activity according to the Disease Activity Score in 28 joints (triangle upDAS28). Clinical factors such as age, sex, and concomitant medications were tested as possible correlates of treatment response. Thirty-one single-nucleotide polymorphisms (SNPs) associated with the risk of RA were genotyped and tested for any association with treatment response, using univariate and multivariate logistic regression models., Results: Of the 31 RA-associated risk alleles, a SNP at the PTPRC (also known as CD45) gene locus (rs10919563) was associated with the primary end point, a EULAR good response versus no response (odds ratio [OR] 0.55, P = 0.0001 in the multivariate model). Similar results were obtained using the secondary end point, the triangle upDAS28 (P = 0.0002). There was suggestive evidence of a stronger association in autoantibody-positive patients with RA (OR 0.55, 95% confidence interval [95% CI] 0.39-0.76) as compared with autoantibody-negative patients (OR 0.90, 95% CI 0.41-1.99)., Conclusion: Statistically significant associations were observed between the response to anti-TNF therapy and an RA risk allele at the PTPRC gene locus. Additional studies will be required to replicate this finding in additional patient collections.

Background: In Northern European descended populations, genetic susceptibility for multiple sclerosis (MS) is associated with alleles of the human leukocyte antigen (HLA) Class II gene DRB1. Whether other major histocompatibility complex (MHC) genes contribute to MS susceptibility is controversial., Methodology/principal Findings: A case control analysis was performed using 958 single nucleotide polymorphisms (SNPs) spanning the MHC assayed in two independent datasets. The discovery dataset consisted of 1,018 cases and 1,795 controls and the replication dataset was composed of 1,343 cases and 1,379 controls. The most significantly MS-associated SNP in the discovery dataset was rs3135391, a Class II SNP known to tag the HLA-DRB1*15:01 allele, the primary MS susceptibility allele in the MHC (O.R. = 3.04, p < 1 x 10(-78)). To control for the effects of the HLA-DRB1*15:01 haplotype, case control analysis was performed adjusting for this HLA-DRB1*15:01 tagging SNP. After correction for multiple comparisons (false discovery rate = .05) 52 SNPs in the Class I, II and III regions were significantly associated with MS susceptibility in both datasets using the Cochran Armitage trend test. The discovery and replication datasets were merged and subjects carrying the HLA-DRB1*15:01 tagging SNP were excluded. Association tests showed that 48 of the 52 replicated SNPs retained significant associations with MS susceptibility independently of the HLA-DRB1*15:01 as defined by the tagging SNP. 20 Class I SNPs were associated with MS susceptibility with p-values < or = 1 x 10(-8). The most significantly associated SNP was rs4959039, a SNP in the downstream un-translated region of the non-classical HLA-G gene (Odds ratio 1.59, 95% CI 1.40, 1.81, p = 8.45 x 10(-13)) and is in linkage disequilibrium with several nearby SNPs. Logistic regression modeling showed that this SNP's contribution to MS susceptibility was independent of the Class II and Class III SNPs identified in this screen., Conclusions: A MHC Class I locus contributes to MS susceptibility independently of the HLA-DRB1*15:01 haplotype.

To identify new genetic risk factors for rheumatoid arthritis, we conducted a genome-wide association study meta-analysis of 5,539 autoantibody-positive individuals with rheumatoid arthritis (cases) and 20,169 controls of European descent, followed by replication in an independent set of 6,768 rheumatoid arthritis cases and 8,806 controls. Of 34 SNPs selected for replication, 7 new rheumatoid arthritis risk alleles were identified at genome-wide significance (P < 5 x 10(-8)) in an analysis of all 41,282 samples. The associated SNPs are near genes of known immune function, including IL6ST, SPRED2, RBPJ, CCR6, IRF5 and PXK. We also refined associations at two established rheumatoid arthritis risk loci (IL2RA and CCL21) and confirmed the association at AFF3. These new associations bring the total number of confirmed rheumatoid arthritis risk loci to 31 among individuals of European ancestry. An additional 11 SNPs replicated at P < 0.05, many of which are validated autoimmune risk alleles, suggesting that most represent genuine rheumatoid arthritis risk alleles.

Systemic sclerosis (SSc) is an autoimmune disease characterized by fibrosis of the skin and internal organs that leads to profound disability and premature death. To identify new SSc susceptibility loci, we conducted the first genome-wide association study in a population of European ancestry including a total of 2,296 individuals with SSc and 5,171 controls. Analysis of 279,621 autosomal SNPs followed by replication testing in an independent case-control set of European ancestry (2,753 individuals with SSc (cases) and 4,569 controls) identified a new susceptibility locus for systemic sclerosis at CD247 (1q22-23, rs2056626, P = 2.09 x 10(-7) in the discovery samples, P = 3.39 x 10(-9) in the combined analysis). Additionally, we confirm and firmly establish the role of the MHC (P = 2.31 x 10(-18)), IRF5 (P = 1.86 x 10(-13)) and STAT4 (P = 3.37 x 10(-9)) gene regions as SSc genetic risk factors.

Investigating genetic interactions (epistasis) has proven difficult despite the recent advances of both laboratory methods and statistical developments. With no 'best' statistical approach available, combining several analytical methods may be optimal for detecting epistatic interactions. Using a multi-stage analysis that incorporated supervised machine learning and methods of association testing, we investigated epistatic interactions with a well-established genetic factor (PTPN22 1858T) in a complex autoimmune disease (rheumatoid arthritis (RA)). Our analysis consisted of four principal stages: Stage I (data reduction)-identifying candidate chromosomal regions in 292 affected sibling pairs, by predicting PTPN22 concordance using multipoint identity-by-descent probabilities and a supervised machine learning algorithm (Random Forests); Stage II (extension analysis)-testing detailed genetic data within candidate chromosomal regions for epistasis with PTPN22 1858T in 677 cases and 750 controls using logistic regression; Stage III (replication analysis)-confirmation of epistatic interactions in 947 cases and 1756 controls; Stage IV (combined analysis)-a pooled analysis including all 1624 RA cases and 2506 control subjects for final estimates of effect size. A total of seven replicating epistatic interactions were identified. SNP variants within CDH13, MYO3A, CEP72 and near WFDC1 showed significant evidence for interaction with PTPN22, affecting susceptibility to RA.

Introduction: A deletion polymorphism in glutathione S-transferase Mu-1 (GSTM1-null) has previously been implicated to play a role in rheumatoid arthritis (RA) risk and progression, although no prior investigations have examined its associations with anticitrullinated protein antibody (ACPA) positivity. The purpose of this study was to examine the associations of GSTM1-null with ACPA positivity in RA and to assess for evidence of interaction between GSTM1 and HLA-DRB1 shared epitope (SE)., Methods: Associations of GSTM1-null with ACPA positivity were examined separately in two RA cohorts, the Veterans Affairs Rheumatoid Arthritis (VARA) registry (n = 703) and the Study of New-Onset RA (SONORA; n = 610). Interactions were examined by calculating an attributable proportion (AP) due to interaction., Results: A majority of patients in the VARA registry (76%) and SONORA (69%) were positive for ACPA with a similar frequency of GSTM1-null (53% and 52%, respectively) and HLA-DRB1 SE positivity (76% and 71%, respectively). The parameter of patients who had ever smoked was more common in the VARA registry (80%) than in SONORA (65%). GSTM1-null was significantly associated with ACPA positivity in the VARA registry (odds ratio (OR), 1.45; 95% confidence interval (CI), 1.02 to 2.05), but not in SONORA (OR, 1.00; 95% CI, 0.71 to 1.42). There were significant additive interactions between GSTM1 and HLA-DRB1 SE in the VARA registry (AP, 0.49; 95% CI, 0.21 to 0.77; P < 0.001) in ACPA positivity, an interaction replicated in SONORA (AP, 0.38; 95% CI, 0.00 to 0.76; P = 0.050)., Conclusions: This study is the first to show that the GSTM1-null genotype, a common genetic variant, exerts significant additive interaction with HLA-DRB1 SE on the risk of ACPA positivity in RA. Since GSTM1 has known antioxidant functions, these data suggest that oxidative stress may be important in the development of RA-specific autoimmunity in genetically susceptible individuals.

Introduction: Selective IgA deficiency (IgAD; serum IgA concentration of <0.07 g/l) is the most common primary immunodeficiency in Caucasians with an estimated prevalence of 1/600. The frequency of the extended major histocompatibility complex haplotype HLA A1, B8, DR3, DQ2 (the "8.1" haplotype) is increased among patients with IgAD., Materials and Methods: We carried out a direct measurement of the relative risk of homozygosity of the 8.1 haplotype for IgA deficiency in a population-based sample of 117 B8, DR3 homozygous individuals., Results and Discussion: IgA deficiency was found to be present in 2 of 117 (1.7%) of these subjects, a figure that is concordant with estimates of relative risk from large case-control studies in the Swedish population. These data are consistent with a multiplicative model for the 8.1 haplotype contribution to IgA deficiency and contrasts with prior studies, suggesting a much higher risk for 8.1 homozygosity. Using a dense single nucleotide polymorphism marker analysis of the MHC region in HLA B8, DR3, DQ2 homozygous individuals, we did not observe consistent differences between cases (n = 26) and controls (n = 24). Overall, our results do not support the hypothesis that IgA deficiency is associated with a distinct subgroup of 8.1 related haplotypes, but rather indicate that risk is conferred by the common 8.1 haplotype acting in multiplicative manner.

In addition to the HLA locus, over 30 genetic loci have been convincingly associated with risk for rheumatoid arthritis, and the majority of these associations have been identified in the last four years. Although this is a remarkable accomplishment, the majority of the genetic risk for RA still remains to be identified. Some of this "missing heritability" will likely be due to rare genetic variation, and will require extensive resequencing of the genomes of patients with RA. In addition, with few exceptions, the function and role in disease pathogenesis of the newly defined risk genes is unknown. Thus, the initial harvest of RA loci will catalyze new lines of hypothesis driven research to determine their role in disease pathogenesis. In addition, the rapidly advancing genetic technologies should lead to a more complete definition of the genetic underpinnings or RA in the next few years.

For Genetic Analysis Workshop 16 Problem 1, we provided data for genome-wide association analysis of rheumatoid arthritis. Single-nucleotide polymorphism (SNP) genotype data were provided for 868 cases and 1194 controls that had been assayed using an Illumina 550 k platform. In addition, phenotypic data were provided from genotyping DRB1 alleles, which were classified according to the rheumatoid arthritis shared epitope, levels of anti-cyclic citrullinated peptide, and levels of rheumatoid factor IgM. Several questions could be addressed using the data, including analysis of genetic associations using single SNPs or haplotypes, as well as gene-gene and genetic analysis of SNPs for qualitative and quantitative factors.

Objective: To describe a large, multicenter prospective cohort study of first-degree relatives (FDRs) of probands with rheumatoid arthritis (RA), and outline the use of such a study in investigating the natural history of RA development., Methods: A total of 1,058 FDRs, none of whom met the American College of Rheumatology criteria for RA, were enrolled in a prospective study investigating genetic and environmental influences on the development of RA-related autoimmunity. Demographic, epidemiologic, genetic, autoantibody, and physical examination data from the initial study enrollment visit were described for these FDRs, and the relationship was examined between genetic factors, autoantibodies, inflammation, and joint disease., Results: Fifty-five percent of the FDRs had > or =1 copy of the shared epitope, 20% had > or =1 copy of the PTPN22 polymorphism, and approximately 16% were positive for rheumatoid factor (RF; including isotypes) and/or anti-cyclic citrullinated peptide antibody. IgM-RF positivity is associated with > or =1 tender joint on examination (odds ratio [OR] 2.50, 95% confidence interval [95% CI] 1.27-4.89; P < 0.01) and elevated C-reactive protein (CRP) levels (OR 5.31, 95% CI 1.45-19.52; P = 0.01)., Conclusion: FDRs without RA demonstrate high prevalences of genetic risk factors and RA-related autoantibodies. Additionally, an RF association with tender joints and elevated CRP levels suggests that autoantibodies are a valid intermediate marker of RA-related autoimmunity in this cohort. This prospective FDR cohort will be a valuable resource for evaluating the relationship between genetic and epidemiologic factors and the development of RA-related autoimmunity.

Biomarker development for prediction of patient response to therapy is one of the goals of molecular profiling of human tissues. Due to the large number of transcripts, relatively limited number of samples, and high variability of data, identification of predictive biomarkers is a challenge for data analysis. Furthermore, many genes may be responsible for drug response differences, but often only a few are sufficient for accurate prediction. Here we present an analysis approach, the Convergent Random Forest (CRF) method, for the identification of highly predictive biomarkers. The aim is to select from genome-wide expression data a small number of non-redundant biomarkers that could be developed into a simple and robust diagnostic tool. Our method combines the Random Forest classifier and gene expression clustering to rank and select a small number of predictive genes. We evaluated the CRF approach by analyzing four different data sets. The first set contains transcript profiles of whole blood from rheumatoid arthritis patients, collected before anti-TNF treatment, and their subsequent response to the therapy. In this set, CRF identified 8 transcripts predicting response to therapy with 89% accuracy. We also applied the CRF to the analysis of three previously published expression data sets. For all sets, we have compared the CRF and recursive support vector machines (RSVM) approaches to feature selection and classification. In all cases the CRF selects much smaller number of features, five to eight genes, while achieving similar or better performance on both training and independent testing sets of data. For both methods performance estimates using cross-validation is similar to performance on independent samples. The method has been implemented in R and is available from the authors upon request: Jadwiga.Bienkowska@biogenidec.com.