Objective

Systemic lupus erythematosus (SLE; OMIM 152700) is a chronic autoimmune disease for which the aetiology includes genetic and environmental factors. ITGAM, integrin αΜ (complement component 3 receptor 3 subunit) encoding a ligand for intracellular adhesion molecule (ICAM) proteins, is an established SLE susceptibility locus. This study aimed to evaluate the independent and joint effects of genetic variations in the genes that encode ITGAM and ICAM.

Methods

The authors examined several markers in the ICAM1–ICAM4–ICAM5 locus on chromosome 19p13 and the single ITGAM polymorphism (rs1143679) using a large-scale case–control study of 17 481 unrelated participants from four ancestry populations. The single marker association and gene–gene interaction were analysed for each ancestry, and a meta-analysis across the four ancestries was performed.

Results

The A-allele of ICAM1–ICAM4–ICAM5 rs3093030, associated with elevated plasma levels of soluble ICAM1, and the A-allele of ITGAM rs1143679 showed the strongest association with increased SLE susceptibility in each of the ancestry populations and the trans-ancestry meta-analysis (ORmeta=1.16, 95% CI 1.11 to 1.22; p=4.88×10−10 and ORmeta=1.67, 95% CI 1.55 to 1.79; p=3.32×10−46, respectively). The effect of the ICAM single-nucleotide polymorphisms (SNPs) was independent of the effect of the ITGAM SNP rs1143679, and carriers of both ICAM rs3093030-AA and ITGAM rs1143679-AA had an OR of 4.08 compared with those with no risk allele in either SNP (95% CI 2.09 to 7.98; p=3.91×10−5).

Conclusion

These findings are the first to suggest that an ICAM–integrin-mediated pathway contributes to susceptibility to SLE.

IFNα is a potent activator of innate and adaptive immunity, and its administration to pre-autoimmune (NZBxNZW)F1 mice promotes virulent systemic lupus erythematosus (SLE) disease. Given the known contributions of B cells and BAFF to SLE, we evaluated the ability of IFNα administration to induce disease in wild-type (WT), B cell-deficient, and BAFF-deficient NZM 2328 mice. Whereas WT mice rapidly developed proliferative glomerulonephritis (GN), marked proteinuria, and increased mortality in response to IFNα administration, B cell-deficient mice developed neither renal pathology nor clinical disease. Moreover, BAFF-deficient mice, despite developing limited glomerular IgG and C3 deposition, also remained free of histological GN and clinical disease. Strikingly, similar T cell expansion and serum IgG responses were observed in Adv-IFN-treated WT and BAFF-deficient mice despite their disparate pathological and clinical responses, whereas numbers of activated B cells increased in WT mice but not in BAFF-deficient mice. Nonetheless, B cell, plasma cell, and T cell infiltration of the kidneys in Adv-IFN-treated WT mice was similar to that in WT mice treated with Adv-control. Its ability to promote SLE disease in WT mice notwithstanding, IFNα administration failed to drive the preferential expansion of CD4+ memory T cells that occurs during the natural course of disease, and glomerular infiltration of macrophages failed to associate with development of disease. These results collectively suggest that therapeutic targeting in SLE of BAFF and/or B cells in SLE could be successful even in states of IFNα overexpression. Moreover, our results document important biological differences between IFNα-driven and spontaneous “natural” SLE disease.

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) 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 pheterogeneity<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.

Author Summary

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease that can involve virtually any organ system. SLE patients produce antibodies that bind to their own cells and proteins (autoantibodies) which can cause irreversible organ damage. One particular SLE–related autoantibody directed at double-stranded DNA (anti–dsDNA) is associated with kidney involvement and more severe disease. Previous genome-wide association studies (GWAS) in SLE have studied SLE itself, not particular SLE manifestations. Therefore, we conducted this GWAS of anti–dsDNA autoantibody production to identify genetic associations with this clinically important autoantibody. We found that many previously identified SLE–associated genes are more strongly associated with anti–dsDNA autoantibody production than SLE itself, and they may be more accurately described as autoantibody propensity genes. No strong genetic associations were observed for SLE patients who do not produce anti–dsDNA autoantibodies, suggesting that other factors may have more influence in developing this type of SLE. Further investigation of these autoantibody propensity genes may lead to greater insight into the causes of autoantibody production and organ damage in SLE.

