Conditional knockout of oestrogen receptor alpha in CD11c+ cells impacts female survival and inflammatory cytokine profile in murine lupus.

Oestrogen and oestrogen receptor alpha (ERα) have been implicated in systemic lupus erythematosus pathogenesis. ERα signalling influences dendritic cell (DC) development and function, as well as inflammation and downstream immune responses. We previously reported that ERα modulates multiple Toll-like receptor-stimulated pathways in both conventional and plasmacytoid DCs in lupus-prone mice. For example, CD11chi MHCII+ cell numbers are reduced in mice with global ERα deficiency or when expressing a short variant of ERα. Herein, RNA-seq analysis of CD11chi cells from bone marrow of NZM2410 mice expressing WT ERα versus ERα short versus ERα null revealed differentially expressed complement genes, interferon-related genes and cytokine signalling (e.g., IL-17 and Th17 pathways). To better understand the role of ERα in CD11c+ cells, lupus prone NZM2410 mice with selective deletion of the Esr1 gene in CD11c+ cells were generated. Phenotype and survival of these mice were similar with the exception of Cre positive (CrePos) female mice. CrePos females, but not males, all died unexpectedly prior to 35 weeks. DC subsets were not significantly different between groups. Since ERα is necessary for robust development of DCs, this result suggests that DC fate was determined prior to CD11c expression and subsequent ERα deletion (i.e., proximally in DC ontogeny). Overall, findings point to a clear functional role for ERα in regulating cytokine signalling and inflammation, suggesting that further study into ERα-mediated regulatory mechanisms in DCs and other immune cell types is warranted.

Human SLE variant NCF1-R90H promotes kidney damage and murine lupus through enhanced Tfh2 responses induced by defective efferocytosis of macrophages.

OBJECTIVES: We previously identified a hypomorphic variant, p.Arg90His (p.R90H) of neutrophil cytosolic factor 1 (NCF1, a regulatory subunit of phagocyte nicotinamide adenine dinucleotide phosphate (NADPH) oxidase 2 complex), as an putative causal variant for systemic lupus erythematosus (SLE), and established a knock-in (KI) H90 variant in the C57BL/6 background to study how this variant promotes lupus development. METHODS: Wild type (WT) and KI littermates were assessed for immune profiles and lupus-like features. Disease activity and renal damage of patients with SLE were assessed by systemic lupus erythematosus disease activity index (SLEDAI) and renal items of systemic lupus international collaborating clinics (SLICC), respectively. RESULTS: Compared with WT littermates, 5-week-old homozygous KI mice had reduced oxidative burst, splenomegaly, elevated type I interferon (IFN-I) scores, increased ratios of splenic follicular T helper 2 (Tfh2) to either T follicular regulatory (Tfr) or Tfh1 cells, increased ANA+ follicular, germinal centre and plasma cells without spontaneous kidney disease up to 1 year of age. Pristane treatment exacerbated the immune dysregulation and induced IFN-I-dependent kidney disease in 36-week-old H90 KI female mice. Decreased efferocytosis of macrophages derived from KI mice and patients with homozygous H90 SLE promoted elevated ratios of Tfh2/Tfr and Tfh2/Tfh1 as well as dysregulated humoral responses due to reduced voltage-gated proton channel 1 (Hv1)-dependent acidification of phagosome pH to neutralise the decreased electrogenic effect of the H90 variant, resulting in impaired maturation and phagosome proteolysis, and increased autoantibody production and kidney damage in mice and patients with SLE of multiple ancestries. CONCLUSIONS: A lupus causal variant, NCF1-H90, reduces macrophage efferocytosis, enhances Tfh2 responses and promotes autoantibody production and kidney damage in both mice and patients with SLE.

Loss-of-function variants in SAT1 cause X-linked childhood-onset systemic lupus erythematosus.

OBJECTIVES: Families that contain multiple siblings affected with childhood onset of systemic lupus erythematosus (SLE) likely have strong genetic predispositions. We performed whole exome sequencing (WES) to identify familial rare risk variants and to assess their effects in lupus. METHODS: Sanger sequencing validated the two ultra-rare, predicted pathogenic risk variants discovered by WES and identified additional variants in 562 additional patients with SLE. Effects of a splice site variant and a frameshift variant were assessed using a Minigene assay and CRISPR/Cas9-mediated knock-in (KI) mice, respectively. RESULTS: The two familial ultra-rare, predicted loss-of-function (LOF) SAT1 variants exhibited X-linked recessive Mendelian inheritance in two unrelated African-American families. Each LOF variant was transmitted from the heterozygous unaffected mother to her two sons with childhood-onset SLE. The p.Asp40Tyr variant affected a splice donor site causing deleterious transcripts. The young hemizygous male and homozygous female Sat1 p.Glu92Leufs*6 KI mice spontaneously developed splenomegaly, enlarged glomeruli with leucocyte infiltration, proteinuria and elevated expression of type I interferon-inducible genes. SAT1 is highly expressed in neutrophils and encodes spermidine/spermine-N1-acetyltransferase 1 (SSAT1), a rate-limiting enzyme in polyamine catabolism. Young male KI mice exhibited neutrophil defects and decreased proportions of Foxp3 +CD4+ T-cell subsets. Circulating neutrophil counts and proportions of Foxp3 +CD4+ T cells correlated with decreased plasma levels of spermine in treatment-naive, incipient SLE patients. CONCLUSIONS: We identified two novel SAT1 LOF variants, showed the ability of the frameshift variant to confer murine lupus, highlighted the pathogenic role of dysregulated polyamine catabolism and identified SAT1 LOF variants as new monogenic causes for SLE.

