Genetic variation was first shown to be part of the cause of systemic lupus erythematosus (SLE or lupus) in the 1970s with associations in the human leukocyte antigen (HLA) region. Almost four decades later, and with the help of increasingly powerful genetic approaches, more than 25 genes are now known to contribute to the mechanisms that predispose individuals to lupus. Over half of these loci have been discovered in the past two years, underscoring the extraordinary success of recent genome-wide association approaches in SLE. The now well established genetic risk factors include alleles in the MHC region (multiple genes), IRF5, ITGAM, STAT4, BLK, BANK1, PDCD1, PTPN22, TNFSF4, TNFAIP3, SPP1, ATG5, XKR6, PXK, some of the Fcγ receptors, and deficiencies in several complement components, including C1q, C4, and C2. As reviewed here, many of these genes fall into key pathways that are consistent with previous studies implicating immune complexes, host immune signal transduction, and interferon pathways in the pathogenesis of SLE. Other genetic loci have no known function or apparent immunological role and have the potential to reveal novel disease mechanisms. Certainly, as our understanding of the genetic etiology of SLE continues to mature, important new opportunities will emerge for developing more targeted and effective diagnostic and clinical management tools for this complex autoimmune disease.

Many autoimmune diseases (ADs) share similar underlying pathology and have a tendency to cluster within families, supporting the involvement of shared susceptibility genes. To date, most of the genetic variants associated with systemic lupus erythematosus (SLE) susceptibility also show association with others ADs. ITGAM and its associated ‘predisposing’ variant (rs1143679, Arg77His), predicted to alter the tertiary structures of the ligand-binding domain of ITGAM, may play a key role for SLE pathogenesis. The aim of this study is to examine whether the ITGAM variant is also associated with other ADs. We evaluated case-control association between rs1143679 and ADs (N=18,457) including primary Sjögren’s syndrome, systemic sclerosis, multiple sclerosis, rheumatoid arthritis, juvenile idiopathic arthritis, celiac disease, and type-1 diabetes. We also performed meta-analyses using our data in addition to available published data. Although the risk allele ‘A’ is relatively more frequent among cases for each disease, it was not significantly associated with any other ADs tested in this study. However, the meta-analysis for systemic sclerosis was associated with rs1143679 (pmeta=0.008). In summary, this study explored the role of ITGAM in general autoimmunity in seven non-lupus ADs, and only found association for systemic sclerosis when our results were combined with published results. Thus ITGAM may not be a general autoimmunity gene but this variant may be specifically associated with SLE and systemic sclerosis.

Recent advances in the field of genetics have dramatically changed our understanding of autoimmune disease. Candidate gene and, more recently, genome-wide association (GWA) studies have led to an explosion in the number of loci and pathways known to contribute to autoimmune phenotypes. Since the 1970s, researchers have known that several alleles in the MHC region play a role in the pathogenesis of many autoimmune diseases. More recent work has identified numerous risk loci involving both the innate and adaptive immune responses. However, much remains to be learned about the heritability of autoimmune conditions. Most regions found through GWA scans have yet to isolate the association to the causal allele(s) responsible for conferring disease risk. A role for rare variants (allele frequencies of <1%) has begun to emerge. Future research will use next generation sequencing (NGS) technology to comprehensively evaluate the human genome for risk variants. Whole transcriptome sequencing is now possible, which will provide much more detailed gene expression data. The dramatic drop in the cost and time required to sequence the entire human genome will ultimately make it possible for this technology to be used as a clinical diagnostic tool.

Purpose

While the key inciting events that drive the progression from autoantibodies to clinical disease remain to be clarified, new light has been shed on the factors contributing to disease susceptibility and the role of genetic factors in determining Sjogren's syndrome (SS) disease phenotypes. The purpose of this review is to provide an update on the role of genetic markers in the susceptibility to and pathogenesis of Sjogren's syndrome. This paper also discusses how genomic and proteomic technology can help in the design of specific therapeutics.

Key Findings

Recent evidence suggests that inflammatory genes associated with interferon pathways, and specific regulatory genes that control the maturation and proliferation of B cells, contribute to the pathogenesis of Sjogren's syndrome. Both gene expression profiling technology and gene association studies have been used to identify these key biologic pathways. Molecularly defined subsets of pSS patients are also being revealed by these studies. Previously identified gene loci which predispose to multiple autoimmune disorders have been confirmed supporting the paradigm of “general” autoimmune disease genes. Association of SS with many additional susceptibility loci are likely to be established through ongoing genome-wide association scans (GWAS). Clues from genetic studies suggest that targeting B cells will prove to be an effective way of reducing the systemic manifestations of pSS and are supported by early clinical trials.

Summary

