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.

Systemic lupus erythematosus (SLE) is an inflammatory autoimmune disease with a strong genetic component. African-Americans (AA) are at increased risk of SLE, but the genetic basis of this risk is largely unknown. To identify causal variants in SLE loci in AA, we performed admixture mapping followed by fine mapping in AA and European-Americans (EA). Through genome-wide admixture mapping in AA, we identified a strong SLE susceptibility locus at 2q22–24 (LOD = 6.28), and the admixture signal is associated with the European ancestry (ancestry risk ratio ∼1.5). Large-scale genotypic analysis on 19,726 individuals of African and European ancestry revealed three independently associated variants in the IFIH1 gene: an intronic variant, rs13023380 [Pmeta = 5.20×10−14; odds ratio, 95% confidence interval = 0.82 (0.78–0.87)], and two missense variants, rs1990760 (Ala946Thr) [Pmeta = 3.08×10−7; 0.88 (0.84–0.93)] and rs10930046 (Arg460His) [Pdom = 1.16×10−8; 0.70 (0.62–0.79)]. Both missense variants produced dramatic phenotypic changes in apoptosis and inflammation-related gene expression. We experimentally validated function of the intronic SNP by DNA electrophoresis, protein identification, and in vitro protein binding assays. DNA carrying the intronic risk allele rs13023380 showed reduced binding efficiency to a cellular protein complex including nucleolin and lupus autoantigen Ku70/80, and showed reduced transcriptional activity in vivo. Thus, in SLE patients, genetic susceptibility could create a biochemical imbalance that dysregulates nucleolin, Ku70/80, or other nucleic acid regulatory proteins. This could promote antibody hypermutation and auto-antibody generation, further destabilizing the cellular network. Together with molecular modeling, our results establish a distinct role for IFIH1 in apoptosis, inflammation, and autoantibody production, and explain the molecular basis of these three risk alleles for SLE pathogenesis.

Author Summary

African-Americans (AA) are at increased risk of systemic lupus erythematosus (SLE), but the genetic basis of this risk increase is largely unknown. We used admixture mapping to localize disease-causing genetic variants that differ in frequency across populations. This approach is advantageous for localizing susceptibility genes in recently admixed populations like AA. Our genome-wide admixture scan identified seven admixture signals, and we followed the best signal at 2q22–24 with fine-mapping, imputation-based association analysis and experimental validation. We identified two independent coding variants and a non-coding variant within the IFIH1 gene associated with SLE. Together with molecular modeling, our results establish a distinct role for IFIH1 in apoptosis, inflammation, and autoantibody production, and explain the molecular basis of these three risk alleles for SLE pathogenesis.

Objective

There is increasing evidence of the impact of SLE on employment, but few studies have had sufficient sample size and longitudinal follow-up to estimate the impact of specific manifestations or of increasing disease activity on employment.

Methods

Data derive from the UCSF Lupus Outcomes Study, a longitudinal cohort of 1204 persons with SLE sampled between 2002 and 2009. Of the 1204, 484 were working at baseline and had at least one follow-up interview. We used the Kaplan-Meier method to estimate the time between onset of thrombotic, neuropsychiatric, or musculoskeletal manifestations or of increased disease activity and work loss, and Cox proportional hazards regression to estimate the risk of work loss associated with the onset of specific manifestations, the number of manifestations, and increased activity, with and without adjustment for sociodemographic, employment, and SLE duration.

Results

By four years of follow-up, 57%, 34%, and 38% of those with thrombotic, musculoskeletal, and neuropsychiatric manifestations, respectively, had stopped working as had 42% of those with increased disease activity. On a bivariable basis, the risk of work loss was significantly higher among persons 55–64, and those with increased disease activity and each kind of manifestation. In multivariable analysis, older age, shorter job tenure, thrombotic and musculoskeletal manifestations, greater number of manifestations, and high levels of activity increased the risk of work loss.

Conclusions

Incident thrombosis and musculoskeletal manifestations, multiple manifestations, and increased disease activity are associated with the risk of work loss in SLE.

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.

