Objectives

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.

Methods

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.

Results

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

Conclusion

These data suggest contributions of both IRAK1 and MECP2 to SLE susceptibility.

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.

Objectives

Altered signaling in B-cells is a predominant feature of systemic lupus erythematosus (SLE). The genes BANK1 and BLK were recently described as associated with SLE. BANK1 codes for a B-cell-specific cytoplasmic protein involved in B-cell receptor signaling and BLK codes for an Src tyrosine kinase with important roles in B-cell development. To characterize the role of BANK1 and BLK in SLE, we performed a genetic interaction analysis hypothesizing that genetic interactions could reveal functional pathways relevant to disease pathogenesis.

Methods

We Used the method GPAT16 to analyze the gene-gene interactions of BANK1 and BLK. Confocal microscopy was used to investigate co-localization, and immunoprecipitation was used to verify the physical interaction of BANK1 and BLK.

Results

Epistatic interactions between BANK1 and BLK polymorphisms associated with SLE were observed in a discovery set of 279 patients and 515 controls from Northern Europe. A meta-analysis with 4399 European individuals confirmed the genetic interactions between BANK1 and BLK.

As BANK1 was identified as a binding partner of the Src tyrosine kinase LYN, we tested the possibility that BANK1 and BLK could also show a protein-protein interaction. We demonstrated co-immunoprecipitation and co-localization of BLK and BANK1. In a Daudi cell line and primary naïve B-cells the endogenous binding was enhanced upon B-cell receptor stimulation using anti-IgM antibodies.

Conclusions

Here, we show a genetic interaction between BANK1 and BLK, and demonstrate that these molecules interact physically. Our results have important consequences for the understanding of SLE and other autoimmune diseases and identify a potential new signaling pathway.

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.

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.

Objectives

To perform fine mapping of the autoimmunity susceptibility gene BLK and identify functional variants involved in SLE.

Methods

Genotyping of 1163 European SLE patients and 1482 controls and imputation were performed covering the BLK gene with 158 SNPs. Logistic regression analysis was done using PLINK and conditional analyses using GENABEL’s test score. Transfections of BLK constructs on HEK293 cells containing the novel mutation or the wild-type form were made to analyse their effect on protein half-life using a protein stability assay, cycloheximide and Western blot. CHiP-qPCR for NFkB binding.

Results

Fine mapping of BLK identified two independent genetic effects with functional consequences: one represented by two tightly linked associated haplotype blocks significantly enriched for NFκB-binding sites and numerous putative regulatory variants whose risk alleles correlated with low BLK mRNA levels. Binding of NFkBp50 and p65 to an associated 1.2Kb haplotype segment was confirmed. A second independent genetic effect was represented by an Ala71Thr, low-frequency missense substitution with an OR = 2.31 (95% c.i. 1.38–3.86). The 71Thr decreased BLK protein half-life.

Conclusions

Our results show that rare and common regulatory variants in BLK are involved in disease susceptibility and both, albeit independently lead to reduced levels of BLK protein.

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.

Argentine population genetic structure was examined using a set of 78 ancestry informative markers (AIMs) to assess the contributions of European, Amerindian, and African ancestry in 94 individuals members of this population. Using the Bayesian clustering algorithm STRUCTURE, the mean European contribution was 78%, the Amerindian contribution was 19.4%, and the African contribution was 2.5%. Similar results were found using weighted least mean square method: European, 80.2%; Amerindian, 18.1%; and African, 1.7%. Consistent with previous studies the current results showed very few individuals (four of 94) with greater than 10% African admixture. Notably, when individual admixture was examined, the Amerindian and European admixture showed a very large variance and individual Amerindian contribution ranged from 1.5 to 84.5% in the 94 individual Argentine subjects. These results indicate that admixture must be considered when clinical epidemiology or case control genetic analyses are studied in this population. Moreover, the current study provides a set of informative SNPs that can be used to ascertain or control for this potentially hidden stratification. In addition, the large variance in admixture proportions in individual Argentine subjects shown by this study suggests that this population is appropriate for future admixture mapping studies.

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, 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.

SUMMARY