Genome-wide technologies are likely to identify new genes and molecular pathways in the pathogenesis of SS that will be useful not only to identify patients at risk for SS, but also to identify subsets of patients at risk for variable levels of disease severity. In the future, these studies could identify novel biomarkers that will lead to significant advances in management by providing the means to tailor therapeutic strategies to individual patients.

Sjögren’s syndrome (SS) is a common chronic autoimmune disease characterized by lymphocytic infiltration of exocrine glands. Affected cases commonly present with oral and ocular dryness, thought to be the result of inflammatory cell-mediated gland dysfunction. To identify important molecular pathways involved in SS, we used high-density microarrays to define global gene expression profiles in peripheral blood. We first analyzed 21 SS cases and 23 controls and identified a prominent pattern of overexpressed genes that are inducible by interferons (IFNs). These results were confirmed by evaluation of a second independent dataset of 17 SS cases and 22 controls. Additional inflammatory and immune-related pathways with altered expression patterns in SS cases included B and T cell receptor, IGF-1, GM-CSF, PPARα/RXRα, and PI3/AKT signaling. Exploration of these data for relationships to clinical features of disease revealed that expression levels for most IFN-inducible genes were positively correlated with titers of anti-Ro/SSA (P<0.001) and anti-La/SSB (P<0.001) autoantibodies. Diagnostic and therapeutic approaches targeting IFN signaling pathway may prove most effective in the subset of SS cases who produce anti-Ro/SSA and anti-La/SSB autoantibodies. Our results strongly support innate and adaptive immune processes in the pathogenesis of SS and provide numerous candidate disease markers for further study.

Sjoögren's syndrome (SS) is a chronic, progressive exocrinopathy characterized by infiltration and proliferation of lymphocytes into affected glands. Although patients are clinically identified through oral and ocular features, the full spectrum of disease encompasses a complex and myriad systemic symptoms. The primary pathophysiology includes concurrent mechanisms of dysregulated innate immunity and adaptive autoimmunity involving cell-mediated and humoral disease processes. Etiology involves environmental and genetic factors; however, large-scale genetic studies have not yet been conducted in SS and the genetic basis for SS is largely unexplored.

Although few genetic studies have been completed to date in SS, the overall evidence to support a genetic basis for SS continues to grow. Current data strongly suggest SS is a complex, polygenic disorder likely sharing common genetic determinants with related autoimmune diseases, such as systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA). Recent advances in SLE and RA provide valuable insight into the potential genetic complexity of SS. This article reviews association studies in various candidate genes for SS completed to date and highlights insights from SS mouse models. Advanced genetic and genomic technologies now are available for assaying gene expression and genetic associations across the entire genome, providing important opportunities to conduct unbiased interrogation of essentially every gene for a role in SS.

The ATP-binding cassette transporter (TAP) proteins are functionally relevant candidates for predisposition to Systemic Lupus Erythematosus (SLE) by virtue of their role in autoantigen presentation and location in the MHC. We tested if variation in the TAP genes (TAP1 and TAP2) is associated with SLE. We genotyped tag single nucleotide polymorphisms (SNPs) and performed family-based association analysis on 390 Caucasian pedigrees. We found significant evidence of association between TAP2 and SLE (rs241453, P = 1.33 × 10-6). Conditional logistic regression analysis suggests that this TAP2 effect is separate from the HLA-DRB1 alleles. Our analyses show that both rs241453 (P = 1.6 × 10-4) and HLA-DRB1*03xx (P = 2.3 × 10-4) have significant autonomous effects not due to linkage disequilibrium. Moreover, these loci exhibit a significant statistical interaction (P < 1.0 × 10-6), demonstrated by an increase in the odds ratio for the TAP2 association from OR = 2.00 (CI=1.17-3.42) in HLA-DRB1*03xx-negative subjects to OR = 4.29 (CI=1.88-9.76) in the subjects with at least one HLA-DRB1*03xx allele group. We report the largest association study of the TAP genes with SLE to date, and the first to test for its separate effect and interaction with the HLA alleles consistently associated with SLE.

Complement cascade plasma proteins have a complex role in the etiopathogenesis of SLE. Hereditary C1q deficiency has been strongly related to SLE; however, there are very few published SLE studies that evaluate the polymorphisms of the genes encoding for C1q (A, B, and C). In this study, we evaluated 17 single nucleotide polymorphisms (SNPs) across 37 kb of C1QA, B and C in a lupus cohort of peoples of African-American and Hispanic origin. In a case only analysis, significant association at multiple SNPs in the C1QA gene was detected in African-Americans with kidney nephritis (best p=4.91 × 10−6). In addition, C1QA was associated with SLE in African-Americans with a lack of nephritis and accompanying photosensitivity when compared to normal controls (p=6.80 × 10−6). A similar trend was observed in the Hispanic subjects (p=0.003). Quantitative analysis demonstrates that some SNPs in the C1q genes might be correlated with C3 complement levels in an additive model among African-Americans (best p=0.0001). The CIQA gene is associated with subphenotypes of lupus in African-American and Hispanic subjects. Further studies with higher SNP densities in this region and other complement components are necessary to elucidate the complex genetics and phenotypic interactions between complement components and SLE.

