Evaluation of RNase therapy in systemic lupus erythematosus: a randomised phase 2a clinical trial of RSLV-132.

BACKGROUND: Circulating, extracellular RNA is the primary trigger of type I interferon in systemic lupus erythematosus (SLE), and interferon is known to play a central pathogenic role in the disease. RSLV-132 is a catalytically active human RNase molecule fused to human IgG1 Fc designed to digest RNA and thereby decrease the chronic inflammation associated with SLE. The drug was evaluated in a cohort of patients with SLE with moderate-severe cutaneous disease activity and the presence of RNA immune complexes. The primary objective of the study was the assessment of the impact of 13 doses of 10 mg/kg RSLV-132 over 6 months on the mean Cutaneous Lupus Erythematosus Disease Area and Severity Index (CLASI) score. METHODS: Sixty-five patients meeting the entry criteria of a baseline CLASI score of 10 or greater and positivity of at least one of five autoantibodies to RNA-binding proteins (SM/RNP, SSA/Ro, SSB/La, Sm, RNP) were randomly assigned (2:1) to receive 13 doses of RSLV-132 10 mg/kg or placebo, respectively. Participants received study drug for 24 weeks on days 1, 8, 15, 29, 43, 57, 71, 85, 99, 113, 127, 141 and 155 with an end-of-treatment visit on day 169 and a follow-up visit at the end of the study on day 215. The primary objective was assessed on days 85 and 169. Secondary objectives included assessment of systemic disease activity using the Systemic Lupus Erythematosus Disease Activity Index 2000 (SLEDAI-2K), the British Isles Lupus Assessment Group 2004 Index and the Physician's Global Assessment. Data from these instruments were used to calculate the SLE Responder Index 4 (SRI-4) and the British Isles Lupus Assessment Group-based Composite Lupus Assessment (BICLA) scores. RESULTS: The mean CLASI score change from baseline at day 169 was -5.7 (±7.0) in the placebo group and -6.2 (±8.5) in the RSLV-132 group. A subgroup of participants with moderate-severe systemic disease activity and high baseline SLEDAI scores (≥9) were analysed with respect to BICLA and SRI-4 responses. The RSLV-132 treated participants in the high SLEDAI subgroup had a greater percentage of BICLA responses (62% vs 44%) and SRI-4 responses (23% vs 11%) as compared with placebo. A second subgroup of participants with high baseline CLASI scores (≥21) were analysed with respect to BICLA and SRI-4 responses. The RSLV-132 treated participants in the high CLASI subgroup had a greater percentage of BICLA responses (28% vs 8%) and SRI-4 responses (39% vs 8%) as compared with placebo. CONCLUSIONS: Six months of RSLV-132 therapy consisting of a weekly loading dose of RSLV-132 for 1 month, followed by 5 months of biweekly administrations did not significantly improve the mean CLASI score relative to placebo in this cohort of patients with SLE. The study entry criteria selected patients with moderate-severe cutaneous disease activity and no minimum SLEDAI score, which resulted in a wide range of systemic disease activity from inactive to severe as measured by SLEDAI. When the participants with higher SLEDAI and CLASI scores were analysed, a trend towards clinical improvement favouring RSLV-132 was observed. The results warrant further evaluation of RSLV-132 in SLE and suggest that patients with more active systemic disease are most likely to benefit from RNase therapy.

Regulatory Architecture of the RCA Gene Cluster Captures an Intragenic TAD Boundary, CTCF-Mediated Chromatin Looping and a Long-Range Intergenic Enhancer.

The Regulators of Complement Activation (RCA) gene cluster comprises several tandemly arranged genes with shared functions within the immune system. RCA members, such as complement receptor 2 (CR2), are well-established susceptibility genes in complex autoimmune diseases. Altered expression of RCA genes has been demonstrated at both the functional and genetic level, but the mechanisms underlying their regulation are not fully characterised. We aimed to investigate the structural organisation of the RCA gene cluster to identify key regulatory elements that influence the expression of CR2 and other genes in this immunomodulatory region. Using 4C, we captured extensive CTCF-mediated chromatin looping across the RCA gene cluster in B cells and showed these were organised into two topologically associated domains (TADs). Interestingly, an inter-TAD boundary was located within the CR1 gene at a well-characterised segmental duplication. Additionally, we mapped numerous gene-gene and gene-enhancer interactions across the region, revealing extensive co-regulation. Importantly, we identified an intergenic enhancer and functionally demonstrated this element upregulates two RCA members (CR2 and CD55) in B cells. We have uncovered novel, long-range mechanisms whereby autoimmune disease susceptibility may be influenced by genetic variants, thus highlighting the important contribution of chromatin topology to gene regulation and complex genetic disease.