Type I interferons are potent regulators of innate and adaptive immunity and are implicated in the pathogenesis of systemic lupus erythematosus. Here we report that clinical and pathological lupus nephritis and serum anti-nuclear antibody levels are greatly attenuated in NZM 2328 mice deficient in type I IFN receptors (IFNAR). To determine if the inflammatory environment in NZM 2328 mice leads to IFNAR-regulated changes in dendritic cells (DC), the number, activation, and function of DC subsets were compared in 2 and 5 month-old (clinically healthy) female NZM and NZM-IFNAR-/- mice. Numbers of activated CD40hi plasmacytoid DC (pDC) were significantly increased in renal lymph nodes of 2 month-old NZM but not NZM-IFNAR-/- mice, suggesting an early IFNAR-dependent expansion and activation of pDC at disease sites. Relative to NZM spleens, NZM-IFNAR-/- spleens in 5 month-old mice were significantly decreased in size and contained reduced numbers of conventional DC (cDC) subsets, but not pDC. Splenic and renal lymph node NZM-IFNAR-/- DC analyzed directly ex vivo expressed significantly less CD40, CD86 and PDL1 than NZM DC. Upon activation with synthetic TLR9 ligands in vitro, splenic NZM-IFNAR-/- DC produced less IL-12p40/70 and TNFα than NZM DC. The limited IFNAR-/- DC response to endogenous activating stimuli correlated with reduced numbers of splenic activated memory CD4+ T cells and CD19+ B cells in older mice. Thus, IFNAR signaling significantly increases DC numbers, acquisition of antigen presentation competence, and pro-inflammatory function prior to onset of clinically apparent lupus disease.

Constitutive overexpression of B cell-activating factor belonging to the TNF family (BAFF) promotes development of systemic lupus erythematosus (SLE), and treatment of SLE mice with BAFF antagonists ameliorates disease. To determine whether SLE can develop de novo in BAFF-deficient hosts, BAFF-deficient New Zealand Mixed (NZM) 2328 (NZM.Baff−/−) mice were generated. In NZM.Baff−/− mice, spleen B cells (including CD5+ B1a and CD5− B1b B cells), germinal centers, Ig-secreting cells, and T cells were reduced in comparison to NZM.Baff+/+ mice. Serum total Ig and autoantibody levels were reduced at 4–6 mo but approached wild-type levels with increasing age, indicating that autoreactive B cells can survive and secrete autoantibodies despite the complete absence of BAFF. At least some of these autoantibodies are nephrophilic in that glomerular deposition of total IgG and IgG1 (but not of IgG2a, IgG2b, or C3) was substantial in NZM.Baff−/− mice by 12–13 mo of age. Despite proliferative glomerulonephritis, highlighted by widespread glomerular hyaline thrombi, being common among NZM.Baff−/− mice by 6–7 mo of age, severe proteinuria and mortality were greatly attenuated. These results demonstrate that the lifelong absence of BAFF does not protect NZM 2328 mice from serological autoimmunity and renal pathology. Nevertheless, the character of the renal pathology is altered, and the mice are largely spared from clinically overt disease (severe proteinuria and premature death). These observations may have profound ramifications for the use of BAFF antagonists in human SLE and related diseases.