We targeted LYN, a src-tyosine kinase involved in B cell activation, in case-control association studies using populations of European American, African American and Korean subjects. Our combined European-derived population, consisting of 2463 independent cases and 3131 unrelated controls, demonstrates significant association with rs6983130 in a female-only analysis with 2254 cases and 2228 controls (p=1.1 × 10−4, OR=0.81 (95% CI: 0.73 – 0.90)). This SNP is located in the 5′ UTR within the first intron near the transcription initiation site of LYN. Additional SNPs upstream of the first exon also show weak and sporadic association in subsets of the total European American population. Multivariate logistic regression analysis implicates rs6983130 as a protective factor for SLE susceptibility when anti-dsDNA, anti-chromatin, anti-52 kDa Ro or anti-Sm autoantibody status were used as covariates. Subset analysis of the European American female cases by ACR classification criteria reveals a reduction in the risk of hematologic disorder with rs6983130 compared to cases without hematologic disorders (p=1.5 × 10−3, OR=0.75 (95% C.I.=0.62-0.89)). None of the 90 SNPs tested demonstrate significant association with SLE in the African American or Korean populations. These results support an association of LYN with European-derived individuals with SLE, especially within autoantibody or clinical subsets.

Systemic lupus erythematosus (SLE) is an autoimmune disease with highly variable clinical presentation. Patients suffer from immunological abnormalities that target T cell, B cell and accessory cell functions. B cells are hyperactive in SLE patients. An adaptor protein expressed in B cells called BANK1 (B-cell scaffold protein with ankyrin repeats) was reported in a previous study to be associated with SLE in a European population. The objective of this study is to assess the BANK1 genotype-phenotype association in an independent replication sample. We genotyped 38 single nucleotide polymorphisms (SNPs) in BANK1 on 1892 European-derived SLE patients and 2652 European-derived controls. The strongest associations with SLE and BANK1 were at rs17266594 (corrected p-value=1.97 × 10−5, OR=1.22, 95% C.I.(1.12–1.34)) and rs10516487 (corrected p-value=2.59 × 10−5, OR=1.22, 95% C.I.(1.11–1.34)). Our findings suggest that the association is explained by these two SNPs, confirming previous reports that these polymorphisms contribute to the risk of developing lupus. Analysis of patient subsets enriched for hematological, immunological and renal ACR criteria or the levels of autoantibodies, such as anti-RNP A and anti-SmRNP, uncovers additional BANK1 associations. Our results suggest that BANK1 polymorphisms alter immune system development and function to increase the risk for developing lupus.

TNFAIP3 encodes the ubiquitin modifying enzyme, A20, a key regulator of inflammatory signaling pathways. We previously reported association between TNFAIP3 variants and systemic lupus erythematosus (SLE). In order to further localize the risk variant(s), we performed a meta-analysis using genetic data available from two Caucasian case/control datasets (1453 total cases, 3381 total controls) and 713 SLE trio families. The best result was found at rs5029939 (P = 1.67 × 10−14, OR = 2.09, 95% CI 1.68–2.60). We then imputed SNPs from the CEU Phase II HapMap using genotypes from 431 SLE cases and 2155 controls. Imputation identified eleven SNPs in addition to three observed SNPs, which together, defined a 109 kb SLE risk segment surrounding TNFAIP3. When evaluating whether the rs5029939 risk allele was associated with SLE clinical manifestations, we observed that heterozygous carriers of the TNFAIP3 risk allele at rs5029939 have a two-fold increased risk of developing renal or hematologic manifestations compared to homozygous non-risk subjects. In summary, our study strengthens the genetic evidence that variants in the region of TNFAIP3 influence risk for SLE, particularly in patients with renal and hematologic manifestations, and narrows the risk effect to a 109 kb DNA segment that spans the TNFAIP3 gene.

Background

Current statistical methods for sib-pair linkage analysis of complex diseases include linear models, generalized linear models, and novel data mining techniques. The purpose of this study was to further investigate the utility and properties of a novel pattern recognition technique (step-wise discriminant analysis) using the chromosome 10 linkage data from the Framingham Heart Study and by comparing it with step-wise logistic regression and linear regression.

Results

