ßIV Spectrinopathies Cause Profound Intellectual Disability, Congenital Hypotonia, and Motor Axonal Neuropathy.

ßIV spectrin links ankyrinG (AnkG) and clustered ion channels at axon initial segments (AISs) and nodes of Ranvier to the axonal cytoskeleton. Here, we report bi-allelic pathogenic SPTBN4 variants (three homozygous and two compound heterozygous) that cause a severe neurological syndrome that includes congenital hypotonia, intellectual disability, and motor axonal and auditory neuropathy. We introduced these variants into ßIV spectrin, expressed these in neurons, and found that 5/7 were loss-of-function variants disrupting AIS localization or abolishing phosphoinositide binding. Nerve biopsies from an individual with a loss-of-function variant had reduced nodal Na+ channels and no nodal KCNQ2 K+ channels. Modeling the disease in mice revealed that although ankyrinR (AnkR) and ßI spectrin can cluster Na+ channels and partially compensate for the loss of AnkG and ßIV spectrin at nodes of Ranvier, AnkR and ßI spectrin cannot cluster KCNQ2- and KCNQ3-subunit-containing K+ channels. Our findings define a class of spectrinopathies and reveal the molecular pathologies causing nervous-system dysfunction.

Genetic fine mapping of systemic lupus erythematosus MHC associations in Europeans and African Americans.

Genetic variation within the major histocompatibility complex (MHC) contributes substantial risk for systemic lupus erythematosus, but high gene density, extreme polymorphism and extensive linkage disequilibrium (LD) have made fine mapping challenging. To address the problem, we compared two association techniques in two ancestrally diverse populations, African Americans (AAs) and Europeans (EURs). We observed a greater number of Human Leucocyte Antigen (HLA) alleles in AA consistent with the elevated level of recombination in this population. In EUR we observed 50 different A-C-B-DRB1-DQA-DQB multilocus haplotype sequences per hundred individuals; in the AA sample, these multilocus haplotypes were twice as common compared to Europeans. We also observed a strong narrow class II signal in AA as opposed to the long-range LD observed in EUR that includes class I alleles. We performed a Bayesian model choice of the classical HLA alleles and a frequentist analysis that combined both single nucleotide polymorphisms (SNPs) and classical HLA alleles. Both analyses converged on a similar subset of risk HLA alleles: in EUR HLA- B*08:01 + B*18:01 + (DRB1*15:01 frequentist only) + DQA*01:02 + DQB*02:01 + DRB3*02 and in AA HLA-C*17:01 + B*08:01 + DRB1*15:03 + (DQA*01:02 frequentist only) + DQA*02:01 + DQA*05:01+ DQA*05:05 + DQB*03:19 + DQB*02:02. We observed two additional independent SNP associations in both populations: EUR rs146903072 and rs501480; AA rs389883 and rs114118665. The DR2 serotype was best explained by DRB1*15:03 + DQA*01:02 in AA and by DRB1*15:01 + DQA*01:02 in EUR. The DR3 serotype was best explained by DQA*05:01 in AA and by DQB*02:01 in EUR. Despite some differences in underlying HLA allele risk models in EUR and AA, SNP signals across the extended MHC showed remarkable similarity and significant concordance in direction of effect for risk-associated variants.

Identification of a Sjögren's syndrome susceptibility locus at OAS1 that influences isoform switching, protein expression, and responsiveness to type I interferons.