Elevated Detection of Dual Antibody B Cells Identifies Lupus Patients With B Cell-Reactive VH4-34 Autoantibodies.

About 5% of B cells in healthy mice and humans are allelically or isotypically included and hence co-express two different antibodies. In mice, dual antibody B cells (B2R) expand with systemic autoimmunity, co-express autoreactive and non-autoreactive antibodies, and participate in immune responses, but this phenomenon is strain dependent. This study was developed with two goals: 1) to establish the contribution of TLR and IFN receptor signaling to the development of germinal center B cells that express two antibodies in MRL/lpr mice; and 2) to determine whether B2R B cells are increased and particularly activated in a subset of adult patients diagnosed with systemic lupus erythematosus (SLE). Results from the MRL/lpr studies indicate that the enhanced differentiation of dual-κ B cells into germinal center B cells is due to a heightened response to TLR7 and TLR9 signaling, further fueled by an increased response to type II IFN. To understand the clinical and translational implications of our observations in mouse B2R B cells, cohorts of SLE patients and healthy controls were recruited and evaluated for expression of dual BCRs. Results from flow cytometry and microscopy revealed supraphysiological frequencies of κ+λ+ B2R cells in one fourth of the SLE patients. Abnormal numbers of κ+λ+ B cells correlated with higher frequencies of activated naïve B cells and age-associated B cells, and a lower proportion of "B cells that are naïve IgD+" (BND). However, results from single cell V(D)J sequencing demonstrated that these high κ+λ+ SLE patients harbored normal frequencies of κ+λ+ and other B2R B cells. and we further show that their B cells were instead decorated by κ and λ VH4-34 autoantibodies. Thus, our findings indicate that elevated flow cytometric detection of isotypically-included B cells can identify patients with high titers of B cell-reactive VH4-34 autoantibodies and abnormal distribution of B cell subsets relevant to autoimmunity.

Central human B cell tolerance manifests with a distinctive cell phenotype and is enforced via CXCR4 signaling in hu-mice.

Central B cell tolerance, the process restricting the development of many newly generated autoreactive B cells, has been intensely investigated in mouse cells while studies in humans have been hampered by the inability to phenotypically distinguish autoreactive and nonautoreactive immature B cell clones and the difficulty in accessing fresh human bone marrow samples. Using a human immune system mouse model in which all human Igκ+ B cells undergo central tolerance, we discovered that human autoreactive immature B cells exhibit a distinctive phenotype that includes lower activation of ERK and differential expression of CD69, CD81, CXCR4, and other glycoproteins. Human B cells exhibiting these characteristics were observed in fresh human bone marrow tissue biopsy specimens, although differences in marker expression were smaller than in the humanized mouse model. Furthermore, the expression of these markers was slightly altered in autoreactive B cells of humanized mice engrafted with some human immune systems genetically predisposed to autoimmunity. Finally, by treating mice and human immune system mice with a pharmacologic antagonist, we show that signaling by CXCR4 is necessary to prevent both human and mouse autoreactive B cell clones from egressing the bone marrow, indicating that CXCR4 functionally contributes to central B cell tolerance.

Phase II Randomized Trial of Rituximab Plus Cyclophosphamide Followed by Belimumab for the Treatment of Lupus Nephritis.

