The TNFAIP3 locus at 6q23, encoding A20, has been associated with multiple autoimmune diseases (AIDs). Here, we sequence the coding portions of the gene in order to identify contributing causal polymorphisms that may explain some of the observed associations. A collection of 123 individuals from the Multiple Autoimmune Disease Genetics Consortium (MADGC) collection, each with multiple AIDs (mean=2.2 confirmed diagnoses), and 397 unrelated healthy controls were used for initial sequencing. Thirty-two polymorphisms were identified in the sequencing experiments, including 16 novel and 11 coding variants. Association testing in the entire MADGC collection (1,008 Caucasians with one or more AIDs and 770 unaffected family controls) revealed association of a novel intronic insertion/deletion polymorphism with rheumatoid arthritis (OR = 2.48, p = 0.041). Genotyping of the most common coding polymorphism, rs2230926, in the MADGC collection and additional control individuals revealed significant association with Sjögren’s syndrome (OR = 3.38, p = 0.038), Crohn’s disease (OR = 2.25, p = 0.041), psoriasis (OR = 0.037, p = 0.036) and rheumatoid arthritis (OR = 1.9, p = 0.025). Lastly, haplotype and additional testing of polymorphisms revealed that cases were enriched for 5’ and 3’ untranslated region (UTR) variants (one-sided p-value=0.04), but not specifically for common (>2% MAF), rare, exonic, intronic, non-synonymous, or synonymous variants.

Systemic lupus erythematosus (SLE) is a clinically heterogeneous, systemic autoimmune disease characterized by autoantibody formation. Previously published genome-wide association studies (GWAS) have investigated SLE as a single phenotype. Therefore, we conducted a GWAS to identify genetic factors associated with anti–dsDNA autoantibody production, a SLE–related autoantibody with diagnostic and clinical importance. Using two independent datasets, over 400,000 single nucleotide polymorphisms (SNPs) were studied in a total of 1,717 SLE cases and 4,813 healthy controls. Anti–dsDNA autoantibody positive (anti–dsDNA +, n = 811) and anti–dsDNA autoantibody negative (anti–dsDNA –, n = 906) SLE cases were compared to healthy controls and to each other to identify SNPs associated specifically with these SLE subtypes. SNPs in the previously identified SLE susceptibility loci STAT4, IRF5, ITGAM, and the major histocompatibility complex were strongly associated with anti–dsDNA + SLE. Far fewer and weaker associations were observed for anti–dsDNA – SLE. For example, rs7574865 in STAT4 had an OR for anti–dsDNA + SLE of 1.77 (95% CI 1.57–1.99, p = 2.0E-20) compared to an OR for anti–dsDNA – SLE of 1.26 (95% CI 1.12–1.41, p = 2.4E-04), with pheterogeneity<0.0005. SNPs in the SLE susceptibility loci BANK1, KIAA1542, and UBE2L3 showed evidence of association with anti–dsDNA + SLE and were not associated with anti–dsDNA – SLE. In conclusion, we identified differential genetic associations with SLE based on anti–dsDNA autoantibody production. Many previously identified SLE susceptibility loci may confer disease risk through their role in autoantibody production and be more accurately described as autoantibody propensity loci. Lack of strong SNP associations may suggest that other types of genetic variation or non-genetic factors such as environmental exposures have a greater impact on susceptibility to anti–dsDNA – SLE.

Author Summary

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease that can involve virtually any organ system. SLE patients produce antibodies that bind to their own cells and proteins (autoantibodies) which can cause irreversible organ damage. One particular SLE–related autoantibody directed at double-stranded DNA (anti–dsDNA) is associated with kidney involvement and more severe disease. Previous genome-wide association studies (GWAS) in SLE have studied SLE itself, not particular SLE manifestations. Therefore, we conducted this GWAS of anti–dsDNA autoantibody production to identify genetic associations with this clinically important autoantibody. We found that many previously identified SLE–associated genes are more strongly associated with anti–dsDNA autoantibody production than SLE itself, and they may be more accurately described as autoantibody propensity genes. No strong genetic associations were observed for SLE patients who do not produce anti–dsDNA autoantibodies, suggesting that other factors may have more influence in developing this type of SLE. Further investigation of these autoantibody propensity genes may lead to greater insight into the causes of autoantibody production and organ damage in SLE.