Sjögren's syndrome (SS) is a common, autoimmune exocrinopathy distinguished by keratoconjunctivitis sicca and xerostomia. Patients frequently develop serious complications including lymphoma, pulmonary dysfunction, neuropathy, vasculitis, and debilitating fatigue. Dysregulation of type I interferon (IFN) pathway is a prominent feature of SS and is correlated with increased autoantibody titers and disease severity. To identify genetic determinants of IFN pathway dysregulation in SS, we performed cis-expression quantitative trait locus (eQTL) analyses focusing on differentially expressed type I IFN-inducible transcripts identified through a transcriptome profiling study. Multiple cis-eQTLs were associated with transcript levels of 2'-5'-oligoadenylate synthetase 1 (OAS1) peaking at rs10774671 (PeQTL = 6.05 × 10-14). Association of rs10774671 with SS susceptibility was identified and confirmed through meta-analysis of two independent cohorts (Pmeta = 2.59 × 10-9; odds ratio = 0.75; 95% confidence interval = 0.66-0.86). The risk allele of rs10774671 shifts splicing of OAS1 from production of the p46 isoform to multiple alternative transcripts, including p42, p48, and p44. We found that the isoforms were differentially expressed within each genotype in controls and patients with and without autoantibodies. Furthermore, our results showed that the three alternatively spliced isoforms lacked translational response to type I IFN stimulation. The p48 and p44 isoforms also had impaired protein expression governed by the 3' end of the transcripts. The SS risk allele of rs10774671 has been shown by others to be associated with reduced OAS1 enzymatic activity and ability to clear viral infections, as well as reduced responsiveness to IFN treatment. Our results establish OAS1 as a risk locus for SS and support a potential role for defective viral clearance due to altered IFN response as a genetic pathophysiological basis of this complex autoimmune disease.

Genetic variants at the 16p13 locus confer risk for eosinophilic esophagitis.

Eosinophilic esophagitis (EoE) is a chronic inflammatory disease of the esophagus triggered by immune hypersensitivity to food. Herein, we tested whether genetic risk factors for known, non-allergic, immune-mediated diseases, particularly those involving autoimmunity, were associated with EoE risk. We used the high-density Immunochip platform, encoding 200,000 genetic variants for major auto-immune disease. Accordingly, 1214 subjects with EoE of European ancestry and 3734 population controls were genotyped and assessed using data directly generated or imputed from the previously published GWAS. We found lack of association of EoE with the genetic variants in the major histocompatibility complex (MHC) class I, II, and III genes and nearly all other loci using a highly powered study design with dense genotyping throughout the locus. Importantly, we identified an EoE risk locus at 16p13 with genome-wide significance (Pcombined=2.05 × 10-9, odds ratio = 0.76-0.81). This region is known to encode for the genes CLEC16A, DEXI, and CIITI, which are expressed in immune cells and esophageal epithelial cells. Suggestive EoE risk were also seen 5q23 (intergenic) and 7p15 (JAZF1). Overall, we have identified an additional EoE risk locus at 16p13 and highlight a shared and unique genetic etiology of EoE with a spectrum of immune-associated diseases.

Association of Epstein-Barr virus serological reactivation with transitioning to systemic lupus erythematosus in at-risk individuals.

OBJECTIVE: Systemic lupus erythematosus (SLE) is a systemic autoimmune disease with unknown aetiology. Epstein-Barr virus (EBV) is an environmental factor associated with SLE. EBV maintains latency in B cells with frequent reactivation measured by antibodies against viral capsid antigen (VCA) and early antigen (EA). In this study, we determined whether EBV reactivation and single nucleotide polymorphisms (SNPs) in EBV-associated host genes are associated with SLE transition. METHODS: SLE patient relatives (n=436) who did not have SLE at baseline were recontacted after 6.3 (±3.9) years and evaluated for interim transitioning to SLE (≥4 cumulative American College of Rheumatology criteria); 56 (13%) transitioned to SLE prior to the follow-up visit. At both visits, detailed demographic, environmental, clinical information and blood samples were obtained. Antibodies against viral antigens were measured by ELISA. SNPs in IL10, CR2, TNFAIP3 and CD40 genes were typed by ImmunoChip. Generalised estimating equations were used to test associations between viral antibody levels and transitioning to SLE. RESULTS: Mean baseline VCA IgG (4.879±1.797 vs 3.866±1.795, p=0.0003) and EA IgG (1.192±1.113 vs 0.7774±0.8484, p=0.0236) levels were higher in transitioned compared with autoantibody negative non-transitioned relatives. Increased VCA IgG and EA IgG were associated with transitioning to SLE (OR 1.28 95% CI 1.07 to 1.53, p=0.007, OR 1.43 95% CI 1.06 to 1.93, p=0.02, respectively). Significant interactions were observed between CD40 variant rs48100485 and VCA IgG levels and IL10 variant rs3024493 and VCA IgA levels in transitioning to SLE. CONCLUSION: Heightened serologic reactivation of EBV increases the probability of transitioning to SLE in unaffected SLE relatives.

