Genome-wide association meta-analysis in Chinese and European individuals identifies ten new loci associated with systemic lupus erythematosus

High-throughput identification of functional regulatory SNPs in systemic lupus erythematosus

Genome-wide association studies implicate multiple loci in risk for systemic lupus erythematosus (SLE), but few contain exonic variants, rendering systematic identification of non-coding variants essential to decoding SLE genetics. We utilized SNP-seq and bioinformatic enrichment to interrogate 2180 single-nucleotide polymorphisms (SNPs) from 87 SLE risk loci for potential binding of transcription factors and related proteins from B cells. 52 SNPs that passed initial screening were tested by electrophoretic mobility shift and luciferase reporter assays. To validate the approach, we studied rs2297550 in detail, finding that the risk allele enhanced binding to the transcription factor Ikaros (encoded by IKZF1), thereby modulating expression of IKBKE. Correspondingly, primary cells from genotyped healthy donors bearing the risk allele expressed higher levels of the interferon / NF-κB regulator IKKε. Together, these findings define a set of likely functional non-coding lupus risk variants and identify a regulatory pathway involving rs2297550, Ikaros, and IKKε implicated by human genetics in risk for SLE.

SPATS2L is a positive feedback regulator of the type I interferon signaling pathway and plays a vital role in lupus

Through genome-wide association studies (GWAS) and integrated expression quantitative trait locus (eQTL) analyses, numerous susceptibility genes ("eGenes", whose expressions are significantly associated with common variants) associated with systemic lupus erythematosus (SLE) have been identified. Notably, a subset of these eGenes is correlated with disease activity. However, the precise mechanisms through which these genes contribute to the initiation and progression of the disease remain to be fully elucidated. In this investigation, we initially identify SPATS2L as an SLE eGene correlated with disease activity. eSignaling and transcriptomic analyses suggest its involvement in the type I interferon (IFN) pathway. We observe a significant increase in SPATS2L expression following type I IFN stimulation, and the expression levels are dependent on both the concentration and duration of stimulation. Furthermore, through dual-luciferase reporter assays, western blot analysis, and imaging flow cytometry, we confirm that SPATS2L positively modulates the type I IFN pathway, acting as a positive feedback regulator. Notably, siRNA-mediated intervention targeting SPATS2L, an interferon-inducible gene, in peripheral blood mononuclear cells (PBMCs) from patients with SLE reverses the activation of the interferon pathway. In conclusion, our research highlights the pivotal role of SPATS2L as a positive-feedback regulatory molecule within the type I IFN pathway. Our findings suggest that SPATS2L plays a critical role in the onset and progression of SLE and may serve as a promising target for disease activity assessment and intervention strategies.

BACH2: The Future of Induced T-Regulatory Cell Therapies

BACH2 (BTB Domain and CNC Homolog 2) is a transcription factor that serves as a central regulator of immune cell differentiation and function, particularly in T and B lymphocytes. A picture is emerging that BACH2 may function as a master regulator of cell fate that is exquisitely sensitive to cell activation status. In particular, BACH2 plays a key role in stabilizing the phenotype and suppressive function of transforming growth factor-beta (TGF-β)-derived human forkhead box protein P3 (FOXP3)+ inducible regulatory T cells (iTregs), a cell type that holds great clinical potential as a cell therapeutic for diverse inflammatory conditions. As such, BACH2 potentially could be targeted to overcome the instability of the iTreg phenotype and suppressive function that has hampered their clinical application. In this review, we focus on the role of BACH2 in T cell fate and iTreg function and stability. We suggest approaches to modulate BACH2 function that may lead to more stable and efficacious Treg cell therapies.

Polygenic Risk Score Predicts Earlier-Onset Adult Systemic Lupus Erythematosus and First-Year Renal Diseases in a Taiwanese Cohort

OBJECTIVES

This study aimed to develop a predictive model using polygenic risk score (PRS) to forecast renal outcomes for adult systemic lupus erythematosus (SLE) in a Taiwanese population.

METHODS

Patients with SLE (n=2782) and matched non-SLE controls (n=11 128) were genotyped using Genome-Wide TWB 2.0 single-nucleotide polymorphism (SNP) array. PRS models (C+T, LDpred2, Lassosum, PRSice-2, PRS-continuous shrinkage (CS)) were constructed for predicting SLE susceptibility. Logistic regression was assessed for C+T-based PRS association with renal involvement in patients with SLE.

RESULTS

In the training set, C+T-based SLE-PRS, only incorporating 27 SNPs, outperformed other models with area under the curve (AUC) values of 0.629, surpassing Lassosum (AUC=0.621), PRSice-2 (AUC=0.615), LDpred2 (AUC=0.609) and PRS-CS (AUC=0.602). Additionally, C+T-based SLE-PRS demonstrated consistent predictive capacity in the testing set (AUC=0.620). Individuals in the highest quartile exhibited earlier SLE onset (39.06 vs 44.22 years, p<0.01), higher Systemic Lupus Erythematosus Disease Activity Index scores (3.00 vs 2.37, p=0.04), elevated risks of renal involvement within the first year of SLE diagnosis, including WHO class III-IV lupus nephritis (OR 2.36, 95% CI 1.47 to 3.80, p<0.01), estimated glomerular filtration rate <60 mL/min/1.73m2 (OR 1.49, 95% CI 1.18 to 1.89, p<0.01) and urine protein-to-creatinine ratio >150 mg/day (OR 2.07, 95% CI 1.49 to 2.89, p<0.01), along with increased seropositivity risks, compared with those in the lowest quartile. Furthermore, among patients with SLE with onset before 50 years, the highest PRS quartile was significantly associated with more serious renal diseases within the first year of SLE diagnosis.

CONCLUSIONS

PRS of SLE is associated with earlier onset, renal involvement within the first year of SLE diagnosis and seropositivity in Taiwanese patients. Integrating PRS with clinical decision-making may enhance lupus nephritis screening and early treatment to improve renal outcomes in patients with SLE.

High-throughput identification of functional regulatory SNPs in systemic lupus erythematosus

Genome-wide association studies implicate multiple loci in risk for systemic lupus erythematosus (SLE), but few contain exonic variants, rendering systematic identification of non-coding variants essential to decoding SLE genetics. We utilized SNP-seq and bioinformatic enrichment to interrogate 2180 single-nucleotide polymorphisms (SNPs) from 87 SLE risk loci for potential binding of transcription factors and related proteins from B cells. 52 SNPs that passed initial screening were tested by electrophoretic mobility shift and luciferase reporter assays. To validate the approach, we studied rs2297550 in detail, finding that the risk allele enhanced binding to the transcription factor Ikaros (IKZF1), thereby modulating expression of IKBKE. Correspondingly, primary cells from genotyped healthy donors bearing the risk allele expressed higher levels of the interferon / NF-κB regulator IKKϵ. Together, these findings define a set of likely functional non-coding lupus risk variants and identify a new regulatory pathway involving rs2297550, Ikaros, and IKKϵ implicated by human genetics in risk for SLE.

Massively parallel reporter assay confirms regulatory potential of hQTLs and reveals important variants in lupus and other autoimmune diseases

We designed a massively parallel reporter assay (MPRA) in an Epstein-Barr virus transformed B cell line to directly characterize the potential for histone post-translational modifications, i.e., histone quantitative trait loci (hQTLs), expression QTLs (eQTLs), and variants on systemic lupus erythematosus (SLE) and autoimmune (AI) disease risk haplotypes to modulate regulatory activity in an allele-dependent manner. Our study demonstrates that hQTLs, as a group, are more likely to modulate regulatory activity in an MPRA compared with other variant classes tested, including a set of eQTLs previously shown to interact with hQTLs and tested AI risk variants. In addition, we nominate 17 variants (including 11 previously unreported) as putative causal variants for SLE and another 14 for various other AI diseases, prioritizing these variants for future functional studies in primary and immortalized B cells. Thus, we uncover important insights into the mechanistic relationships among genotype, epigenetics, and gene expression in SLE and AI disease phenotypes.

Rare SH2B3 coding variants in lupus patients impair B cell tolerance and predispose to autoimmunity