TNF-α has both proinflammatory and immunoregulatory functions. Whereas a protective role for TNF administration in systemic lupus erythematosus (SLE)-prone (New Zealand Black × New Zealand White)F1 mice has been established, it remains uncertain whether this effect segregates at the individual TNFR. We generated SLE-prone New Zealand Mixed 2328 mice genetically deficient in TNFR1, in TNFR2, or in both receptors. Doubly-deficient mice developed accelerated pathological and clinical nephritis with elevated levels of circulating IgG anti-dsDNA autoantibodies and increased numbers of CD4+ T lymphocytes, especially activated memory (CD44highCD62Llow) CD4+ T cells. We show that these cells expressed a Th17 gene profile, were positive for IL-17 intracellular staining by FACS, and produced exogenous IL-17 in culture. In contrast, immunological, pathological, and clinical profiles of mice deficient in either TNFR alone did not differ from those in each other or from those in wild-type controls. Thus, total ablation of TNF-α-mediated signaling was highly deleterious to the host in the New Zealand Mixed 2328 SLE model. These observations may have profound ramifications for the use of TNF and TNFR antagonists in human SLE and related autoimmune disorders, as well as demonstrate, for the first time, the association of the Th17 pathway with an animal model of 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.

Systemic Lupus Erythematosus (SLE) disproportionately affects minorities, such as Hispanic-Americans. Prevalence of SLE is 3–5 times higher in Hispanic Americans (HA) than European derived populations, and have more active disease at the time of diagnosis, with more serious organ system involvement. HA is an admixed population, it is possible that there is an effect of admixture on the relative risk of disease. This admixture can create substantial increase of linkage disequilibrium (LD) in both magnitude and range, which can provide a unique opportunity for admixture mapping. Main objectives of this study are to (a) estimate hidden population structure in HA individuals; (b) estimate individual ancestry proportions and its impact on SLE risk; (c) assess impact of admixture on ITGAM association, a recently identified SLE susceptibility gene; and (d) estimate power of admixture mapping in HA. Our dataset contained 1,125 individuals, of whom 884 (657 SLE cases and 227 controls) were self classified as HA. Using 107 unlinked ancestry informative markers (AIMs) we estimated hidden population structure and individual ancestry in HA. Out of 5,671 possible pair-wise LD, 54% were statistically significant, indicating recent population admixture. The best fitted model for HA was a four population model with average ancestry of European (48%), American-Indian (40%), African (8%) and a fourth population (4%) with unknown ancestry. We also identified significant higher risk associated with American-Indian ancestry (OR=4.84, P=0.0001, 95%CI=2.14—10.95) on overall SLE. We showed that ITGAM is associated as a risk factor for SLE (OR= 2.06, P=8.74×10−5, 95%CI=1.44–2.97). This association is not affected by population substructure or admixture. We have demonstrated that HA have great potential and are an 3 appropriate population for admixture mapping. As expected, the case-only design is more powerful than case-control design, for any given admixture proportion or ancestry risk ratio.

Background

Many common disorders have multiple genetic components which convey increased susceptibility. SNPs have been used to identify genetic components which are associated with a disease. Unfortunately, many studies using these methods suffer from low reproducibility due to lack of power.

Results

We present a set of programs which implement a novel method for searching for disease-associated genes using prior information to select and order genes from publicly available databases by their prior likelihood of association with the disease. These programs were used in a published study of childhood-onset SLE which yielded novel associations with modest sample size.

Conclusion

Using prior information to decrease the size of the problem space to an amount commensurate with available samples and resources while maintaining appropriate power enables researchers to increase their likelihood of discovering reproducible associations.

Mechanisms that initiate lupus nephritis and cause progression to end-stage renal disease remain poorly understood. In this study, we show that lupus-prone New Zealand Mixed 2410 mice that develop a severe glomerulosclerosis and rapidly progressive renal disease overexpress IL-4 in vivo. In these mice, STAT6 deficiency or anti-IL-4 Ab treatment decreases type 2 cytokine responses and ameliorates kidney disease, particularly glomerulosclerosis, despite the presence of high levels of IgG anti-dsDNA Abs. STAT4 deficiency, however, decreases type 1 and increases type 2 cytokine responses, and accelerates nephritis, in the absence of high levels of IgG anti-dsDNA Abs. Thus, STAT6 and IL-4 may selectively contribute to the development of glomerulosclerosis, whereas STAT4 may play a role in autoantibody production.