The three step-wise approaches were compared in terms of statistical significance and gene localization. Step-wise discriminant linkage analysis approach performed best; next was step-wise logistic regression; and step-wise linear regression was the least efficient because it ignored the categorical nature of disease phenotypes. Nevertheless, all three methods successfully identified the previously reported chromosomal region linked to human hypertension, marker GATA64A09. We also explored the possibility of using the discriminant analysis to detect gene × gene and gene × environment interactions. There was evidence to suggest the existence of gene × environment interactions between markers GATA64A09 or GATA115E01 and hypertension treatment and gene × gene interactions between markers GATA64A09 and GATA115E01. Finally, we answered the theoretical question "Is a trichotomous phenotype more efficient than a binary?" Unlike logistic regression, discriminant sib-pair linkage analysis might have more power to detect linkage to a binary phenotype than a trichotomous one.

Conclusion

We confirmed our previous speculation that step-wise discriminant analysis is useful for genetic mapping of complex diseases. This analysis also supported the possibility of the pattern recognition technique for investigating gene × gene or gene × environment interactions.

Background

Longitudinal data often have multiple (repeated) measures recorded along a time trajectory. For example, the two cohorts from the Framingham Heart Study (GAW13 Problem 1) contain 21 and 5 repeated measures for hypertension phenotypes as well as epidemiological risk factors, respectively. Direct modelling of a large number of serially and biologically correlated traits in the context of linkage analysis can be prohibitively complex. Alternatively, we may consider using univariate transformation for linkage analysis of longitudinal repeated measures.

Results

We evaluated the utility of three conventional summary measures (mean, slope, and principal components) for genetic linkage analysis of longitudinal phenotypes by analyzing the chromosome 10 data of the Framingham Heart Study. Except for the temporal slope, all of the summary methods and the multivariate analysis identified the previously reported region, marker GATA64A09, for systolic blood pressure or high blood pressure. Further analysis revealed that this region may harbor gene(s) affecting human blood pressure at multiple stages of life.

Conclusion

We conclude that mean and principal components are feasible alternatives for genetic linkage analysis of longitudinal phenotypes, but the slope might have a separate genetic basis from that of the original longitudinal phenotypes.

Objective

High serum interferon α (IFNα) activity is a heritable risk factor for systemic lupus erythematosus (SLE). Auto-antibodies found in SLE form immune complexes which can stimulate IFNα production by activating endosomal Toll-like receptors and interferon regulatory factors (IRFs), including IRF5. Genetic variation in IRF5 is associated with SLE susceptibility; however, it is unclear how IRF5 functional genetic elements contribute to human disease.

Methods

1034 patients with SLE and 989 controls of European ancestry, 555 patients with SLE and 679 controls of African–American ancestry, and 73 patients with SLE of South African ancestry were genotyped at IRF5 polymorphisms, which define major haplotypes. Serum IFNα activity was measured using a functional assay.

Results

In European ancestry subjects, anti-double-stranded DNA (dsDNA) and anti-Ro antibodies were each associated with different haplotypes characterised by a different combination of functional genetic elements (OR > 2.56, p >003C; 1.9×10−14 for both). These IRF5 haplotype-auto-antibody associations strongly predicted higher serum IFNα in patients with SLE and explained > 70% of the genetic risk of SLE due to IRF5. In African–American patients with SLE a similar relationship between serology and IFNα was observed, although the previously described European ancestry-risk haplotype was present at admixture proportions in African–American subjects and absent in African patients with SLE.

Conclusions

The authors define a novel risk haplotype of IRF5 that is associated with anti-dsDNA antibodies and show that risk of SLE due to IRF5 genotype is largely dependent upon particular auto-antibodies. This suggests that auto-antibodies are directly pathogenic in human SLE, resulting in increased IFNα in cooperation with particular combinations of IRF5 functional genetic elements.

SLE is a systemic autoimmune disorder affecting multiple organ systems including the skin, musculoskeletal, renal and haematopoietic systems. Humoral autoimmunity is a hallmark of SLE, and patients frequently have circulating auto-antibodies directed against dsDNA, as well as RNA binding proteins (RBP). Anti-RBP autoantibodies include antibodies which recognize Ro, La, Smith (anti-Sm), and ribonucleoprotein (anti-nRNP), collectively referred to as anti-retinol-binding protein). Anti-retinol-binding protein and anti-dsDNA auto-antibodies are rare in the healthy population.1 These auto-antibodies can be present in sera for years preceding the onset of clinical SLE illness2 and are likely pathogenic in SLE.34