OBJECTIVE: To assess the safety, mechanism of action, and preliminary efficacy of rituximab followed by belimumab in the treatment of refractory lupus nephritis (LN). METHODS: In a multicenter, randomized, open-label clinical trial, 43 patients with recurrent or refractory LN were treated with rituximab, cyclophosphamide (CYC), and glucocorticoids followed by weekly belimumab infusions until week 48 (RCB group), or treated with rituximab and CYC but no belimumab infusions (RC group). Patients were followed up until week 96. Percentages of total and autoreactive B cell subsets in the patients' peripheral blood were analyzed by flow cytometry. RESULTS: Treatment with belimumab did not increase the incidence of adverse events in patients with refractory LN. At week 48, a complete or partial renal response occurred in 11 (52%) of 21 patients receiving belimumab, compared to 9 (41%) of 22 patients in the RC group who did not receive belimumab (P = 0.452). Lack of improvement in or worsening of LN was the major reason for treatment failure. B cell depletion occurred in both groups, but the percentage of B cells remained lower in those receiving belimumab (geometric mean number of B cells at week 60, 53 cells/µl in the RCB group versus 11 cells/µl in the RC group; P = 0.0012). Percentages of total and autoreactive transitional B cells increased from baseline to week 48 in both groups. However, percentages of total and autoreactive naive B cells decreased from baseline to week 48 in the belimumab group compared to the no belimumab RC group (P = 0.0349), a finding that is consistent with the observed impaired maturation of naive B cells and enhanced censoring of autoreactive B cells. CONCLUSION: The addition of belimumab to a treatment regimen with rituximab and CYC was safe in patients with refractory LN. This regimen diminished maturation of transitional to naive B cells during B cell reconstitution, and enhanced the negative selection of autoreactive B cells. Clinical efficacy was not improved with rituximab and CYC in combination with belimumab when compared to a therapeutic strategy of B cell depletion alone in patients with LN.

A plausibly causal functional lupus-associated risk variant in the STAT1-STAT4 locus.

Systemic lupus erythematosus (SLE or lupus) (OMIM: 152700) is a chronic autoimmune disease with debilitating inflammation that affects multiple organ systems. The STAT1-STAT4 locus is one of the first and most highly replicated genetic loci associated with lupus risk. We performed a fine-mapping study to identify plausible causal variants within the STAT1-STAT4 locus associated with increased lupus disease risk. Using complementary frequentist and Bayesian approaches in trans-ancestral Discovery and Replication cohorts, we found one variant whose association with lupus risk is supported across ancestries in both the Discovery and Replication cohorts: rs11889341. In B cell lines from patients with lupus and healthy controls, the lupus risk allele of rs11889341 was associated with increased STAT1 expression. We demonstrated that the transcription factor HMGA1, a member of the HMG transcription factor family with an AT-hook DNA-binding domain, has enriched binding to the risk allele compared with the non-risk allele of rs11889341. We identified a genotype-dependent repressive element in the DNA within the intron of STAT4 surrounding rs11889341. Consistent with expression quantitative trait locus (eQTL) analysis, the lupus risk allele of rs11889341 decreased the activity of this putative repressor. Altogether, we present a plausible molecular mechanism for increased lupus risk at the STAT1-STAT4 locus in which the risk allele of rs11889341, the most probable causal variant, leads to elevated STAT1 expression in B cells due to decreased repressor activity mediated by increased binding of HMGA1.

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.

Focused transcription from the human CR2/CD21 core promoter is regulated by synergistic activity of TATA and Initiator elements in mature B cells.

Complement receptor 2 (CR2/CD21) is predominantly expressed on the surface of mature B cells where it forms part of a coreceptor complex that functions, in part, to modulate B-cell receptor signal strength. CR2/CD21 expression is tightly regulated throughout B-cell development such that CR2/CD21 cannot be detected on pre-B or terminally differentiated plasma cells. CR2/CD21 expression is upregulated at B-cell maturation and can be induced by IL-4 and CD40 signaling pathways. We have previously characterized elements in the proximal promoter and first intron of CR2/CD21 that are involved in regulating basal and tissue-specific expression. We now extend these analyses to the CR2/CD21 core promoter. We show that in mature B cells, CR2/CD21 transcription proceeds from a focused TSS regulated by a non-consensus TATA box, an initiator element and a downstream promoter element. Furthermore, occupancy of the general transcriptional machinery in pre-B versus mature B-cell lines correlate with CR2/CD21 expression level and indicate that promoter accessibility must switch from inactive to active during the transitional B-cell window.

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.