Objective

To determine whether genetic substructure in European-derived populations is associated with specific manifestations of systemic lupus erythematosus (SLE), including mucocutaneous phenotypes, autoantibody production, and renal disease.

Methods

SLE patients of European descent (n=1754) from 8 case collections were genotyped for over 1,400 ancestry informative markers that define a north/south gradient of European substructure. Based on these genetic markers, we used the STRUCTURE program to characterize each SLE patient in terms of percent northern (vs. southern) European ancestry. Non-parametric methods, including tests of trend, were used to identify associations between northern European ancestry and specific SLE manifestations.

Results

In multivariate analyses, increasing levels of northern European ancestry were significantly associated with photosensitivity (ptrend=0.0021, OR for highest quartile of northern European ancestry compared to lowest quartile 1.64, 95% CI 1.13–2.35) and discoid rash (ptrend=0.014, ORhigh-low 1.93, 95% CI 0.98–3.83). In contrast, northern European ancestry was protective for anticardiolipin (ptrend=1.6 × 10−4, ORhigh-low 0.46, 95% CI 0.30–0.69) and anti-dsDNA (ptrend=0.017, ORhigh-low 0.67, 95% CI 0.46–0.96) autoantibody production.

Conclusions

This study demonstrates that specific SLE manifestations vary according to northern vs. southern European ancestry. Thus, genetic ancestry may contribute to the clinical heterogeneity and variation in disease outcomes among SLE patients of European descent. Moreover, these results suggest that genetic studies of SLE subphenotypes will need to carefully address issues of population substructure due to genetic ancestry.

Objective

To define the contribution of polymorphisms in genes encoding tumor necrosis factor (TNF), mannose-binding lectin (MBL), and Fcγ receptor IIa (FCGR2A) as well as clinical factors, to the development of pneumonia in patients with systemic lupus erythematosus (SLE).

Methods

We studied 282 SLE patients from a multiethnic cohort. Pneumonia events and clinical risk factors for pneumonia were identified through medical record review. Genotyping was performed for MBL (+223, +230, and +239), TNF (−308, −238, and +488), and FCGR2A (−131H/R) polymorphisms. Univariate analyses were performed to identify clinical and genetic risk factors for pneumonia. Covariates for multivariate analysis included sex, ethnicity, treatment with immunomodulators, and leukopenia.

Results

Forty-two patients (15%) had at least 1 episode of pneumonia. Polymorphism of the TNF gene, particularly the −238A allele and a related haplotype, revealed the most striking and consistent association with pneumonia in univariate analyses. Results of multivariate analyses indicated an odds ratio (OR) for the TNF −238A allele of 3.5 (P = 0.007) and an OR for the related haplotype of 5.4 (P = 0.001). Male sex, treatment with immunomodulators, and leukopenia also influenced the risk of pneumonia.

Conclusion

These findings suggest that specific TNF variants may identify SLE patients who are at particularly high risk of developing pneumonia. Given the prevalence and excessive morbidity associated with pneumonia in SLE, these findings have clinical relevance and provide insight into the pathogenesis.

A substantial genetic contribution to systemic lupus erythematosus (SLE) risk is conferred by major histocompatibility complex (MHC) gene(s) on chromosome 6p21. Previous studies in SLE have lacked statistical power and genetic resolution to fully define MHC influences. We characterized 1,610 Caucasian SLE cases and 1,470 parents for 1,974 MHC SNPs, the highly polymorphic HLA-DRB1 locus, and a panel of ancestry informative markers. Single-marker analyses revealed strong signals for SNPs within several MHC regions, as well as with HLA-DRB1 (global p = 9.99×10−16). The most strongly associated DRB1 alleles were: *0301 (odds ratio, OR = 2.21, p = 2.53×10−12), *1401 (OR = 0.50, p = 0.0002), and *1501 (OR = 1.39, p = 0.0032). The MHC region SNP demonstrating the strongest evidence of association with SLE was rs3117103, with OR = 2.44 and p = 2.80×10−13. Conditional haplotype and stepwise logistic regression analyses identified strong evidence for association between SLE and the extended class I, class I, class III, class II, and the extended class II MHC regions. Sequential removal of SLE–associated DRB1 haplotypes revealed independent effects due to variation within OR2H2 (extended class I, rs362521, p = 0.006), CREBL1 (class III, rs8283, p = 0.01), and DQB2 (class II, rs7769979, p = 0.003, and rs10947345, p = 0.0004). Further, conditional haplotype analyses demonstrated that variation within MICB (class I, rs3828903, p = 0.006) also contributes to SLE risk independent of HLA-DRB1*0301. Our results for the first time delineate with high resolution several MHC regions with independent contributions to SLE risk. We provide a list of candidate variants based on biologic and functional considerations that may be causally related to SLE risk and warrant further investigation.