We previously reported that the G allele of rs3853839 at 3′untranslated region (UTR) of Toll-like receptor 7 (TLR7) was associated with elevated transcript expression and increased risk for systemic lupus erythematosus (SLE) in 9,274 Eastern Asians [P = 6.5×10−10, odds ratio (OR) (95%CI) = 1.27 (1.17–1.36)]. Here, we conducted trans-ancestral fine-mapping in 13,339 subjects including European Americans, African Americans, and Amerindian/Hispanics and confirmed rs3853839 as the only variant within the TLR7-TLR8 region exhibiting consistent and independent association with SLE (Pmeta = 7.5×10−11, OR = 1.24 [1.18–1.34]). The risk G allele was associated with significantly increased levels of TLR7 mRNA and protein in peripheral blood mononuclear cells (PBMCs) and elevated luciferase activity of reporter gene in transfected cells. TLR7 3′UTR sequence bearing the non-risk C allele of rs3853839 matches a predicted binding site of microRNA-3148 (miR-3148), suggesting that this microRNA may regulate TLR7 expression. Indeed, miR-3148 levels were inversely correlated with TLR7 transcript levels in PBMCs from SLE patients and controls (R2 = 0.255, P = 0.001). Overexpression of miR-3148 in HEK-293 cells led to significant dose-dependent decrease in luciferase activity for construct driven by TLR7 3′UTR segment bearing the C allele (P = 0.0003). Compared with the G-allele construct, the C-allele construct showed greater than two-fold reduction of luciferase activity in the presence of miR-3148. Reduced modulation by miR-3148 conferred slower degradation of the risk G-allele containing TLR7 transcripts, resulting in elevated levels of gene products. These data establish rs3853839 of TLR7 as a shared risk variant of SLE in 22,613 subjects of Asian, EA, AA, and Amerindian/Hispanic ancestries (Pmeta = 2.0×10−19, OR = 1.25 [1.20–1.32]), which confers allelic effect on transcript turnover via differential binding to the epigenetic factor miR-3148.

Author Summary

Systemic lupus erythematosus (SLE) is a debilitating autoimmune disease contributed to by excessive innate immune activation involving toll-like receptors (TLRs, particularly TLR7/8/9) and type I interferon (IFN) signaling pathways. TLR7 responds against RNA–containing nuclear antigens and activates IFN-α pathway, playing a pivotal role in the development of SLE. While a genomic duplication of Tlr7 promotes lupus-like disease in the Y-linked autoimmune accelerator (Yaa) murine model, the lack of common copy number variations at TLR7 in humans led us to identify a functional single nucleotide polymorphism (SNP), rs3853839 at 3′ UTR of the TLR7 gene, associated with SLE susceptibility in Eastern Asians. In this study, we fine-mapped the TLR7-TLR8 region and confirmed rs3853839 exhibiting the strongest association with SLE in European Americans, African Americans, and Amerindian/Hispanics. Individuals carrying the risk G allele of rs3853839 exhibited increased TLR7 expression at the both mRNA and protein level and decreased transcript degradation. MicroRNA-3148 (miR-3148) downregulated the expression of non-risk allele (C) containing transcripts preferentially, suggesting a likely mechanism for increased TLR7 levels in risk-allele carriers. This trans-ancestral mapping provides evidence for the global association with SLE risk at rs3853839, which resides in a microRNA–gene regulatory site affecting TLR7 expression.

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.