Most of the recently identified autoimmunity loci are shared among multiple autoimmune diseases. The pattern of genetic association with autoimmune phenotypes varies, suggesting that certain subgroups of autoimmune diseases are likely to share etiological similarities and underlying mechanisms of disease. In this review, we summarize the major findings from recent studies that have sought to refine genotype-phenotype associations in autoimmune disease by identifying both shared and distinct autoimmunity loci. More specifically, we focus on information from recent genome-wide association studies of rheumatoid arthritis, ankylosing spondylitis, celiac disease, multiple sclerosis, systemic lupus erythematosus, type 1 diabetes and inflammatory bowel disease. Additional work in this area is warranted given both the opportunity it provides to elucidate pathogenic mechanisms in autoimmunity and its potential to inform the development of improved diagnostic and therapeutic tools for this group on complex human disorders.

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.

Background

Obesity is associated with high rates of disability in the general population. The nature of the relationship between obesity and disability in lupus, a condition with a high background rate of disability, is unknown.

Methods

Data were from two interviews, 4 years apart, of a longitudinal cohort of individuals with SLE (n=716 women). Body mass index (BMI) was calculated from self-reported height and weight; obesity was classified by usual (BMI≥30) and revised (BMI≥26.8) definitions. Three measures of functioning were examined: SF-36 Physical Function (PF) subscale, Valued Life Activities (VLA) Disability questionnaire, and employment. Multivariate analyses controlled for demographics, SLE duration and disease activity, glucocorticoid use, depression, and comorbidities. Prospective analyses also controlled for baseline function.

Results

At BMI≥30, 27.8% were obese; at BMI≥26.8, 40.6% were obese. Regardless of obesity definition, obese women exhibited poorer baseline function, with decrements ranging from 20%-33%, depending on the functional measure and obesity definition. With BMI≥26.8, adjusted SF-36 PF scores were 4.3 points lower for obese women (p<.0001), VLA difficulty was 0.09 higher (p=.01), and odds of employment were 80% of non-obese women (OR=0.8 [95% CI 0.5, 1.1]). At 4-year follow-up, women who were obese at baseline had poorer function and experienced greater functional declines.

Conclusion

Obesity was associated with clinically significant negative effects on function, both concurrently and prospectively. This negative impact occurred at a lower BMI than is often considered problematic clinically. Because of the high rate of SLE-related disability, addressing preventable risk factors such as obesity may improve long-term SLE outcomes.

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.

The major histocompatibility complex (MHC) class II transactivator gene (CIITA) encodes an important transcription factor regulating genes required for human leukocyte antigen (HLA) class II MHC-restricted antigen presentation. Major histocompatibility complex (MHC) genes, particularly HLA class II, are strongly associated with risk of developing rheumatoid arthritis (RA). Given the strong biological relationship between CIITA and HLA class II genes, a comprehensive investigation of CIITA variation in RA was conducted. This study tested 31 CIITA SNPs in 2542 RA cases and 3690 controls (N = 6232). All individuals were of European ancestry, as determined by ancestry informative genetic markers. No evidence for association between CIITA variation and RA was observed after a correction for multiple testing was applied. This is the largest study to fully characterize common genetic variation in CIITA, including an assessment of haplotypes. Results exclude even a modest role for common CIITA polymorphisms in susceptibility to RA.

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.

Dendritic cells (DCs), known to support immune activation during infections, may also regulate immune homeostasis in resting animals. Here we show that mice lacking A20 specifically in DCs spontaneously exhibited DC activation and expansion of activated T cells. DC-specific epistasis experiments using A20fl/fl Myd88fl/fl Cd11c-Cre compound mice revealed that A20 restricts both MyD88-independent signals, which drive DC and T cell activation, and MyD88-dependent signals, which drive T cell expansion. In addition, A20fl/fl Cd11c-Cre mice spontaneously developed lymphocyte-dependent colitis, sero-negative ankylosing arthritis and enthesitis, conditions stereotypical for human inflammatory bowel disease (IBD). These findings indicate that DCs require A20 to preserve immune quiescence and suggest that A20-dependent DC functions may underlie IBD and IBD-associated arthritides.