Genetic associations of leptin-related polymorphisms with systemic lupus erythematosus.

Leptin is abnormally elevated in the plasma of patients with systemic lupus erythematosus (SLE), where it is thought to promote and/or sustain pro-inflammatory responses. Whether this association could reflect an increased genetic susceptibility to develop SLE is not known, and studies of genetic associations with leptin-related polymorphisms in SLE patients have been so far inconclusive. Here we genotyped DNA samples from 15,706 SLE patients and healthy matched controls from four different ancestral groups, to correlate polymorphisms of genes of the leptin pathway to risk for SLE. It was found that although several SNPs showed weak associations, those associations did not remain significant after correction for multiple testing. These data do not support associations between defined leptin-related polymorphisms and increased susceptibility to develop SLE.

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.

Two functional lupus-associated BLK promoter variants control cell-type- and developmental-stage-specific transcription.

Efforts to identify lupus-associated causal variants in the FAM167A/BLK locus on 8p21 are hampered by highly associated noncausal variants. In this report, we used a trans-population mapping and sequencing strategy to identify a common variant (rs922483) in the proximal BLK promoter and a tri-allelic variant (rs1382568) in the upstream alternative BLK promoter as putative causal variants for association with systemic lupus erythematosus. The risk allele (T) at rs922483 reduced proximal promoter activity and modulated alternative promoter usage. Allelic differences at rs1382568 resulted in altered promoter activity in B progenitor cell lines. Thus, our results demonstrated that both lupus-associated functional variants contribute to the autoimmune disease association by modulating transcription of BLK in B cells and thus potentially altering immune responses.

Lupus risk variants in the PXK locus alter B-cell receptor internalization.

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.

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.

Advances in lupus genetics.

PURPOSE OF REVIEW: The purpose of this review is to highlight recent progress in defining the effects of variation in the germline DNA sequence on the susceptibility to systemic lupus erythematosus. RECENT FINDINGS: The genetic contributions to lupus susceptibility have been increasingly defined in recent years, although a large proportion of the genetic heritability of this disease remains undiscovered. Various approaches have been taken to determine this 'missing heritability' and to further characterize the known genetic effects in order to more fully elucidate the cause and pathogenesis of this complex and heterogenous disease. SUMMARY: Characterizing the genetic contributions to lupus and their functional effects will advance our understanding of the biology of this disease and aid in the development of new therapeutics.

ABIN1 dysfunction as a genetic basis for lupus nephritis.

The genetic factors underlying the pathogenesis of lupus nephritis associated with systemic lupus erythematosus are largely unknown, although animal studies indicate that nuclear factor (NF)-κB is involved. We reported previously that a knockin mouse expressing an inactive form of ABIN1 (ABIN1[D485N]) develops lupus-like autoimmune disease and demonstrates enhanced activation of NF-κB and mitogen-activated protein kinases in immune cells after toll-like receptor stimulation. In the current study, we show that ABIN1[D485N] mice develop progressive GN similar to class III and IV lupus nephritis in humans. To investigate the clinical relevance of ABIN1 dysfunction, we genotyped five single-nucleotide polymorphisms in the gene encoding ABIN1, TNIP1, in samples from European-American, African American, Asian, Gullah, and Hispanic participants in the Large Lupus Association Study 2. Comparing cases of systemic lupus erythematosus with nephritis and cases of systemic lupus erythematosus without nephritis revealed strong associations with lupus nephritis at rs7708392 in European Americans and rs4958881 in African Americans. Comparing cases of systemic lupus erythematosus with nephritis and healthy controls revealed a stronger association at rs7708392 in European Americans but not at rs4958881 in African Americans. Our data suggest that variants in the TNIP1 gene are associated with the risk for lupus nephritis and could be mechanistically involved in disease development via aberrant regulation of NF-κB and mitogen-activated protein kinase activity.

PTPN22 association in systemic lupus erythematosus (SLE) with respect to individual ancestry and clinical sub-phenotypes.

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.

Trans-ancestral studies fine map the SLE-susceptibility locus TNFSF4.

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.