Systemic lupus erythematosus (SLE) is an autoimmune disease with diverse clinical manifestations characterized by the development of pathogenic autoantibodies manifesting in inflammation of target organs such as the kidneys, skin and joints. Genome-wide association studies have identified genetic variants in the UBE2L3 region that are associated with SLE in subjects of European and Asian ancestry. UBE2L3 encodes an ubiquitin-conjugating enzyme, UBCH7, involved in cell proliferation and immune function. In this study, we sought to further characterize the genetic association in the region of UBE2L3 and use molecular methods to determine the functional effect of the risk haplotype. We identified significant associations between variants in the region of UBE2L3 and SLE in individuals of European and Asian ancestry that exceeded a Bonferroni corrected threshold (P < 1 × 10−4). A single risk haplotype was observed in all associated populations. Individuals harboring the risk haplotype display a significant increase in both UBE2L3 mRNA expression (P = 0.0004) and UBCH7 protein expression (P = 0.0068). The results suggest that variants carried on the SLE associated UBE2L3 risk haplotype influence autoimmunity by modulating UBCH7 expression.

Objective

Candidate gene and genome-wide association studies have identified several disease susceptibility loci in lupus patients. These studies have been largely performed in European-derived and Asian lupus patients. In this study, we examine if some of these same susceptibility loci increase lupus risk in African-American individuals.

Methods

Single nucleotide polymorphisms tagging 15 independent lupus susceptibility loci were genotyped in a set of 1,724 lupus patients and 2,024 normal healthy controls of African-American descent. The loci examined included: PTPN22, FCGR2A, TNFSF4, STAT4, CTLA4, PDCD1, PXK, BANK1, MSH5 (HLA region), CFB (HLA region), C8orf13-BLK region, MBL2, KIAA1542, ITGAM, and MECP2/IRAK1.

Results

We provide the first evidence for genetic association between lupus and five susceptibility loci in African-American patients (C8orf13-BLK, BANK1, TNFSF4, KIAA1542 andCTLA4; P values= 8.0 × 10−6, 1.9 × 10−5, 5.7 × 10−5, 0.00099, 0.0045, respectively). Further, we confirm the genetic association between lupus and five additional lupus susceptibility loci (ITGAM, MSH5, CFB, STAT4, and FCGR2A; P values= 7.5 × 10−11, 5.2 × 10−8, 8.7 × 10−7, 0.0058, and 0.0070, respectively), and provide evidence for a genome-wide significance for the association between ITGAM and MSH5 (HLA region) for the first time in African-American lupus patients.

Conclusion

These findings provide evidence for novel genetic susceptibility loci for lupus in African-Americans and demonstrate that the majority of lupus susceptibility loci examined confer lupus risk across multiple ethnicities.

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.

Objective

T cells from patients with SLE express increased amounts of PP2Ac which contribute to decreased production of IL-2. Because IL-2 is important in the regulation of several aspects of the immune response, it has been proposed that PP2Ac contributes to the expression of SLE. This study was designed to determine whether genetic variants of PPP2AC are linked to the expression of SLE and specific clinical manifestations and account for the increased expression of PP2Ac.

Methods

We conducted a trans-ethnic study consisting of 8,695 SLE cases and 7,308 controls from four different ancestries. Eighteen single-nucleotide polymorphisms (SNPs) across the PPP2CA were genotyped using an Illumina custom array. PPP2CA expression in SLE and control T cells was analyzed by real-time PCR.

Results

A 32-kb haplotype comprised of multiple SNPs of PPP2CA showed significant association with SLE in Hispanic Americans (HA), European Americans (EA) and Asians but not in African-Americans (AA). Conditional analyses revealed that SNP rs7704116 in intron 1 showed consistently strong association with SLE across Asian, EA and HA populations (pmeta=3.8×10−7, OR=1.3[1.14–1.31]). In EA, the largest ethnic dataset, the risk A allele of rs7704116 was associated with the presence of renal disease, anti-dsDNA and anti-RNP antibodies. PPP2CA expression was approximately 2-fold higher in SLE patients carrying the rs7704116 AG genotype than those carrying GG genotype (p = 0.008).

Conclusion

Our data provide the first evidence for an association between PPP2CA polymorphisms and elevated PP2Ac transcript levels in T cells, which implicates a new molecular pathway for SLE susceptibility in EA, HA and Asians.

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.