Systemic lupus erythematosus (SLE) is an inflammatory autoimmune disease with a strong genetic component. African-Americans (AA) are at increased risk of SLE, but the genetic basis of this risk is largely unknown. To identify causal variants in SLE loci in AA, we performed admixture mapping followed by fine mapping in AA and European-Americans (EA). Through genome-wide admixture mapping in AA, we identified a strong SLE susceptibility locus at 2q22–24 (LOD = 6.28), and the admixture signal is associated with the European ancestry (ancestry risk ratio ∼1.5). Large-scale genotypic analysis on 19,726 individuals of African and European ancestry revealed three independently associated variants in the IFIH1 gene: an intronic variant, rs13023380 [Pmeta = 5.20×10−14; odds ratio, 95% confidence interval = 0.82 (0.78–0.87)], and two missense variants, rs1990760 (Ala946Thr) [Pmeta = 3.08×10−7; 0.88 (0.84–0.93)] and rs10930046 (Arg460His) [Pdom = 1.16×10−8; 0.70 (0.62–0.79)]. Both missense variants produced dramatic phenotypic changes in apoptosis and inflammation-related gene expression. We experimentally validated function of the intronic SNP by DNA electrophoresis, protein identification, and in vitro protein binding assays. DNA carrying the intronic risk allele rs13023380 showed reduced binding efficiency to a cellular protein complex including nucleolin and lupus autoantigen Ku70/80, and showed reduced transcriptional activity in vivo. Thus, in SLE patients, genetic susceptibility could create a biochemical imbalance that dysregulates nucleolin, Ku70/80, or other nucleic acid regulatory proteins. This could promote antibody hypermutation and auto-antibody generation, further destabilizing the cellular network. Together with molecular modeling, our results establish a distinct role for IFIH1 in apoptosis, inflammation, and autoantibody production, and explain the molecular basis of these three risk alleles for SLE pathogenesis.

Author Summary

African-Americans (AA) are at increased risk of systemic lupus erythematosus (SLE), but the genetic basis of this risk increase is largely unknown. We used admixture mapping to localize disease-causing genetic variants that differ in frequency across populations. This approach is advantageous for localizing susceptibility genes in recently admixed populations like AA. Our genome-wide admixture scan identified seven admixture signals, and we followed the best signal at 2q22–24 with fine-mapping, imputation-based association analysis and experimental validation. We identified two independent coding variants and a non-coding variant within the IFIH1 gene associated with SLE. Together with molecular modeling, our results establish a distinct role for IFIH1 in apoptosis, inflammation, and autoantibody production, and explain the molecular basis of these three risk alleles for SLE pathogenesis.

Objective

Several confirmed genetic susceptibility loci for lupus have been described. To date, no clear evidence for genetic epistasis is established in lupus. We test for gene-gene interactions in a number of known lupus susceptibility loci.

Methods

Eighteen SNPs tagging independent and confirmed lupus susceptibility loci were genotyped in a set of 4,248 lupus patients and 3,818 normal healthy controls of European descent. Epistasis was tested using a 2-step approach utilizing both parametric and non-parametric methods. The false discovery rate (FDR) method was used to correct for multiple testing.

Results

We detected and confirmed gene-gene interactions between the HLA region and CTLA4, IRF5, and ITGAM, and between PDCD1 and IL21 in lupus patients. The most significant interaction detected by parametric analysis was between rs3131379 in the HLA region and rs231775 in CTLA4 (Interaction odds ratio=1.19, z-score= 3.95, P= 7.8×10−5 (FDR≤0.05), PMDR= 5.9×10−45). Importantly, our data suggest that in lupus patients the presence of the HLA lupus-risk alleles in rs1270942 and rs3131379 increases the odds of also carrying the lupus-risk allele in IRF5 (rs2070197) by 17% and 16%, respectively (P= 0.0028 and 0.0047).

Conclusion

We provide evidence for gene-gene epistasis in systemic lupus erythematosus. These findings support a role for genetic interaction contributing to the complexity of lupus heritability.

Objective

Systemic lupus erythematosus is a clinically heterogeneous autoimmune disease. A number of genetic loci that increase lupus susceptibility have been established. This study examines if these genetic loci also contribute to the clinical heterogeneity in lupus.

Materials and methods

4001 European-derived, 1547 Hispanic, 1590 African-American and 1191 Asian lupus patients were genotyped for 16 confirmed lupus susceptibility loci. Ancestry informative markers were genotyped to calculate and adjust for admixture. The association between the risk allele in each locus was determined and compared in patients with and without the various clinical manifestations included in the ACR criteria.

Results

