The majority of the genetic variance of systemic lupus erythematosus (SLE) remains unexplained by the common disease-common variant hypothesis. Rare variants, which are not detectable by genome-wide association studies because of their low frequencies, are predicted to explain part of this ”missing heritability.” However, recent studies identifying rare variants within known disease-susceptibility loci have failed to show genetic associations because of their extremely low frequencies, leading to the questioning of the contribution of rare variants to disease susceptibility. A common (minor allele frequency = 17.4% in cases) nonsynonymous coding variant rs1143679 (R77H) in ITGAM (CD11b), which forms half of the heterodimeric integrin receptor, complement receptor 3 (CR3), is robustly associated with SLE and has been shown to impair CR3-mediated phagocytosis.
We resequenced ITGAM in 73 SLE cases and identified two previously unidentified, case-specific nonsynonymous variants, F941V and G1145S. Both variants were genotyped in 2,107 and 949 additional SLE cases, respectively, to estimate their frequencies in a disease population. An in vitro model was used to assess the impact of F941V and G1145S, together with two nonsynonymous ITGAM polymorphisms, A858V (rs1143683) and M441T (rs11861251), on CR3-mediated phagocytosis. A paired two-tailed t test was used to compare the phagocytic capabilities of each variant with that of wild-type CR3.
Both rare variants, F941V and G1145S, significantly impair CR3-mediated phagocytosis in an in vitro model (61% reduction, P = 0.006; 26% reduction, P = 0.0232). However, neither of the common variants, M441T and A858V, had an effect on phagocytosis. Neither rare variant was observed again in the genotyping of additional SLE cases, suggesting that there frequencies are extremely low.
Our results add further evidence to the functional importance of ITGAM in SLE pathogenesis through impaired phagocytosis. Additionally, this study provides a new example of the identification of rare variants in common-allele-associated loci, which, because of their extremely low frequencies, are not statistically associated. However, the demonstration of their functional effects adds support to their contribution to disease risk, and questions the current notion of dismissing the contribution of very rare variants on purely statistical analyses.
Naturally occurring variation in gene copy number is increasingly recognized as a heritable source of susceptibility to genetically complex diseases. Here we report strong association between FCGR3B copy number and risk of systemic lupus erythematosus (P = 2.7 × 10-8), microscopic polyangiitis (P = 2.9 × 10-4) and Wegener’s granulomatosis in two independent cohorts from the UK (P = 3 × 10-3) and France (P = 1.1 × 10-4). We did not observe this association in the organ-specific Graves’ disease or Addison’s disease. Our findings suggest that low FCGR3B copy number, and in particular complete FCGR3B deficiency, has a key role in the development of systemic autoimmunity.
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.
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.
Endogenous retroviruses are implicated in murine lupus nephritis. They provide a source of nephritogenic retroviral gp70-anti-gp70 immune complexes through the production of serum gp70 protein and anti-gp70 autoantibodies as a result of the activation of TLR7. The Sgp (serum gp70 production) loci identified in lupus-prone mice play distinct roles for the expression of different classes of endogenous retroviruses, as Sgp3 regulates the transcription of xenotropic, polytropic and modified polytropic (mPT) viruses, and Sgp4 the transcription of only xenotropic viruses. In the present study, we extended these analyses to a third locus, Sgp5, using BALB/c mice congenic for the NZW-derived Sgp5 allele and also explored the possible interaction of Sgp3 and Sgp4 loci to promote the expression of endogenous retroviruses and serum gp70. The analysis of Sgp5 BALB/c congenic mice demonstrated that the Sgp5 locus enhanced the expression of xenotropic and mPT viruses, thereby upregulating the production of serum gp70. These data indicate a distinct action of the Sgp5 locus on the expression of endogenous retroviruses, as compared with two other Sgp loci. Moreover, comparative analysis of C57BL/6 double congenic mice for Sgp3 and Sgp4 loci with single congenic mice revealed that Sgp3 and Sgp4 acted synergistically to elevate the transcription of the potentially replicationcompetent Xmv18 provirus and the production of serum gp70. This indicates that the combined effect of three different Sgp loci markedly enhance the expression of endogenous retroviruses and their gene product, serum gp70, thereby contributing to the formation of nephritogenic gp70-anti-gp70 immune complexes in murine lupus.
Systemic lupus erythematosus; Endogenous retrovirus; gp70
Sjögren’s syndrome is a common autoimmune disease (~0.7% of European Americans) typically presenting as keratoconjunctivitis sicca and xerostomia. In addition to strong association within the HLA region at 6p21 (Pmeta=7.65×10−114), we establish associations with IRF5-TNPO3 (Pmeta=2.73×10−19), STAT4 (Pmeta=6.80×10−15), IL12A (Pmeta =1.17×10−10), FAM167A-BLK (Pmeta=4.97×10−10), DDX6-CXCR5 (Pmeta=1.10×10−8), and TNIP1 (Pmeta=3.30×10−8). Suggestive associations with Pmeta<5×10−5 were observed with 29 regions including TNFAIP3, PTTG1, PRDM1, DGKQ, FCGR2A, IRAK1BP1, ITSN2, and PHIP amongst others. These results highlight the importance of genes involved in both innate and adaptive immunity in Sjögren’s syndrome.
Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by autoantibody production and altered type I interferon expression. Genetic surveys and genome-wide association studies have identified more than 30 SLE susceptibility genes. One of these genes, TNIP1, encodes the ABIN1 protein. ABIN1 functions in the immune system by restricting the NF-κB signaling. In order to better understand the genetic factors that influence association with SLE in genes that regulate the NF-κB pathway, we analyzed a dense set of genetic markers spanning TNIP1 and TAX1BP1, as well as the TNIP1 homolog, TNIP2, in case-control sets of diverse ethnic origins.
We fine-mapped TNIP1, TNIP2, and TAX1BP1 in a total of 8372 SLE cases and 7492 healthy controls from European-ancestry, African-American, Hispanic, East Asian, and African-American Gullah populations. Levels of TNIP1 mRNA and ABIN1 protein were analyzed using quantitative RT-PCR and Western blotting, respectively, in EBV-transformed human B cell lines.
We found significant associations between genetic variants within TNIP1 and SLE but not in TNIP2 or TAX1BP1. After resequencing and imputation, we identified two independent risk haplotypes within TNIP1 in individuals of European-ancestry that were also present in African-American and Hispanic populations. These risk haplotypes produced lower levels of TNIP1 mRNA and ABIN1 protein suggesting they harbor hypomorphic functional variants that influence susceptibility to SLE by restricting ABIN1 expression.
Our results confirmed the association signals between SLE and TNIP1 variants in multiple populations and provide new insight into the mechanism by which TNIP1 variants may contribute to SLE pathogenesis.
Amerindian-Europeans, Asians and African-Americans have an excess morbidity from SLE and higher prevalence of lupus nephritis than Caucasians. The aim of this study was to analyze the relationship between genetic ancestry and socio-demographic characteristics and clinical features in a large cohort of Amerindian-European SLE patients.
A total of 2116 SLE patients of Amerindian-European origin and 4001 SLE patients of European descent with clinical data were used in the study. Genotyping of 253 continental ancestry informative markers was performed on the Illumina platform. The STRUCTURE and ADMIXTURE software were used to determine genetic ancestry of each individual. Correlation between ancestry and socio-demographic and clinical data were analyzed using logistic regression.
The average Amerindian genetic ancestry of 2116 SLE patients was 40.7%. There was an increased risk of having renal involvement (P<0.0001, OR= 3.50 95%CI 2.63-4.63) and an early age of onset with the presence of Amerindian genetic ancestry (P<0.0001). Amerindian ancestry protected against photosensitivity (P<0.0001, OR= 0.58 95%CI 0.44-0.76), oral ulcers (P<0.0001, OR= 0.55 95%CI 0.42-0.72), and serositis (P<0.0001, OR= 0.56 95%CI 0.41-0.75) after adjustment by age, gender and age of onset. However, gender and age of onset had stronger effects on malar rash, discoid rash, arthritis and neurological involvement than genetic ancestry.
In general, genetic Amerindian ancestry correlates with lower socio-demographic status and increases the risk for developing renal involvement and SLE at an earlier age of onset.
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.
Little is known about the genetic etiology of systemic lupus erythematosus (SLE) in individuals of African ancestry, despite its higher prevalence and greater disease severity. Overproduction of nitric oxide (NO) and reactive oxygen species are implicated in the pathogenesis and severity of SLE, making NO synthases and other reactive intermediate related genes biological candidates for disease susceptibility. This study analyzed variation in reactive intermediate genes for association with SLE in two populations with African ancestry.
A total of 244 SNPs from 53 regions were analyzed in non-Gullah African Americans (AA; 1432 cases and 1687 controls) and the genetically more homogeneous Gullah of the Sea Islands of South Carolina (133 cases and 112 controls) and. Single-marker, haplotype, and two-locus interaction tests were computed for these populations.
The glutathione reductase gene GSR (rs2253409, P=0.0014, OR [95% CI]=1.26 [1.09–1.44]) was the most significant single-SNP association in AA. In the Gullah, the NADH dehydrogenase NDUFS4 (rs381575, P=0.0065, OR [95%CI]=2.10 [1.23–3.59]) and nitric oxide synthase gene NOS1 (rs561712, P=0.0072, OR [95%CI]=0.62 [0.44–0.88]) were most strongly associated with SLE. When both populations were analyzed together, GSR remained the most significant effect (rs2253409, P=0.00072, OR [95%CI]=1.26 [1.10–1.44]). Haplotype and two-locus interaction analyses also uncovered different loci in each population.
These results suggest distinct patterns of association with SLE in African-derived populations; specific loci may be more strongly associated within select population groups.
systemic lupus erythematosus; African Americans; genetic association studies; oxygen compounds; single nucleotide polymorphism
C-src tyrosine kinase, Csk, physically interacts with the intracellular phosphatase Lyp (PTPN22) and can modify the activation state of downstream Src kinases, such as Lyn, in lymphocytes. We identified an association of Csk with systemic lupus erythematosus (SLE) and refined its location to an intronic polymorphism rs34933034 (OR 1.32, p = 1.04 × 10−9). The risk allele is associated with increased CSK expression and augments inhibitory phosphorylation of Lyn. In carriers of the risk allele, B cell receptor (BCR)-mediated activation of mature B cells, as well as plasma IgM, are increased. Moreover, the fraction of transitional B cells is doubled in the cord blood of carriers of the risk allele compared to non-risk haplotypes due to an expansion of the late transitional cells, a stage targeted by selection mechanisms. This suggests that the Lyp-Csk complex increases susceptibility to lupus at multiple maturation and activation points of B cells.
Systemic lupus erythematosus (SLE) is a common systemic autoimmune disease with complex etiology but strong clustering in families (λS = ~30). We performed a genome-wide association scan using 317,501 SNPs in 720 women of European ancestry with SLE and in 2,337 controls, and we genotyped consistently associated SNPs in two additional independent sample sets totaling 1,846 affected women and 1,825 controls. Aside from the expected strong association between SLE and the HLA region on chromosome 6p21 and the previously confirmed non-HLA locus IRF5 on chromosome 7q32, we found evidence of association with replication (1.1 × 10−7 < Poverall < 1.6 × 10−23; odds ratio 0.82–1.62)in four regions: 16p11.2 (ITGAM), 11p15.5 (KIAA1542), 3p14.3 (PXK) and 1q25.1 (rs10798269). We also found evidence for association (P < 1 × 10−5) at FCGR2A, PTPN22 and STAT4, regions previously associated with SLE and other autoimmune diseases, as well as at ≥9 other loci (P < 2 × 10−7). Our results show that numerous genes, some with known immune-related functions, predispose to SLE.
Relative deficiency of pentraxin proteins is implicated in the pathogenesis of systemic lupus erythematosus. The C-reactive protein (CRP) response is defective in patients with acute flares of disease, and mice with targeted deletions of the serum amyloid P component gene (Sap) develop a lupus-like illness. In humans, the genes for CRP (CRP) and SAP (APCS) map to 1q23.2 within an interval linked with SLE. We have investigated the candidate genes CRP and APCS in two cohorts totalling 586 UK simplex SLE families. The inheritance of an intronic dinucleotide repeat and seven single nucleotide polymorphisms in the CRP and APCS genes was examined by application of family-based tests of association and linkage. Basal levels of CRP were influenced independently by two polymorphisms at the CRP locus, CRP 2 and CRP 4. Furthermore, the latter polymorphism was linked/associated with SLE and antinuclear autoantibody production. Thus, the polymorphism associated with reduced basal CRP was also associated with the development of SLE. These data support the hypothesis that defective disposal of potentially immunogenic material is a contributory factor in lupus pathogenesis. The identification of polymorphisms that determine basal CRP levels has implications in ischaemic heart disease, where CRP level is an important predictor of risk.
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.
Systemic lupus erythematosus (SLE) is a heterogeneous autoimmune disease with significant immune system aberrations resulting from complex heritable genetics as well as environmental factors. TRAF6 is a candidate gene for SLE, which has a major role in several signaling pathways that are important for immunity and organ development.
Fifteen single-nucleotide polymorphisms (SNPs), across TRAF6 were evaluated in 7,490 SLE and 6,780 control subjects from different ancestries. Population-based case-control association analyses and meta-analyses were performed. P values, false discovery rate q values, and odds ratios with 95% confidence intervals were calculated.
Evidence of associations in multiple SNPs was detected. The best overall p values were obtained for SNPs rs5030437 and rs4755453 (p=7.85×10−5 and p=4.73×10−5, respectively) without significant heterogeneity among populations (p=0.67 and p=0.50 in Q-statistic). In addition, rs540386 previously reported to be associated with RA was found to be in LD with these two SNPs (r2= 0.95) and demonstrated evidence of association with SLE in the same direction (meta-analysis p=9.15×10−4, OR=0.89, 95%CI=0.83–0.95). Thrombocytopenia improved the overall results in different populations (meta-analysis p=1.99×10−6, OR=0.57, 95%CI=0.45–0.72, for rs5030470). Finally evidence of family based association in 34 African-American pedigrees with the presence of thrombocytopenia were detected in one available SNP rs5030437 with Z score magnitude of 2.28 (p=0.02) under a dominant model.
Our data indicate the presence of association of TRAF6 with SLE in agreement with the previous report of association with RA. These data provide further support for the involvement of TRAF6 in the pathogenesis of autoimmunity.
TRAF6; polymorphism; systemic lupus erythematosus
Systemic lupus erythematosus (SLE) is a sexually dimorphic autoimmune disease which is more common in women, but affected men often experience a more severe disease. The genetic basis of sexual dimorphism in SLE is not clearly defined. A study was undertaken to examine sex-specific genetic effects among SLE susceptibility loci.
A total of 18 autosomal genetic susceptibility loci for SLE were genotyped in a large set of patients with SLE and controls of European descent, consisting of 5932 female and 1495 male samples. Sex-specific genetic association analyses were performed. The sex–gene interaction was further validated using parametric and nonparametric methods. Aggregate differences in sex-specific genetic risk were examined by calculating a cumulative genetic risk score for SLE in each individual and comparing the average genetic risk between male and female patients.
A significantly higher cumulative genetic risk for SLE was observed in men than in women. (P = 4.52×10−8) A significant sex–gene interaction was seen primarily in the human leucocyte antigen (HLA) region but also in IRF5, whereby men with SLE possess a significantly higher frequency of risk alleles than women. The genetic effect observed in KIAA1542 is specific to women with SLE and does not seem to have a role in men.
The data indicate that men require a higher cumulative genetic load than women to develop SLE. These observations suggest that sex bias in autoimmunity could be influenced by autosomal genetic susceptibility loci.
An increase in leucocyte apoptosis and impaired clearance of apoptotic cells has been observed in patients with systemic lupus erythematosus (SLE). Apoptotic cells are likely to be a key source of autoantigens in SLE as they express many of the nuclear autoantigens (in surface blebs and apoptotic bodies) that are relevant to this disease. The clearance of apoptotic cells is usually a rapid process, such that few cells are usually seen in the extracellular environment in vivo. We report a case in which multiple apoptotic bodies were observed in the bone marrow of a patient with SLE that was complicated by an immune‐mediated pancytopenia. We have subsequently examined the frequency of apoptotic cells, identified morphologically, and by caspase‐3 staining in bone‐marrow trephine samples taken from patients with SLE over a 10‐year period of follow‐up. A high proportion of bone marrows contained apoptotic debris. The novel demonstration of apoptotic bodies in vivo in patients with SLE is unusual and supports the notion that the marrow may be a target organ in the disease. Their abundance is also consistent with the hypothesis that normal clearance mechanisms are defective and/or overwhelmed in SLE.
apoptosis; bone marrow; pancytopenia; SLE
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
The Xq28 region containing IRAK1 and MECP2 has been identified as a risk locus for systemic lupus erythematosus (SLE) in previous genetic association studies. However, due to the strong linkage disequilibrium between IRAK1 and MECP2, it remains unclear which gene is affected by the underlying causal variant(s) conferring risk of SLE.
We fine-mapped ≥136 SNPs in a ~227kb region on Xq28, containing IRAK1, MECP2 and 7 adjacent genes (L1CAM, AVPR2, ARHGAP4, NAA10, RENBP, HCFC1 and TMEM187), for association with SLE in 15,783 case-control subjects derived from 4 different ancestral groups.
Multiple SNPs showed strong association with SLE in European Americans, Asians and Hispanics at P<5×10−8 with consistent association in subjects with African ancestry. Of these, 6 SNPs located in the TMEM187-IRAK1-MECP2 region captured the underlying causal variant(s) residing in a common risk haplotype shared by all 4 ancestral groups. Among them, rs1059702 best explained the Xq28 association signals in conditional testings and exhibited the strongest P value in trans-ancestral meta-analysis (Pmeta=1.3×10−27, OR=1.43), and thus was considered to be the most-likely causal variant. The risk allele of rs1059702 results in the amino acid substitution S196F in IRAK1 and had previously been shown to increase NF-κB activity in vitro. We also found that the homozygous risk genotype of rs1059702 was associated with lower mRNA levels of MECP2, but not IRAK1, in SLE patients (P=0.0012) and healthy controls (P=0.0064).
These data suggest contributions of both IRAK1 and MECP2 to SLE susceptibility.
Systemic Lupus Erythematosus; Gene Polymorphism; Xq28; IRAK1; MECP2
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.
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.
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.
Several members of the GIMAP gene family have been suggested as being involved in different aspects of the immune system in different species. Recently, a mutation in the GIMAP5 gene was shown to cause lymphopenia in a rat model of autoimmune insulin‐dependent diabetes. Thus it was hypothesised that genetic variation in GIMAP5 may be involved in susceptibility to other autoimmune disorders where lymphopenia is a key feature, such as systemic lupus erythematosus (SLE).
Material and methods
To investigate this, seven single nucleotide polymorphisms in GIMAP5 were analysed in five independent sets of family‐based SLE collections, containing more than 2000 samples.
A significant increase in SLE risk associated with the most common GIMAP5 haplotype was found (OR 1.26, 95% CI 1.02 to 1.54, p = 0.0033). In families with probands diagnosed with trombocytopenia, the risk was increased (OR 2.11, 95% CI 1.09 to 4.09, p = 0.0153). The risk haplotype bears a polymorphic polyadenylation signal which alters the 3′ part of GIMAP5 mRNA by producing an inefficient polyadenylation signal. This results in higher proportion of non‐terminated mRNA for homozygous individuals (p<0.005), a mechanism shown to be causal in thalassaemias. To further assess the functional effect of the polymorphic polyadenylation signal in the risk haplotype, monocytes were treated with several cytokines affecting apoptosis. All the apoptotic cytokines induced GIMAP5 expression in two monocyte cell lines (1.5–6 times, p<0.0001 for all tests).
Taken together, the data suggest the role of GIMAP5 in the pathogenesis of SLE.
genetic association; autoimmune; apoptosis; susceptibility gene
The complement C4 locus is in the class III region of the MHC, and exhibits copy number variation. Complement C4 null alleles have shown association with a number of diseases including systemic lupus erythematosus (SLE). However, most studies to date have used protein immunophenotyping and not direct interrogation of the genome to determine C4 null allele status. Moreover, a lack of accurate C4 gene copy number (GCN) estimation and tight linkage disequilibrium across the disease-associated MHC haplotypes has confounded attempts to establish whether or not these associations are causal. We have therefore developed a high through-put paralog ratio test (PRT) in association with two restriction enzyme digest variant ratio tests (REDVRs) to determine total C4 GCN, C4A GCN, and C4B GCN. In the densely genotyped CEU cohort we show that this method is accurate and reproducible when compared to gold standard Southern blot copy number estimation with a discrepancy rate of 9%. We find a broad range of C4 GCNs in the CEU and the 1958 British Birth Cohort populations under study. In addition, SNP-C4 CNV analyses show only moderate levels of correlation and therefore do not support the use of SNP genotypes as proxies for complement C4 GCN.
complement C4; CNV; lupus; paralog ratio test
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.
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.
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).
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.
Systemic Lupus Erythematosus; UBE2L3; Multi Ethnic Association Study; UBCH7 Expression