Staphylococcus aureus peptidoglycan (PGN) induces pathogenic autoantibody production via autoreactive B cell receptor clonal selection, implications in systemic lupus erythematosus.

OBJECTIVES: There is an intricate interplay between the microbiome and the immune response impacting development of normal immunity and autoimmunity. However, we do not fully understand how the microbiome affects production of natural-like and pathogenic autoantibodies. Peptidoglycan (PGN) is a component of the bacterial cell wall which is highly antigenic. PGNs from different bacteria can differ in their immune regulatory activities. METHODS: C57BL/6 and MRL/lpr mice were intraperitoneally injected with saline or PGN from Staphylococcus aureus or Bacillus subtilis. Spleen anti-double-stranded DNA (dsDNA) IgG + B cells were sorted for B-cell receptor sequencing. Serum autoantibody levels and kidney damage were analyzed. Further, the association between plasma S. aureus translocation and systemic lupus erythematosus (SLE) pathogenesis was assessed in women. RESULTS: Administration of B. subtilis PGN induced natural-like anti-dsDNA autoantibodies (e.g., IgM, short lived IgG response, and no tissue damage), whereas S. aureus PGN induced pathogenic anti-dsDNA autoantibodies (e.g., prolonged IgG production, low IgM, autoantibody-mediated kidney damage) in C57BL/6 and/or MRL/lpr mice. However, serum total IgG did not differ. S. aureus PGN induced antibodies with reduced clonality and greater hypermutation of IGHV3-74 in splenic anti-dsDNA IgG + B cells from C57BL/6 mice. Further, S. aureus PGN promoted IgG class switch recombination via toll-like receptor 2. Plasma S. aureus DNA levels were increased in women with SLE versus control women and correlated with levels of lupus-related autoantibodies and renal involvement. CONCLUSIONS: S. aureus PGN induces pathogenic autoantibody production, whereas B. subtilis PGN drives production of natural nonpathogenic autoantibodies.

Genetic landscape of Gullah African Americans.

OBJECTIVES: Gullah African Americans are descendants of formerly enslaved Africans living in the Sea Islands along the coast of the southeastern U.S., from North Carolina to Florida. Their relatively high numbers and geographic isolation were conducive to the development and preservation of a unique culture that retains deep African features. Although historical evidence supports a West-Central African ancestry for the Gullah, linguistic and cultural evidence of a connection to Sierra Leone has led to the suggestion of this country/region as their ancestral home. This study sought to elucidate the genetic structure and ancestry of the Gullah. MATERIALS AND METHODS: We leveraged whole-genome genotype data from Gullah, African Americans from Jackson, Mississippi, African populations from Sierra Leone, and population reference panels from Africa and Europe to infer population structure, ancestry proportions, and global estimates of admixture. RESULTS: Relative to non-Gullah African Americans from the Southeast US, the Gullah exhibited higher mean African ancestry, lower European admixture, a similarly small Native American contribution, and increased male-biased European admixture. A slightly tighter bottleneck in the Gullah 13 generations ago suggests a largely shared demographic history with non-Gullah African Americans. Despite a slightly higher relatedness to populations from Sierra Leone, our data demonstrate that the Gullah are genetically related to many West African populations. DISCUSSION: This study confirms that subtle differences in African American population structure exist at finer regional levels. Such observations can help to inform medical genetics research in African Americans, and guide the interpretation of genetic data used by African Americans seeking to explore ancestral identities.

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.

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.

Novel mechanism for estrogen receptor alpha modulation of murine lupus.

Female sex is a risk factor for lupus. Sex hormones, sex chromosomes and hormone receptors are implicated in the pathogenic pathways in lupus. Estrogen receptor alpha (ERα) knockout (KO) mice are used for defining hormone receptor effects in lupus. Prior studies of ERα KO in lupus have conflicting results, likely due to sex hormone levels, different lupus strains and different ERα KO constructs. Our objective was to compare a complete KO of ERα vs. the original functional KO of ERα (expressing a short ERα) on disease expression and immune phenotype, while controlling sex hormone levels. We studied female lupus prone NZM2410 WT and ERα mutant mice. All mice (n = 44) were ovariectomized (OVX) for hormonal control. Groups of each genotype were estrogen (E2)-repleted after OVX. We found that OVXed NZM mice expressing the truncated ERα (ERα short) had significantly reduced nephritis and prolonged survival compared to both wildtype and the complete ERαKO (ERα null) mice, but surprisingly only if E2-repleted. ERα null mice were not protected regardless of E2 status. We observed significant differences in splenic B cells and dendritic cells and a decrease in cDC2 (CD11b+CD8-) dendritic cells, without a concomitant decrease in cDC1 (CD11b-CD8a+) cells comparing ERα short to ERα null or WT mice. Our data support a protective role for the ERα short protein. ERα short is similar to an endogenously expressed ERα variant (ERα46). Modulating its expression/activity represents a potential approach for treating female-predominant autoimmune diseases.

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.

