A classic T-cell phenotype in Systemic lupus erythematosus (SLE) is the downregulation and replacement of the CD3ζ chain that alters TCR signaling. However, genetic associations with SLE in the human CD247 locus that encodes CD3ζ are not well established and require replication in independent cohorts. Our aim was therefore to examine, localize and validate CD247-SLE association in a large multi-ethnic population. We typed 44 contiguous CD247 SNPs in 8 922 SLE patients and 8 077 controls from four ethnically distinct populations. The strongest associations were found in the Asian population (11 SNPs in intron 1, 4.99×10−4
Exploiting genotyping, DNA sequencing, imputation and trans-ancestral mapping, we used Bayesian and frequentist approaches to model the IRF5–TNPO3 locus association, now implicated in two immunotherapies and seven autoimmune diseases. Specifically, in systemic lupus erythematosus (SLE), we resolved separate associations in the IRF5 promoter (all ancestries) and with an extended European haplotype. We captured 3230 IRF5–TNPO3 high-quality, common variants across 5 ethnicities in 8395 SLE cases and 7367 controls. The genetic effect from the IRF5 promoter can be explained by any one of four variants in 5.7 kb (P-valuemeta = 6 × 10−49; OR = 1.38–1.97). The second genetic effect spanned an 85.5-kb, 24-variant haplotype that included the genes IRF5 and TNPO3 (P-valuesEU = 10−27–10−32, OR = 1.7–1.81). Many variants at the IRF5 locus with previously assigned biological function are not members of either final credible set of potential causal variants identified herein. In addition to the known biologically functional variants, we demonstrated that the risk allele of rs4728142, a variant in the promoter among the lowest frequentist probability and highest Bayesian posterior probability, was correlated with IRF5 expression and differentially binds the transcription factor ZBTB3. Our analytical strategy provides a novel framework for future studies aimed at dissecting etiological genetic effects. Finally, both SLE elements of the statistical model appear to operate in Sjögren's syndrome and systemic sclerosis whereas only the IRF5–TNPO3 gene-spanning haplotype is associated with primary biliary cirrhosis, demonstrating the nuance of similarity and difference in autoimmune disease risk mechanisms at IRF5–TNPO3.
Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by complex interplay among immune cell types. SLE activity is experimentally assessed by several blood tests, including gene expression profiling of heterogeneous populations of cells in peripheral blood. To better understand the contribution of different cell types in SLE pathogenesis, we applied the two methods in cell-type-specific differential expression analysis, csSAM and DSection, to identify cell-type-specific gene expression differences in heterogeneous gene expression measures obtained using RNA-seq technology. We identified B-cell-, monocyte-, and neutrophil-specific gene expression differences. Immunoglobulin-coding gene expression was altered in B-cells, while a ribosomal signature was prominent in monocytes. On the contrary, genes differentially expressed in the heterogeneous mixture of cells did not show any functional enrichment. Our results identify antigen binding and structural constituents of ribosomes as functions altered by B-cell- and monocyte-specific gene expression differences, respectively. Finally, these results position both csSAM and DSection methods as viable techniques for cell-type-specific differential expression analysis, which may help uncover pathogenic, cell-type-specific processes in SLE.
cell-type-specific; deconvolution; RNA-seq; SLE; csSAM; DSection
Both genetic and environmental interactions affect systemic lupus erythematosus (SLE) development and pathogenesis. One known genetic factor associated with lupus is a haplotype of the interferon regulatory factor 5 (IRF5) gene. Analysis of global gene expression microarray data using gene set enrichment analysis identified multiple interferon- and inflammation-related gene sets significantly overrepresented in cells with the risk haplotype. Pathway analysis using expressed genes from the significant gene sets impacted by the IRF5 risk haplotype confirmed significant correlation with the interferon pathway, Toll-like receptor pathway, and the B-cell receptor pathway. SLE patients with the IRF5 risk haplotype have a heightened interferon signature, even in an unstimulated state (P = 0.011), while patients with the IRF5 protective haplotype have a B cell interferon signature similar to that of controls. These results identify multiple genes in functionally significant pathways which are affected by IRF5 genotype. They also establish the IRF5 risk haplotype as a key determinant of not only the interferon response, but also other B-cell pathways involved in SLE.
Systemic lupus erythematosus (SLE) is a multifaceted disease characterized by immune dysregulation and unpredictable disease activity. This study evaluated changes in plasma concentrations of soluble mediators preceding clinically-defined disease flares.
Soluble mediators (n=52) were examined, including cytokines, chemokines, and soluble receptors, using validated multiplex bead-based or enzyme-linked immunosorbent assays in plasma from European American SLE patients who developed disease flare 6 or 12 weeks after baseline assessment were compared to 28 matched SLE patients without impending flare and 28 matched healthy controls (n=84). For a subset, mediators within samples preceding SLE disease flare and during a clinically stable period from the same individual were compared.
Compared to clinically stable patients, patients with impending flare had significant (p≤0.01) alterations in 27 soluble mediators at baseline with significantly higher levels of pro-inflammatory mediators, including Th1, Th2, and Th17-type cytokines, several weeks before clinical flare. Baseline levels of regulatory cytokines, including IL-10 and TGF-β, were higher in non-flare SLE patients, while baseline levels of soluble TNFRI, TNFRII, Fas, FasL, and CD40L were significantly greater in pre-flare patients (p≤0.002). A normalized and weighted combined soluble mediator score was significantly higher in pre-flare SLE patients versus those with stable disease (p≤0.0002).
Pro-inflammatory adaptive cytokines and shed TNF receptors, are elevated prior to disease flare, while regulatory mediators are elevated during periods of stable disease. Alterations in the balance between inflammatory and regulatory mediators may help identify patients at risk of disease flare and help decipher SLE pathogenic mechanisms.
SLE; disease flare; cytokines
We targeted LYN, a src-tyosine kinase involved in B cell activation, in case-control association studies using populations of European American, African American and Korean subjects. Our combined European-derived population, consisting of 2463 independent cases and 3131 unrelated controls, demonstrates significant association with rs6983130 in a female-only analysis with 2254 cases and 2228 controls (p=1.1 × 10−4, OR=0.81 (95% CI: 0.73 – 0.90)). This SNP is located in the 5′ UTR within the first intron near the transcription initiation site of LYN. Additional SNPs upstream of the first exon also show weak and sporadic association in subsets of the total European American population. Multivariate logistic regression analysis implicates rs6983130 as a protective factor for SLE susceptibility when anti-dsDNA, anti-chromatin, anti-52 kDa Ro or anti-Sm autoantibody status were used as covariates. Subset analysis of the European American female cases by ACR classification criteria reveals a reduction in the risk of hematologic disorder with rs6983130 compared to cases without hematologic disorders (p=1.5 × 10−3, OR=0.75 (95% C.I.=0.62-0.89)). None of the 90 SNPs tested demonstrate significant association with SLE in the African American or Korean populations. These results support an association of LYN with European-derived individuals with SLE, especially within autoantibody or clinical subsets.
systemic lupus erythematosus; association; LYN; SNP
Systemic lupus erythematosus (SLE) is an autoimmune disease with highly variable clinical presentation. Patients suffer from immunological abnormalities that target T cell, B cell and accessory cell functions. B cells are hyperactive in SLE patients. An adaptor protein expressed in B cells called BANK1 (B-cell scaffold protein with ankyrin repeats) was reported in a previous study to be associated with SLE in a European population. The objective of this study is to assess the BANK1 genotype-phenotype association in an independent replication sample. We genotyped 38 single nucleotide polymorphisms (SNPs) in BANK1 on 1892 European-derived SLE patients and 2652 European-derived controls. The strongest associations with SLE and BANK1 were at rs17266594 (corrected p-value=1.97 × 10−5, OR=1.22, 95% C.I.(1.12–1.34)) and rs10516487 (corrected p-value=2.59 × 10−5, OR=1.22, 95% C.I.(1.11–1.34)). Our findings suggest that the association is explained by these two SNPs, confirming previous reports that these polymorphisms contribute to the risk of developing lupus. Analysis of patient subsets enriched for hematological, immunological and renal ACR criteria or the levels of autoantibodies, such as anti-RNP A and anti-SmRNP, uncovers additional BANK1 associations. Our results suggest that BANK1 polymorphisms alter immune system development and function to increase the risk for developing lupus.
systemic lupus erythematosus; replication; association; European; BANK1
To examine whether smoking is associated with autoantibody production in systemic lupus erythematosus (SLE) patients, unaffected first-degree relatives (FDR) of individuals with SLE - a group at increased risk of developing SLE, or unaffected, unrelated controls.
Detailed demographic, environmental, clinical, and therapeutic information was collected by questionnaire on 1,242 SLE patients, 981 FDRs, and 946 controls in the Lupus Family Registry and Repository; a blood sample was obtained. All sera were tested for multiple lupus autoantibodies by immunofluorescence and luminex bead-based assays. Generalized estimating equations, adjusting for age, gender, and ethnicity and accounting for correlation within families, were used to assess smoking status with the dichotomous outcome variables of positivity for SLE status, positivity of ANA by immunofluorescence (≥ 1:120), positivity for ≥ 1 autoantibody by the luminex assay, and positivity for each of the 11 autoantibodies.
Current smoking was associated with being positive for ≥ 1 autoantibody (excluding ANA) (adjusted OR=1.53, 95% CI 1.04–2.24) in our subjects with SLE. No association was observed in unaffected FDRs or healthy controls. Former smoking was associated with anti-Ro/SS-A60 in our unaffected FDRs. There was an increased association with anti-nRNP A seropositivity, as well as a decreased association with anti-nRNP 68 positivity, in current smokers in SLE subjects.
No clear association between smoking status and individual autoantibodies was detected in SLE patients, unaffected FDRs, nor healthy controls within this collection. The association of smoking with SLE may therefore manifest its risk through mechanisms outside of autoantibody production, at least for the specificities tested.
Smoking; autoantibodies; systemic lupus erythematosus
Recent reports have associated NCF2, encoding a core component of the multi-protein NADPH oxidase (NADPHO), with systemic lupus erythematosus (SLE) susceptibility in individuals of European ancestry. To identify ethnicity-specific and -robust variants within NCF2, we assessed 145 SNPs in and around the NCF2 gene in 5325 cases and 21 866 controls of European-American (EA), African-American (AA), Hispanic (HS) and Korean (KR) ancestry. Subsequent imputation, conditional, haplotype and bioinformatic analyses identified seven potentially functional SLE-predisposing variants. Association with non-synonymous rs17849502, previously reported in EA, was detected in EA, HS and AA (PEA = 1.01 × 10−54, PHS = 3.68 × 10−10, PAA = 0.03); synonymous rs17849501 was similarly significant. These SNPs were monomorphic in KR. Novel associations were detected with coding variants at rs35937854 in AA (PAA = 1.49 × 10−9), and rs13306575 in HS and KR (PHS = 7.04 × 10−7, PKR = 3.30 × 10−3). In KR, a 3-SNP haplotype was significantly associated (P = 4.20 × 10−7), implying that SLE predisposing variants were tagged. Significant SNP–SNP interaction (P = 0.02) was detected between rs13306575 and rs17849502 in HS, and a dramatically increased risk (OR = 6.55) with a risk allele at each locus. Molecular modeling predicts that these non-synonymous mutations could disrupt NADPHO complex assembly. The risk allele of rs17849501, located in a conserved transcriptional regulatory region, increased reporter gene activity, suggesting in vivo enhancer function. Our results not only establish allelic heterogeneity within NCF2 associated with SLE, but also emphasize the utility of multi-ethnic cohorts to identify predisposing variants explaining additional phenotypic variance (‘missing heritability’) of complex diseases like SLE.
Genome wide association studies have identified variants in PXK that confer risk for humoral autoimmune diseases, including systemic lupus erythematosus (SLE or lupus), rheumatoid arthritis and more recently systemic sclerosis. While PXK is involved in trafficking of epidermal growth factor Receptor (EGFR) in COS-7 cells, mechanisms linking PXK to lupus pathophysiology have remained undefined. In an effort to uncover the mechanism at this locus that increases lupus-risk, we undertook a fine-mapping analysis in a large multi-ancestral study of lupus patients and controls. We define a large (257kb) common haplotype marking a single causal variant that confers lupus risk detected only in European ancestral populations and spans the promoter through the 3′ UTR of PXK. The strongest association was found at rs6445972 with P < 4.62 × 10−10, OR 0.81 (0.75–0.86). Using stepwise logistic regression analysis, we demonstrate that one signal drives the genetic association in the region. Bayesian analysis confirms our results, identifying a 95% credible set consisting of 172 variants spanning 202 kb. Functionally, we found that PXK operates on the B-cell antigen receptor (BCR); we confirmed that PXK influenced the rate of BCR internalization. Furthermore, we demonstrate that individuals carrying the risk haplotype exhibited a decreased rate of BCR internalization, a process known to impact B cell survival and cell fate. Taken together, these data define a new candidate mechanism for the genetic association of variants around PXK with lupus risk and highlight the regulation of intracellular trafficking as a genetically regulated pathway mediating human autoimmunity.
lupus; PXK; fine-mapping; B cells; BCR
To explore the potential influence of the polymorphic 8p23.1 inversion on known autoimmune susceptibility risk at or near BLK locus, we validated a new bioinformatics method that utilizes SNP data to enable accurate, high-throughput genotyping of the 8p23.1 inversion in a Caucasian population. Methods: Principal components analysis (PCA) was performed using markers inside the inversion territory followed by k-means cluster analyses on 7416 European derived and 267 HapMaP CEU and TSI samples. A logistic regression conditional analysis was performed. Results: Three subgroups have been identified; inversion homozygous, heterozygous and non-inversion homozygous. The status of inversion was further validated using HapMap samples that had previously undergone Fluorescence in situ hybridization (FISH) assays with a concordance rate of above 98%. Conditional analyses based on the status of inversion were performed. We found that overall association signals in the BLK region remain significant after controlling for inversion status. The proportion of lupus cases and controls (cases/controls) in each subgroup was determined to be 0.97 for the inverted homozygous group (1067 cases and 1095 controls), 1.12 for the inverted heterozygous group (1935 cases 1717 controls) and 1.36 for non-inverted subgroups (924 cases and 678 controls). After calculating the linkage disequilibrium between inversion status and lupus risk haplotype we found that the lupus risk haplotype tends to reside on non-inversion background. As a result, a new association effect between non-inversion status and lupus phenotype has been identified ((p = 8.18×10−7, OR = 1.18, 95%CI = 1.10–1.26). Conclusion: Our results demonstrate that both known lupus risk haplotype and inversion status act additively in the pathogenesis of lupus. Since inversion regulates expression of many genes in its territory, altered expression of other genes might also be involved in the development of lupus.
Patients with systemic lupus erythematosus (SLE) are at increased risk of herpes zoster (HZ). Although a vaccine for HZ has been FDA approved, its use in immunocompromised individuals remains controversial because it is a live-attenuated virus vaccine. We performed a pilot study of the immunogenicity of Zostavax® in SLE patients.
Ten SLE patients and 10 controls ≥50 years old participated in this open label vaccination study. All were seropositive for varicella zoster virus (VZV). SLE patients were excluded for SLEDAI>4, use of mycophenolatemofetil, cyclophosphamide, biologics, or >10 mg prednisone daily. Follow-up visits occurred at 2, 6, and 12 weeks. Clinical outcomes included the development of adverse events, particularly HZ or vesicular lesions, and SLE flare. Immunogenicity was assessed with VZV-specific IFN-γ producing ELISPOT assays and with antibody concentrations.
All subjects were women. SLE patients were slightly older than controls (60.5 vs. 55.3 years, p<0.05) Median baseline SLEDAI was 0 (range 0–2) for SLE patients. No episodes of HZ, vesicular rash, serious adverse events, or SLE flares occurred. Three injection site reactions occurred in each group: mild erythema or tenderness. The proportion of subjects with a >50% increase in ELISPOT results following vaccination was comparable between both groups, although absolute SLE responses were lower than controls. Antibody titers increased only among controls following vaccination (p<0.05).
Zostavax vaccination yielded a measurable immuneresponse in this cohort of mild SLE patients on mild-moderate immunosuppressive medications. No herpetiform lesions or lupus flares were seen in this small cohort of patients.
Systemic lupus erythematosus; herpes zoster; vaccine; Zostavax; infection; clinical trial
Sjögren’s syndrome is a common autoimmune disease (~0.7% of European Americans) typically presenting as keratoconjunctivitis sicca and xerostomia. In addition to strong association within the HLA region at 6p21 (Pmeta=7.65×10−114), we establish associations with IRF5-TNPO3 (Pmeta=2.73×10−19), STAT4 (Pmeta=6.80×10−15), IL12A (Pmeta =1.17×10−10), FAM167A-BLK (Pmeta=4.97×10−10), DDX6-CXCR5 (Pmeta=1.10×10−8), and TNIP1 (Pmeta=3.30×10−8). Suggestive associations with Pmeta<5×10−5 were observed with 29 regions including TNFAIP3, PTTG1, PRDM1, DGKQ, FCGR2A, IRAK1BP1, ITSN2, and PHIP amongst others. These results highlight the importance of genes involved in both innate and adaptive immunity in Sjögren’s syndrome.
In recent years, vitamin D has been shown to possess a wide range of immunomodulatory effects. Although there is extensive amount of research on vitamin D, we lack a comprehensive understanding of the prevalence of vitamin D deficiency or the mechanism by which vitamin D regulates the human immune system. This study examined the prevalence and correlates of vitamin D deficiency and the relationship between vitamin D and the immune system in healthy individuals.
Healthy individuals (n = 774) comprised of European-Americans (EA, n = 470), African–Americans (AA, n = 125), and Native Americans (NA, n = 179) were screened for 25-hydroxyvitamin D [25(OH)D] levels by ELISA. To identify the most noticeable effects of vitamin D on the immune system, 20 EA individuals with severely deficient (<11.3 ng/mL) and sufficient (>24.8 ng/mL) vitamin D levels were matched and selected for further analysis. Serum cytokine level measurement, immune cell phenotyping, and phosphoflow cytometry were performed.
Vitamin D sufficiency was observed in 37.5% of the study cohort. By multivariate analysis, AA, NA, and females with a high body mass index (BMI, >30) demonstrate higher rates of vitamin D deficiency (p<0.05). Individuals with vitamin D deficiency had significantly higher levels of serum GM-CSF (p = 0.04), decreased circulating activated CD4+ (p = 0.04) and CD8+ T (p = 0.04) cell frequencies than individuals with sufficient vitamin D levels.
A large portion of healthy individuals have vitamin D deficiency. These individuals have altered T and B cell responses, indicating that the absence of sufficient vitamin D levels could result in undesirable cellular and molecular alterations ultimately contributing to immune dysregulation.
Genetic polymorphism in MECP2/IRAK1 on chromosome Xq28 is a confirmed and replicated susceptibility locus for lupus. High linkage disequilibrium in this locus suggests that both MECP2 and IRAK1 are candidate genes for the disease. DNA methylation changes in lupus T cells play a central role in the pathogenesis of lupus, and MeCp-2 (encoded by MECP2) is a master regulator of gene expression and is also known to recruit DNA methyltransferase 1 (DNMT1) during DNA synthesis. Using human T cells from normal individuals with either the lupus risk or the lupus protective haplotype in MECP2/IRAK1, we demonstrate that polymorphism in this locus increases MECP2 isoform 2 mRNA expression in stimulated but not unstimulated T cells. By assessing DNA methylation levels across over 485,000 methylation sites across the entire genome, we further demonstrate that the lupus risk variant in this locus is associated with significant DNA methylation changes, including in the HLA-DR and HLA-DQ loci, as well as interferon-related genes such as IFI6, IRF6, and BST2. Further, using a human MECP2 transgenic mouse, we show that overexpression of MECP2 alters gene expression in stimulated T cells. This includes overexpression of Eif2c2 that regulates the expression of multiple microRNAs (such as miR-21), and the histone demethylase Jhdm1d. In addition, we show that MECP2 transgenic mice develop antinuclear antibodies. Our data suggest that the lupus associated variant in the MECP2/IRAK1 locus has the potential to affect all 3 epigenetic mechanisms: DNA methylation, microRNA expression, and histone modification. Importantly, these data support the notion that variants within the MECP2 gene can alter DNA methylation in other genetic loci including the HLA and interferon-regulated genes, thereby providing evidence for genetic-epigenetic interaction in lupus.
MECP2; IRAK1; lupus; epigenetics; polymorphism; DNA methylation; T cells; transgenic mouse
Analysis of the ImmunoChip single nucleotide polymorphism (SNP) array in 2816 individuals, comprising the most common subtypes (oligoarticular and RF negative polyarticular) of juvenile idiopathic arthritis (JIA) and 13056 controls strengthens the evidence for association to three known JIA-risk loci (HLA, PTPN22 and PTPN2) and has identified fourteen risk loci reaching genome-wide significance (p < 5 × 10-8) for the first time. Eleven additional novel regions showed suggestive evidence for association with JIA (p < 1 × 10-6). Dense-mapping of loci along with bioinformatic analysis has refined the association to one gene for eight regions, highlighting crucial pathways, including the IL-2 pathway, in JIA disease pathogenesis. The entire ImmunoChip loci, HLA region and the top 27 loci (p < 1 × 10-6) explain an estimated 18%, 13% and 6% risk of JIA, respectively. Analysis of the ImmunoChip dataset, the largest cohort of JIA cases investigated to date, provides new insight in understanding the genetic basis for this childhood autoimmune disease.
Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by autoantibody production and altered type I interferon expression. Genetic surveys and genome-wide association studies have identified more than 30 SLE susceptibility genes. One of these genes, TNIP1, encodes the ABIN1 protein. ABIN1 functions in the immune system by restricting the NF-κB signaling. In order to better understand the genetic factors that influence association with SLE in genes that regulate the NF-κB pathway, we analyzed a dense set of genetic markers spanning TNIP1 and TAX1BP1, as well as the TNIP1 homolog, TNIP2, in case-control sets of diverse ethnic origins.
We fine-mapped TNIP1, TNIP2, and TAX1BP1 in a total of 8372 SLE cases and 7492 healthy controls from European-ancestry, African-American, Hispanic, East Asian, and African-American Gullah populations. Levels of TNIP1 mRNA and ABIN1 protein were analyzed using quantitative RT-PCR and Western blotting, respectively, in EBV-transformed human B cell lines.
We found significant associations between genetic variants within TNIP1 and SLE but not in TNIP2 or TAX1BP1. After resequencing and imputation, we identified two independent risk haplotypes within TNIP1 in individuals of European-ancestry that were also present in African-American and Hispanic populations. These risk haplotypes produced lower levels of TNIP1 mRNA and ABIN1 protein suggesting they harbor hypomorphic functional variants that influence susceptibility to SLE by restricting ABIN1 expression.
Our results confirmed the association signals between SLE and TNIP1 variants in multiple populations and provide new insight into the mechanism by which TNIP1 variants may contribute to SLE pathogenesis.
The replacement of standard immunofluorescence anti-nuclear antibody (ANA) methods with bead-based assays is a new clinical option. A large, multi-racial cohort of SLE patients, blood relatives and unaffected control individuals was evaluated for familial aggregation and subset clustering of autoantibodies by high-throughput serum screening technology and traditional methods.
Serum samples (1,540 SLE patients, 1,127 unaffected relatives, and 906 healthy, population-based controls) were analyzed for SLE autoantibodies using a bead-based assay, immunofluorescence, and immunodiffusion. Autoantibody prevalence, disease sensitivity, clustering, and association with standard immunodiffusion results were evaluated.
ANA frequency in SLE patient sera were 89%, 73%, and 67% by BioPlex 2200 and 94%, 84%, and 86% by immunofluorescence in African-American, Hispanic, and European-American patients respectively. 60kD Ro, La, Sm, nRNP A, and ribosomal P prevalence were compared across assays, with sensitivities ranging from 0.92 to 0.83 and specificities ranging from 0.90 to 0.79. Cluster autoantibody analysis showed association of three subsets: 1) 60kD Ro, 52kD Ro and La, 2) spliceosomal proteins, and 3) dsDNA, chromatin, and ribosomal P. Familial aggregation of Sm/RNP, ribosomal P, and 60kD Ro in SLE patient sibling pairs was observed (p ≤ 0.004). Simplex pedigree patients had a greater prevalence for dsDNA (p=0.0003) and chromatin (p=0.005) autoantibodies than multiplex patients.
ANA frequencies detected by a bead-based assay are lower in European-American SLE patients compared to immunofluorescence. These assays have strong positive predictive values across racial groups, provide useful information for clinical care, and provide unique insights to familial aggregation and autoantibody clustering.
systemic lupus erythematosus; autoantibodies; ancestry
Recent application of gene expression profiling to the immune system has shown a great potential for characterization of complex regulatory processes. It is becoming increasingly important to characterize functional systems through multigene interactions to provide valuable insights into differences between healthy controls and autoimmune patients. Here we apply an original systematic approach to the analysis of changes in regulatory gene interconnections between in Epstein-Barr virus transformed hyperresponsive B cells from SLE patients and normal control B cells. Both traditional analysis of differential gene expression and analysis of the dynamics of gene expression variations were performed in combination to establish model networks of functional gene expression. This Pathway Dysregulation Analysis identified known transcription factors and transcriptional regulators activated uniquely in stimulated B cells from SLE patients.
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.
Little is known about the genetic etiology of systemic lupus erythematosus (SLE) in individuals of African ancestry, despite its higher prevalence and greater disease severity. Overproduction of nitric oxide (NO) and reactive oxygen species are implicated in the pathogenesis and severity of SLE, making NO synthases and other reactive intermediate related genes biological candidates for disease susceptibility. This study analyzed variation in reactive intermediate genes for association with SLE in two populations with African ancestry.
A total of 244 SNPs from 53 regions were analyzed in non-Gullah African Americans (AA; 1432 cases and 1687 controls) and the genetically more homogeneous Gullah of the Sea Islands of South Carolina (133 cases and 112 controls) and. Single-marker, haplotype, and two-locus interaction tests were computed for these populations.
The glutathione reductase gene GSR (rs2253409, P=0.0014, OR [95% CI]=1.26 [1.09–1.44]) was the most significant single-SNP association in AA. In the Gullah, the NADH dehydrogenase NDUFS4 (rs381575, P=0.0065, OR [95%CI]=2.10 [1.23–3.59]) and nitric oxide synthase gene NOS1 (rs561712, P=0.0072, OR [95%CI]=0.62 [0.44–0.88]) were most strongly associated with SLE. When both populations were analyzed together, GSR remained the most significant effect (rs2253409, P=0.00072, OR [95%CI]=1.26 [1.10–1.44]). Haplotype and two-locus interaction analyses also uncovered different loci in each population.
These results suggest distinct patterns of association with SLE in African-derived populations; specific loci may be more strongly associated within select population groups.
systemic lupus erythematosus; African Americans; genetic association studies; oxygen compounds; single nucleotide polymorphism
We previously established an 80 kb haplotype upstream of TNFSF4 as a susceptibility locus in the autoimmune disease SLE. SLE-associated alleles at this locus are associated with inflammatory disorders, including atherosclerosis and ischaemic stroke. In Europeans, the TNFSF4 causal variants have remained elusive due to strong linkage disequilibrium exhibited by alleles spanning the region. Using a trans-ancestral approach to fine-map the locus, utilising 17,900 SLE and control subjects including Amerindian/Hispanics (1348 cases, 717 controls), African-Americans (AA) (1529, 2048) and better powered cohorts of Europeans and East Asians, we find strong association of risk alleles in all ethnicities; the AA association replicates in African-American Gullah (152,122). The best evidence of association comes from two adjacent markers: rs2205960-T (P = 1.71×10−34, OR = 1.43[1.26–1.60]) and rs1234317-T (P = 1.16×10−28, OR = 1.38[1.24–1.54]). Inference of fine-scale recombination rates for all populations tested finds the 80 kb risk and non-risk haplotypes in all except African-Americans. In this population the decay of recombination equates to an 11 kb risk haplotype, anchored in the 5′ region proximal to TNFSF4 and tagged by rs2205960-T after 1000 Genomes phase 1 (v3) imputation. Conditional regression analyses delineate the 5′ risk signal to rs2205960-T and the independent non-risk signal to rs1234314-C. Our case-only and SLE-control cohorts demonstrate robust association of rs2205960-T with autoantibody production. The rs2205960-T is predicted to form part of a decameric motif which binds NF-κBp65 with increased affinity compared to rs2205960-G. ChIP-seq data also indicate NF-κB interaction with the DNA sequence at this position in LCL cells. Our research suggests association of rs2205960-T with SLE across multiple groups and an independent non-risk signal at rs1234314-C. rs2205960-T is associated with autoantibody production and lymphopenia. Our data confirm a global signal at TNFSF4 and a role for the expressed product at multiple stages of lymphocyte dysregulation during SLE pathogenesis. We confirm the validity of trans-ancestral mapping in a complex trait.
Systemic lupus erythematosus (SLE/lupus) is a complex disease in which the body's immune cells cause inflammation in one or more systems to cause the associated morbidity. Hormones, the environment and genes are all causal contributors to SLE and over the past several years the genetic component of SLE has been firmly established. Several genes which are regulators of the immune system are associated with disease risk. We have established one of these, the tumour-necrosis family superfamily member 4 (TNFSF4) gene, as a lupus susceptibility gene in Northern Europeans. A major obstacle in pinpointing the marker(s) at TNFSF4 which best explain the risk of SLE has been the strong correlation (linkage disequilibrium, LD) between adjacent markers across the TNFSF4 region in this population. To address this, we have typed polymorphisms in several populations in addition to the European groups. The mixed ancestry of these populations gives a different LD pattern than that found in Europeans, presenting a method of pinpointing the section of the TNFSF4 region which results in SLE susceptibility. The Non-European populations have allowed identification of a polymorphism likely to regulate expression of TNFSF4 to increase susceptibility to SLE.
Systemic lupus erythematosus (SLE) is a common systemic autoimmune disease with complex etiology but strong clustering in families (λS = ~30). We performed a genome-wide association scan using 317,501 SNPs in 720 women of European ancestry with SLE and in 2,337 controls, and we genotyped consistently associated SNPs in two additional independent sample sets totaling 1,846 affected women and 1,825 controls. Aside from the expected strong association between SLE and the HLA region on chromosome 6p21 and the previously confirmed non-HLA locus IRF5 on chromosome 7q32, we found evidence of association with replication (1.1 × 10−7 < Poverall < 1.6 × 10−23; odds ratio 0.82–1.62)in four regions: 16p11.2 (ITGAM), 11p15.5 (KIAA1542), 3p14.3 (PXK) and 1q25.1 (rs10798269). We also found evidence for association (P < 1 × 10−5) at FCGR2A, PTPN22 and STAT4, regions previously associated with SLE and other autoimmune diseases, as well as at ≥9 other loci (P < 2 × 10−7). Our results show that numerous genes, some with known immune-related functions, predispose to SLE.
In order to identify the combination of antibody-mediated mechanisms of neutralization that result from vaccination with anthrax vaccine adsorbed (AVA), we isolated antibody secreting cells from a single donor seven days after booster vaccination with AVA and generated nine fully human monoclonal antibodies (hmAb) with high specificity for protective antigen (PA). Two of the antibodies were able to neutralize lethal toxin in vitro at low concentrations (IC50: p6C01, 0.12 µg/ml and p6F01, 0.45 µg/ml). Passive transfer of either of these hmAbs to A/J mice prior to challenge with lethal toxin conferred 80–90% protection. We demonstrate that hmAb p6C01 is neutralizing by preventing furin cleavage of PA in a dose-dependent manner, but the mechanism of p6F01 is unclear. Three additional antibodies were found to bind to domain 3 of PA and prevent oligomerization, although they did not confer significant protection in vivo and showed a significant prozone-like effect in vitro. These fully human antibodies provide insight into the neutralizing response to AVA for future subunit vaccine and passive immunotherapeutic cocktail design.
anthrax; Anthrax Vaccine Adsorbed; human monoclonal antibodies; passive immunotherapeutics; protective antigen
Systemic lupus erythematosus (SLE) is a heterogeneous autoimmune disease with significant immune system aberrations resulting from complex heritable genetics as well as environmental factors. TRAF6 is a candidate gene for SLE, which has a major role in several signaling pathways that are important for immunity and organ development.
Fifteen single-nucleotide polymorphisms (SNPs), across TRAF6 were evaluated in 7,490 SLE and 6,780 control subjects from different ancestries. Population-based case-control association analyses and meta-analyses were performed. P values, false discovery rate q values, and odds ratios with 95% confidence intervals were calculated.
Evidence of associations in multiple SNPs was detected. The best overall p values were obtained for SNPs rs5030437 and rs4755453 (p=7.85×10−5 and p=4.73×10−5, respectively) without significant heterogeneity among populations (p=0.67 and p=0.50 in Q-statistic). In addition, rs540386 previously reported to be associated with RA was found to be in LD with these two SNPs (r2= 0.95) and demonstrated evidence of association with SLE in the same direction (meta-analysis p=9.15×10−4, OR=0.89, 95%CI=0.83–0.95). Thrombocytopenia improved the overall results in different populations (meta-analysis p=1.99×10−6, OR=0.57, 95%CI=0.45–0.72, for rs5030470). Finally evidence of family based association in 34 African-American pedigrees with the presence of thrombocytopenia were detected in one available SNP rs5030437 with Z score magnitude of 2.28 (p=0.02) under a dominant model.
Our data indicate the presence of association of TRAF6 with SLE in agreement with the previous report of association with RA. These data provide further support for the involvement of TRAF6 in the pathogenesis of autoimmunity.
TRAF6; polymorphism; systemic lupus erythematosus