Systemic lupus erythematosus (SLE) is a heterogeneous autoimmune disease with a clear genetic component. While most SLE patients carry rare gene variants in lupus risk genes, little is known about their contribution to disease pathogenesis. Amongst them, SH2B3-a negative regulator of cytokine and growth factor receptor signaling-harbors rare coding variants in over 5% of SLE patients. Here, we show that unlike the variant found exclusively in healthy controls, SH2B3 rare variants found in lupus patients are predominantly hypomorphic alleles, failing to suppress IFNGR signaling via JAK2-STAT1. The generation of two mouse lines carrying patients' variants revealed that SH2B3 is important in limiting the number of immature and transitional B cells. Furthermore, hypomorphic SH2B3 was shown to impair the negative selection of immature/transitional self-reactive B cells and accelerate autoimmunity in sensitized mice, at least in part due to increased IL-4R signaling and BAFF-R expression. This work identifies a previously unappreciated role for SH2B3 in human B cell tolerance and lupus risk.

Bayesian approach to assessing population differences in genetic risk of disease with application to prostate cancer

Population differences in risk of disease are common, but the potential genetic basis for these differences is not well understood. A standard approach is to compare genetic risk across populations by testing for mean differences in polygenic scores, but existing studies that use this approach do not account for statistical noise in effect estimates (i.e., the GWAS betas) that arise due to the finite sample size of GWAS training data. Here, we show using Bayesian polygenic score methods that the level of uncertainty in estimates of genetic risk differences across populations is highly dependent on the GWAS training sample size, the polygenicity (number of causal variants), and genetic distance (FST) between the populations considered. We derive a Wald test for formally assessing the difference in genetic risk across populations, which we show to have calibrated type 1 error rates under a simplified assumption that all SNPs are independent, which we achieve in practise using linkage disequilibrium (LD) pruning. We further provide closed-form expressions for assessing the uncertainty in estimates of relative genetic risk across populations under the special case of an infinitesimal genetic architecture. We suggest that for many complex traits and diseases, particularly those with more polygenic architectures, current GWAS sample sizes are insufficient to detect moderate differences in genetic risk across populations, though more substantial differences in relative genetic risk (relative risk > 1.5) can be detected. We show that conventional approaches that do not account for sampling error from the training sample, such as using a simple t-test, have very high type 1 error rates. When applying our approach to prostate cancer, we demonstrate a higher genetic risk in African Ancestry men, with lower risk in men of European followed by East Asian ancestry.

Non-coding DNA variants for risk in lupus

Systemic Lupus Erythematosus (SLE) is a multifactorial autoimmune disease that arises from a dynamic interplay between genetics and environmental triggers. The advent of sophisticated genomics technology has catalyzed a shift in our understanding of disease etiology, spotlighting the pivotal role of non-coding DNA variants in SLE pathogenesis. In this review, we present a comprehensive examination of the non-coding variants associated with SLE, shedding light on their role in influencing disease risk and progression. We discuss the latest methodological advancements that have been instrumental in the identification and functional characterization of these genomic elements, with a special focus on the transformative power of CRISPR-based gene-editing technologies. Additionally, the review probes into the therapeutic opportunities that arise from modulating non-coding regions associated with SLE. Through an exploration of the complex network of non-coding DNA, this review aspires to decode the genetic puzzle of SLE and set the stage for groundbreaking gene-based therapeutic interventions and the advancement of precision medicine strategies tailored to SLE management.

Relationship between CD274 gene polymorphism and systemic lupus erythematosus risk in a Chinese Han population

OBJECTIVE

Relationship between surface antigen differentiation cluster 274 (CD274) gene polymorphism and systemic lupus erythematosus (SLE) risk is limited. This study aims to discuss whether in a Chinese Han population, CD274 gene polymorphisms may relate to SLE susceptibility.

METHODS

Three hundred and ten SLE patients and 390 healthy controls were included in this case-control study. Using the Kompetitive Allele-Specific PCR (KASP) approach, five single nucleotide polymorphisms (SNPs), including rs2890658, rs4143815, rs822339, rs2282055, and rs2297137, were genotyped for CD274 gene polymorphisms. Correlation between the polymorphisms and clinical, laboratory features in SLE patients were discussed.

RESULTS

Frequency of C allele was substantially lower in SLE patients than in healthy controls (p = .015), and CC genotype was significantly negatively related to developing SLE at locus rs4143815 (p = .013). At locus rs822339, frequency of GA genotype was higher than that of the healthy controls (p = .006). At locus rs2282055, frequency of GG genotype was lower than that of healthy controls (p = .024). According to subgroup analysis, the CD274 gene polymorphisms rs2890658, rs4143815, rs822339, rs2282055, and rs2297137 were partly linked to some clinical symptoms of SLE patients, such as Complement 4 (C4), C-reactive protein (CRP), and erythrocyte sedimentation rate (ESR).

CONCLUSION

CD274 gene polymorphisms may be susceptible to SLE in the Chinese Han people.

Association between systemic lupus erythematosus and inflammatory bowel disease in European and East Asian populations: a two-sample Mendelian randomization study

Background

Previous studies have shown a coexistence phenomenon between systemic lupus erythematosus (SLE) and inflammatory bowel disease (IBD), but the causal relationship between them is still unclear. Therefore, we conducted a two-sample Mendelian randomization (MR) analysis using publicly available summary statistics data to evaluate whether there was a causal relationship between the two diseases.

Methods

Summary statistics for SLE and IBD were downloaded from the Open Genome-Wide Association Study and the International Inflammatory Bowel Disease Genetics Consortium. European and East Asian populations were included in this MR work. We adopted a series of methods to select instrumental variables that are closely related to SLE and IBD. To make the conclusion more reliable, we applied a variety of different analysis methods, among which the inverse variance-weighted (IVW) method was the main method. In addition, heterogeneity, pleiotropy, and sensitivity were assessed to make the conclusions more convincing.

Results

In the European population, a negative causal relationship was observed between SLE and overall IBD (OR = 0.94; 95% CI = 0.90, 0.98; P < 0.004) and ulcerative colitis (UC) (OR = 0.93; 95% CI = 0.88, 0.98; P = 0.006). After removing outliers with Mendelian Randomization Pleiotropy RESidual Sum and Outlier (MR-PRESSO), the results remained consistent with IVW. However, there was no causal relationship between SLE and Crohn's disease. In the East Asian population, no causal relationship was found between SLE and IBD.

Conclusion

Our results found that genetic susceptibility to SLE was associated with lower overall IBD risk and UC risk in European populations. In contrast, no association between SLE and IBD was found in East Asian populations. This work might enrich the previous research results, and it may provide some references for research in the future.

Genetic examination of hematological parameters in SARS-CoV-2 infection and COVID-19

People hospitalized with COVID-19 often exhibit altered hematological traits associated with disease prognosis (e.g., lower lymphocyte and platelet counts). We investigated whether inter-individual variability in baseline hematological traits influences risk of acute SARS-CoV-2 infection or progression to severe COVID-19. We report inconsistent associations between blood cell traits with incident SARS-CoV-2 infection and severe COVID-19 in UK Biobank and the Vanderbilt University Medical Center Synthetic Derivative (VUMC SD). Since genetically determined blood cell measures better represent cell abundance across the lifecourse, we also assessed the shared genetic architecture of baseline blood cell traits on COVID-19 related outcomes by Mendelian randomization (MR) analyses. We found significant relationships between COVID-19 severity and mean sphered cell volume after adjusting for multiple testing. However, MR results differed significantly across different freezes of COVID-19 summary statistics and genetic correlation between these traits was modest (0.1), decreasing our confidence in these results. We observed overlapping genetic association signals between other hematological and COVID-19 traits at specific loci such as MAPT and TYK2. In conclusion, we did not find convincing evidence of relationships between the genetic architecture of blood cell traits and either SARS-CoV-2 infection or COVID-19 hospitalization, though we do see evidence of shared signals at specific loci.

An update on the roles of transcription factor Ets1 in autoimmune diseases

