Recently an association was demonstrated between the single nucleotide polymorphism (SNP), rs10516487, within the B-cell gene BANK1 and systemic lupus erythematosus (SLE) as a consequence of a genome wide association study of this disease in European and Argentinean populations. In a bid for replication, we examined the effects of the R61H non-synonymous variant with respect to SLE in our genotyped American cohorts of European and African ancestry. Utilizing data from our ongoing genome-wide association study in our cohort of 178 Caucasian SLE cases and 1808 Caucasian population-based controls plus 148 African American (AA) SLE cases and 1894 AA population-based controls we investigated the association of the previously described non-synonymous SNP at the BANK1 locus with the disease in the two ethnicities separately. Using a Fisher’s exact test, the minor allele frequency (MAF) of rs10516487 in the Caucasian cases was 22.6% while it was 31.2% in Caucasian controls, yielding a protective odds ratio (OR) of 0.64 (95% CI 0.49-0.85; one-sided p = 7.07 × 10−4). Furthermore, the MAF of rs10516487 in the AA cases was 18.7% while it was 23.3% in AA controls, yielding a protective OR of 0.75 (95% CI 0.55–1.034; one-sided p = 0.039). The OR of the BANK1 variant in our study cohorts is highly comparable with that reported previously in a South American/European SLE case-control cohort (OR = 0.72). As such, R61H in the BANK1 gene confers a similar magnitude of SLE protection, not only in European Americans, but also in African Americans.
systemic lupus erythematosus; African Americans; European Americans; BANK1 gene
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.
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.
Alternative splicing; systemic lupus erythematosus; complement receptors; single-nucleotide polymorphisms; B cells; follicular dendritic cells
Recent genome-wide association studies demonstrated association of single nucleotide polymorphisms (SNPs) in the TNFAIP3 region at 6q23 with systemic lupus erythematosus (SLE) in European-American populations. In this study, we investigated whether SNPs in the TNFAIP3 region are associated with SLE also in a Japanese population. A case-control association study was performed on the SNPs rs13192841, rs2230926, and rs6922466 in 318 Japanese SLE patients and 444 healthy controls. Association of rs2230926 G allele with SLE was replicated in Japanese (allelic association P = .033, odds ratio [OR] 1.47, recessive model P = .023, OR 8.52). The association was preferentially observed in the SLE patients with nephritis. When the TNFAIP3 mRNA levels of the HapMap samples were examined using GENEVAR database, the presence of TNFAIP3 rs2230926 G allele was associated with lower mRNA expression of TNFAIP3 (P = .013). These results indicated that TNFAIP3 is a susceptibility gene to SLE both in the Caucasian and Asian populations.
The autoimmune disease systemic lupus erythematosus (SLE) is characterized by loss of tolerance to nuclear antigens and a heightened inflammatory environment, which together result in end organ damage. Lyn-deficient mice, a model of SLE, lack an inhibitor of B- and myeloid cell activation. This results in B-cell hyperresponsiveness, plasma cell (PC) accumulation, autoantibodies, and glomerulonephritis (GN). IL-21 is associated with autoimmunity in mice and humans and promotes B-cell differentiation and class switching. Here, we explore the role of IL-21 in the autoimmune phenotypes of lyn−/− mice. We find that IL-21 mRNA is reduced in the spleens of lyn−/−IL-6−/− and lyn−/−Btklo mice, neither of which produce pathogenic autoantibodies or develop significant GN. While IL-21 is dispensable for PC accumulation and IgM autoantibodies in lyn−/− mice, it is required for anti-DNA IgG antibodies and some aspects of T cell activation. Surprisingly, GN still develops in lyn−/−IL-21−/− mice. This likely results from the presence of IgG autoantibodies against a limited set of non-DNA antigens. These studies identify a specific role for IL-21 in the class switching of anti-DNA B-cells and demonstrate that neither IL-21 nor anti-DNA IgG is required for kidney damage in lyn−/− mice.
autoimmunity; Lyn; IL-21; IL-6; autoantibody; glomerulonephritis
The major histocompatibility complex (MHC) on chromosome 6p21 is a key contributor to the genetic basis of systemic lupus erythemathosus (SLE). Although SLE affects African Americans disproportionately compared to European Americans, there has been no comprehensive analysis of the MHC region in relationship to SLE in African Americans. We conducted a screening of the MHC region for 1,536 single nucleotide polymorphisms (SNPs) and the deletion of the C4A gene in a SLE case-control study (380 cases, 765 age-matched controls) nested within the prospective Black Women’s Health Study. We also genotyped 1,509 ancestral informative markers throughout the genome to estimate European ancestry in order to control for population stratification due to population admixture. The most strongly associated SNP with SLE was the rs9271366 (odds ratio, OR = 1.70, p = 5.6×10−5) near the HLA-DRB1 gene. Conditional haplotype analysis revealed three other SNPs, rs204890 (OR = 1.86, p = 1.2×10−4), rs2071349 (OR = 1.53, p = 1.0×10−3), and rs2844580 (OR = 1.43, p = 1.3×10−3) to be associated with SLE independent of the rs9271366 SNP. In univariate analysis, the OR for the C4A deletion was 1.38, p = 0.075, but after simultaneous adjustment for the other four SNPs the odds ratio was 1.01, p = 0.98. A genotype score combining the four newly identified SNPs showed an additive risk according to the number of high-risk alleles (OR = 1.67 per high-risk allele, p< 0.0001). Our strongest signal, the rs9271366 SNP, was also associated with higher risk of SLE in a previous Chinese genome-wide association study (GWAS). In addition, two SNPs found in a GWAS of European ancestry women were confirmed in our study, indicating that African Americans share some genetic risk factors for SLE with European and Chinese subjects. In summary, we found four independent signals in the MHC region associated with risk of SLE in African American women.
systemic lupus erythemathosus; African Americans; major histocompatibility complex; single nucleotide polymorphisms
The increased risk of thrombosis in systemic lupus erythematosus (SLE) may be partially explained by interrelated genetic pathways for thrombosis and SLE. In a case-control analysis, we investigated whether 33 established and novel single nucleotide polymorphisms (SNP) in 20 genes involved in hemostasis pathways that have been associated with deep venous thrombosis in the general population were risk factors for SLE development among Asians.
Patients in the discovery cohort were enrolled in one of two North American SLE cohorts. Patients in the replication cohort were enrolled in one of four Asian or two North American cohorts. SLE cases met American College of Rheumatology classification criteria. We first genotyped 263 Asian SLE and 357 healthy Asian control individuals for 33 SNPs using Luminex multiplex technology in the discovery phase, and then used Taqman and Immunochip assays to examine 5 SNPs in up to an additional 1496 cases and 993 controls in the Replication phase. SLE patients were compared to healthy controls for association with minor alleles in allelic models. Principal components analysis was used to control for intra-Asian ancestry in an analysis of the replication cohort.
Two genetic variants in the gene VKORC1, rs9934438 and rs9923231, were highly significant in both the discovery and replication cohorts: OR(disc) = 2.45 (p=2×10−9), OR(rep) = 1.53 (p=5×10−6) and OR(disc) = 2.40 (p=6×10−9), OR(rep) = 1.53 (p=5×10−6), respectively. These associations were significant in the replication cohort after adjustment for intra-Asian ancestry: rs9934438 OR(adj) = 1.34 (p=0.0029) and rs9923231 OR(adj) = 1.34 (p=0.0032).
Genetic variants in VKORC1, involved in vitamin K reduction and associated with DVT, are associated with SLE development in Asians. These results suggest intersecting genetic pathways for the development of SLE and thrombosis.
systemic lupus erythematosus; single nucleotide polymorphisms; genetic risk factors
The Toll-like receptor 7 (TLR7) gene, encoded on human chromosome Xp22.3, is crucial for type I interferon production. A recent multicenter study in East Asian populations, comprising Chinese, Korean and Japanese participants, identified an association of a TLR7 single-nucleotide polymorphism (SNP) located in the 3' untranslated region (3' UTR), rs3853839, with systemic lupus erythematosus (SLE), especially in males, although some difference was observed among the tested populations. To test whether additional polymorphisms contribute to SLE in Japanese, we systematically analyzed the association of TLR7 with SLE in a Japanese female population.
A case-control association study was conducted on eight tag SNPs in the TLR7 region, including rs3853839, in 344 Japanese females with SLE and 274 healthy female controls.
In addition to rs3853839, two SNPs in intron 2, rs179019 and rs179010, which were in moderate linkage disequilibrium with each other (r2 = 0.53), showed an association with SLE (rs179019: P = 0.016, odds ratio (OR) 2.02, 95% confidence interval (95% CI) 1.15 to 3.54; rs179010: P = 0.018, OR 1.75, 95% CI 1.10 to 2.80 (both under the recessive model)). Conditional logistic regression analysis revealed that the association of the intronic SNPs and the 3' UTR SNP remained significant after we adjusted them for each other. When only the patients and controls carrying the risk genotypes at the 3' UTR SNPpositionwere analyzed, the risk of SLE was significantly increased when the individuals also carried the risk genotypes at both of the intronic SNPs (P = 0.0043, OR 2.45, 95% CI 1.31 to 4.60). Furthermore, the haplotype containing the intronic risk alleles in addition to the 3' UTR risk allele was associated with SLE under the recessive model (P = 0.016, OR 2.37, 95% CI 1.17 to 4.80), but other haplotypes were not associated with SLE.
The TLR7 intronic SNPs rs179019 and rs179010 are associated with SLE independently of the 3' UTR SNP rs3853839 in Japanese women. Our findings support a role of TLR7 in predisposition for SLE in Asian populations.
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.
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 systemic autoimmune disease characterized by increased type I interferons (IFNs) and multiorgan inflammation frequently targeting the skin. IFN-kappa is a type I IFN expressed in skin. A pooled genome-wide scan implicated the IFNK locus in SLE susceptibility. We studied IFNK single nucleotide polymorphisms (SNPs) in 3982 SLE cases and 4275 controls, composed of European (EA), African-American (AA), and Asian ancestry. rs12553951C was associated with SLE in EA males (odds ratio = 1.93, P = 2.5 × 10−4), but not females. Suggestive associations with skin phenotypes in EA and AA females were found, and these were also sex-specific. IFNK SNPs were associated with increased serum type I IFN in EA and AA SLE patients. Our data suggest a sex-dependent association between IFNK SNPs and SLE and skin phenotypes. The serum IFN association suggests that IFNK variants could influence type I IFN producing plasmacytoid dendritic cells in affected skin.
African Americans (AA) disproportionately develop lupus nephritis (LN) relative to European Americans and familial clustering supports causative genes. Since MYH9 underlies approximately 40% of end-stage renal disease (ESRD) in AA, we tested for genetic association with LN.
Seven MYH9 single nucleotide polymorphisms (SNPs) and the E1 risk haplotype were tested for association with LN in three cohorts of AA.
A preliminary analysis revealed that the MYH9 E1 risk haplotype was associated with ESRD in 25 cases with presumed systemic lupus erythematosus (SLE)-associated ESRD, compared to 735 non-SLE controls (odds ratio 3.1; p = 0.010 recessive). Replication analyses were performed in 583 AA with SLE in the PROFILE cohort (318 with LN; 265 with SLE but without nephropathy) and 60 AA from the NIH (39 with LN; 21 with SLE but without nephropathy). Analysis of the NIH and larger PROFILE cohorts, as well as a combined analysis, did not support this association.
These results suggest that AA with ESRD and coincident SLE who were recruited from dialysis clinics more likely have kidney diseases in the MYH9-associated spectrum of focal segmental glomerulosclerosis. PROFILE and NIH participants, recruited from rheumatology practices, demonstrate that MYH9 does not contribute substantially to the development of LN in AA.
African Americans; Genetics; Lupus nephritis; Kidney; MYH9; Systemic lupus erythematosus
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.
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
ITGAM was recently found to be associated with systemic lupus erythematosus (SLE) in populations of not only European ancestry, but also in Hispanic- and African-Americans, Mexicans and Colombians. The risk alleles in the gene, however, were found to be monomorphic in two Asian populations examined: Japanese and Korean. In this study, using a collection of 910 SLE patients and 2360 controls from Chinese living in Hong Kong, analyzed by both genome-wide association and direct sequencing, we confirmed the association of the same risk alleles in ITGAM with the disease. These findings were further replicated in the Thai population with 278 patients and 383 ethnicity- and geography-matched controls. Subphenotype stratification analyses showed significantly more involvement of the gene in patients with renal nephritis and neurological disorders. Although our results support a pivotal role by rs1143679 (R77H) in disease association, our data also suggests an additional contribution from rs1143683, another non-synonymous polymorphism in this gene (A858V). Therefore, despite the low-allele frequencies of the risk alleles of the gene in our two Asian populations, ITGAM was confirmed to be a risk factor related to disease susceptibility and probably severe manifestations of SLE.
Systemic Lupus Erythematosus (SLE) shows a spectrum of clinical manifestations that complicate its diagnosis, treatment and research. This variability is likely related with environmental exposures and genetic factors among which known SLE susceptibility loci are prime candidates. The first published analyses seem to indicate that this is the case for some of them, but results are still inconclusive and we aimed to further explore this question.
European SLE patients, 1444, recruited at 17 centres from 10 countries were analyzed. Genotypes for 26 SLE associated SNPs were compared between patients with and without each of 11 clinical features: ten of the American College of Rheumatology (ACR) classification criteria (except ANAs) and age of disease onset. These analyses were adjusted for centre of recruitment, top ancestry informative markers, gender and time of follow-up. Overlap of samples with previous studies was excluded for assessing replication.
There were three new associations: the SNPs in XKR6 and in FAM167A-BLK were associated with lupus nephritis (OR = 0.76 and 1.30, Pcorr = 0.007 and 0.03, respectively) and the SNP of MECP2, which is in chromosome X, with earlier age of disease onset in men. The previously reported association of STAT4 with early age of disease onset was replicated. Some other results were suggestive of the presence of additional associations. Together, the association signals provided support to some previous findings and to the characterization of lupus nephritis, autoantibodies and age of disease onset as the clinical features more associated with SLE loci.
Some of the SLE loci shape the disease phenotype in addition to increase susceptibility to SLE. This influence is more prominent for some clinical features than for others. However, results are only partially consistent between studies and subphenotype specific GWAS are needed to unravel their genetic component.
Genome-wide association studies (GWAS) have been shown to be a powerful way of identifying novel susceptibility genes in systemic lupus erythematosus (SLE), as demonstrated by a series of publications in the past year. Lupus has been a late-comer to the GWAS community, being preceded by success stories for the GWAS approach in other autoimmune diseases, including type I diabetes, ankylosing spondylitis, rheumatoid arthritis, Crohn's disease and ulcerative colitis. The paper by Suarez-Gestal and colleagues seeks to exploit the wealth of data available from a total of four GWAS in SLE, three in European-American populations and one in a Swedish population. The authors describe replication of ten lupus susceptibility alleles in a Spanish SLE case-control study.
African Americans, East Asians, and Hispanics with systemic lupus erythematosus (SLE) are more likely to develop renal disease than SLE patients of European descent. We investigated whether European genetic ancestry protects against the development of lupus nephritis and explored genetic and socioeconomic factors that might explain this effect.
This was a cross-sectional study of 1906 adults with SLE. Participants were genotyped for 126 single nucleotide polymorphisms (SNPs) informative for ancestry. A subset of participants was also genotyped for 80 SNPs in 14 candidate genes for renal disease in SLE. We used logistic regression to test the association between European ancestry and renal disease. Analyses adjusted for continental ancestries, socioeconomic status, and candidate genes.
Participants (n=1906) had on average 62.4% European, 15.8% African, 11.5% East Asian, 6.5% Amerindian, and 3.8% South Asian ancestry. Among participants, 34% (n=656) had renal disease. A 10% increase in European ancestry was associated with a 15% reduction in the odds of having renal disease after adjustment for disease duration and sex (OR 0.85, 95% CI 0.82-0.87, p=1.9 × 10−30). Adjusting for other genetic ancestries, measures of socioeconomic status, or SNPs in genes most associated with renal disease (IRF5 (rs4728142), BLK (rs2736340), STAT4 (rs3024912), ITGAM (rs9937837) and HLA-DRB1*0301 and DRB1*1501, p<0.05) did not substantively alter this relationship.
European ancestry is protective against the development of renal disease in SLE, an effect independent of other genetic ancestries, common risk alleles, and socioeconomic status.
In systemic lupus erythematosus (SLE) self-reactive antibodies can target the kidney (lupus nephritis) leading to functional failure and possible mortality. We report that activation of basophils by autoreactive IgE, causes their homing to lymph nodes, promoting TH2 cell differentiation, and enhancing the production of self-reactive antibodies that cause lupus-like nephritis in Lyn−/− mice. SLE patients also have elevated serum IgE, self-reactive IgE's, and activated basophils that express CD62L and the MHC Class II molecule, HLA-DR; parameters that were found to be associated with increased disease activity and active lupus nephritis. Basophils were also present in the lymph nodes and spleen of SLE patients. Thus, in Lyn−/− mice, basophils and IgE autoantibodies amplify autoantibody production that leads to lupus nephritis, and in SLE patients, the presence of IgE autoantibodies and activated basophils are factors associated with disease activity and nephritis.
A recent study in the North American White population has documented the association of a common STAT4 haplotype (tagged by rs7574865) with risk for rheumatoid arthritis (RA) and systemic lupus erythematosus. To replicate this finding in the Korean population, we performed a case-control association study. We genotyped 67 single nucleotide polymorphisms (SNPs) within the STAT1 and STAT4 regions in 1123 Korean patients with RA and 1008 ethnicity-matched controls. The most significant four risk SNPs (rs11889341, rs7574865, rs8179673, and rs10181656 located within the third intron of STAT4) among 67 SNPs are identical with those in the North American study. All four SNPs have modest risk for RA susceptibility (odds ratio 1.21–1.27). A common haplotype defined by these markers (TTCG) carries significant risk for RA in Koreans [34 percent versus 28 percent, P = 0.0027, OR (95 percent CI) = 1.33 (1.10–1.60)]. By logistic regression analysis, this haplotype is an independent risk factor in addition to the classical shared epitope alleles at the HLA-DRB1 locus. There were no significant associations with age of disease onset, radiographic progression, or serologic status using either allelic or haplotypic analysis. Unlike several other risk genes for RA such as PTPN22, PADI4, and FCRL3, a haplotype of the STAT4 gene shows consistent association with RA susceptibility across Whites and Asians, suggesting that this risk haplotype predates the divergence of the major racial groups.
Lyn, a Src-family protein tyrosine kinase expressed in B lymphocytes, contributes to initiation of BCR signaling and is also responsible for feedback inhibition of BCR signaling. Lyn-deficient mice have a decreased number of follicular B cells and also spontaneously develop a lupus-like autoimmunity. We used flow cytometric analysis, BrdU labeling, and our mathematical models of B cell population dynamics, to analyze how Lyn deficiency impacts B-cell maturation and survival. We found that Lyn-deficient transitional 1 (T1) cells develop normally, but T2 cells develop primarily from the T1 subset in the spleen and fail to also develop directly from bone marrow immature B cells. Lyn-deficient T2 cells either mature to the follicular B-cell type at a close to normal rate, or die in this compartment rather than access the T3 anergic subset. The ~40% of wild type follicular cells that were short-lived exited primarily by joining the T3 anergic subset, whereas the ~15% Lyn−/− follicular cells that were not long-lived had a high death rate and died in this compartment rather than entering the T3 subset. We hypothesize that exaggerated BCR signaling resulting from weak interactions with self-antigens is largely responsible for these alterations in Lyn-deficient B cells.
B lymphocytes; Lyn; BrdU; mathematical modeling
UBE2L3 is associated with susceptibility to systemic lupus erythematosus (SLE) and rheumatoid arthritis in European ancestry populations, and this locus has not been investigated fully in non-European populations. We studied the UBE2L3 risk allele for association with SLE, interferon-α (IFN-α), and autoantibodies in a predominantly African American SLE cohort.
We studied 395 patients with SLE and 344 controls. The UBE2L3 rs5754217 polymorphism was genotyped using Taqman primer-probe sets, and IFN-α was measured using a reporter cell assay.
The UBE2L3 rs5754217 T allele was strongly enriched in African American patients with anti-La antibodies as compared to controls, and a recessive model was the best fit for this association (OR 2.55, p = 0.0061). Serum IFN-α also demonstrated a recessive association with the rs5754217 genotype in African American patients, and the TT/anti-La-positive patients formed a significantly high IFN-α subgroup (p = 0.0040). Similar nonstatistically significant patterns of association were observed in the European American patients with SLE. Case-control analysis did not show large allele frequency differences, supporting the idea that this allele is most strongly associated with anti-La-positive patients.
This pattern of recessive influence within a subgroup of patients may explain why this allele does not produce a strong signal in standard case-control studies, and subphenotypes should be included in future studies of UBE2L3. The interaction we observed between UBE2L3 genotype and autoantibodies upon serum IFN-α suggests a biological role for this locus in patients with SLE in vivo.
SYSTEMIC LUPUS ERYTHEMATOSUS; GENETICS INTERFERON-α; AUTOANTIBODIES; UBE2L3 GENOTYPE
Interferon-α (IFNα) is a heritable risk factor for systemic lupus erythematosus (SLE). Genetic variation near IRF7 is implicated in SLE susceptibility. SLE-associated autoantibodies can stimulate IFNα production through the Toll-like receptor/IRF7 pathway. This study was undertaken to determine whether variants of IRF7 act as risk factors for SLE by increasing IFNα production and whether autoantibodies are important to this phenomenon.
We studied 492 patients with SLE (236 African American, 162 European American, and 94 Hispanic American subjects). Serum levels of IFNα were measured using a reporter cell assay, and single-nucleotide polymorphisms (SNPs) in the IRF7/PHRF1 locus were genotyped.
In a joint analysis of European American and Hispanic American subjects, the rs702966 C allele was associated with the presence of anti–double-stranded DNA (anti-dsDNA) antibodies (odds ratio [OR] 1.83, P = 0.0069). The rs702966 CC genotype was only associated with higher serum levels of IFNα in European American and Hispanic American patients with anti-dsDNA antibodies (joint analysis P = 4.1 × 10−5 in anti-dsDNA–positive patients and P = 0.99 in anti-dsDNA–negative patients). In African American subjects, anti-Sm antibodies were associated with the rs4963128 SNP near IRF7 (OR 1.95, P = 0.0017). The rs4963128 CT and TT genotypes were associated with higher serum levels of IFNα only in African American patients with anti-Sm antibodies (P = 0.0012). In African American patients lacking anti-Sm antibodies, an effect of anti-dsDNA–rs702966 C allele interaction on serum levels of IFNα was observed, similar to the other patient groups (overall joint analysis P = 1.0 × 10−6). In European American and Hispanic American patients, the IRF5 SLE risk haplotype showed an additive effect with the rs702966 C allele on IFNα level in anti-dsDNA–positive patients.
Our findings indicate that IRF7/PHRF1 variants in combination with SLE-associated autoantibodies result in higher serum levels of IFNα, providing a biologic relevance for this locus at the protein level in human SLE in vivo.
Previous genome wide association study conducted in a population of European ancestry identified rs4963128, a KIAA1542 SNP 23kb telomeric to IRF7, in strong association with SLE. This study was undertaken to investigate whether genetic polymorphism within IRF7 is a risk factor for the development of SLE.
We genotyped one KIAA1542 SNP rs4963128 and one IRF7 SNP rs1131665 (Q412R) in an Asian population (cases vs. controls: 1302 vs.1479) to assess their association with SLE using custom-designed Beadstation Infinium II platform (Illumina). Subsequently, rs1131665 was further genotyped in independent panels of Chinese (528 vs.527), European American (EA) (446 vs.461) and African American (AA) (159 vs.115) by Taqman genotyping assay to seek confirmation of association in various ethnic groups. Luciferase reporter assay was used to assess the effect of Q412R polymorphism on the activation of IRF7.
Consistent association of rs1131665 (Q412R) with SLE was identified in Asian, EA and AA populations (case vs. control: 2435 vs. 2582; Pmeta = 6.18×10−6, OR = 1.42[1.22–1.65]). Expression of IRF7 412Q risk allele resulted in a 2-fold increase in ISRE transcriptional activity compared with expression of IRF7 412R (P = 0.0003), suggesting IRF7 412Q confers elevated IRF7 activity and may therefore affect downstream IFN pathway.
We showed that the major allele of a nonsynonymous SNP rs1131665 (412Q) in IRF7 confers elevated IRF7 activation and predisposes to the development of SLE in multiple ethnic groups. This result provides direct genetic evidence supporting IRF7 may be a risk gene for human SLE.
Systemic lupus erythematosus (SLE) is a chronic multisystem genetically complex autoimmune disease characterised by the production of autoantibodies to nuclear and cellular antigens, tissue inflammation and organ damage. Genome-wide association studies have shown that variants within the major histocompatibility complex (MHC) region on chromosome 6 confer the greatest genetic risk for SLE in European and Chinese populations. However, the causal variants remain elusive due to tight linkage disequilibrium across disease-associated MHC haplotypes, the highly polymorphic nature of many MHC genes and the heterogeneity of the SLE phenotype.
A high-density case-control single nucleotide polymorphism (SNP) study of the MHC region was undertaken in SLE cohorts of Spanish and Filipino ancestry using a custom Illumina chip in order to fine-map association signals in these haplotypically diverse populations. In addition, comparative analyses were performed between these two datasets and a northern European UK SLE cohort. A total of 1433 cases and 1458 matched controls were examined.
Using this transancestral SNP mapping approach, novel independent loci were identified within the MHC region in UK, Spanish and Filipino patients with SLE with some evidence of interaction. These loci include HLA-DPB1, HLA-G and MSH5 which are independent of each other and HLA-DRB1 alleles. Furthermore, the established SLE-associated HLA-DRB1*15 signal was refined to an interval encompassing HLA-DRB1 and HLA-DQA1. Increased frequencies of MHC region risk alleles and haplotypes were found in the Filipino population compared with Europeans, suggesting that the greater disease burden in non-European SLE may be due in part to this phenomenon.
These data highlight the usefulness of mapping disease susceptibility loci using a transancestral approach, particularly in a region as complex as the MHC, and offer a springboard for further fine-mapping, resequencing and transcriptomic analysis.