Plasmacytoid dendritic cell distribution and maturation are altered in lupus prone mice prior to the onset of clinical disease.

Plasmacytoid dendritic cells (pDCs) and their production of type I interferons (IFN) are key pathogenic mediators of systemic lupus erythematosus (SLE). Despite the key role of pDCs in SLE, the mechanism by which pDCs promote disease is not well understood. The first objective for this study was to assess the number and maturation state of pDCs in pre-disease NZM2410 lupus prone mice compared to control mice. Second, we sought to identify mechanisms responsible for the alteration in pDCs in NZM mice prior to onset of clinical disease. We compared the number and percent of pDCs in the spleens and bone marrow (BM) of pre-disease NZM24010 (NZM) mice to C57BL/6 (B6) control mice. In the spleens of pre-disease NZM mice, pDC percent and number were increased. This increase occurs in parallel with a decrease in BM pDC number and percent in the NZM mice. The decrease in BM pDC number suggests the increase in spleen pDCs is a result of altered pDC distribution and not increased production of pDCs in the BM. To determine if pDC developmental potential is altered in lupus prone mice, we cultured BM from NZM and B6 mice in vitro. We found a reduced percentage/number of pDCs developing from the BM of NZM mice compared to B6 mice, which further supports that the increase in pDC number is a result of altered pDC distribution rather than increased pDC production. To better characterize the pDC population, we compared the percentage of mature pDCs in the spleens and BM of NZM mice to controls. In the NZM mice, there is a dramatic reduction in the number of mature pDCs in the BM of NZM mice, suggesting that mature pDCs exit the BM at a higher rate/earlier maturation time compared to healthy mice. We conclude that pDCs contribution to disease pathogenesis in NZM mice may include the alteration of pDC distribution to increase the number of pDCs in the spleen prior to disease onset.

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.

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.

Estrogen Receptor α Deficiency Modulates TLR Ligand-Mediated PDC-TREM Expression in Plasmacytoid Dendritic Cells in Lupus-Prone Mice.

Female lupus-prone NZM2410 estrogen receptor α (ERα)-deficient mice are protected from renal disease and have prolonged survival compared with wild-type littermates; however, the mechanism of protection is unknown. Plasmacytoid dendritic cells (pDCs) and type I IFN drive lupus pathogenesis. Estrogen acting via ERα enhances both pDC development and IFN production. The objectives for this study were to determine if ERα modulates pDC function and IFN activity in predisease NZM2410 mice as a possible protective mechanism of ERα deficiency in lupus-prone mice. We measured the effect of ERα deficiency on spleen pDC frequency, number, maturation, and activation state. ERα deficiency reduced type I IFN activity and the frequency of MHC class II(+) pDCs in the spleen without altering overall pDC frequency, number, or maturation state. Additionally, ERα-deficient NZM2410 mice had a significantly decreased frequency of pDCs expressing PDC-TREM, a modulator of TLR-mediated IFN production. After in vitro TLR9 stimulation, ERα deficiency significantly reduced the expression of PDC-TREM on pDCs from both NZM2410 and C57BL/6 mice. Thus, we have identified a significant effect of ERα deficiency on pDCs in predisease NZM2410 mice, which may represent a mechanism by which ERα deficiency protects NZM2410 mice from lupuslike disease.

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.

Differential efficacy of human mesenchymal stem cells based on source of origin.

Mesenchymal stem cells (MSCs) are useful in tissue repair but also possess immunomodulatory properties. Murine and uncontrolled human trials suggest efficacy of MSCs in treating lupus. Autologous cells are preferable; however, recent studies suggest that lupus-derived MSCs lack efficacy in treating disease. Thus, the optimum derivation of MSCs for use in lupus is unknown. It is also unknown which in vitro assays of MSC function predict in vivo efficacy. The objectives for this study were to provide insight into the optimum source of MSCs and to identify in vitro assays that predict in vivo efficacy. We derived MSCs from four umbilical cords, four healthy bone marrows (BMs), and four lupus BMs. In diseased MRL/lpr mice, MSCs from healthy BM and umbilical cords significantly decreased renal disease, whereas lupus BM MSCs only delayed disease. Current in vitro assays did not differentiate efficacy of the different MSCs. However, differences in MSC efficacy were observed in B cell proliferation assays. Our results suggest that autologous MSCs from lupus patients are not effective in treating disease. Furthermore, standard in vitro assays for MSC licensing are not predictive of in vivo efficacy, whereas inhibiting B cell proliferation appears to differentiate effective MSCs from ineffective MSCs.

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.

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.

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.

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.