Transcription factors are crucial to regulate gene expression in immune cells and in other cell types. In lymphocytes, there are a large number of different transcription factors that are known to contribute to cell differentiation and the balance between quiescence and activation. One such transcription factor is E26 oncogene homolog 1 (Ets1). Ets1 expression is high in quiescent B and T lymphocytes and its levels are decreased upon activation. The human ETS1 gene has been identified as a susceptibility locus for many autoimmune and inflammatory diseases. In accord with this, gene knockout of Ets1 in mice leads to development of a lupus-like autoimmune disease, with enhanced activation and differentiation of both B cells and T cells. Prior reviews have summarized functional roles for Ets1 based on studies of Ets1 knockout mice. In recent years, numerous additional studies have been published that further validate ETS1 as a susceptibility locus for human diseases where immune dysregulation plays a causative role. In this update, new information that further links Ets1 to human autoimmune diseases is organized and collated to serve as a resource. This update also describes recent studies that seek to understand molecularly how Ets1 regulates immune cell activation, either using human cells and tissues or mouse models. This resource is expected to be useful to investigators seeking to understand how Ets1 may regulate the human immune response, particularly in terms of its roles in autoimmunity and inflammation. This article is categorized under: Immune System Diseases > Genetics/Genomics/Epigenetics Immune System Diseases > Molecular and Cellular Physiology.

Systemic lupus erythematosus is associated with lower risk of hepatitis B virus infection: A multivariable Mendelian randomization study in East Asian population

The relationship between systemic lupus erythematosus (SLE) and hepatitis B virus (HBV) infection is still unclear. We conducted a two-sample Mendelian randomization (MR) analysis using summary statistics from genome-wide association studies for SLE and HBV infection in individuals of East Asian ancestry. The inverse-variance weighted (IVW) method, weighted median (WM) method, and MR-Egger method were used to estimate the causal effect of SLE on HBV infection. Additionally, we performed a multivariable MR analysis adjusting for the effects of body mass index and rheumatoid arthritis. This MR study included a total of 225 106 individuals of East Asian ancestry, comprising 5616 cases and 219 490 controls. The IVW method (OR: 0.79, p = 3.34E-08) and the WM method (OR: 0.79, p = 9.09E-06) revealed a causal relationship between genetically predicted SLE and a low risk of HBV infection. The multivariable MR analysis still suggested a low risk of HBV infection associated with SLE (OR: 0.83, p = 2.89E-06). Our MR analysis supports a causal relationship between SLE and a low risk of HBV infection in individuals of East Asian ancestry.

SOCS-JAK-STAT inhibitors and SOCS mimetics as treatment options for autoimmune uveitis, psoriasis, lupus, and autoimmune encephalitis

Autoimmune diseases arise from atypical immune responses that attack self-tissue epitopes, and their development is intricately connected to the disruption of the JAK-STAT signaling pathway, where SOCS proteins play crucial roles. Conditions such as autoimmune uveitis, psoriasis, lupus, and autoimmune encephalitis exhibit immune system dysfunctions associated with JAK-STAT signaling dysregulation. Emerging therapeutic strategies utilize JAK-STAT inhibitors and SOCS mimetics to modulate immune responses and alleviate autoimmune manifestations. Although more research and clinical studies are required to assess their effectiveness, safety profiles, and potential for personalized therapeutic approaches in autoimmune conditions, JAK-STAT inhibitors and SOCS mimetics show promise as potential treatment options. This review explores the action, effectiveness, safety profiles, and future prospects of JAK inhibitors and SOCS mimetics as therapeutic agents for psoriasis, autoimmune uveitis, systemic lupus erythematosus, and autoimmune encephalitis. The findings underscore the importance of investigating these targeted therapies to advance treatment options for individuals suffering from autoimmune diseases.

B cell polygenic risk scores associate with anti-dsDNA antibodies and nephritis in systemic lupus erythematosus

OBJECTIVE

B cell function and autoantibodies are important in SLE pathogenesis. In this work, we aimed to investigate the impact of cumulative SLE B cell genetics on SLE subphenotype and autoantibody profile.

METHODS

Female patients with SLE (n=1248) and healthy controls (n=400) were genotyped using Illumina's Global Screening Array. Two polygenic risk scores (PRSs), one representing B cell genes and the other B cell activation genes, were calculated for each individual using risk loci for SLE in genes assigned to B cell-related pathways according to the Kyoto Encyclopedia of Genes and Genomes, Gene Ontology and Reactome Databases.

RESULTS

