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.

Objective

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

Methods

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

Results

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

Conclusion

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

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

Objective

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.

Defective clearance of apoptotic cells has been shown in systemic lupus erythematosus (SLE) and is postulated to enhance autoimmune responses by increasing access to intracellular autoantigens. Until now, research has emphasized inherited rather than acquired impairment of apoptotic cell engulfment in the pathogenesis of SLE. Here, we confirm previous results that efficient removal of apoptotic cells (efferocytosis) is bolstered in the presence of wild type mouse serum, through the C3 deposition on the apoptotic cell surface. In contrast, sera from three mouse models of SLE, MerKD, MRLlpr and NZBWF1, did not support and in fact actively inhibited apoptotic cell uptake. IgG autoantibodies were responsible for the inhibition, through the blockade of C3 recognition by macrophages. Consistent with this, IgG removal reversed the inhibitory activity within autoimmune serum and purified autoimmune IgG blocked both the detection of C3 on apoptotic cells and C3-dependent efferocytosis. Sera from SLE patients demonstrated elevated anti-C3b IgG that blocked detection of C3 on apoptotic cells, activity that was not found in healthy controls or patients with rheumatoid arthritis, nor in mice prior to the onset of autoimmunity. We propose that the suppression of apoptotic cell disposal by antibodies against deposited C3 may contribute to increasing severity and/or exacerbations in SLE.

Objective

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

Methods

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

Results

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

Conclusion

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

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.

The Sle1c subinterval on distal murine chromosome 1 confers loss of tolerance to chromatin. Cr2, which encodes complement receptors 1 and 2 (CR1/CR2; CD35/CD21), is a strong candidate gene for lupus susceptibility within this interval based on structural and functional alterations in its protein products. Complement receptor 1 related protein/gene Y (Crry) lies 10 kb from Cr2 and encodes a ubiquitously expressed complement regulatory protein that could also play a role in the pathogenesis of systemic lupus erythematosus (SLE). Crry derived from B6.Sle1c congenic mice migrated at a higher molecular weight by SDS-PAGE compared to B6 Crry as a result of differential glycosylation. A single-nucleotide polymorphism (SNP) in the first short consensus repeat (SCR) of Sle1c Crry introduced a novel N-linked glycosylation site likely responsible for this structural alteration. Five additional SNPs in the signal peptide and SCR1 of Sle1c Crry were identified. However, the cellular expression of B6 and B6.Sle1c Crry and their ability to regulate the classical pathway of complement were not significantly different, and although soluble Sle1c Crry regulated the alternative pathway of complement more efficiently than B6 Crry, as a membrane protein, it regulated the alternative pathway equivalently to B6 Crry. These data fail to provide evidence for a functional effect of the structural alterations in Sle1c Crry, and suggest that the role of Cr2 in the Sle1c autoimmune phenotypes can be isolated in recombinant congenic mice containing both genes.

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.

Genetic factors influence susceptibility to systemic lupus erythematosus (SLE). A recent family-based analysis in Caucasian and Chinese populations provided evidence for association of single-nucleotide polymorphisms (SNPs) in the complement receptor 2 (CR2/CD21) gene with SLE. Here we confirmed this result in a case-control analysis of an independent European-derived population including 2084 patients with SLE and 2853 healthy controls. A haplotype formed by the minor alleles of three CR2 SNPs (rs1048971, rs17615, rs4308977) showed significant association with decreased risk of SLE (30.4% in cases vs. 32.6% in controls, P = 0.016, OR = 0.90 [0.82-0.98]). Two of these SNPs are in exon 10, directly 5′ of an alternatively spliced exon preferentially expressed in follicular dendritic cells (FDC), and the third is in the alternatively spliced exon. Effects of these SNPs as well as a fourth SNP in exon 11 (rs17616) on alternative splicing were evaluated. We found that the minor alleles of these SNPs decreased splicing efficiency of exon 11 both in vitro and ex vivo. These findings further implicate CR2 in the pathogenesis of SLE and suggest that CR2 variants alter the maintenance of tolerance and autoantibody production in the secondary lymphoid tissues where B cells and FDCs interact.

Rationale: Severe pulmonary arterial hypertension (SPH) is a frequently lethal condition characterized by pulmonary vascular remodeling and right heart strain or failure. SPH is also often associated with autoimmune and collagen vascular disorders.

Objectives: To study the effects of T cells on the development of experimental SPH.

Methods: Athymic nude rats lacking T cells were treated with a single subcutaneous injection of vascular endothelial growth factor (VEGF) receptor blocker SU5416 (20 mg/kg) to induce pulmonary vascular endothelial cell apoptosis. Immunohistochemical analysis and IL-4 levels of the lung tissue were performed. Cell death and proliferation were assessed by Western blot and immunohistochemistry.

Measurements and Main Results: In contrast to SU5416-treated euthymic rats that develop SPH only in combination with chronic hypoxia, athymic nude rats developed SPH and vascular remodeling (similar to clinical SPH) at normoxic conditions as demonstrated by measurements of pulmonary artery pressure and right ventricle hypertrophy. Pulmonary arterioles became occluded with proliferating endothelial cells and were surrounded by mast cells, B cells, and macrophages. IL-4, proliferating cell nuclear antigen, and collagen type I levels were markedly increased in SU5416-treated athymic rat lungs. Antibody deposition was noted along the vascular endothelium in rats with SPH. Finally, protection from SPH was conferred by immune challenge with spleen cells from euthymic nude rats.

Conclusions: These studies demonstrate the importance of a complete, intact immune system in protecting against pulmonary angioproliferation in this new model of SPH as well as the importance of intact VEGF receptor signaling for lung endothelial cell homeostasis.