Lupus Risk Variant Increases pSTAT1 Binding and Decreases ETS1 Expression.

Genetic variants at chromosomal region 11q23.3, near the gene ETS1, have been associated with systemic lupus erythematosus (SLE), or lupus, in independent cohorts of Asian ancestry. Several recent studies have implicated ETS1 as a critical driver of immune cell function and differentiation, and mice deficient in ETS1 develop an SLE-like autoimmunity. We performed a fine-mapping study of 14,551 subjects from multi-ancestral cohorts by starting with genotyped variants and imputing to all common variants spanning ETS1. By constructing genetic models via frequentist and Bayesian association methods, we identified 16 variants that are statistically likely to be causal. We functionally assessed each of these variants on the basis of their likelihood of affecting transcription factor binding, miRNA binding, or chromatin state. Of the four variants that we experimentally examined, only rs6590330 differentially binds lysate from B cells. Using mass spectrometry, we found more binding of the transcription factor signal transducer and activator of transcription 1 (STAT1) to DNA near the risk allele of rs6590330 than near the non-risk allele. Immunoblot analysis and chromatin immunoprecipitation of pSTAT1 in B cells heterozygous for rs6590330 confirmed that the risk allele increased binding to the active form of STAT1. Analysis with expression quantitative trait loci indicated that the risk allele of rs6590330 is associated with decreased ETS1 expression in Han Chinese, but not other ancestral cohorts. We propose a model in which the risk allele of rs6590330 is associated with decreased ETS1 expression and increases SLE risk by enhancing the binding of pSTAT1.

ABIN1 Determines Severity of Glomerulonephritis via Activation of Intrinsic Glomerular Inflammation.

Transcription factor NF-κB regulates expression of numerous genes that control inflammation and is activated in glomerular cells in glomerulonephritis (GN). We previously identified genetic variants for a NF-κB regulatory, ubiquitin-binding protein ABIN1 as risk factors for GN in systemic autoimmunity. The goal was to define glomerular inflammatory events controlled by ABIN1 function in GN. Nephrotoxic serum nephritis was induced in wild-type (WT) and ubiquitin-binding deficient ABIN1[D485N] mice, and renal pathophysiology and glomerular inflammatory phenotypes were assessed. Proteinuria was also measured in ABIN1[D485N] mice transplanted with WT mouse bone marrow. Inflammatory activation of ABIN1[D472N] (D485N homolog) cultured human-derived podocytes, and interaction with primary human neutrophils were also assessed. Disruption of ABIN1 function exacerbated proteinuria, podocyte injury, glomerular NF-κB activity, glomerular expression of inflammatory mediators, and glomerular recruitment and retention of neutrophils in antibody-mediated nephritis. Transplantation of WT bone marrow did not prevent the increased proteinuria in ABIN1[D845N] mice. Tumor necrosis factor-stimulated enhanced expression and secretion of NF-κB-targeted proinflammatory mediators in ABIN1[D472N] cultured podocytes compared with WT cells. Supernatants from ABIN1[D472N] podocytes accelerated chemotaxis of human neutrophils, and ABIN1[D472N] podocytes displayed a greater susceptibility to injurious morphologic findings induced by neutrophil granule contents. These studies define a novel role for ABIN1 dysfunction and NF-κB in mediating GN through proinflammatory activation of podocytes.

Functional activation of PPARγ in human upper aerodigestive cancer cell lines.