We have previously developed and validated a self-administered questionnaire, modeled after the Systemic Lupus International Collaborating Clinics Damage Index (SDI), the Lupus Damage Index Questionnaire (LDIQ), which may allow the ascertainment of this construct in Systemic Lupus Erythematosus (SLE) patients followed in the community and thus expand observations made about damage. We have now translated, back-translated and adapted the LDIQ to Spanish, Portuguese and French and applied it to patients followed at academic and non-academic centers in North and South America, Portugal and Spain while their physicians scored the SDI. A total of 887 patients (659 Spanish-, 140 Portuguese- and 80 French-speaking) and 40 physicians participated. Overall patients scored higher than their physicians (total score and all domains) for all versions of the LDIQ. Infrequent manifestations had less than optional clinimetric properties but overall agreement was over 95% for the majority of items. The larger sample size may explain the higher correlations observed among the Spanish-speaking patients. Pending minor adjustments, the LDIQ may prove to be useful in community-based studies. The relationship between the LDIQ and other outcome parameters is currently being investigated in a different patient sample.

A gain-of-function R620W polymorphism in the PTPN22 gene, encoding the lymphoid tyrosine phosphatase LYP, has recently emerged as an important risk factor for human autoimmunity. Here we report that another missense substitution (R263Q) within the catalytic domain of LYP leads to reduced phosphatase activity. High-resolution structural analysis revealed the molecular basis for this loss of function. Furthermore, the Q263 variant conferred protection against human systemic lupus erythematosus, reinforcing the proposal that inhibition of LYP activity could be beneficial in human autoimmunity.

Immune-mediated nephritis contributes to disease in systemic lupus erythematosus, Goodpasture syndrome (caused by antibodies specific for glomerular basement membrane [anti-GBM antibodies]), and spontaneous lupus nephritis. Inbred mouse strains differ in susceptibility to anti-GBM antibody–induced and spontaneous lupus nephritis. This study sought to clarify the genetic and molecular factors that may be responsible for enhanced immune-mediated renal disease in these models. When the kidneys of 3 mouse strains sensitive to anti-GBM antibody–induced nephritis were compared with those of 2 control strains using microarray analysis, one-fifth of the underexpressed genes belonged to the kallikrein gene family, which encodes serine esterases. Mouse strains that upregulated renal and urinary kallikreins exhibited less evidence of disease. Antagonizing the kallikrein pathway augmented disease, while agonists dampened the severity of anti-GBM antibody–induced nephritis. In addition, nephritis-sensitive mouse strains had kallikrein haplotypes that were distinct from those of control strains, including several regulatory polymorphisms, some of which were associated with functional consequences. Indeed, increased susceptibility to anti-GBM antibody–induced nephritis and spontaneous lupus nephritis was achieved by breeding mice with a genetic interval harboring the kallikrein genes onto a disease-resistant background. Finally, both human SLE and spontaneous lupus nephritis were found to be associated with kallikrein genes, particularly KLK1 and the KLK3 promoter, when DNA SNPs from independent cohorts of SLE patients and controls were compared. Collectively, these studies suggest that kallikreins are protective disease-associated genes in anti-GBM antibody–induced nephritis and lupus.

Objectives

To perform fine mapping of the autoimmunity susceptibility gene BLK and identify functional variants involved in systemic lupus erythematosus (SLE).

Methods

Genotyping of 1163 European SLE patients and 1482 controls and imputation were performed covering the BLK gene with 158 single-nucleotide polymorphisms. Logistic regression analysis was done using PLINK and conditional analyses using GENABEL's test score. Transfections of BLK constructs on HEK293 cells containing the novel mutation or the wild type form were analysed for their effect on protein half-life using a protein stability assay, cycloheximide and western blot. CHiP-qPCR for detection of nuclear factor κ B (NFkB) binding.

Results

Fine mapping of BLK identified two independent genetic effects with functional consequences: one represented by two tightly linked associated haplotype blocks significantly enriched for NFκB-binding sites and numerous putative regulatory variants whose risk alleles correlated with low BLK mRNA levels. Binding of NFkBp50 and p65 to an associated 1.2 Kb haplotype segment was confirmed. A second independent genetic effect was represented by an Ala71Thr, low-frequency missense substitution with an OR=2.31 (95% CI 1.38 to 3.86). The 71Thr decreased BLK protein half-life.

Conclusions

These results show that rare and common regulatory variants in BLK are involved in disease susceptibility and both, albeit independently, lead to reduced levels of BLK protein.