Short- and mid-term outcomes in systemic lupus erythematosus patients presenting with disease exacerbation after SARS-CoV-2 mRNA vaccination: A cohort study from Puerto Rico.

OBJECTIVE: To determine if SARS-CoV-2 mRNA vaccination has an impact on the clinical course of systemic lupus erythematosus (SLE). METHODS: Puerto Ricans with SLE who received mRNA COVID-19 vaccines were studied. Demographic parameters, clinical manifestations, disease activity (per Systemic Lupus Erythematosus Disease Activity Index (SLEDAI), disease damage (per Systemic Lupus International Collaborating Clinics/American College of Rheumatology Damage Index), emergency room visits, hospitalizations, and pharmacologic therapy were determined. Baseline variables (prior to vaccination) were compared between patients with and without exacerbation after SARS-CoV-2 vaccination. Among those with exacerbation, clinical outcomes were determined up to 1 year after vaccination. RESULTS: Of the entire cohort (n = 247), 14 (5.7%) had post-vaccination exacerbations. Photosensitivity, oral ulcers, anti-Ro antibodies, higher SLEDAI score, and corticosteroids exposure were associated with post-vaccination flares. Among those with post-vaccination flares, 10 (71.4%) had major organ involvement. No significant differences were observed for mean SLEDAI scores, emergency room visits, hospitalizations, disease damage, and exposure to immunosuppressive drugs before and after SARS-CoV-2 mRNA vaccination. At 12 months of follow-up, all patients were fully controlled without evidence of active disease. CONCLUSION: In our group of SLE patients, 5.7% had a disease flare after SARS-CoV-2 mRNA vaccination. Most had exacerbations involving major organs/systems. Mucocutaneous manifestations, anti-Ro antibodies, disease activity, and corticosteroids were associated with flares. Awareness of these factors and the possibility of a major lupus flare after vaccination with COVD-19 vaccines is critical to provide timely and effective therapy.

Does hydroxychloroquine improve patient reported outcomes in patients with lupus?

BACKGROUND: Hydroxychloroquine (HCQ) use is associated with less disease activity, flares, damage and improved survival in Systemic Lupus Erythematosus (SLE). However, its effect on patient reported health outcomes (PROs) such as quality of life (QOL) is not known. METHODS: International data from Study on Outcomes of Lupus (SOUL) from 2,161 SLE patients were compared by HCQ use. Disease activity and damage were assessed using SELENA-SLEDAI and SLICC-ACR/SDI. QOL was evaluated using LupusPRO and Lupus Impact Tracker (LIT). Linear regression analyses were performed with LupusPRO summary scores health related HRQOL, non-health related NHRQOL and LIT as dependent and HCQ use as independent variable. Analyses were undertaken to test mediation of effects of HCQ use on QOL through disease activity. RESULTS: Mean age was 40.5 ± 12.8 years, 93% were women. Sixty-three (1363/2161) percent were on HCQ. On univariate analysis, HCQ use was associated with (a) better QOL (LupusPRO-HRQOL: ß 6.19, 95% CI 4.15, 8.24, P ≤ 0.001, LupusPRO NHRQOL: ß 5.83, 95% CI 4.02, 7.64, P ≤ 0.001) and less impact on daily life (LIT: ß -9.37, 95% CI -12.24, -6.50, P ≤ 0.001). On multivariate and mediational analyses, the effects of HCQ on QOL were indirectly and completely mediated through disease activity. CONCLUSIONS: HCQ use in SLE is associated with better patient reported health outcomes (LupusPRO-HRQOL and NHRQOL and impact on daily life), and the effects are mediated through disease activity. This information can facilitate patients and physician's communication with decision-making regarding the use of HCQ for SLE management.

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.

Factors Associated With Disease Damage in Puerto Ricans With Primary Sjögren Syndrome.

OBJECTIVE: The aim of this study was to determine the association of demographic parameters, clinical manifestations, disease activity, and pharmacologic therapy with disease damage in a group of Puerto Ricans with primary Sjögren syndrome (pSS). METHODS: A cross-sectional study was conducted in 100 Hispanics of Puerto Rico with pSS. Patients were 21 years or older and fulfilled the 2012 American College of Rheumatology classification criteria for pSS. Demographic factors, lifestyle behaviors, extraglandular manifestations, serologic tests, comorbidities, pharmacologic therapy, disease activity (per European League Against Rheumatism Sjögren Syndrome Disease Activity Index), and disease damage (per Sjögren Syndrome Disease Damage Index [SSDDI]) were assessed. Patients with disease damage (SSDDI ≥1) and without damage (SSDDI = 0) were compared using bivariate analysis and multivariate regression analysis adjusted for age, sex, and disease duration. RESULTS: The mean age of patients was 52.8 years; 94% were women. The mean disease duration was 5.9 years. Thirty-nine patients had disease damage. Disease damage was mainly attributed to pulmonary fibrosis and peripheral neuropathy. In the bivariate analysis, disease damage was associated with low C3 and C4, coronary artery disease, infections, and higher activity index and was more frequently treated corticosteroids and azathioprine. In the multivariate analysis, low C3, disease activity, and corticosteroid exposure retained significance. CONCLUSIONS: In this population of Puerto Ricans with pSS, C3 and C4 hypocomplementemia, coronary artery disease, infections, and exposure to corticosteroids and azathioprine were associated with damage accrual. Clinicians should be aware of these factors to identify those who may require close follow-up and early therapeutic intervention.

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.

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.

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.