The TNF α-induced protein 3 (TNFAIP3) is an ubiquitin-modifying enzyme and an essential negative regulator of inflammation. Genome-wide association studies have implicated the TNFAIP3 locus in susceptibility to autoimmune disorders in European cohorts, including rheumatoid arthritis, coronary artery disease, psoriasis, celiac disease, type 1 diabetes, inflammatory bowel disease, and systemic lupus erythematosus (SLE). There are two nonsynonymous coding polymorphisms in the deubiquitinating (DUB) domain of TNFAIP3: F127C, which is in high-linkage disequilibrium with reported SLE-risk variants, and A125V, which has not been previously studied. We conducted a case–control study in African-American SLE patients using these coding variants, along with tagging polymorphisms in TNFAIP3, and identified a novel African-derived risk haplotype that is distinct from previously reported risk variants (odds ratio = 1.6, p = 0.006). In addition, a rare protective haplotype was defined by A125V (odds ratio = 0.31, p = 0.027). Although A125V was associated with protection from SLE, surprisingly the same allele was associated with increased risk of inflammatory bowel disease. We tested the functional activity of nonsynonymous coding polymorphisms within TNFAIP3, and found that the A125V coding-change variant alters the DUB activity of the protein. Finally, we used computer modeling to depict how the A125V amino acid change in TNFAIP3 may affect the three-dimensional structure of the DUB domain to a greater extent than F127C. This is the first report of an association between TNFAIP3 polymorphisms and autoimmunity in African-Americans.
Functional characterization of causal variants present on risk haplotypes identified through genome-wide association studies (GWAS) is a primary objective of human genetics. In this report, we evaluate the function of a pair of tandem polymorphic dinucleotides, 42 kb downstream of the promoter of TNFAIP3, (rs148314165, rs200820567, collectively referred to as TT>A) recently nominated as causal variants responsible for genetic association of systemic lupus erythematosus (SLE) with tumor necrosis factor alpha inducible protein 3 (TNFAIP3). TNFAIP3 encodes the ubiquitin-editing enzyme, A20, a key negative regulator of NF-κB signaling. A20 expression is reduced in subjects carrying the TT>A risk alleles; however, the underlying functional mechanism by which this occurs is unclear. We used a combination of electrophoretic mobility shift assays (EMSA), mass spectrometry (MS), reporter assays, chromatin immunoprecipitation-PCR (ChIP-PCR) and chromosome conformation capture (3C) EBV transformed lymphoblastoid cell lines (LCL) from individuals carrying risk and non-risk TNFAIP3 haplotypes to characterize the effect of TT>A on A20 expression. Our results demonstrate that the TT>A variants reside in an enhancer element that binds NF-κB and SATB1 enabling physical interaction of the enhancer with the TNFAIP3 promoter through long-range DNA looping. Impaired binding of NF-κB to the TT>A risk alleles or knockdown of SATB1 expression by shRNA, inhibits the looping interaction resulting in reduced A20 expression. Together, these data reveal a novel mechanism of TNFAIP3 transcriptional regulation and establish the functional basis by which the TT>A risk variants attenuate A20 expression through inefficient delivery of NF-κB to the TNFAIP3 promoter. These results provide critical functional evidence supporting a direct causal role for TT>A in the genetic predisposition to SLE.
A key objective of human genetics is the identification and characterization of variants responsible for association with complex diseases. A pair of single nucleotide polymorphisms (rs148314165, rs200820567) 42 kb downstream from the promoter of TNFAIP3, have been proposed as the variants responsible for association with systemic lupus erythematosus based on comprehensive genetic and bioinformatic analyses. TNFAIP3 encodes for the ubiquitin-editing enzyme, A20, which plays a central role in maintaining immune system homeostasis through restriction of NF-κB signaling. Cells that carry this risk haplotype express low levels of TNFAIP3 compared to cells carrying the nonrisk haplotype. How the risk alleles of rs148314165 and rs200820567 might influence low TNFAIP3 expression is unknown. In this paper, we demonstrate that these variants reside in an enhancer element that binds NF-κB and SATB1 enabling the interaction of the enhancer with the TNFAIP3 promoter through long-range DNA looping. Impaired binding of NF-κB directly to the risk alleles or shRNA-mediated knockdown of SATB1 inhibits interaction of the enhancer with the TNFAIP3 promoter resulting in reduced A20 expression. These results clarify the functional mechanism by which rs148314165 and rs200820567 attenuate A20 expression and support a causal role for these variants in the predisposition to autoimmune disease.
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.
Systemic Lupus Erythematosus; UBE2L3; Multi Ethnic Association Study; UBCH7 Expression
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.
systemic lupus erythematosus; TNFAIP3; imputation; meta-analysis
Systemic lupus erythematosus (SLE) is a chronic multisystem genetically complex autoimmune disease characterised by the production of autoantibodies to nuclear and cellular antigens, tissue inflammation and organ damage. Genome-wide association studies have shown that variants within the major histocompatibility complex (MHC) region on chromosome 6 confer the greatest genetic risk for SLE in European and Chinese populations. However, the causal variants remain elusive due to tight linkage disequilibrium across disease-associated MHC haplotypes, the highly polymorphic nature of many MHC genes and the heterogeneity of the SLE phenotype.
A high-density case-control single nucleotide polymorphism (SNP) study of the MHC region was undertaken in SLE cohorts of Spanish and Filipino ancestry using a custom Illumina chip in order to fine-map association signals in these haplotypically diverse populations. In addition, comparative analyses were performed between these two datasets and a northern European UK SLE cohort. A total of 1433 cases and 1458 matched controls were examined.
Using this transancestral SNP mapping approach, novel independent loci were identified within the MHC region in UK, Spanish and Filipino patients with SLE with some evidence of interaction. These loci include HLA-DPB1, HLA-G and MSH5 which are independent of each other and HLA-DRB1 alleles. Furthermore, the established SLE-associated HLA-DRB1*15 signal was refined to an interval encompassing HLA-DRB1 and HLA-DQA1. Increased frequencies of MHC region risk alleles and haplotypes were found in the Filipino population compared with Europeans, suggesting that the greater disease burden in non-European SLE may be due in part to this phenomenon.
These data highlight the usefulness of mapping disease susceptibility loci using a transancestral approach, particularly in a region as complex as the MHC, and offer a springboard for further fine-mapping, resequencing and transcriptomic analysis.
Genome-wide association studies of rheumatoid arthritis (RA) have identified an association of the disease with a 6q23 region devoid of genes. TNFAIP3, an RA candidate gene, flanks this region, and polymorphisms in both the TNFAIP3 gene and the intergenic region are associated with systemic lupus erythematosus. We hypothesized that there is a similar association with RA, including polymorphisms in TNFAIP3 and the intergenic region.
To test this hypothesis, we selected tag-single nucleotide polymorphisms (SNPs) in both loci. They were analyzed in 1,651 patients with RA and 1,619 control individuals of Spanish ancestry.
Weak evidence of association was found both in the 6q23 intergenic region and in the TNFAIP3 locus. The rs582757 SNP and a common haplotype in the TNFAIP3 locus exhibited association with RA. In the intergenic region, two SNPs were associated, namely rs609438 and rs13207033. The latter was only associated in patients with anti-citrullinated peptide antibodies. Overall, statistical association was best explained by the interdependent contribution of SNPs from the two loci TNFAIP3 and the 6q23 intergenic region.
Our data are consistent with the hypothesis that several RA genetic factors exist in the 6q23 region, including polymorphisms in the TNFAIP3 gene, like that previously described for systemic lupus erythematosus.
Recent genome-wide association studies demonstrated association of single nucleotide polymorphisms (SNPs) in the TNFAIP3 region at 6q23 with systemic lupus erythematosus (SLE) in European-American populations. In this study, we investigated whether SNPs in the TNFAIP3 region are associated with SLE also in a Japanese population. A case-control association study was performed on the SNPs rs13192841, rs2230926, and rs6922466 in 318 Japanese SLE patients and 444 healthy controls. Association of rs2230926 G allele with SLE was replicated in Japanese (allelic association P = .033, odds ratio [OR] 1.47, recessive model P = .023, OR 8.52). The association was preferentially observed in the SLE patients with nephritis. When the TNFAIP3 mRNA levels of the HapMap samples were examined using GENEVAR database, the presence of TNFAIP3 rs2230926 G allele was associated with lower mRNA expression of TNFAIP3 (P = .013). These results indicated that TNFAIP3 is a susceptibility gene to SLE both in the Caucasian and Asian populations.
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.
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.
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.
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.
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).
These findings are the first to suggest that an ICAM–integrin-mediated pathway contributes to susceptibility to SLE.
Systemic lupus erythematosus (SLE) is a systemic multisystem autoimmune disorder influenced by genetic background and environmental factors. Our aim here was to replicate findings of associations between 7 of the implicated single nucleotide polymorphisms (SNPs) in IRF5, BLK, STAT4, TNFAIP3, SPP1, TNIP1 and ETS1 genes with susceptibility to childhood-onset SLE in the Japanese population. In particular, we focused on gender differences in allelic frequencies.
The 7 SNPs were genotyped using TaqMan assays in 75 patients with childhood-onset SLE and in 190 healthy controls. The relationship between the cumulative number of risk alleles and SLE manifestations was explored in childhood-onset SLE. Logistic regression was used to test the effect of each polymorphism on susceptibility to SLE, and Wilcoxon rank sum testing was used for comparison of total risk alleles. Data on rs7574865 in the STAT4 gene and rs9138 in SPP1 were replicated for associations with SLE when comparing cases and controls (corrected P values ranging from 0.0043 to 0.027). The rs2230926 allele of TNFAIP3 was associated with susceptibility to SLE in males, but after Bonferroni correction there were no significant associations with any of the other four SNPs in IRF5, BLK, TNIP1 and ETS1 genes. The cumulative number of risk alleles was significantly increased in childhood-onset SLE relative to healthy controls (P = 0.0000041). Male SLE patients had a slightly but significantly higher frequency of the TNFAIP3 (rs2230926G) risk allele than female patients (odds ratio [OR] = 4.05, 95% confidence interval [95%CI] = 1.46–11.2 P<0.05).
Associations of polymorphisms in STAT4 and SPP1 with childhood-onset SLE were confirmed in a Japanese population. Although these are preliminary results for a limited number of cases, TNFAIP3 rs2230926G may be an important predictor of disease onset in males. We also replicated findings that the cumulative number of risk alleles was significantly increased in childhood-onset SLE.
Genetic variation in interferon regulatory factor 5 (IRF5) has been associated with risk of developing systemic lupus erythematosus (SLE), and this association is largely dependent upon anti-Ro autoantibodies. We studied a unique cohort of anti-Ro positive individuals with diverse diagnoses to determine if IRF5 genotype associated with maternal diagnosis or progression of autoimmunity.
We genotyped haplotype-tagging polymorphisms in IRF5 in 93 European ancestry subjects recruited to the Research Registry for Neonatal Lupus who all had high titer anti-Ro autoantibodies and a child with neonatal lupus (NL), and allele frequencies were compared to non-autoimmune controls. The mothers diagnoses included SLE, Sjogren’s syndrome (SS), undifferentiated autoimmune syndrome (UAS), and asymptomatic.
The SLE-risk haplotype of IRF5 was enriched in all anti-Ro positive subjects except those with SS (OR = 2.55, p=8.8×10−4). Even asymptomatic individuals with anti-Ro antibodies were enriched for the SLE-risk haplotype (OR=2.69, p=0.019). The same haplotype was more prevalent in subjects who were initially asymptomatic, but developed symptomatic SLE during follow up (OR=5.83, p=0.0024). Interestingly, SS was associated with two minor IRF5 haplotypes, and these same haplotypes were decreased in frequency in those with SLE and UAS.
The IRF5 SLE-risk haplotype was associated with anti-Ro antibodies in asymptomatic individuals as well as progression to SLE in asymptomatic anti-Ro positive individuals. SS in NL mothers was associated with different IRF5 haplotypes. These data suggest that IRF5 polymorphisms play a role in serologic autoimmunity in humans and may promote the progression to clinical autoimmunity.
systemic lupus erythematosus; interferon; autoantibodies; neonatal lupus; Sjogren’s syndrome
A haplotype of the interferon regulatory factor 5 (IRF5) gene has been associated with the risk of developing systemic lupus erythematosus (SLE), and our previous studies have demonstrated that high levels of serum interferon-α (IFNα) activity are a heritable risk factor for SLE. The aim of this study was to determine whether the IRF5 SLE risk haplotype mediates the risk of SLE by predisposing patients to the development of high levels of serum IFNα activity.
IFNα levels in 199 SLE patients of European and Hispanic ancestry were measured with a sensitive functional reporter cell assay. The rs2004640, rs3807306, rs10488631, and rs2280714 single-nucleotide polymorphisms (SNPs) in IRF5 were genotyped in these patients. Haplotypes were categorized as SLE risk, neutral, or protective based on published data.
SLE patients with risk/risk and risk/neutral IRF5 genotypes had higher serum IFNα activity than did those with protective/protective and neutral/protective genotypes (P = 0.025). This differential effect of IRF5 genotype on serum IFNα levels was driven largely by SLE patients who were positive for either anti–RNA binding protein (anti-RBP) or anti–double-stranded DNA (anti-dsDNA) autoantibodies (P = 0.012 for risk/risk or risk/neutral versus protective/protective or neutral/protective). The rs3807306 genotype was independently associated with high serum IFNα in this autoantibody group. We found no difference in IFNα activity according to IRF5 genotype in patients lacking either type of autoantibody or in patients positive for both classes of autoantibody.
The IRF5 SLE risk haplotype is associated with higher serum IFNα activity in SLE patients, and this effect is most prominent in patients positive for either anti-RBP or anti-dsDNA autoantibodies. This study demonstrates the biologic relevance of the SLE risk haplotype of IRF5 at the protein level.
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.
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.
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.
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 a multisystem complex autoimmune disease of uncertain etiology (OMIM 152700). Over recent years a genetic component to SLE susceptibility has been established1–3. Recent successes with association studies in SLE have identified genes including IRF5 (refs. 4,5) and FCGR3B6. Two tumor necrosis factor (TNF) superfamily members located within intervals showing genetic linkage with SLE are TNFSF4 (also known as OX40L; 1q25), which is expressed on activated antigen-presenting cells (APCs)7,8 and vascular endothelial cells9, and also its unique receptor, TNFRSF4 (also known as OX40; 1p36), which is primarily expressed on activated CD4+ T cells10. TNFSF4 produces a potent co-stimulatory signal for activated CD4+ T cells after engagement of TNFRSF4 (ref. 11). Using both a family-based and a case-control study design, we show that the upstream region of TNFSF4 contains a single risk haplotype for SLE, which is correlated with increased expression of both cell-surface TNFSF4 and the TNFSF4 transcript. We hypothesize that increased expression of TNFSF4 predisposes to SLE either by quantitatively augmenting T cell–APC interaction or by influencing the functional consequences of T cell activation via TNFRSF4.
Immunoregulatory cytokine interleukin-10 (IL-10) is elevated in sera from patients with systemic lupus erythematosus (SLE) correlating with disease activity. The established association of IL10 with SLE and other autoimmune diseases led us to fine map causal variant(s) and to explore underlying mechanisms. We assessed 19 tag SNPs, covering the IL10 gene cluster including IL19, IL20 and IL24, for association with SLE in 15,533 case and control subjects from four ancestries. The previously reported IL10 variant, rs3024505 located at 1 kb downstream of IL10, exhibited the strongest association signal and was confirmed for association with SLE in European American (EA) (P = 2.7×10−8, OR = 1.30), but not in non-EA ancestries. SNP imputation conducted in EA dataset identified three additional SLE-associated SNPs tagged by rs3024505 (rs3122605, rs3024493 and rs3024495 located at 9.2 kb upstream, intron 3 and 4 of IL10, respectively), and SLE-risk alleles of these SNPs were dose-dependently associated with elevated levels of IL10 mRNA in PBMCs and circulating IL-10 protein in SLE patients and controls. Using nuclear extracts of peripheral blood cells from SLE patients for electrophoretic mobility shift assays, we identified specific binding of transcription factor Elk-1 to oligodeoxynucleotides containing the risk (G) allele of rs3122605, suggesting rs3122605 as the most likely causal variant regulating IL10 expression. Elk-1 is known to be activated by phosphorylation and nuclear localization to induce transcription. Of interest, phosphorylated Elk-1 (p-Elk-1) detected only in nuclear extracts of SLE PBMCs appeared to increase with disease activity. Co-expression levels of p-Elk-1 and IL-10 were elevated in SLE T, B cells and monocytes, associated with increased disease activity in SLE B cells, and were best downregulated by ERK inhibitor. Taken together, our data suggest that preferential binding of activated Elk-1 to the IL10 rs3122605-G allele upregulates IL10 expression and confers increased risk for SLE in European Americans.
Systemic lupus erythematosus (SLE), a debilitating autoimmune disease characterized by the production of pathogenic autoantibodies, has a strong genetic basis. Variants of the IL10 gene, which encodes cytokine interleukin-10 (IL-10) with known function of promoting B cell hyperactivity and autoantibody production, are associated with SLE and other autoimmune diseases, and serum IL-10 levels are elevated in SLE patients correlating with increased disease activity. In this study, to discover SLE-predisposing causal variant(s), we assessed variants within the genomic region containing IL10 and its gene family member IL19, IL20 and IL24 for association with SLE in case and control subjects from diverse ancestries. We identified SLE-associated SNP rs3122605 located at 9.2 kb upstream of IL10 as the most likely causal variant in subjects of European ancestry. The SLE-risk allele of rs3122605 was dose-dependently associated with elevated IL10 expression at both mRNA and protein levels in peripheral blood samples from SLE patients and controls, which could be explained, at least in part, by its preferential binding to Elk-1, a transcription factor activated in B cells during active disease of SLE patients. Elk-1-mediated IL-10 overexpression could be downregulated by inhibiting activation of mitogen-activated protein kinases, suggesting a potential therapeutic target for SLE.
Protein tyrosine phosphatase non-receptor type 22 (PTPN22) is a negative regulator of T-cell activation associated with several autoimmune diseases, including systemic lupus erythematosus (SLE). Missense rs2476601 is associated with SLE in individuals with European ancestry. Since the rs2476601 risk allele frequency differs dramatically across ethnicities, we assessed robustness of PTPN22 association with SLE and its clinical sub-phenotypes across four ethnically diverse populations. Ten SNPs were genotyped in 8220 SLE cases and 7369 controls from in European-Americans (EA), African-Americans (AA), Asians (AS), and Hispanics (HS). We performed imputation-based association followed by conditional analysis to identify independent associations. Significantly associated SNPs were tested for association with SLE clinical sub-phenotypes, including autoantibody profiles. Multiple testing was accounted for by using false discovery rate. We successfully imputed and tested allelic association for 107 SNPs within the PTPN22 region and detected evidence of ethnic-specific associations from EA and HS. In EA, the strongest association was at rs2476601 (P = 4.7×10−9, OR = 1.40 (95% CI = 1.25–1.56)). Independent association with rs1217414 was also observed in EA, and both SNPs are correlated with increased European ancestry. For HS imputed intronic SNP, rs3765598, predicted to be a cis-eQTL, was associated (P = 0.007, OR = 0.79 and 95% CI = 0.67–0.94). No significant associations were observed in AA or AS. Case-only analysis using lupus-related clinical criteria revealed differences between EA SLE patients positive for moderate to high titers of IgG anti-cardiolipin (aCL IgG >20) versus negative aCL IgG at rs2476601 (P = 0.012, OR = 1.65). Association was reinforced when these cases were compared to controls (P = 2.7×10−5, OR = 2.11). Our results validate that rs2476601 is the most significantly associated SNP in individuals with European ancestry. Additionally, rs1217414 and rs3765598 may be associated with SLE. Further studies are required to confirm the involvement of rs2476601 with aCL IgG.
Recent genome-wide and candidate gene association studies in large numbers of systemic lupus erythematosus (SLE) patients have suggested approximately 30 susceptibility genes. These genes are involved in three types of biological processes, including immune complex processing, toll-like receptor function and type I interferon production, and immune signal transduction in lymphocytes, and they may contribute to the pathogenesis of SLE. To better understand the genetic risk factors of SLE, we investigated the associations of seven SLE susceptibility genes in a Chinese population, including FCGR3A, FCGR2A, TNFAIP3, TLR9, TREX1, ETS1 and TNIP1.
A total of 20 SNPs spanning the seven SLE susceptibility genes were genotyped in a sample of 564 unrelated SLE patients and 504 unrelated healthy controls recruited from Yunnan, southwestern China. The associations of SNPs with SLE were assessed by statistical analysis.
Five SNPs in two genes (TNFAIP3 and ETS1) were significantly associated with SLE (corrected P values ranging from 0.03 to 5.5 × 10-7). Through stratified analysis, TNFAIP3 and ETS1 showed significant associations with multiple SLE subphenotypes (such as malar rash, arthritis, hematologic disorder and antinuclear antibody) while TNIP1 just showed relatively weak association with onset age. The associations of the SNPs in the other four genes were not replicated.
The replication analysis indicates that TNFAIP3, ETS1 and TNIP1 are probably common susceptibility genes for SLE in Chinese populations, and they may contribute to the pathogenesis of multiple SLE subphenotypes.
A20 is a ubiquitin modifying enzyme that restricts NF-κB signals and protects cells against tumor necrosis factor (TNF) induced programmed cell death. Given recent data linking A20 (TNFAIP3) with human B cell lymphomas and systemic lupus erythematosus (SLE), we have generated mice bearing a floxed allele of Tnfaip3 to interrogate A20’s roles in regulating B cell functions. A20-deficient B cells are hyper-responsive to multiple stimuli and display exaggerated NF-κB responses to CD40 induced signals. Mice expressing absent or hypomorphic amounts of A20 in B cells possess elevated numbers of germinal center B cells, autoantibodies, and glomerular immunoglobulin deposits. A20 deficient B cells are resistant to Fas mediated cell death, likely due to increased expression of NF-κB-dependent anti-apoptotic proteins such as Bcl-x. These findings show that A20 can restrict B cell survival, while A20 protects other cells from TNF induced cell death. Our studies demonstrate how reduced A20 expression predisposes to autoimmunity.
A20; ubiquitin; haploinsufficiency; NFkB; B cell selection; autoimmunity; programmed cell death; SLE
TNFAIP3 interacting protein 1, TNIP1 (ABIN-1) is involved in inhibition of nuclear factor-κB (NF-κB) activation by interacting with TNF alpha-induced protein 3, A20 (TNFAIP3), an established susceptibility gene to systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA). Recent genome-wide association studies revealed association of TNIP1 with SLE in the Caucasian and Chinese populations. In this study, we investigated whether the association of TNIP1 with SLE was replicated in a Japanese population. In addition, association of TNIP1 with RA was also examined.
A case-control association study was conducted on the TNIP1 single nucleotide polymorphism (SNP) rs7708392 in 364 Japanese SLE patients, 553 RA patients and 513 healthy controls.
Association of TNIP1 rs7708392C was replicated in Japanese SLE (allele frequency in SLE: 76.5%, control: 69.9%, P = 0.0022, odds ratio [OR] 1.40, 95% confidence interval [CI] 1.13-1.74). Notably, the risk allele frequency in the healthy controls was considerably greater in Japanese (69.9%) than in Caucasians (24.3%). A tendency of stronger association was observed in the SLE patients with renal disorder (P = 0.00065, OR 1.60 [95%CI 1.22-2.10]) than in all SLE patients (P = 0.0022, OR 1.40 [95%CI 1.13-1.74]). Significant association with RA was not observed, regardless of the carriage of human leukocyte antigen DR β1 (HLA-DRB1) shared epitope. Significant gene-gene interaction between TNIP1 and TNFAIP3 was detected neither in SLE nor RA.
Association of TNIP1 with SLE was confirmed in a Japanese population. TNIP1 is a shared SLE susceptibility gene in the Caucasian and Asian populations, but the genetic contribution appeared to be greater in the Japanese and Chinese populations because of the higher risk allele frequency. Taken together with the association of TNFAIP3, these observations underscore the crucial role of NF-κB regulation in the pathogenesis of SLE.
Variants of the osteopontin (OPN) gene have been associated with systemic lupus erythematosus (SLE) susceptibility and cytokine profiles in SLE patients. It is not known whether these alleles are associated with specific clinical phenotypes in SLE. We studied 252 well-characterized SLE patients from a multiethnic cohort, genotyping the rs11730582, rs28357094, rs6532040, and rs9138 SNPs in the OPN gene. Ancestry informative markers were used to control for genetic ancestry. The SLE-risk allele rs9138C in the 3′ UTR region was associated with photosensitivity in lupus patients across all ancestral backgrounds (meta-analysis OR = 3.2, 95% CI = 1.6–6.5, P = 1.0 × 10−3). Additionally, the promoter variant rs11730582C demonstrated suggestive evidence for association with two hematologic traits: thrombocytopenia (OR = 2.1, P = 0.023) and hemolytic anemia (OR = 2.6, P = 0.036). These clinical associations with SNPs in the promoter and 3′ UTR regions align with previously reported SLE-susceptibility SNPs in OPN and suggest potential roles for these variants in antibody-mediated cytopenias and skin inflammation in SLE.
Our understanding of the genetic basis of systemic lupus erythematosus (SLE) has been rapidly advanced using large-scale, case–control, candidate gene studies as well as genome-wide association studies during the past 3 years. These techniques have identified more than 30 robust genetic associations with SLE including genetic variants of HLA and Fcγ receptor genes, IRF5, STAT4, PTPN22, TNFAIP3, BLK, BANK1, TNFSF4 and ITGAM. Most SLE-associated gene products participate in key pathogenic pathways, including Toll-like receptor and type I interferon signaling pathways, immune regulation pathways and those that control the clearance of immune complexes. Disease-associated loci that have not yet been demonstrated to have important functions in the immune system might provide new clues to the underlying molecular mechanisms that contribute to the pathogenesis or progression of SLE. Of note, genetic risk factors that are shared between SLE and other immune-related diseases highlight common pathways in the pathophysiology of these diseases, and might provide innovative molecular targets for therapeutic interventions.
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.
The TNFAIP3 gene has been associated with psoriasis, rheumatoid arthritis, type 1 diabetes mellitus, systemic lupus erythematosus and celiac disease. TNFAIP3 encodes A20, a TNF-α-inducible zinc finger protein thought to limit NF-κB mediated immune responses. In this study we report association of response of psoriasis to TNF blockers with two TNFAIP3 SNPs (rs2230926 in exon 3 and rs610604 in intron 6) and their haplotypes. Treatment response was self-evaluated using a 0–5 visual analog scale in 433 psoriasis patients who received TNF blockers. Confirmation was sought in 199 psoriasis and psoriatic arthritis patients from Toronto who were followed prospectively. Response variables were dichotomized separately in the two cohorts, yielding similar proportions of good responses. While significant associations were observed only for the Michigan cohort, fixed-effects meta-analysis retained significant association between dosage of the G allele of rs610604 and good combined response to all TNF blockers (OR = 1.50, pcorr = 0.050) and etanercept (OR = 1.64, pcorr = 0.016). The rs2230926 T–rs610604 G haplotype was similarly associated. By demonstrating an association with therapeutic response, these results provide a clinically relevant functional correlate to the recently described genetic association between psoriasis and TNFAIP3.
psoriasis; TNFAIP3; tumor necrosis factor; pharmacogenetics
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.
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.
We aimed to investigate whether the effect size of the systemic lupus erythematosus (SLE) risk alleles varies across European subpopulations.
European SLE patients (n = 1,742) and ethnically matched healthy controls (n = 2,101) were recruited at 17 centres from 10 different countries. Only individuals with self-reported ancestry from the country of origin were included. In addition, participants were genotyped for top ancestry informative markers and for 25 SLE associated SNPs. The results were used to compare effect sizes between the Central Eureopan and Southern European subgroups.
Twenty of the 25 SNPs showed independent association with SLE, These SNPs showed a significant bias to larger effect sizes in the Southern subgroup, with 15/20 showing this trend (P = 0.019) and a larger mean odds ratio of the 20 SNPs (1.46 vs. 1.34, P = 0.02) as well as a larger difference in the number of risk alleles (2.06 vs. 1.63, P = 0.027) between SLE patients and controls than for Central Europeans. This bias was reflected in a very significant difference in the cumulative genetic risk score (4.31 vs. 3.48, P = 1.8 × 10-32). Effect size bias was accompanied by a lower number of SLE risk alleles in the Southern subjects, both patients and controls, the difference being more marked between the controls (P = 1.1 × 10-8) than between the Southern and Central European patients (P = 0.016). Seven of these SNPs showed significant allele frequency clines.
Our findings showed a bias to larger effect sizes of SLE loci in the Southern Europeans relative to the Central Europeans together with clines of SLE risk allele frequencies. These results indicate the need to study risk allele clines and the implications of the polygenic model of inheritance in SLE.