Faculty Bibliography

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 <or= 5 x 10(-7)). Here we show an association with genomic regions containing several genes controlling the activation and proliferation of regulatory T cells (T(reg) cells), cytotoxic T lymphocyte-associated antigen 4 (CTLA4), interleukin (IL)-2/IL-21, IL-2 receptor A (IL-2RA; CD25) and Eos (also known as Ikaros family zinc finger 4; IKZF4), as well as the human leukocyte antigen (HLA) region. We also find association evidence for regions containing genes expressed in the hair follicle itself (PRDX5 and STX17). A region of strong association resides within the ULBP (cytomegalovirus UL16-binding protein) gene cluster on chromosome 6q25.1, encoding activating ligands of the natural killer cell receptor NKG2D that have not previously been implicated in an autoimmune disease. By probing the role of ULBP3 in disease pathogenesis, we also show that its expression in lesional scalp from patients with AA is markedly upregulated in the hair follicle dermal sheath during active disease. This study provides evidence for the involvement of both innate and acquired immunity in the pathogenesis of AA. We have defined the genetic underpinnings of AA, placing it within the context of shared pathways among autoimmune diseases, and implicating a novel disease mechanism, the upregulation of ULBP ligands, in triggering autoimmunity

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: 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

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

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