Systemic lupus erythematosus (SLE) is a complex autoimmune disease characterized by autoantibody production and organ damage. Lupus nephritis (LN) is one of the most severe manifestations of SLE. Multiple studies reported associations between renal diseases and variants in the non-muscle myosin heavy chain 9 (MYH9) and the neighboring apolipoprotein L 1 (APOL1) genes. We evaluated 167 variants spanning MYH9 for association with LN in a multiethnic sample. The two previously identified risk variants in APOL1 were also tested for association with LN in European-Americans (EAs) (N = 579) and African-Americans (AAs) (N = 407). Multiple peaks of association exceeding a Bonferroni corrected p-value of p < 2.03 × 10−3 were observed between LN and MYH9 in EAs (N=4620), with the most pronounced association at rs2157257 (p = 4.7 × 10−4; odds ratio [OR]=1.205). A modest effect with MYH9 was also detected in Gullah (rs8136069, p = 0.0019, OR = 2.304). No association between LN and MYH9 was found in AAs, Asians, Amerindians or Hispanics. This study provides the first investigation of MYH9 in LN in non-Africans and of APOL1 in LN in any population and presents novel insight into the potential role of MYH9 in LN in EAs.

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, OMIM 152700) 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 resequencing in ethnically diverse populations we fully characterized the TNFAIP3 risk haplotype and isolated a novel TT>A polymorphic dinucleotide associated with SLE in subjects of European (P = 1.58 × 10−8; odds ratio (OR) = 1.70) and Korean (P = 8.33 × 10−10; OR = 2.54) ancestry. This variant, located in a region of high conservation and regulatory potential, bound a nuclear protein complex comprised of NF-κB subunits with reduced avidity. Furthermore, 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), 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.

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×10−8) 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×10−128. 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 (ORhigh-low = 2.36, p = 9e−9), the immunologic criterion (ORhigh-low = 2.23, p = 3e−7), and age at diagnosis (ORhigh-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.

Author Summary

Systemic lupus erythematosus is a chronic disabling autoimmune disease, most commonly striking women in their thirties or forties. It can cause a wide variety of clinical manifestations, including kidney disease, arthritis, and skin disorders. Prognosis varies greatly depending on these clinical features, with kidney disease and related characteristics leading to greater morbidity and mortality. It is also complex genetically; while lupus runs in families, genes increase one's risk for lupus but do not fully determine the outcome. The interactions of multiple genes and/or interactions between genes and environmental factors may cause lupus, but the causes and disease pathways of this very heterogeneous disease are not well understood. By examining relationships between the presence of multiple lupus risk genes, lupus susceptibility, and clinical manifestations, we hope to better understand how lupus is triggered and by what biological pathways it progresses. We show in this work that certain clinical manifestations of lupus are highly associated with cumulative genetic variations, i.e. multiple risk alleles, while others are associated with a single variation or none at all.

Objective

Interleukin (IL) 21 is a member of the type I cytokine superfamily that exerts a variety of effects on the immune system including B cell activation, plasma cell differentiation, and immunoglobulin production. The expression of IL21R is reduced in B cells from lupus patients, while IL21 serum levels are increased in both lupus patients and some lupus-murine models. We recently reported that polymorphisms within the IL21 gene are associated with increased susceptibility to lupus. Herein, we examined the genetic association between SNPs within IL21R and lupus.

Methods

We genotyped 17 SNPs in the IL21R gene in two large cohorts of lupus patients and ethnically-matched healthy controls. Genotyping was performed with the Illumina BeadStation 500GX instrument using Illumina Infinum II genotyping assays.

Results

We identified and confirmed the association between rs3093301 within the IL21R gene and lupus in two independent European-derived and Hispanic cohorts (meta analysis odds ratio= 1.16, 95% CI= 1.08-1.25, meta analysis p=1.0×10-4).

Conclusion

We identified IL21R as a novel susceptibility gene for lupus.