Renal disorder was significantly correlated with the lupus risk allele in ITGAM (p=5.0×10−6, OR 1.25, 95% CI 1.12 to 1.35) and in TNFSF4 (p=0.0013, OR 1.14, 95% CI 1.07 to 1.25). Other significant findings include the association between risk alleles in FCGR2A and malar rash (p=0.0031, OR 1.11, 95% CI 1.17 to 1.33), ITGAM and discoid rash (p=0.0020, OR 1.20, 95% CI 1.06 to 1.33), STAT4 and protection from oral ulcers (p=0.0027, OR 0.89, 95% CI 0.83 to 0.96) and IL21 and haematological disorder (p=0.0027, OR 1.13, 95% CI 1.04 to 1.22). All these associations are significant with a false discovery rate of <0.05 and pass the significance threshold using Bonferroni correction for multiple testing.

Conclusion

Significant associations were found between lupus clinical manifestations and the FCGR2A, ITGAM, STAT4, TNSF4 and IL21 genes. The findings suggest that genetic profiling might be a useful tool to predict disease manifestations in lupus patients in the future.

Sarcoidosis is a systemic inflammatory disease characterized by the formation of granulomas in affected organs. Genome-wide association studies (GWASs) of this disease have been conducted only in European population. We present the first sarcoidosis GWAS in African Americans (AAs, 818 cases and 1,088 related controls) followed by replication in independent sets of AAs (455 cases and 557 controls) and European Americans (EAs, 442 cases and 2,284 controls). We evaluated >6 million SNPs either genotyped using the Illumina Omni1-Quad array or imputed from the 1000 Genomes Project data. We identified a novel sarcoidosis-associated locus, NOTCH4, that reached genome-wide significance in the combined AA samples (rs715299, PAA-meta = 6.51×10−10) and demonstrated the independence of this locus from others in the MHC region in the same sample. We replicated previous European GWAS associations within HLA-DRA, HLA-DRB5, HLA-DRB1, BTNL2, and ANXA11 in both our AA and EA datasets. We also confirmed significant associations to the previously reported HLA-C and HLA-B regions in the EA but not AA samples. We further identified suggestive associations with several other genes previously reported in lung or inflammatory diseases.

The major histocompatibility complex (MHC) class II transactivator gene (CIITA) encodes an important transcription factor required for HLA class II MHC-restricted antigen presentation. MHC genes, including the HLA class II DRB1*03:01 allele, are strongly associated with systemic lupus erythematosus (SLE). Recently the rs4774 CIITA missense variant (+1632G/C) was reported to be associated with susceptibility to multiple sclerosis. In the current study, we investigated CIITA, DRB1*03:01 and risk of SLE using a multi-stage analysis. In stage 1, 9 CIITA variants were tested in 658 cases and 1,363 controls (N = 2,021). In stage 2, rs4774 was tested in 684 cases and 2,938 controls (N = 3,622). We also performed a meta-analysis of the pooled 1,342 cases and 4,301 controls (N = 5,643). In stage 1, rs4774*C was associated with SLE (odds ratio [OR] = 1.24, 95% confidence interval [95% CI] = 1.07–1.44, P = 4.2 × 10−3). Similar results were observed in stage 2 (OR = 1.16, 95% CI = 1.02–1.33, P = 8.5×10−3) and the meta-analysis of the combined dataset (OR = 1.20, 95% CI = 1.09–1.33, Pmeta = 2.5×10−4). In all three analyses, the strongest evidence for association between rs4774*C and SLE was present in individuals who carried at least one copy of DRB1*03:01 (Pmeta= 1.9×10−3). Results support a role for CIITA in SLE, which appears to be stronger in the presence of DRB1*03:01.

Objective

The overexpression of interferon (IFN)-inducible genes is a prominent feature of SLE, serves as a marker for active and more severe disease, and is also observed in other autoimmune and inflammatory conditions. The genetic variations responsible for sustained activation of IFN responsive genes are unknown.

Methods

We systematically evaluated association of SLE with a total of 1,754 IFN-pathway related genes, including IFN-inducible genes known to be differentially expressed in SLE patients and their direct regulators. We performed a three-stage design where two cohorts (total n=939 SLE cases, 3,398 controls) were analyzed independently and jointly for association with SLE, and the results were adjusted for the number of comparisons.

Results

A total of 16,137 SNPs passed all quality control filters of which 316 demonstrated replicated association with SLE in both cohorts. Nine variants were further genotyped for confirmation in an average of 1,316 independent SLE cases and 3,215 independent controls. Association with SLE was confirmed for several genes, including the transmembrane receptor CD44 (rs507230, P = 3.98×10−12), cytokine pleiotrophin (PTN) (rs919581, P = 5.38×10−04), the heat-shock DNAJA1 (rs10971259, P = 6.31×10−03), and the nuclear import protein karyopherin alpha 1 (KPNA1) (rs6810306, P = 4.91×10−02).

Conclusion

This study expands the number of candidate genes associated with SLE and highlights the potential of pathway-based approaches for gene discovery. Identification of the causal alleles will help elucidate the molecular mechanisms responsible for activation of the IFN system in SLE.