Author Summary

Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by autoantibody production and involvement of multiple organ systems. Although the cause of SLE remains unknown, several lines of evidence underscore the importance of genetic factors. As is true for most autoimmune diseases, a substantial genetic contribution to disease risk is conferred by major histocompatibility complex (MHC) gene(s) on chromosome 6. This region of the genome contains a large number of genes that participate in the immune response. However, the full contribution of this genomic region to SLE risk has not yet been defined. In the current study we characterize a large number of SLE patients and family members for approximately 2,000 MHC region variants to identify the specific genes that influence disease risk. Our results, for the first time, implicate four different MHC regions in SLE risk. We provide a list of candidate variants based on biologic and functional considerations that may be causally related to SLE risk and warrant further investigation.

Systemic lupus erythematosus (SLE) is a genetically complex disease with heterogeneous clinical manifestations. A polymorphism in the STAT4 gene has recently been established as a risk factor for SLE, but the relationship with specific SLE subphenotypes has not been studied. We studied 137 SNPs in the STAT4 region genotyped in 4 independent SLE case series (total n = 1398) and 2560 healthy controls, along with clinical data for the cases. Using conditional testing, we confirmed the most significant STAT4 haplotype for SLE risk. We then studied a SNP marking this haplotype for association with specific SLE subphenotypes, including autoantibody production, nephritis, arthritis, mucocutaneous manifestations, and age at diagnosis. To prevent possible type-I errors from population stratification, we reanalyzed the data using a subset of subjects determined to be most homogeneous based on principal components analysis of genome-wide data. We confirmed that four SNPs in very high LD (r2 = 0.94 to 0.99) were most strongly associated with SLE, and there was no compelling evidence for additional SLE risk loci in the STAT4 region. SNP rs7574865 marking this haplotype had a minor allele frequency (MAF) = 31.1% in SLE cases compared with 22.5% in controls (OR = 1.56, p = 10−16). This SNP was more strongly associated with SLE characterized by double-stranded DNA autoantibodies (MAF = 35.1%, OR = 1.86, p<10−19), nephritis (MAF = 34.3%, OR = 1.80, p<10−11), and age at diagnosis<30 years (MAF = 33.8%, OR = 1.77, p<10−13). An association with severe nephritis was even more striking (MAF = 39.2%, OR = 2.35, p<10−4 in the homogeneous subset of subjects). In contrast, STAT4 was less strongly associated with oral ulcers, a manifestation associated with milder disease. We conclude that this common polymorphism of STAT4 contributes to the phenotypic heterogeneity of SLE, predisposing specifically to more severe disease.

Author Summary

Systemic lupus erythematosus is a chronic disabling autoimmune disease, most commonly striking women in their thirties or forties. It can cause a wide variety of clinical manifestations, including kidney disease, arthritis, and skin disorders. Prognosis varies greatly depending on these clinical features, with kidney disease and related characteristics leading to greater morbidity and mortality. It is also complex genetically; while lupus runs in families, genes increase one’s risk for lupus but do not fully determine the outcome. It is thought that the interactions of multiple genes and/or interactions between genes and environmental factors may cause lupus, but the causes and disease pathways of this very heterogeneous disease are not well understood. By examining relationships between subtypes of lupus and specific genes, we hope to better understand how lupus is triggered and by what biological pathways it progresses. We show in this work that the STAT4 gene, very recently identified as a lupus risk gene, predisposes specifically to severe manifestations of lupus, including kidney disease.