Upper aerodigestive cancer is an aggressive malignancy with relatively stagnant long-term survival rates over 20 yr. Recent studies have demonstrated that exploitation of PPARγ pathways may be a novel therapy for cancer and its prevention. We tested whether PPARγ is expressed and inducible in aerodigestive carcinoma cells and whether it is present in human upper aerodigestive tumors. Human oral cancer CA-9-22 and NA cell lines were treated with the PPAR activators eicosatetraynoic acid (ETYA), 15-deoxy-δ- 12,14-prostaglandin J2 (PG-J2), and the thiazolidinedione, ciglitazone, and evaluated for their ability to functionally activate PPARγ luciferase reporter gene constructs. Cellular proliferation and clonogenic potential after PPARγ ligand treatment were also evaluated. Aerodigestive cancer specimens and normal tissues were evaluated for PPARγ expression on gene expression profiling and immunoblotting. Functional activation of PPARγ reporter gene constructs and increases in PPARγ protein were confirmed in the nuclear compartment after PPARγ ligand treatment. Significant decreases in cell proliferation and clonogenic potential resulted from treatment. Lipid accumulation was induced by PPARγ activator treatment. 75% of tumor specimens and 100% of normal control tissues expressed PPARγ RNA, and PPARγ protein was confirmed in 66% of tumor specimens analyzed by immunoblotting. We conclude PPARγ can be functionally activated in upper aerodigestive cancer and that its activation downregulates several features of the neoplastic phenotype. PPARγ expression in human upper aerodigestive tract tumors and normal cells potentially legitimizes it as a novel intervention target in this disease. © 2016 Wiley Periodicals, Inc.

Combined role of vitamin D status and CYP24A1 in the transition to systemic lupus erythematosus.

OBJECTIVE: We examined whether measures of vitamin D were associated with transitioning to systemic lupus erythematosus (SLE) in individuals at risk for SLE. METHODS: 436 individuals who reported having a relative with SLE but who did not have SLE themselves were evaluated at baseline and again an average of 6.3 (±3.9) years later. Fifty-six individuals transitioned to SLE (≥4 cumulative American College of Rheumatology criteria). 25-Hydroxyvitamin D (25[OH]D) levels were measured by ELISA. Six single-nucleotide polymorphisms in four vitamin D genes were genotyped. Generalised estimating equations, adjusting for correlation within families, were used to test associations between the vitamin D variables and the outcome of transitioning to SLE. RESULTS: Mean baseline 25[OH]D levels (p=0.42) and vitamin D supplementation (p=0.65) were not different between those who did and did not transition to SLE. Vitamin D deficiency (25[OH]D <20 ng/mL) was greater in those who transitioned compared with those who did not transition to SLE (46% vs 33%, p=0.05). The association between 25[OH]D and SLE was modified by CYP24A1 rs4809959, where for each additional minor allele increased 25[OH]D was associated with decreased SLE risk: zero minor alleles (adjusted OR: 1.03, CI 0.98 to 1.09), one minor allele (adjusted OR: 1.01, CI 0.97 to 1.05) and two minor alleles (adjusted OR: 0.91, CI 0.84 to 0.98). Similarly, vitamin D deficiency significantly increased the risk of transitioning to SLE in those with two minor alleles at rs4809959 (adjusted OR: 4.90, CI 1.33 to 18.04). CONCLUSIONS: Vitamin D status and CYP24A1 may have a combined role in the transition to SLE in individuals at increased genetic risk for SLE.

Activating mutations in STIM1 and ORAI1 cause overlapping syndromes of tubular myopathy and congenital miosis.

Signaling through the store-operated Ca(2+) release-activated Ca(2+) (CRAC) channel regulates critical cellular functions, including gene expression, cell growth and differentiation, and Ca(2+) homeostasis. Loss-of-function mutations in the CRAC channel pore-forming protein ORAI1 or the Ca(2+) sensing protein stromal interaction molecule 1 (STIM1) result in severe immune dysfunction and nonprogressive myopathy. Here, we identify gain-of-function mutations in the cytoplasmic domain of STIM1 (p.R304W) associated with thrombocytopenia, bleeding diathesis, miosis, and tubular myopathy in patients with Stormorken syndrome, and in ORAI1 (p.P245L), associated with a Stormorken-like syndrome of congenital miosis and tubular aggregate myopathy but without hematological abnormalities. Heterologous expression of STIM1 p.R304W results in constitutive activation of the CRAC channel in vitro, and spontaneous bleeding accompanied by reduced numbers of thrombocytes in zebrafish embryos, recapitulating key aspects of Stormorken syndrome. p.P245L in ORAI1 does not make a constitutively active CRAC channel, but suppresses the slow Ca(2+)-dependent inactivation of the CRAC channel, thus also functioning as a gain-of-function mutation. These data expand our understanding of the phenotypic spectrum of dysregulated CRAC channel signaling, advance our knowledge of the molecular function of the CRAC channel, and suggest new therapies aiming at attenuating store-operated Ca(2+) entry in the treatment of patients with Stormorken syndrome and related pathologic conditions.