Objective

Genetic association of the IL2/IL21 region at 4q27 has been previously reported in lupus and a number of autoimmune and inflammatory diseases. Herein, using a very large cohort of lupus patients and controls, we localize this genetic effect to the IL21 gene.

Methods

We genotyped 45 tag SNPs across the IL2/IL21 locus in two large independent lupus sample sets. We studied a European-derived set consisting of 4,248 lupus patients and 3,818 healthy controls, and an African-American set of 1,569 patients and 1,893 healthy controls. Imputation in 3,004 WTCCC additional control individuals was also performed. Genetic association between the genotyped markers was determined, and pair-wise conditional analysis was performed to localize the independent genetic effect in the IL2/IL21 locus in lupus.

Results

We established and confirmed the genetic association between IL2/IL21 and lupus. Using conditional analysis and trans-ethnic mapping, we localized the genetic effect in this locus to two SNPs in high linkage disequilibrium; rs907715 located within IL21 (OR=1.16 (1.10–1.22), P= 2.17 ×10−8), and rs6835457 located in the 3’-UTR flanking region of IL21 (OR= 1.11 (1.05–1.17), P= 9.35×10−5).

Conclusion

We have established the genetic association between lupus and IL2/IL21 with a genome-wide level of significance. Further, we localized this genetic association within the IL2/IL21 linkage disequilibrium block to IL21. If other autoimmune IL2/IL21 genetic associations are similarly localized, then the IL21 risk alleles would be predicted to operate in a fundamental mechanism that influences the course of a number of autoimmune disease processes.

In spite of the well-known clustering of multiple autoimmune disorders in families, analyses of specific shared genes and polymorphisms between systemic lupus erythematosus (SLE) and other autoimmune diseases (ADs) have been limited. Therefore, we comprehensively tested autoimmune variants for association with SLE, aiming to identify pleiotropic genetic associations between these diseases. We compiled a list of 446 non–Major Histocompatibility Complex (MHC) variants identified in genome-wide association studies (GWAS) of populations of European ancestry across 17 ADs. We then tested these variants in our combined Caucasian SLE cohorts of 1,500 cases and 5,706 controls. We tested a subset of these polymorphisms in an independent Caucasian replication cohort of 2,085 SLE cases and 2,854 controls, allowing the computation of a meta-analysis between all cohorts. We have uncovered novel shared SLE loci that passed multiple comparisons adjustment, including the VTCN1 (rs12046117, P = 2.02×10−06) region. We observed that the loci shared among the most ADs include IL23R, OLIG3/TNFAIP3, and IL2RA. Given the lack of a universal autoimmune risk locus outside of the MHC and variable specificities for different diseases, our data suggests partial pleiotropy among ADs. Hierarchical clustering of ADs suggested that the most genetically related ADs appear to be type 1 diabetes with rheumatoid arthritis and Crohn's disease with ulcerative colitis. These findings support a relatively distinct genetic susceptibility for SLE. For many of the shared GWAS autoimmune loci, we found no evidence for association with SLE, including IL23R. Also, several established SLE loci are apparently not associated with other ADs, including the ITGAM-ITGAX and TNFSF4 regions. This study represents the most comprehensive evaluation of shared autoimmune loci to date, supports a relatively distinct non–MHC genetic susceptibility for SLE, provides further evidence for previously and newly identified shared genes in SLE, and highlights the value of studies of potentially pleiotropic genes in autoimmune diseases.

Author Summary

It is well known that multiple autoimmune disorders cluster in families. However, all of the genetic variants that explain this clustering have not been discovered, and the specific genetic variants shared between systemic lupus erythematosus (SLE) and other autoimmune diseases (ADs) are not known. In order to better understand the genetic factors that explain this predisposition to autoimmunity, we performed a comprehensive evaluation of shared autoimmune genetic variants. First we considered results from 17 ADs and compiled a list with 446 significant genetic variants from these studies. We identified some genetic variants extensively shared between ADs, as well as the ADs that share the most variants. The genetic overlap between SLE and other ADs was modest. Next we tested how important all the 446 genetic variants were in our collection with a minimum of 1,500 SLE patients. Among the most significant variants in SLE, the majority had already been identified in previous studies, but we also discovered variants in two important immune genes. In summary, our data identified diseases with common genetic risk factors and novel SLE effects, and this supports a relatively distinct genetic susceptibility for SLE. This study helps delineate the genetic architecture of ADs.