Double-stranded DNA (dsDNA) antibodies were more prevalent among patients with a high compared with a low SLE B cell PRS (OR 1.47 (1.07 to 2.01), p=0.018), and effect sizes were augmented in patients with human leucocyte antigen (HLA) risk haplotypes HLA-DRB1*03:01 and HLA-DRB1*15:01 (DRB1*03/15 -/- (OR 0.99 (0.56 to 1.77), p=0.98; DRB1*03/15 +/- or -/+ (OR 1.64 (1.06 to 2.54), p=0.028; and DRB1*03/15 +/+ (OR 4.47 (1.21 to 16.47), p=0.024). Further, a high compared with a low B cell PRS was associated with low complement levels in DRB1*03/15 +/+ patients (OR 3.92 (1.22 to 12.64), p=0.022). The prevalence of lupus nephritis (LN) was higher in patients with a B cell activation PRS above the third quartile compared with patients below (OR 1.32 (1.00 to 1.74), p=0.048).

CONCLUSIONS

High genetic burden related to B cell function is associated with dsDNA antibody development and LN. Assessing B cell PRSs may be important in order to determine immunological pathways influencing SLE and to predict clinical phenotype.

A multilayered post-GWAS analysis pipeline defines functional variants and target genes for systemic lupus erythematosus (SLE)

Objectives

Systemic lupus erythematosus (SLE), an autoimmune disease with incompletely understood etiology, has a strong genetic component. Although genome-wide association studies (GWAS) have revealed multiple SLE susceptibility loci and associated single nucleotide polymorphisms (SNPs), the precise causal variants, target genes, cell types, tissues, and mechanisms of action remain largely unknown.

Methods

Here, we report a comprehensive post-GWAS analysis using extensive bioinformatics, molecular modeling, and integrative functional genomic and epigenomic analyses to optimize fine-mapping. We compile and cross-reference immune cell-specific expression quantitative trait loci ( cis - and trans -eQTLs) with promoter-capture Hi-C, allele-specific chromatin accessibility, and massively parallel reporter assay data to define predisposing variants and target genes. We experimentally validate a predicted locus using CRISPR/Cas9 genome editing, qPCR, and Western blot.

Results

Anchoring on 452 index SNPs, we selected 9,931 high-linkage disequilibrium (r 2 >0.8) SNPs and defined 182 independent non-HLA SLE loci. 3,746 SNPs from 143 loci were identified as regulating 564 unique genes. Target genes are enriched in lupus-related tissues and associated with other autoimmune diseases. Of these, 329 SNPs (106 loci) showed significant allele-specific chromatin accessibility and/or enhancer activity, indicating regulatory potential. Using CRISPR/Cas9, we validated rs57668933 as a functional variant regulating multiple targets, including SLE risk gene ELF1 , in B-cells.

Conclusion

We demonstrate and validate post-GWAS strategies for utilizing multi-dimensional data to prioritize likely causal variants with cognate gene targets underlying SLE pathogenesis. Our results provide a catalog of significantly SLE-associated SNPs and loci, target genes, and likely biochemical mechanisms, to guide experimental characterization.

Adaptive MCMC for Bayesian Variable Selection in Generalised Linear Models and Survival Models

Developing an efficient computational scheme for high-dimensional Bayesian variable selection in generalised linear models and survival models has always been a challenging problem due to the absence of closed-form solutions to the marginal likelihood. The Reversible Jump Markov Chain Monte Carlo (RJMCMC) approach can be employed to jointly sample models and coefficients, but the effective design of the trans-dimensional jumps of RJMCMC can be challenging, making it hard to implement. Alternatively, the marginal likelihood can be derived conditional on latent variables using a data-augmentation scheme (e.g., Pólya-gamma data augmentation for logistic regression) or using other estimation methods. However, suitable data-augmentation schemes are not available for every generalised linear model and survival model, and estimating the marginal likelihood using a Laplace approximation or a correlated pseudo-marginal method can be computationally expensive. In this paper, three main contributions are presented. Firstly, we present an extended Point-wise implementation of Adaptive Random Neighbourhood Informed proposal (PARNI) to efficiently sample models directly from the marginal posterior distributions of generalised linear models and survival models. Secondly, in light of the recently proposed approximate Laplace approximation, we describe an efficient and accurate estimation method for marginal likelihood that involves adaptive parameters. Additionally, we describe a new method to adapt the algorithmic tuning parameters of the PARNI proposal by replacing Rao-Blackwellised estimates with the combination of a warm-start estimate and the ergodic average. We present numerous numerical results from simulated data and eight high-dimensional genetic mapping data-sets to showcase the efficiency of the novel PARNI proposal compared with the baseline add-delete-swap proposal.

Antibody Cross-Reactivity in Auto-Immune Diseases

Autoimmunity is defined by the presence of antibodies and/or T cells directed against self-components. Although of unknown etiology, autoimmunity commonly is associated with environmental factors such as infections, which have been reported to increase the risk of developing autoimmune diseases. Occasionally, similarities between infectious non-self and self-tissue antigens may contribute to immunological cross-reactivity in autoimmune diseases. These reactions may be interpreted as molecular mimicry, which describes cross-reactivity between foreign pathogens and self-antigens that have been reported to cause tissue damage and to contribute to the development of autoimmunity. By focusing on the nature of antibodies, cross-reactivity in general, and antibody-antigen interactions, this review aims to characterize the nature of potential cross-reactive immune reactions between infectious non-self and self-tissue antigens which may be associated with autoimmunity but may not actually be the cause of disease onset.

Phenotype Risk Score but Not Genetic Risk Score Aids in Identifying Individuals With Systemic Lupus Erythematosus in the Electronic Health Record

OBJECTIVE

Systemic lupus erythematosus (SLE) poses diagnostic challenges. We undertook this study to evaluate the utility of a phenotype risk score (PheRS) and a genetic risk score (GRS) to identify SLE individuals in a real-world setting.

METHODS

Using a de-identified electronic health record (EHR) database with an associated DNA biobank, we identified 789 SLE cases and 2,261 controls with available MEGAEX genotyping. A PheRS for SLE was developed using billing codes that captured American College of Rheumatology SLE criteria. We developed a GRS with 58 SLE risk single-nucleotide polymorphisms (SNPs).

RESULTS

SLE cases had a significantly higher PheRS (mean ± SD 7.7 ± 8.0 versus 0.8 ± 2.0 in controls; P < 0.001) and GRS (mean ± SD 12.2 ± 2.3 versus 11.0 ± 2.0 in controls; P < 0.001). Black individuals with SLE had a higher PheRS compared to White individuals (mean ± SD 10.0 ± 10.1 versus 7.1 ± 7.2, respectively; P = 0.002) but a lower GRS (mean ± SD 9.0 ± 1.4 versus 12.3 ± 1.7, respectively; P < 0.001). Models predicting SLE that used only the PheRS had an area under the curve (AUC) of 0.87. Adding the GRS to the PheRS resulted in a minimal difference with an AUC of 0.89. On chart review, controls with the highest PheRS and GRS had undiagnosed SLE.

CONCLUSION

We developed a SLE PheRS to identify established and undiagnosed SLE individuals. A SLE GRS using known risk SNPs did not add value beyond the PheRS and was of limited utility in Black individuals with SLE. More work is needed to understand the genetic risks of SLE in diverse populations.

Genetics of SLE: mechanistic insights from monogenic disease and disease-associated variants

The past few years have provided important insights into the genetic architecture of systemic autoimmunity through aggregation of findings from genome-wide association studies (GWAS) and whole-exome or whole-genome sequencing studies. In the prototypic systemic autoimmune disease systemic lupus erythematosus (SLE), monogenic disease accounts for a small fraction of cases but has been instrumental in the elucidation of disease mechanisms. Defects in the clearance or digestion of extracellular or intracellular DNA or RNA lead to increased sensing of nucleic acids, which can break B cell tolerance and induce the production of type I interferons leading to tissue damage. Current data suggest that multiple GWAS SLE risk alleles act in concert with rare functional variants to promote SLE development. Moreover, introduction of orthologous variant alleles into mice has revealed that pathogenic X-linked dominant and recessive SLE can be caused by novel variants in TLR7 and SAT1, respectively. Such bespoke models of disease help to unravel pathogenic pathways and can be used to test targeted therapies. Cell type-specific expression data revealed that most GWAS SLE risk genes are highly expressed in age-associated B cells (ABCs), which supports the view that ABCs produce lupus autoantibodies and contribute to end-organ damage by persisting in inflamed tissues, including the kidneys. ABCs have thus emerged as key targets of promising precision therapeutics.

Massively Parallel Reporter Assay Confirms Regulatory Potential of hQTLs and Reveals Important Variants in Lupus and Other Autoimmune Diseases

Objective

To systematically characterize the potential for histone post-translational modifications, i.e., histone quantitative trait loci (hQTLs), expression QTLs (eQTLs), and variants on systemic lupus erythematosus (SLE) and autoimmune (AI) disease risk haplotypes to modulate gene expression in an allele dependent manner.

Methods

We designed a massively parallel reporter assay (MPRA) containing ~32K variants and transfected it into an Epstein-Barr virus transformed B cell line generated from an SLE case.

Results

Our study expands our understanding of hQTLs, illustrating that epigenetic QTLs are more likely to contribute to functional mechanisms than eQTLs and other variant types, and a large proportion of hQTLs overlap transcription start sites (TSS) of noncoding RNAs. In addition, we nominate 17 variants (including 11 novel) as putative causal variants for SLE and another 14 for various other AI diseases, prioritizing these variants for future functional studies primary and immortalized B cells.

Conclusion

We uncover important insights into the mechanistic relationships between genotype, epigenetics, gene expression, and SLE and AI disease phenotypes.

Modulation of the immunity and inflammation by autophagy

Autophagy, a highly conserved cellular self-degradation pathway, has emerged with novel roles in the realms of immunity and inflammation. Genome-wide association studies have unveiled a correlation between genetic variations in autophagy-related genes and heightened susceptibility to autoimmune and inflammatory diseases. Subsequently, substantial progress has been made in unraveling the intricate involvement of autophagy in immunity and inflammation through functional studies. The autophagy pathway plays a crucial role in both innate and adaptive immunity, encompassing various key functions such as pathogen clearance, antigen processing and presentation, cytokine production, and lymphocyte differentiation and survival. Recent research has identified novel approaches in which the autophagy pathway and its associated proteins modulate the immune response, including noncanonical autophagy. This review provides an overview of the latest advancements in understanding the regulation of immunity and inflammation through autophagy. It summarizes the genetic associations between variants in autophagy-related genes and a range of autoimmune and inflammatory diseases, while also examining studies utilizing transgenic animal models to uncover the in vivo functions of autophagy. Furthermore, the review delves into the mechanisms by which autophagy dysregulation contributes to the development of three common autoimmune and inflammatory diseases and highlights the potential for autophagy-targeted therapies.

A review of genetic risk in systemic lupus erythematosus

INTRODUCTION

Systemic Lupus Erythematosus (SLE) is a complex multisystem autoimmune disease with a wide range of signs and symptoms in affected individuals. The utilization of genome-wide association study (GWAS) technology has led to an explosion in the number of genetic risk factors mapped for autoimmune diseases, including SLE.

AREAS COVERED

In this review, we summarize the more recent genetic risk loci mapped in SLE, which bring the total number of loci mapped to approximately 200. We review prioritization analyses of the associated variants and experimental validation of the putative causal variants. This includes the implementation of new bioinformatic techniques to align genomic and functional data and the use of transcriptomics with single-cell RNA-sequencing, CRISPR genome editing, and Massive Parallel Reporter Assays to analyze non-coding regulatory genetics.

EXPERT OPINION

Despite progress in identifying more genetic risk loci and variant-gene pairs for SLE, understanding its pathogenesis and applying findings clinically remains challenging. The polygenic risk score (PRS) has been used as an application of SLE genetics, but with limited performance in non-EUR populations. In the next few years, advancements in proteomics, post-translational modification estimation, and whole-genome sequencing will enhance disease understanding.

Study of NLRP3 Single Nucleotide Polymorphisms in Juvenile Systemic Lupus Erythematosus (JSLE)

Juvenile systemic lupus erythematosus (JSLE) is a multifaceted multifactorial disorder with an unclear etiopathogenesis. Environmental factors, genetic factors, and dysregulated and defective immune system responses are known to have a role in JSLE etiopathogenesis. NLRP3 inflammasome, as an important contributor to immune-mediated inflammatory responses, is assumed to be involved in JSLE etiopathogenesis. To determine whether the NLRP3 genetic variants are altered in patients with JSLE. Fifty-three patients diagnosed with JSLE and 56 healthy sex-matched controls were studied. NLRP3 (C/G rs10754558, C/T rs3806265, C/T rs4612666, A/C rs35829419) gene polymorphisms were evaluated using a TaqMan single-nucleotide polymorphism assay. C allele at position rs3806265 was detected in higher frequencies in patients than in the control group (37.74% vs 23.21%, P-value = .028). At the genotype level at the same position, CT has a significantly higher frequency in patients than the healthy subjects (75.47% vs 46.43%, P-value = .003). The NLRP3 rs3806265 CT genotype was detected at a higher frequency in patients with JSLE than in the healthy control group.

Differential impact of environmental factors on systemic and localized autoimmunity

The influence of environmental factors on the development of autoimmune disease is being broadly investigated to better understand the multifactorial nature of autoimmune pathogenesis and to identify potential areas of intervention. Areas of particular interest include the influence of lifestyle, nutrition, and vitamin deficiencies on autoimmunity and chronic inflammation. In this review, we discuss how particular lifestyles and dietary patterns may contribute to or modulate autoimmunity. We explored this concept through a spectrum of several autoimmune diseases including Multiple Sclerosis (MS), Systemic Lupus Erythematosus (SLE) and Alopecia Areata (AA) affecting the central nervous system, whole body, and the hair follicles, respectively. A clear commonality between the autoimmune conditions of interest here is low Vitamin D, a well-researched hormone in the context of autoimmunity with pleiotropic immunomodulatory and anti-inflammatory effects. While low levels are often correlated with disease activity and progression in MS and AA, the relationship is less clear in SLE. Despite strong associations with autoimmunity, we lack conclusive evidence which elucidates its role in contributing to pathogenesis or simply as a result of chronic inflammation. In a similar vein, other vitamins impacting the development and course of these diseases are explored in this review, and overall diet and lifestyle. Recent work exploring the effects of dietary interventions on MS showed that a balanced diet was linked to improvement in clinical parameters, comorbid conditions, and overall quality of life for patients. In patients with MS, SLE and AA, certain diets and supplements are linked to lower incidence and improved symptoms. Conversely, obesity during adolescence was linked with higher incidence of MS while in SLE it was associated with organ damage. Autoimmunity is thought to emerge from the complex interplay between environmental factors and genetic background. Although the scope of this review focuses on environmental factors, it is imperative to elaborate the interaction between genetic susceptibility and environment due to the multifactorial origin of these disease. Here, we offer a comprehensive review about the influence of recent environmental and lifestyle factors on these autoimmune diseases and potential translation into therapeutic interventions.

CLEC16A-An Emerging Master Regulator of Autoimmunity and Neurodegeneration

CLEC16A is emerging as an important genetic risk factor for several autoimmune disorders and for Parkinson disease (PD), opening new avenues for translational research and therapeutic development. While the exact role of CLEC16A in health and disease is still being elucidated, the gene plays a critical role in the regulation of autophagy, mitophagy, endocytosis, intracellular trafficking, immune function, and in biological processes such as insulin secretion and others that are important to cellular homeostasis. As shown in both human and animal modeling studies, CLEC16A hypofunction predisposes to both autoinflammatory phenotype and neurodegeneration. While the two are clearly related, further functional studies are needed to fully understand the mechanisms involved for optimized therapeutic interventions. Based on recent data, mitophagy-inducing drugs may be warranted, and such therapy should be tested in clinical trials as these drugs would tackle the underlying pathogenic mechanism (s) and could treat or prevent symptoms of autoimmunity and neurodegeneration in individuals with CLEC16A risk variants. Accordingly, interventions directed at reversing the dysregulated mitophagy and the consequences of loss of function of CLEC16A without activating other detrimental cellular pathways could present an effective therapy. This review presents the emerging role of CLEC16A in health and disease and provides an update on the disease processes that are attributed to variants located in the CLEC16A gene, which are responsible for autoimmune disorders and neurodegeneration with emphasis on how this information is being translated into practical and effective applications in the clinic.

B lymphocytes in treatment-naive paediatric patients with lupus are epigenetically distinct from healthy children

BACKGROUND

SLE is likely triggered by gene-environment interactions. We have shown that most SLE-associated haplotypes encompass genomic regions enriched for epigenetic marks associated with enhancer function in lymphocytes, suggesting genetic risk is exerted through altered gene regulation. Data remain scarce on how epigenetic variance contributes to disease risk in paediatric SLE (pSLE). We aim to identify differences in epigenetically regulated chromatin architecture in treatment-naive patients with pSLE compared with healthy children.

METHODS

Using the assay for transposase-accessible chromatin with sequencing (ATACseq), we surveyed open chromatin in 10 treatment-naive patients with pSLE, with at least moderate disease severity, and 5 healthy children. We investigated whether regions of open chromatin unique to patients with pSLE demonstrate enrichment for specific transcriptional regulators, using standard computational approaches to identify unique peaks and a false discovery rate of <0.05. Further analyses for histone modification enrichment and variant calling were performed using bioinformatics packages in R and Linux.

RESULTS

We identified 30 139 differentially accessible regions (DAR) unique to pSLE B cells; 64.3% are more accessible in pSLE than healthy children. Many DAR are found in distal, intergenic regions and enriched for enhancer histone marks (p=0.027). B cells from adult patients with SLE contain more regions of inaccessible chromatin than those in pSLE. In pSLE B cells, 65.2% of the DAR are located within or near known SLE haplotypes. Further analysis revealed enrichment of transcription factor binding motifs within these DAR that may regulate genes involved in pro-inflammatory responses and cellular adhesion.

CONCLUSIONS

We demonstrate an epigenetically distinct profile in pSLE B cells when compared with healthy children and adults with lupus, indicating that pSLE B cells are predisposed for disease onset/development. Increased chromatin accessibility in non-coding genomic regions controlling activation of inflammation suggest that transcriptional dysregulation by regulatory elements controlling B cell activation plays an important role in pSLE pathogenesis.

Molecular pathways identified from single nucleotide polymorphisms demonstrate mechanistic differences in systemic lupus erythematosus patients of Asian and European ancestry

Systemic lupus erythematosus (SLE) is a multi-organ autoimmune disorder with a prominent genetic component. Individuals of Asian-Ancestry (AsA) disproportionately experience more severe SLE compared to individuals of European-Ancestry (EA), including increased renal involvement and tissue damage. However, the mechanisms underlying elevated severity in the AsA population remain unclear. Here, we utilized available gene expression data and genotype data based on all non-HLA SNP associations in EA and AsA SLE patients detected using the Immunochip genotyping array. We identified 2778 ancestry-specific and 327 trans-ancestry SLE-risk polymorphisms. Genetic associations were examined using connectivity mapping and gene signatures based on predicted biological pathways and were used to interrogate gene expression datasets. SLE-associated pathways in AsA patients included elevated oxidative stress, altered metabolism and mitochondrial dysfunction, whereas SLE-associated pathways in EA patients included a robust interferon response (type I and II) related to enhanced cytosolic nucleic acid sensing and signaling. An independent dataset derived from summary genome-wide association data in an AsA cohort was interrogated and identified similar molecular pathways. Finally, gene expression data from AsA SLE patients corroborated the molecular pathways predicted by SNP associations. Identifying ancestry-related molecular pathways predicted by genetic SLE risk may help to disentangle the population differences in clinical severity that impact AsA and EA individuals with SLE.

Abnormalities in intron retention characterize patients with systemic lupus erythematosus

Regulation of intron retention (IR), a form of alternative splicing, is a newly recognized checkpoint in gene expression. Since there are numerous abnormalities in gene expression in the prototypic autoimmune disease systemic lupus erythematosus (SLE), we sought to determine whether IR was intact in patients with this disease. We, therefore, studied global gene expression and IR patterns of lymphocytes in SLE patients. We analyzed RNA-seq data from peripheral blood T cell samples from 14 patients suffering from systemic lupus erythematosus (SLE) and 4 healthy controls and a second, independent data set of RNA-seq data from B cells from16 SLE patients and 4 healthy controls. We identified intron retention levels from 26,372 well annotated genes as well as differential gene expression and tested for differences between cases and controls using unbiased hierarchical clustering and principal component analysis. We followed with gene-disease enrichment analysis and gene-ontology enrichment analysis. Finally, we then tested for significant differences in intron retention between cases and controls both globally and with respect to specific genes. Overall decreased IR was found in T cells from one cohort and B cells from another cohort of patients with SLE and was associated with increased expression of numerous genes, including those encoding spliceosome components. Different introns within the same gene displayed both up- and down-regulated retention profiles indicating a complex regulatory mechanism. These results indicate that decreased IR in immune cells is characteristic of patients with active SLE and may contribute to the abnormal expression of specific genes in this autoimmune disease.

Inhibitor of NF-κB Kinase Subunit ε Contributes to Neuropsychiatric Manifestations in Lupus-Prone Mice Through Microglial Activation

OBJECTIVE

Systemic lupus erythematosus (SLE) is a systemic autoimmune disease characterized by multiorgan dysfunction. Neuropsychiatric SLE (NPSLE) occurs in 30-40% of lupus patients and is the most severe presentation of SLE, frequently resulting in limitation of daily life. Recent studies have shown that microglia, tissue-resident macrophages in the central nervous system, are involved in the pathogenesis of NPSLE. This study was undertaken to explore new therapeutic targets for NPSLE focusing on microglia.

METHODS

RNA sequencing of microglia in MRL/lpr, lupus-prone mice, as well as that of microglia cultured in vitro with cytokines were performed. A candidate gene, which could be a therapeutic target for NPSLE, was identified, and its role in microglial activation and phagocytosis was investigated using specific inhibitors and small interfering RNA. The effect of intracerebroventricular administration of the inhibitor on the behavioral abnormalities of MRL/lpr was also evaluated.

RESULTS

Transcriptome analysis revealed the up-regulation of Ikbke, which encodes the inhibitor of NF-κB kinase subunit ɛ (IKBKε) in both microglia from MRL/lpr mice and cytokine-stimulated microglia in vitro. Intracerebroventricular administration of an IKBKε inhibitor ameliorated cognitive function and suppressed microglial activation in MRL/lpr mice. Mechanistically, IKBKε inhibition reduced glycolysis, which dampened microglial activation and phagocytosis.

CONCLUSION

These findings suggest that IKBKε plays a vital role in the pathogenesis of NPSLE via microglial activation, and it could serve as a therapeutic target for NPSLE.

Multitrait genome-wide analyses identify new susceptibility loci and candidate drugs to primary sclerosing cholangitis

Primary sclerosing cholangitis (PSC) is a rare autoimmune bile duct disease that is strongly associated with immune-mediated disorders. In this study, we implemented multitrait joint analyses to genome-wide association summary statistics of PSC and numerous clinical and epidemiological traits to estimate the genetic contribution of each trait and genetic correlations between traits and to identify new lead PSC risk-associated loci. We identified seven new loci that have not been previously reported and one new independent lead variant in the previously reported locus. Functional annotation and fine-mapping nominated several potential susceptibility genes such as MANBA and IRF5. Network-based in silico drug efficacy screening provided candidate agents for further study of pharmacological effect in PSC.

Risk prediction models for incident systemic lupus erythematosus among women in the Nurses' health study cohorts using genetics, family history, and lifestyle and environmental factors

OBJECTIVE

Systemic lupus erythematosus (SLE) is a severe multisystem autoimmune disease that predominantly affects women. Its etiology is complex and multifactorial, with several known genetic and environmental risk factors, but accurate risk prediction models are still lacking. We developed SLE risk prediction models, incorporating known genetic, lifestyle and environmental risk factors, and family history.

METHODS

We performed a nested case-control study within the Nurses' Health Study cohorts (NHS). NHS began in 1976 and enrolled 121,700 registered female nurses ages 30-55 from 11 U.S. states; NHSII began in 1989 and enrolled 116,430 registered female nurses ages 25-42 from 14 U.S. states. Participants were asked about lifestyle, reproductive and environmental exposures, as well as medical information, on biennial questionnaires. Incident SLE cases were self-reported and validated by medical record review (Updated 1997 American College of Rheumatology classification criteria). Those with banked blood samples for genotyping (∼25% of each cohort), were selected and matched by age (± 4 years) and race/ethnicity to women who had donated a blood sample but did not develop SLE. Lifestyle and reproductive variables, including smoking, alcohol use, body mass index, sleep, socioeconomic status, U.S. region, menarche age, oral contraceptive use, menopausal status/postmenopausal hormone use, and family history of SLE or rheumatoid arthritis (RA) were assessed through the questionnaire prior to SLE diagnosis questionnaire cycle (or matched index date). Genome-wide genotyping results were used to calculate a SLE weighted genetic risk score (wGRS) using 86 published single nucleotide polymorphisms (SNPs) and 10 classical HLA alleles associated with SLE. We compared four sequential multivariable logistic regression models of SLE risk prediction, each calculating the area under the receiver operating characteristic curve (AUC): 1) SLE wGRS, 2) SLE/RA family history, 3) lifestyle, environmental and reproductive factors and 4) combining model 1-3 factors. Models were internally validated using a bootstrapped estimate of optimism of the AUC. We also examined similar sequential models to predict anti-dsDNA positive SLE risk.

RESULTS

We identified and matched 138 women who developed incident SLE to 1136 women who did not. Models 1-4 yielded AUCs 0.63 (95%CI 0.58-0.68), 0.64 (95%CI 0.59-0.68), 0.71(95% CI 0.66-0.75), and 0.76 (95% CI 0.72-0.81). Model 4 based on genetics, family history and eight lifestyle and environmental factors had best discrimination, with an optimism-corrected AUC 0.75. AUCs for similar models predicting anti-dsDNA positive SLE risk, were 0.60, 0.63, 0.81 and 0.82, with optimism corrected AUC of 0.79 for model 4.

CONCLUSION

A final model including SLE weighted genetic risk score, family history and eight lifestyle and environmental SLE risk factors accurately classified future SLE risk with optimism corrected AUC of 0.75. To our knowledge, this is the first SLE prediction model based on known risk factors. It might be feasibly employed in at-risk populations as genetic data are increasingly available and the risk factors easily assessed. The NHS cohorts include few non-White women and mean age at incident SLE was early 50s, calling for further research in younger and more diverse cohorts.

Two cases of systemic lupus erythematosus (SLE) after administration of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccine

Although there is great demand for increased COVID-19 vaccination worldwide, rare side effects of the vaccines in susceptible individuals are attracting attention. We recently treated two patients who developed SLE after administration of a SARS-CoV-2 vaccine from Pfizer-BioNTech or Moderna. While causal relationships between vaccination and adverse events are difficult to discern due to both confounding and masking factors, our findings suggest that attention to possible adjuvant-related autoimmune diseases in certain individuals receiving SARS-CoV-2 vaccines is appropriate.

Perturbed autophagy intervenes systemic lupus erythematosus by active ingredients of traditional Chinese medicine

Systemic lupus erythematosus (SLE) is a common multisystem, multiorgan heterozygous autoimmune disease. The main pathological features of the disease are autoantibody production and immune complex deposition. Autophagy is an important mechanism to maintain cell homeostasis. Autophagy functional abnormalities lead to the accumulation of apoptosis and induce the autoantibodies that result in immune disorders. Therefore, improving autophagy may alleviate the development of SLE. For SLE, glucocorticoids or immunosuppressive agents are commonly used in clinical treatment, but long-term use of these drugs causes serious side effects in humans. Immunosuppressive agents are expensive. Traditional Chinese medicines (TCMs) are widely used for immune diseases due to their low toxicity and few side effects. Many recent studies found that TCM and its active ingredients affected the pathological development of SLE by regulating autophagy. This article explains how autophagy interferes with immune system homeostasis and participates in the occurrence and development of SLE. It also summarizes several studies on TCM-regulated autophagy intervention in SLE to generate new ideas for basic research, the development of novel medications, and the clinical treatment of SLE.

Contribution of genetic variants associated with primary immunodeficiencies to childhood-onset systemic lupus erythematous

BACKGROUND

A dysregulated immune response is a hallmark of autoimmune disorders. Evidence suggests that systemic autoimmune diseases and primary immunodeficiency disorders (PIDs) may be similar diseases with different clinical phenotypes.

OBJECTIVE

This study aimed to investigate the burden of PID-associated genetic variants in patients with childhood-onset systemic lupus erythematosus (cSLE).

METHODS

We enrolled 118 cSLE patients regularly followed at Chang Gung Memorial Hospital. Targeted next-generation sequencing identified PID genetic variants in patients versus 1475 unrelated healthy individuals, which were further filtered by allelic frequency and various functional scores. Customized immune assays tested the functions of the identified variants.

RESULTS

On filtration, 36 patients (30.5%) harbored rare variants in PID-associated genes predicted to be damaging. One homozygous TREX1 (c.294dupA) mutation and 4 heterozygous variants with possible dominant PID traits, including BCL11B (c.G1040T), NFKB1 (c.T695G), and NFKB2 (c.G1210A, c.G1651A), were discovered. With recessive traits, variants were found across all PID types; one fifth involved phagocyte number or function defects. Predicted pathogenic PID variants were more predominant in those with a family history of lupus, regardless of infection susceptibility. Moreover, mutation loads were greater among cSLE patients than controls despite sex or age at disease onset. While greater mutation loads were observed among cSLE patients with peripubertal disease onset, no significant differences in sex or phenotype were noted among cSLE patients.

CONCLUSION

cSLE is mostly not monogenic. Gene-specific analysis and mutation load investigations suggested that rare and predicted damaging variants in PID-related genes can potentially contribute to cSLE susceptibility.

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.

Proteomics analysis of lysine crotonylation and 2-hydroxyisobutyrylation reveals significant features of systemic lupus erythematosus

INTRODUCTION/OBJECTIVES

To seek significant features of systemic lupus erythematosus (SLE) by utilizing bioinformatics analysis.

METHOD

Liquid chromatography-tandem mass spectrometry (LC-MS/MS) was used to quantify lysine crotonylation (Kcr) and lysine 2-hydroxyisobutyrylation (Khib) in peripheral blood mononuclear cells (PBMCs) of systemic lupus erythematosus (SLE) patients and normal controls.

RESULTS

Seventy-six differentially modified proteins (DMPs) dually modified by Kcr and Khib were identified between SLE patients and healthy people. GO enrichment analysis prompted significant enrichment of seventy-six DMPs in MHC class II protein complex binding and leukocyte migration. KEGG pathways were enriched in antigen processing and presentation pathway and leukocyte transendothelial migration pathway. Six DMPs (CLTC, HSPA1B, HSPA8, HSP90AB1, HSPD1, and PDIA3) were identified in antigen processing and presentation pathway, of which HSPA8 was the core protein. Significant changes of Kcr and Khib in HSPA8 may increase ATP hydrolysis and promote antigen binding to MHC II molecule. In leukocyte transendothelial migration pathway, 7 DMPs (ACTN1, ACTN4, EZR, MSN, RAC1, RHOA, and VCL) were identified. MSN was the protein with the most modification sites in this pathway. In amino terminal ferm region of MSN, Kcr and Khib expression change may lead to the adhesion between leukocytes and endothelial cells, which was an important step of leukocyte migration.

CONCLUSION

Kcr and Khib may promote the antigen presentation and jointly regulate the tissue damage mediated by leukocyte migration in SLE patients, which may play key roles in the pathogenesis of SLE probably. Key Points • Antigen processing and presentation and leukocyte transendothelial migration may play key roles in the pathogenesis of SLE.

Immunogenetic Profiling of SLE and LN among Jordanian Patients

Systemic Lupus Erythematosus (SLE) is a prolonged inflammatory autoimmune disease, which is characterized by a high titer of serological autoantibodies. Interactions between environmental and genetic factors play a crucial role in the pathogenesis of SLE. Human Leukocyte Antigen (HLA) genes, namely HLA-class II genes, are one of the main candidate genes that increase susceptibility to SLE. The aim of this study was to investigate, for the first time, the association of HLA-DRB1 and HLA-DQB1 genes among Jordanian patients diagnosed with SLE and Lupus Nephritis (LN) using the Polymerase Chain Reaction-Sequence-Specific Primer (PCR-SSP) technique. This study showed that SLE is positively associated with DRB1*0301, DRB1*1101, DRB1*1102 and HLA-DQB1*0601. Furthermore, HLA-DRB1*0301, DRB1*1101, HLA-DRB1*1501 and HLA-DQB1*0601 were found to be linked to SLE patients with LN. In addition, haplotypes HLA-DRB1*0301/DQB1*0201 and HLA-DRB1*1501/DQB1*0601 were found to be linked to SLE and LN. Our findings may serve as possible predictive markers for early screening for LN risk in SLE patients. In light of these results, the role of HLA gene polymorphisms may help in understanding the clinical course, prognosis of the disease and developing better treatment strategies for SLE patients. In addition, it may help in early diagnosis, prevention, intervention and management of the disease.

Mendelian randomization and pathway analysis demonstrate shared genetic associations between lupus and coronary artery disease

Coronary artery disease (CAD) is a leading cause of death in patients with systemic lupus erythematosus (SLE). Despite clinical evidence supporting an association between SLE and CAD, pleiotropy-adjusted genetic association studies are limited and focus on only a few common risk loci. Here, we identify a net positive causal estimate of SLE-associated non-HLA SNPs on CAD by traditional Mendelian randomization (MR) approaches. Pathway analysis using SNP-to-gene mapping followed by unsupervised clustering based on protein-protein interactions (PPIs) identifies biological networks composed of positive and negative causal sets of genes. In addition, we confirm the casual effects of specific SNP-to-gene modules on CAD using only SNP mapping to each PPI-defined functional gene set as instrumental variables. This PPI-based MR approach elucidates various molecular pathways with causal implications between SLE and CAD and identifies biological pathways likely causative of both pathologies, revealing known and novel therapeutic interventions for managing CAD in SLE.

COVID-19 and systemic lupus erythematosus genetics: A balance between autoimmune disease risk and protection against infection

Genome wide association studies show there is a genetic component to severe COVID-19. We find evidence that the genome-wide genetic association signal with severe COVID-19 is correlated with that of systemic lupus erythematosus (SLE), having formally tested this using genetic correlation analysis by LD score regression. To identify the shared associated loci and gain insight into the shared genetic effects, using summary level data we performed meta-analyses, a local genetic correlation analysis and fine-mapping using stepwise regression and functional annotation. This identified multiple loci shared between the two traits, some of which exert opposing effects. The locus with most evidence of shared association is TYK2, a gene critical to the type I interferon pathway, where the local genetic correlation is negative. Another shared locus is CLEC1A, where the direction of effects is aligned, that encodes a lectin involved in cell signaling, and the anti-fungal immune response. Our analyses suggest that several loci with reciprocal effects between the two traits have a role in the defense response pathway, adding to the evidence that SLE risk alleles are protective against infection.

We observed a correlation between the genetic associations with severe COVID-19 and those with systemic lupus erythematosus (SLE, Lupus), and aimed to discover which genetic loci were shared by these diseases and what biological processes were involved. This resulted in the discovery of several genetic loci, some of which had alleles that were risk for both diseases and some of which were risk for severe COVID-19 yet protective for SLE. The locus with most evidence of shared association (TYK2) is involved in interferon production, a process that is important in response to viral infection and known to be dysregulated in SLE patients. Other shared associated loci contained genes also involved in the defense response and the immune system signaling. These results add to the growing evidence that there are alleles in the human genome that provide protection against viral infection yet are risk for autoimmune disease.

A database of 5305 healthy Korean individuals reveals genetic and clinical implications for an East Asian population

Despite substantial advances in disease genetics, studies to date have largely focused on individuals of European descent. This limits further discoveries of novel functional genetic variants in other ethnic groups. To alleviate the paucity of East Asian population genome resources, we established the Korean Variant Archive 2 (KOVA 2), which is composed of 1896 whole-genome sequences and 3409 whole-exome sequences from healthy individuals of Korean ethnicity. This is the largest genome database from the ethnic Korean population to date, surpassing the 1909 Korean individuals deposited in gnomAD. The variants in KOVA 2 displayed all the known genetic features of those from previous genome databases, and we compiled data from Korean-specific runs of homozygosity, positively selected intervals, and structural variants. In doing so, we found loci, such as the loci of ADH1A/1B and UHRF1BP1, that are strongly selected in the Korean population relative to other East Asian populations. Our analysis of allele ages revealed a correlation between variant functionality and evolutionary age. The data can be browsed and downloaded from a public website ( https://www.kobic.re.kr/kova/ ). We anticipate that KOVA 2 will serve as a valuable resource for genetic studies involving East Asian populations.

Longitudinal Immune Cell Profiling in Patients With Early Systemic Lupus Erythematosus

OBJECTIVES

To investigate the immune cell profiles of patients with systemic lupus erythematosus (SLE), and to identify longitudinal changes in those profiles over time.

METHODS

We employed mass cytometry with 3 different panels of 38-39 markers (an immunophenotyping panel, a T cell/monocyte panel, and a B cell panel) in cryopreserved peripheral blood mononuclear cells (PBMCs) from 9 patients with early SLE, 15 patients with established SLE, and 14 controls without autoimmune disease. We used machine learning-driven clustering, flow self-organizing maps, and dimensional reduction with t-distributed stochastic neighbor embedding to identify unique cell populations in early SLE and established SLE. We used mass cytometry data of PBMCs from 19 patients with early rheumatoid arthritis (RA) and 23 controls to compare levels of specific cell populations in early RA and SLE. For the 9 patients with early SLE, longitudinal mass cytometry analysis was applied to PBMCs at enrollment, 6 months after enrollment, and 1 year after enrollment. Serum samples were also assayed for 65 cytokines using Luminex multiplex assay, and associations between cell types and cytokines/chemokines were assessed.

RESULTS

Levels of peripheral helper T cells, follicular helper T (Tfh) cells, and several Ki-67+ proliferating subsets (ICOS+Ki-67+ CD8 T cells, Ki-67+ regulatory T cells, CD19intermediate Ki-67high plasmablasts, and PU.1high Ki-67high monocytes) were increased in patients with early SLE, with more prominent alterations than were seen in patients with early RA. Longitudinal mass cytometry and multiplex serum cytokine assays of samples from patients with early SLE revealed that levels of Tfh cells and CXCL10 had decreased 1 year after enrollment. Levels of CXCL13 were positively correlated with levels of several of the expanded cell populations in early SLE.

CONCLUSION

Two major helper T cell subsets and unique Ki-67+ proliferating immune cell subsets were expanded in patients in the early phase of SLE, and the immunologic features characteristic of early SLE evolved over time.

Integrative Functional Genomics Identifies Systemic Lupus Erythematosus Causal Genetic Variant in the IRF5 Risk Locus

OBJECTIVE

IRF5 plays a crucial role in the development of lupus. Genome-wide association studies have identified several systemic lupus erythematosus (SLE) risk single-nucleotide polymorphisms (SNPs) enriched in the IRF5 locus. However, no comprehensive genome editing-based functional analysis exists to establish a direct link between these variants and altered IRF5 expression, particularly for enhancer variants. This study was undertaken to dissect the regulatory function and mechanisms of SLE IRF5 enhancer risk variants and to explore the utilization of clustered regularly interspaced short palindromic repeat interference (CRISPRi) to regulate the expression of disease risk gene to intervene in the disease.

METHODS

Epigenomic profiles and expression quantitative trait locus analysis were applied to prioritize putative functional variants in the IRF5 locus. CRISPR-mediated deletion, activation, and interference were performed to investigate the genetic function of rs4728142. Allele-specific chromatin immunoprecipitation-quantitative polymerase chain reaction and allele-specific formaldehyde-assisted isolation of regulatory element-quantitative polymerase chain reaction were used to decipher the mechanism of alleles differentially regulating IRF5 expression. The CRISPRi approach was used to evaluate the intervention effect in monocytes from SLE patients.

RESULTS

SLE risk SNP rs4728142 was located in an enhancer region, indicating a disease-related regulatory function, and risk allele rs4728142-A was closely associated with increased IRF5 expression. We demonstrated that an rs4728142-containing region could act as an enhancer to regulate the expression of IRF5. Moreover, rs4728142 affected the binding affinity of zinc finger and BTB domain-containing protein 3 (ZBTB3), a transcription factor involved in regulation. Furthermore, in monocytes from SLE patients, CRISPR-based interference with the regulation of this enhancer attenuated the production of disease-associated cytokines.

CONCLUSION

These results demonstrate that the rs4728142-A allele increases the SLE risk by affecting ZBTB3 binding, chromatin status, and regulating IRF5 expression, establishing a biologic link between genetic variation and lupus pathogenesis.

Hypomethylation of miR-17-92 cluster in lupus T cells and no significant role for genetic factors in the lupus-associated DNA methylation signature

OBJECTIVES

Lupus T cells demonstrate aberrant DNA methylation patterns dominated by hypomethylation of interferon-regulated genes. The objective of this study was to identify additional lupus-associated DNA methylation changes and determine the genetic contribution to epigenetic changes characteristic of lupus.

METHODS

Genome-wide DNA methylation was assessed in naïve CD4+ T cells from 74 patients with lupus and 74 age-matched, sex-matched and race-matched healthy controls. We applied a trend deviation analysis approach, comparing methylation data in our cohort with over 16 500 samples. Methylation quantitative trait loci (meQTL) analysis was performed by integrating methylation profiles with genome-wide genotyping data.

RESULTS

In addition to the previously reported epigenetic signature in interferon-regulated genes, we observed hypomethylation in the promoter region of the miR-17-92 cluster in patients with lupus. Members of this microRNA cluster play an important role in regulating T cell proliferation and differentiation. Expression of two microRNAs in this cluster, miR-19b1 and miR-18a, showed a significant positive correlation with lupus disease activity. Among miR-18a target genes, TNFAIP3, which encodes a negative regulator of nuclear factor kappa B, was downregulated in lupus CD4+ T cells. MeQTL identified in lupus patients showed overlap with genetic risk loci for lupus, including CFB and IRF7. The lupus risk allele in IRF7 (rs1131665) was associated with significant IRF7 hypomethylation. However, <1% of differentially methylated CpG sites in patients with lupus were associated with an meQTL, suggesting minimal genetic contribution to lupus-associated epigenotypes.

CONCLUSION

The lupus defining epigenetic signature, characterised by robust hypomethylation of interferon-regulated genes, does not appear to be determined by genetic factors. Hypomethylation of the miR-17-92 cluster that plays an important role in T cell activation is a novel epigenetic locus for lupus.

Multiethnic polygenic risk prediction in diverse populations through transfer learning

Polygenic risk scores (PRS) leverage the genetic contribution of an individual’s genotype to a complex trait by estimating disease risk. Traditional PRS prediction methods are predominantly for the European population. The accuracy of PRS prediction in non-European populations is diminished due to much smaller sample size of genome-wide association studies (GWAS). In this article, we introduced a novel method to construct PRS for non-European populations, abbreviated as TL-Multi, by conducting a transfer learning framework to learn useful knowledge from the European population to correct the bias for non-European populations. We considered non-European GWAS data as the target data and European GWAS data as the informative auxiliary data. TL-Multi borrows useful information from the auxiliary data to improve the learning accuracy of the target data while preserving the efficiency and accuracy. To demonstrate the practical applicability of the proposed method, we applied TL-Multi to predict the risk of systemic lupus erythematosus (SLE) in the Asian population and the risk of asthma in the Indian population by borrowing information from the European population. TL-Multi achieved better prediction accuracy than the competing methods, including Lassosum and meta-analysis in both simulations and real applications.