Genetic variants near TNFAIP3 on 6q23 are associated with systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is an autoimmune disease influenced by genetic and environmental factors. We carried out a genome-wide association scan and replication study and found an association between SLE and a variant in TNFAIP3 (rs5029939, meta-analysis P = 2.89 x 10(-12), OR = 2.29). We also found evidence of two independent signals near TNFAIP3 associated with SLE, including one previously associated with rheumatoid arthritis (RA). These results establish that variants near TNFAIP3 contribute to differential risk of SLE and RA.

Allelic heterogeneity in NCF2 associated with systemic lupus erythematosus (SLE) susceptibility across four ethnic populations.

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 (P(EA) = 1.01 × 10(-54), PHS = 3.68 × 10(-10), P(AA) = 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 (P(HS) = 7.04 × 10(-7), P(KR) = 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.

Decreased SMG7 expression associates with lupus-risk variants and elevated antinuclear antibody production.

OBJECTIVES: Following up the systemic lupus erythematosus (SLE) genome-wide association studies (GWAS) identification of NMNAT2 at rs2022013, we fine-mapped its 150 kb flanking regions containing NMNAT2 and SMG7 in a 15 292 case-control multi-ancestry population and tested functions of identified variants. METHODS: We performed genotyping using custom array, imputation by IMPUTE 2.1.2 and allele specific functions using quantitative real-time PCR and luciferase reporter transfections. SLE peripheral blood mononuclear cells (PBMCs) were cultured with small interfering RNAs to measure antinuclear antibody (ANA) and cyto/chemokine levels in supernatants using ELISA. RESULTS: We confirmed association at NMNAT2 in European American (EA) and Amerindian/Hispanic ancestries, and identified independent signal at SMG7 tagged by rs2702178 in EA only (p=2.4×10-8, OR=1.23 (95% CI 1.14 to 1.32)). In complete linkage disequilibrium with rs2702178, rs2275675 in the promoter region robustly associated with SMG7 mRNA levels in multiple expression quantitative trait locus (eQTL) datasets. Its risk allele was dose-dependently associated with decreased SMG7 mRNA levels in PBMCs of 86 patients with SLE and 119 controls (p=1.1×10-3 and 6.8×10-8, respectively) and conferred reduced transcription activity in transfected HEK-293 (human embryonic kidney cell line) and Raji cells (p=0.0035 and 0.0037, respectively). As a critical component in the nonsense-mediated mRNA decay pathway, SMG7 could regulate autoantigens including ribonucleoprotein (RNP) and Smith (Sm). We showed SMG7 mRNA levels in PBMCs correlated inversely with ANA titres of patients with SLE (r=-0.31, p=0.01), and SMG7 knockdown increased levels of ANA IgG and chemokine (C-C motif) ligand 19 in SLE PBMCs (p=2.0×10-5 and 2.0×10-4, respectively). CONCLUSION: We confirmed NMNAT2 and identified independent SMG7 association with SLE. The inverse relationship between levels of the risk allele-associated SMG7 mRNAs and ANA suggested the novel contribution of mRNA surveillance pathway to SLE pathogenesis.

Preferential association of a functional variant in complement receptor 2 with antibodies to double-stranded DNA.

OBJECTIVES: Systemic lupus erythematosus (SLE; OMIM 152700) is characterised by the production of antibodies to nuclear antigens. We previously identified variants in complement receptor 2 (CR2/CD21) that were associated with decreased risk of SLE. This study aimed to identify the causal variant for this association. METHODS: Genotyped and imputed genetic variants spanning CR2 were assessed for association with SLE in 15 750 case-control subjects from four ancestral groups. Allele-specific functional effects of associated variants were determined using quantitative real-time PCR, quantitative flow cytometry, electrophoretic mobility shift assay (EMSA) and chromatin immunoprecipitation (ChIP)-PCR. RESULTS: The strongest association signal was detected at rs1876453 in intron 1 of CR2 (pmeta=4.2×10(-4), OR 0.85), specifically when subjects were stratified based on the presence of dsDNA autoantibodies (case-control pmeta=7.6×10(-7), OR 0.71; case-only pmeta=1.9×10(-4), OR 0.75). Although allele-specific effects on B cell CR2 mRNA or protein levels were not identified, levels of complement receptor 1 (CR1/CD35) mRNA and protein were significantly higher on B cells of subjects harbouring the minor allele (p=0.0248 and p=0.0006, respectively). The minor allele altered the formation of several DNA protein complexes by EMSA, including one containing CCCTC-binding factor (CTCF), an effect that was confirmed by ChIP-PCR. CONCLUSIONS: These data suggest that rs1876453 in CR2 has long-range effects on gene regulation that decrease susceptibility to lupus. Since the minor allele at rs1876453 is preferentially associated with reduced risk of the highly specific dsDNA autoantibodies that are present in preclinical, active and severe lupus, understanding its mechanisms will have important therapeutic implications.

Forward genetic screening identifies a small molecule that blocks Toxoplasma gondii growth by inhibiting both host- and parasite-encoded kinases.

The simultaneous targeting of host and pathogen processes represents an untapped approach for the treatment of intracellular infections. Hypoxia-inducible factor-1 (HIF-1) is a host cell transcription factor that is activated by and required for the growth of the intracellular protozoan parasite Toxoplasma gondii at physiological oxygen levels. Parasite activation of HIF-1 is blocked by inhibiting the family of closely related Activin-Like Kinase (ALK) host cell receptors ALK4, ALK5, and ALK7, which was determined in part by use of an ALK4,5,7 inhibitor named SB505124. Besides inhibiting HIF-1 activation, SB505124 also potently blocks parasite replication under normoxic conditions. To determine whether SB505124 inhibition of parasite growth was exclusively due to inhibition of ALK4,5,7 or because the drug inhibited a second kinase, SB505124-resistant parasites were isolated by chemical mutagenesis. Whole-genome sequencing of these mutants revealed mutations in the Toxoplasma MAP kinase, TgMAPK1. Allelic replacement of mutant TgMAPK1 alleles into wild-type parasites was sufficient to confer SB505124 resistance. SB505124 independently impacts TgMAPK1 and ALK4,5,7 signaling since drug resistant parasites could not activate HIF-1 in the presence of SB505124 or grow in HIF-1 deficient cells. In addition, TgMAPK1 kinase activity is inhibited by SB505124. Finally, mice treated with SB505124 had significantly lower tissue burdens following Toxoplasma infection. These data therefore identify SB505124 as a novel small molecule inhibitor that acts by inhibiting two distinct targets, host HIF-1 and TgMAPK1.

Human effector B lymphocytes express ARID3a and secrete interferon alpha.

Previously, we determined that enhanced disease activity in patients with systemic lupus erythematosus (SLE) was associated with dramatic increases in numbers of B lymphocytes expressing the transcription factor ARID3a. Our data now indicate ARID3a is important for interferon alpha (IFNa) expression and show a strong association between ARID3a expression and transcription of genes associated with lupus IFN signatures. Furthermore, both ARID3a and IFNa production were elicited in healthy control B cells upon stimulation with the TLR 9 agonist, CpG. Importantly, secretion of IFNa from ARID3a+ healthy B lymphocytes stimulated increased IFNa production in plasmacytoid dendritic cells. These data identify ARID3a+ B cells as a novel type of effector B cell, and link ARID3a expression in B lymphocytes to IFN-associated inflammatory responses in SLE.

An enhancer element harboring variants associated with systemic lupus erythematosus engages the TNFAIP3 promoter to influence A20 expression.

Functional characterization of causal variants present on risk haplotypes identified through genome-wide association studies (GWAS) is a primary objective of human genetics. In this report, we evaluate the function of a pair of tandem polymorphic dinucleotides, 42 kb downstream of the promoter of TNFAIP3, (rs148314165, rs200820567, collectively referred to as TT>A) recently nominated as causal variants responsible for genetic association of systemic lupus erythematosus (SLE) with tumor necrosis factor alpha inducible protein 3 (TNFAIP3). TNFAIP3 encodes the ubiquitin-editing enzyme, A20, a key negative regulator of NF-κB signaling. A20 expression is reduced in subjects carrying the TT>A risk alleles; however, the underlying functional mechanism by which this occurs is unclear. We used a combination of electrophoretic mobility shift assays (EMSA), mass spectrometry (MS), reporter assays, chromatin immunoprecipitation-PCR (ChIP-PCR) and chromosome conformation capture (3C) EBV transformed lymphoblastoid cell lines (LCL) from individuals carrying risk and non-risk TNFAIP3 haplotypes to characterize the effect of TT>A on A20 expression. Our results demonstrate that the TT>A variants reside in an enhancer element that binds NF-κB and SATB1 enabling physical interaction of the enhancer with the TNFAIP3 promoter through long-range DNA looping. Impaired binding of NF-κB to the TT>A risk alleles or knockdown of SATB1 expression by shRNA, inhibits the looping interaction resulting in reduced A20 expression. Together, these data reveal a novel mechanism of TNFAIP3 transcriptional regulation and establish the functional basis by which the TT>A risk variants attenuate A20 expression through inefficient delivery of NF-κB to the TNFAIP3 promoter. These results provide critical functional evidence supporting a direct causal role for TT>A in the genetic predisposition to SLE.

Admixture mapping in lupus identifies multiple functional variants within IFIH1 associated with apoptosis, inflammation, and autoantibody production.

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 [P(meta) = 5.20×10(-14); odds ratio, 95% confidence interval = 0.82 (0.78-0.87)], and two missense variants, rs1990760 (Ala946Thr) [P(meta) = 3.08×10(-7); 0.88 (0.84-0.93)] and rs10930046 (Arg460His) [P(dom) = 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.

MicroRNA-3148 modulates allelic expression of toll-like receptor 7 variant associated with systemic lupus erythematosus.

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 (R(2) = 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.

Preferential binding to Elk-1 by SLE-associated IL10 risk allele upregulates IL10 expression.

Immunoregulatory cytokine interleukin-10 (IL-10) is elevated in sera from patients with systemic lupus erythematosus (SLE) correlating with disease activity. The established association of IL10 with SLE and other autoimmune diseases led us to fine map causal variant(s) and to explore underlying mechanisms. We assessed 19 tag SNPs, covering the IL10 gene cluster including IL19, IL20 and IL24, for association with SLE in 15,533 case and control subjects from four ancestries. The previously reported IL10 variant, rs3024505 located at 1 kb downstream of IL10, exhibited the strongest association signal and was confirmed for association with SLE in European American (EA) (P = 2.7×10⁻8, OR = 1.30), but not in non-EA ancestries. SNP imputation conducted in EA dataset identified three additional SLE-associated SNPs tagged by rs3024505 (rs3122605, rs3024493 and rs3024495 located at 9.2 kb upstream, intron 3 and 4 of IL10, respectively), and SLE-risk alleles of these SNPs were dose-dependently associated with elevated levels of IL10 mRNA in PBMCs and circulating IL-10 protein in SLE patients and controls. Using nuclear extracts of peripheral blood cells from SLE patients for electrophoretic mobility shift assays, we identified specific binding of transcription factor Elk-1 to oligodeoxynucleotides containing the risk (G) allele of rs3122605, suggesting rs3122605 as the most likely causal variant regulating IL10 expression. Elk-1 is known to be activated by phosphorylation and nuclear localization to induce transcription. Of interest, phosphorylated Elk-1 (p-Elk-1) detected only in nuclear extracts of SLE PBMCs appeared to increase with disease activity. Co-expression levels of p-Elk-1 and IL-10 were elevated in SLE T, B cells and monocytes, associated with increased disease activity in SLE B cells, and were best downregulated by ERK inhibitor. Taken together, our data suggest that preferential binding of activated Elk-1 to the IL10 rs3122605-G allele upregulates IL10 expression and confers increased risk for SLE in European Americans.