Objectives:

We sought to develop and test an interviewer-administered measure of damage in SLE, the Brief Index of Lupus Damage (BILD), for use in epidemiological studies in which administration of the SLICC/ACR Damage Index (SDI) by trained physicians is not possible or feasible. In addition, we compared the BILD to another recently developed patient-reported measure, the Lupus Damage Index Questionnaire (LDIQ), which was designed as a written survey.

Methods:

A sample of 81 patients from two university-affiliated SLE clinics was used to test the criterion validity of the BILD and the LDIQ. A second sample, the Lupus Outcomes Study (LOS, n=728) was used to ascertain the construct validity of the BILD.

Results:

We found good agreement between most BILD items and corresponding SDI items, and moderately high overall Spearman rank correlations for SDI with BILD (0.64 and with LDIQ (0.54). BILD scores were higher among older individuals, those with longer disease duration, and those with higher mean disease activity in the preceding four years. In addition, higher BILD scores were associated with poorer self-rated health and functional status, greater unemployment and work disability, and increased health care utilization.

Conclusions:

We developed and performed a preliminary validation study demonstrating content, criterion and construct validity of a new practical patient-reported instrument of SLE disease damage. Although further studies are needed to examine reliability and to document psychometric properties in other populations, the BILD appears to represent a promising tool for studies of SLE outside the clinical setting.

The major histocompatibility complex (MHC) class II transactivator gene (CIITA) encodes an important transcription factor required for HLA class II MHC-restricted antigen presentation. MHC genes, including the HLA class II DRB1*03:01 allele, are strongly associated with systemic lupus erythematosus (SLE). Recently the rs4774 CIITA missense variant (+1632G/C) was reported to be associated with susceptibility to multiple sclerosis. In the current study, we investigated CIITA, DRB1*03:01 and risk of SLE using a multi-stage analysis. In stage 1, 9 CIITA variants were tested in 658 cases and 1,363 controls (N = 2,021). In stage 2, rs4774 was tested in 684 cases and 2,938 controls (N = 3,622). We also performed a meta-analysis of the pooled 1,342 cases and 4,301 controls (N = 5,643). In stage 1, rs4774*C was associated with SLE (odds ratio [OR] = 1.24, 95% confidence interval [95% CI] = 1.07–1.44, P = 4.2 × 10−3). Similar results were observed in stage 2 (OR = 1.16, 95% CI = 1.02–1.33, P = 8.5×10−3) and the meta-analysis of the combined dataset (OR = 1.20, 95% CI = 1.09–1.33, Pmeta = 2.5×10−4). In all three analyses, the strongest evidence for association between rs4774*C and SLE was present in individuals who carried at least one copy of DRB1*03:01 (Pmeta= 1.9×10−3). Results support a role for CIITA in SLE, which appears to be stronger in the presence of DRB1*03:01.

Objective

The overexpression of interferon (IFN)-inducible genes is a prominent feature of SLE, serves as a marker for active and more severe disease, and is also observed in other autoimmune and inflammatory conditions. The genetic variations responsible for sustained activation of IFN responsive genes are unknown.

Methods

We systematically evaluated association of SLE with a total of 1,754 IFN-pathway related genes, including IFN-inducible genes known to be differentially expressed in SLE patients and their direct regulators. We performed a three-stage design where two cohorts (total n=939 SLE cases, 3,398 controls) were analyzed independently and jointly for association with SLE, and the results were adjusted for the number of comparisons.

Results

A total of 16,137 SNPs passed all quality control filters of which 316 demonstrated replicated association with SLE in both cohorts. Nine variants were further genotyped for confirmation in an average of 1,316 independent SLE cases and 3,215 independent controls. Association with SLE was confirmed for several genes, including the transmembrane receptor CD44 (rs507230, P = 3.98×10−12), cytokine pleiotrophin (PTN) (rs919581, P = 5.38×10−04), the heat-shock DNAJA1 (rs10971259, P = 6.31×10−03), and the nuclear import protein karyopherin alpha 1 (KPNA1) (rs6810306, P = 4.91×10−02).

Conclusion

This study expands the number of candidate genes associated with SLE and highlights the potential of pathway-based approaches for gene discovery. Identification of the causal alleles will help elucidate the molecular mechanisms responsible for activation of the IFN system in SLE.

Objectives

The Sjögren’s International Collaborative Clinical Alliance (SICCA) is an ongoing NIH-funded registry whose cohort ranges from those with symptoms of possible Sjögren’s syndrome (SS) to those with obvious disease. Using this database we examined associations between labial salivary gland (LSG) histopathology and other phenotypic features of SS.

Methods

LSG biopsy specimens from SICCA participants underwent protocol-directed histopathological assessments. Among 1726 LSG specimens exhibiting any pattern of sialadenitis, we compared biopsy diagnoses against concurrent salivary, ocular and serological assessments.

Results

LSG specimens included 61% with focal lymphocytic sialadenitis, (FLS; 66% of which had focus scores [FS] ≥ 1 per 4 mm2) and 38% with non-specific or sclerosing chronic sialadenitis (NS/SCS). FS ≥ 1 was strongly associated with positive serum anti-SS-A/-B, rheumatoid factor and the ocular component of SS, but not with symptoms of dry mouth or eyes. Those with positive anti-SS-A/-B were 9 times more likely to have a FS ≥ 1 (95% CI: 7.4; 11.9) than FS<1 or another pattern, while those with unstimulated whole salivary flow < 0.1 ml/min were only 2 times more likely to have a FS ≥ 1 (95% CI:1.7; 2.8) than FS<1 or another pattern, while controlling for other phenotypic features of SS.

Conclusions

Distinguishing FLS from NS/SCS is essential in assessing LSG biopsies, before determining FS. A diagnosis of FLS with FS ≥ 1 per 4 mm2, as compared to FLS with FS< 1 or with NS/SCS, was strongly associated with the ocular and serological components of SS and reflects SS autoimmunity.

Introduction

Plasmacytoid dendritic cells (pDCs) constitutively express two members of the Toll-like receptor (TLR) family, TLR-9 and TLR-7, through which they can be stimulated to produce high levels of interferon (IFN)-α, a key mediator of the pathogenesis of systemic lupus erythematosus (SLE). Given the known efficacy of hydroxychloroquine (HCQ) in the treatment of SLE, we examined its ability to inhibit such pDC function in vivo.

Methods

Peripheral blood mononuclear cells (PBMCs) from SLE subjects treated or not with HCQ and from healthy controls were stimulated with the TLR-9 agonist, CpG oligodeoxynucleotides (CpG-A ODN)-2216, and the TLR-7 agonist, imiquimod. The proportion of monocytes, B cells, myeloid dendritic cells, pDCs, and natural killer (NK) cells producing IFN-α and tumor necrosis factor alpha (TNF-α) was then analyzed by multiparameter flow cytometry.

Results

After TLR-9/7 stimulation in both SLE and healthy subjects, significant production of IFN-α and TNF-α was only observed in pDCs. TLR-7 and TLR-9 induced IFN-α and TNF-α production by pDCs from subjects with SLE was decreased relative to that found in controls (TLR-9/IFN-α, P < 0.0001; TLR-9/TNF-α P < 0.0001; TLR-7/TNF-α P = 0.01). TLR-9 and TLR-7 induced IFN-α and TNF-α production by pDCs was severely impaired in 36% (TLR-9) and 33% (TLR-7) of SLE subjects. In almost all cases, these subjects were being treated with HCQ (HCQ vs. no HCQ: impaired TLR-9/IFN-α, P = 0.0003; impaired TLR-7/IFN-α, P = 0.07; impaired TLR-9/TNF-α, P < 0.009; impaired TLR-7/TNF-α, P < 0.01).

Conclusions

Treatment with HCQ is associated with impaired ability of pDCs from subjects with SLE to produce IFN-α and TNF-α upon stimulation with TLR-9 and TLR-7 agonists.

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.

Objective

To determine the prevalence of obesity and evaluate how accurately standard anthropometric measures identify obesity among women with SLE.

Methods

Dual-energy x-ray absorptiometry (DEXA), height, weight, and waist and hip circumference were collected from 145 women with SLE. Three anthropometric proxies of obesity (body mass index (BMI) ≥30 kg/m2, waist circumference (WC) ≥88 cm, and waist-hip ratio (WHR) ≥0.85) were compared to a DEXA-based obesity criterion. Correspondence between measures was assessed with Cohen’s kappa. Receiver operating characteristic (ROC) curves determined optimal cut-points for each anthropometric measure, relative to DEXA. Framingham cardiovascular risk scores were compared among women who were classified as not obese by both traditional and revised anthropometric definitions, obese by both definitions, and obese only by the revised definition.

Results

28%, 29%, 41%, and 50% were classified as obese by WC, BMI, WHR, and DEXA, respectively. Correspondence between anthropometric and DEXA-based measures was moderate. Women misclassified by anthropometric measures had less truncal fat and more appendicular lean and fat mass. Cut-points were identified for anthropometric measures to better approximate DEXA estimates of percent body fat: BMI ≥26.8 kg/m2, waist circumference ≥84.75 cm., and waist-hip ratio ≥0.80. Framingham risk scores were significantly higher in women classified as obese by either traditional or revised criteria.

Conclusions

A large percentage of this group of women with SLE was obese. Substantial portions of women were misclassified by anthropometric measures. Utility of revised cut-points compared to traditional cut-points in identifying risk of cardiovascular disease or disability remains to be examined in prospective studies, but results from the Framingham risk score analysis suggest that traditional cut-points exclude a significant number of at-risk women with SLE.

In spite of the well-known clustering of multiple autoimmune disorders in families, analyses of specific shared genes and polymorphisms between systemic lupus erythematosus (SLE) and other autoimmune diseases (ADs) have been limited. Therefore, we comprehensively tested autoimmune variants for association with SLE, aiming to identify pleiotropic genetic associations between these diseases. We compiled a list of 446 non–Major Histocompatibility Complex (MHC) variants identified in genome-wide association studies (GWAS) of populations of European ancestry across 17 ADs. We then tested these variants in our combined Caucasian SLE cohorts of 1,500 cases and 5,706 controls. We tested a subset of these polymorphisms in an independent Caucasian replication cohort of 2,085 SLE cases and 2,854 controls, allowing the computation of a meta-analysis between all cohorts. We have uncovered novel shared SLE loci that passed multiple comparisons adjustment, including the VTCN1 (rs12046117, P = 2.02×10−06) region. We observed that the loci shared among the most ADs include IL23R, OLIG3/TNFAIP3, and IL2RA. Given the lack of a universal autoimmune risk locus outside of the MHC and variable specificities for different diseases, our data suggests partial pleiotropy among ADs. Hierarchical clustering of ADs suggested that the most genetically related ADs appear to be type 1 diabetes with rheumatoid arthritis and Crohn's disease with ulcerative colitis. These findings support a relatively distinct genetic susceptibility for SLE. For many of the shared GWAS autoimmune loci, we found no evidence for association with SLE, including IL23R. Also, several established SLE loci are apparently not associated with other ADs, including the ITGAM-ITGAX and TNFSF4 regions. This study represents the most comprehensive evaluation of shared autoimmune loci to date, supports a relatively distinct non–MHC genetic susceptibility for SLE, provides further evidence for previously and newly identified shared genes in SLE, and highlights the value of studies of potentially pleiotropic genes in autoimmune diseases.

Author Summary

It is well known that multiple autoimmune disorders cluster in families. However, all of the genetic variants that explain this clustering have not been discovered, and the specific genetic variants shared between systemic lupus erythematosus (SLE) and other autoimmune diseases (ADs) are not known. In order to better understand the genetic factors that explain this predisposition to autoimmunity, we performed a comprehensive evaluation of shared autoimmune genetic variants. First we considered results from 17 ADs and compiled a list with 446 significant genetic variants from these studies. We identified some genetic variants extensively shared between ADs, as well as the ADs that share the most variants. The genetic overlap between SLE and other ADs was modest. Next we tested how important all the 446 genetic variants were in our collection with a minimum of 1,500 SLE patients. Among the most significant variants in SLE, the majority had already been identified in previous studies, but we also discovered variants in two important immune genes. In summary, our data identified diseases with common genetic risk factors and novel SLE effects, and this supports a relatively distinct genetic susceptibility for SLE. This study helps delineate the genetic architecture of ADs.

Systemic lupus erythematosus (SLE) is a complex trait characterised by the production of a range of auto-antibodies and a diverse set of clinical phenotypes. Currently, ∼8% of the genetic contribution to SLE in Europeans is known, following publication of several moderate-sized genome-wide (GW) association studies, which identified loci with a strong effect (OR>1.3). In order to identify additional genes contributing to SLE susceptibility, we conducted a replication study in a UK dataset (870 cases, 5,551 controls) of 23 variants that showed moderate-risk for lupus in previous studies. Association analysis in the UK dataset and subsequent meta-analysis with the published data identified five SLE susceptibility genes reaching genome-wide levels of significance (Pcomb<5×10−8): NCF2 (Pcomb = 2.87×10−11), IKZF1 (Pcomb = 2.33×10−9), IRF8 (Pcomb = 1.24×10−8), IFIH1 (Pcomb = 1.63×10−8), and TYK2 (Pcomb = 3.88×10−8). Each of the five new loci identified here can be mapped into interferon signalling pathways, which are known to play a key role in the pathogenesis of SLE. These results increase the number of established susceptibility genes for lupus to ∼30 and validate the importance of using large datasets to confirm associations of loci which moderately increase the risk for disease.

Author Summary

Genome-wide association studies have revolutionised our ability to identify common susceptibility alleles for systemic lupus erythematosus (SLE). In complex diseases such as SLE, where many different genes make a modest contribution to disease susceptibility, it is necessary to perform large-scale association studies to combine results from several datasets, to have sufficient power to identify highly significant novel loci (P<5×10−8). Using a large SLE collection of 870 UK SLE cases and 5,551 UK unaffected individuals, we firstly replicated ten moderate-risk alleles (P<0.05) from a US–Swedish study of 3,273 SLE cases and 12,188 healthy controls. Combining our results with the US-Swedish data identified five new loci, which crossed the level for genome-wide significance: NCF2 (neutrophil cytosolic factor 2), IKZF1 (Ikaros family zinc-finger 1), IRF8 (interferon regulatory factor 8), IFIH1 (interferon-induced helicase C domain-containing protein 1), and TYK2 (tyrosine kinase 2). Each of these five genes regulates a different aspect of the immune response and contributes to the production of type-I and type-II interferons. Although further studies will be required to identify the causal alleles within these loci, the confirmation of five new susceptibility genes for lupus makes a significant step forward in our understanding of the genetic contribution to SLE.

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.

To discover novel RA risk loci, we systematically examined 370 SNPs from 179 independent loci with p<0.001 in a published meta-analysis of RA GWAS of 3,393 cases and 12,462 controls1. We used GRAIL2, a computational method that applies statistical text mining to PubMed abstracts, to score these 179 loci for functional relationships to genes in 16 established RA disease loci1,3-11. We identified 22 loci with a significant degree of functional connectivity. We genotyped 22 representative SNPs in an independent set of 7,957 cases and 11,958 matched controls. Three validate convincingly: CD2/CD58 (rs11586238, p=1×10−6 replication, p=1×10−9 overall), and CD28 (rs1980422, p=5×10−6 replication, p=1×10−9 overall), PRDM1 (rs548234, p=1×10−5 replication, p=2×10−8 overall). An additional four replicate (p<0.0023): TAGAP (rs394581, p=0.0002 replication, p=4×10−7 overall), PTPRC (rs10919563, p=0.0003 replication, p=7×10−7 overall), TRAF6/RAG1 (rs540386, p=0.0008 replication, p=4×10−6 overall), and FCGR2A (rs12746613, p=0.0022 replication, p=2×10−5 overall). Many of these loci are also associated to other immunologic diseases.

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.