Systemic Lupus Erythematosus (SLE) is a prototypic autoimmune disorder with a complex pathogenesis in which genetic, hormonal and environmental factors play a role. Rare mutations in the TREX1 gene, the major mammalian 3′-5′ exonuclease, have been reported in sporadic SLE cases. Some of these mutations have also been identified in a rare pediatric neurologic condition featuring an inflammatory encephalopathy known as Aicardi-Goutières syndrome (AGS). We sought to investigate the frequency of these mutations in a large multi-ancestral cohort of SLE cases and controls.

Methods

Forty single-nucleotide polymorphisms (SNPs), including both common and rare variants, across the TREX1 gene were evaluated in ∼8370 patients with SLE and ∼7490 control subjects. Stringent quality control procedures were applied and principal components and admixture proportions were calculated to identify outliers for removal from analysis. Population-based case-control association analyses were performed. P values, false discovery rate q values, and odds ratios with 95% confidence intervals were calculated.

Results

The estimated frequency of TREX1 mutations in our lupus cohort was 0.5%. Five heterozygous mutations were detected at the Y305C polymorphism in European lupus cases but none were observed in European controls. Five African cases incurred heterozygous mutations at the E266G polymorphism and, again, none were observed in the African controls. A rare homozygous R114H mutation was identified in one Asian SLE patient whereas all genotypes at this mutation in previous reports for SLE were heterozygous. Analysis of common TREX1 SNPs (MAF >10%) revealed a relatively common risk haplotype in European SLE patients with neurologic manifestations, especially seizures, with a frequency of 58% in lupus cases compared to 45% in normal controls (p=0.0008, OR=1.73, 95% CI=1.25-2.39). Finally, the presence or absence of specific autoantibodies in certain populations produced significant genetic associations. For example, a strong association with anti-nRNP was observed in the European cohort at a coding synonymous variant rs56203834 (p=2.99E-13, OR=5.2, 95% CI=3.18-8.56).

Conclusion

Our data confirm and expand previous reports and provide additional support for the involvement of TREX1 in lupus pathogenesis.

Systemic lupus erythematosus (SLE), a complex polygenic autoimmune disease, is associated with increased complement activation. Variants of genes encoding complement regulator factor H (CFH) and five CFH-related proteins (CFHR1-CFHR5) within the chromosome 1q32 locus linked to SLE, have been associated with multiple human diseases and may contribute to dysregulated complement activation predisposing to SLE. We assessed 60 SNPs covering the CFH-CFHRs region for association with SLE in 15,864 case-control subjects derived from four ethnic groups. Significant allelic associations with SLE were detected in European Americans (EA) and African Americans (AA), which could be attributed to an intronic CFH SNP (rs6677604, in intron 11, Pmeta = 6.6×10−8, OR = 1.18) and an intergenic SNP between CFHR1 and CFHR4 (rs16840639, Pmeta = 2.9×10−7, OR = 1.17) rather than to previously identified disease-associated CFH exonic SNPs, including I62V, Y402H, A474A, and D936E. In addition, allelic association of rs6677604 with SLE was subsequently confirmed in Asians (AS). Haplotype analysis revealed that the underlying causal variant, tagged by rs6677604 and rs16840639, was localized to a ∼146 kb block extending from intron 9 of CFH to downstream of CFHR1. Within this block, the deletion of CFHR3 and CFHR1 (CFHR3-1Δ), a likely causal variant measured using multiplex ligation-dependent probe amplification, was tagged by rs6677604 in EA and AS and rs16840639 in AA, respectively. Deduced from genotypic associations of tag SNPs in EA, AA, and AS, homozygous deletion of CFHR3-1Δ (Pmeta = 3.2×10−7, OR = 1.47) conferred a higher risk of SLE than heterozygous deletion (Pmeta = 3.5×10−4, OR = 1.14). These results suggested that the CFHR3-1Δ deletion within the SLE-associated block, but not the previously described exonic SNPs of CFH, might contribute to the development of SLE in EA, AA, and AS, providing new insights into the role of complement regulators in the pathogenesis of SLE.

Author Summary

Systemic lupus erythematosus (SLE) is a complex autoimmune disease, associated with increased complement activation. Previous studies have provided evidence for the presence of SLE susceptibility gene(s) in the chromosome 1q31-32 locus. Within 1q32, genes encoding complement regulator factor H (CFH) and five CFH-related proteins (CFHR1-CFHR5) may contribute to the development of SLE, because genetic variants of these genes impair complement regulation and predispose to various human diseases. In this study, we tested association of genetic variants in the region containing CFH and CFHRs with SLE. We identified genetic variants predisposing to SLE in European American, African American, and Asian populations, which might be attributed to the deletion of CFHR3 and CFHR1 genes but not previously identified disease-associated exonic variants of CFH. This study provides the first evidence for consistent association between CFH/CFHRs and SLE across multi-ancestral SLE datasets, providing new insights into the role of complement regulators in the pathogenesis of SLE.