Identification of 38 novel loci for systemic lupus erythematosus and genetic heterogeneity between ancestral groups

Multiple polygenic risk scores can improve the prediction of systemic lupus erythematosus in Taiwan

OBJECTIVE

To identify new genetic variants associated with SLE in Taiwan and establish polygenic risk score (PRS) models to improve the early diagnostic accuracy of SLE.

METHODS

The study enrolled 2429 patients with SLE and 48 580 controls from China Medical University Hospital in Taiwan. A genome-wide association study (GWAS) and PRS analyses of SLE and other three SLE markers, namely ANA, anti-double-stranded DNA antibody (dsDNA) and anti-Smith antibody (Sm), were conducted.

RESULTS

Genetic variants associated with SLE were identified through GWAS. Some novel genes, which have been previously reported, such as RCC1L and EGLN3, were revealed to be associated with SLE in Taiwan. Multiple PRS models were established, and optimal cut-off points for each PRS were determined using the Youden Index. Combining the PRSs for SLE, ANA, dsDNA and Sm yielded an area under the curve of 0.64 for the optimal cut-off points. An analysis of human leucocyte antigen (HLA) haplotypes in SLE indicated that individuals with HLA-DQA1*01:01 and HLA-DQB1*05:01 were at a higher risk of being classified into the SLE group.

CONCLUSIONS

The use of PRSs to predict SLE enables the identification of high-risk patients before abnormal laboratory data were obtained or symptoms were manifested. Our findings underscore the potential of using PRSs and GWAS in identifying SLE markers, offering promise for early diagnosis and prediction of SLE.

Association between systemic lupus erythematosus and osteoporosis: a mendelian randomization analysis

BACKGROUND

Systemic Lupus Erythematosus (SLE) and Osteoporosis are two prevalent medical conditions. Previous studies have suggested a possible correlation between SLE and osteoporosis, though the underpinning causal relationship remains largely unknown. The current study aimed to elucidate the causal association between SLE and osteoporosis by employing a Mendelian randomization (MR) approach.

METHODS

We performed two-sample MR analysis using the inverse variance-weighted (IVW), weighted median, and MR-Egger methods on publicly available summary statistics datasets using a SLE genome-wide association study (GWAS) as an exposure and osteoporosis GWASs in people with East Asia ancestry as outcomes. The pleiotropy and heterogeneity were examined using a variety of techniques, including the MR-Egger intercept, the MR-PRESSO approach, and the Cochran's Q test.

RESULTS

We selected 26 single-nucleotide polymorphisms from a SLE GWAS as instrumental variables for osteoporosis. The IVW (p < 0.05) method results support a potential association between SLE and osteoporosis. MR-Egger intercept (p = 0.82) and MR-PRESSO global test (p = 0.80) did not suggest evidence of horizontal or directional pleiotropy. Cochran's Q test (p = 0.78) showed that there was no heterogeneity between IVs.

CONCLUSION

The results of MR analysis indicated that SLE is likely associated with an increased risk of osteoporosis incidence. Our findings highlight the need for increased awareness the potential risk of osteoporosis among SLE patients.

Thioredoxin is a metabolic rheostat controlling regulatory B cells

Metabolic programming is important for B cell fate, but the bioenergetic requirement for regulatory B (Breg) cell differentiation and function is unknown. Here we show that Breg cell differentiation, unlike non-Breg cells, relies on mitochondrial electron transport and homeostatic levels of reactive oxygen species (ROS). Single-cell RNA sequencing analysis revealed that TXN, encoding the metabolic redox protein thioredoxin (Trx), is highly expressed by Breg cells, unlike Trx inhibitor TXNIP which was downregulated. Pharmacological inhibition or gene silencing of TXN resulted in mitochondrial membrane depolarization and increased ROS levels, selectively suppressing Breg cell differentiation and function while favoring pro-inflammatory B cell differentiation. Patients with systemic lupus erythematosus (SLE), characterized by Breg cell deficiencies, present with B cell mitochondrial membrane depolarization, elevated ROS and fewer Trx+ B cells. Exogenous Trx stimulation restored Breg cells and mitochondrial membrane polarization in SLE B cells to healthy B cell levels, indicating Trx insufficiency underlies Breg cell impairment in patients with SLE.

Graves' disease and the risk of five autoimmune diseases: A Mendelian randomization and colocalization study

BACKGROUND

Epidemiological studies have consistently demonstrated a high prevalence of concurrent autoimmune diseases in individuals with Graves' disease (GD).

OBJECTIVE

The objective of this study is to establish a causal association between GD and autoimmune diseases.

METHODS

We employed a two-sample Mendelian randomization (MR) to infer a causal association between GD and five autoimmune diseases, namely rheumatoid arthritis (RA), systemic lupus erythematosus (SLE), Crohn's disease (CD), ulcerative colitis (UC), and amyotrophic lateral sclerosis (ALS), in the East Asian and European population. Genetic correlations were explored through linkage disequilibrium score regression analysis (LDSC). Finally, colocalization analyses were performed to investigate possible genetic foundations.

RESULTS

Bidirectional MR analysis indicated that genetically predicted GD increased the risk of RA (Odds Ratio (OR): 1.34, 95 % Confidence Interval (CI): 1.21 to 1.47, P < 0.001) and SLE (OR: 1.21, 95%CI: 1.08 to 1.35, P < 0.001) in the East Asian population. In contrast, we found that genetically predicted RA (OR: 1.14, 95%CI: 1.05 to 1.24, P = 0.002) and SLE (OR: 1.10, 95%CI: 1.03 to 1.17, P = 0.003) were associated with a higher risk of GD. The results have been partially validated in European cohorts. Colocalization analysis suggested the potential existence of shared causal variants between GD and other autoimmune diseases. In particular, gene ARID5B may play an important role in the incidence of autoimmune diseases.

CONCLUSION

This study has confirmed that GD was associated with RA and SLE and found a possible key gene ARID5B. It may be necessary to strengthen detection to prevent the occurrence of comorbidities in clinical practice.

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.

Shared genetic architecture between autoimmune disorders and B-cell acute lymphoblastic leukemia: insights from large-scale genome-wide cross-trait analysis

BACKGROUND

To study the shared genetic structure between autoimmune diseases and B-cell acute lymphoblastic leukemia (B-ALL) and identify the shared risk loci and genes and genetic mechanisms involved.

METHODS

Based on large-scale genome-wide association study (GWAS) summary-level data sets, we observed genetic overlaps between autoimmune diseases and B-ALL, and cross-trait pleiotropic analysis was performed to detect shared pleiotropic loci and genes. A series of functional annotation and tissue-specific analysis were performed to determine the influence of pleiotropic genes. The heritability enrichment analysis was used to detect crucial immune cells and tissues. Finally, bidirectional Mendelian randomization (MR) methods were utilized to investigate the casual associations.

RESULTS

Our research highlighted shared genetic mechanisms between seven autoimmune disorders and B-ALL. A total of 73 pleiotropic loci were identified at the genome-wide significance level (P < 5 × 10-8), 16 of which had strong evidence of colocalization. We demonstrated that several loci have been previously reported (e.g., 17q21) and discovered some novel loci (e.g., 10p12, 5p13). Further gene-level identified 194 unique pleiotropic genes, for example IKZF1, GATA3, IKZF3, GSDMB, and ORMDL3. Pathway analysis determined the key role of cellular response to cytokine stimulus, B cell activation, and JAK-STAT signaling pathways. SNP-level and gene-level tissue enrichment suggested that crucial role pleiotropic mechanisms involved in the spleen, whole blood, and EBV-transformed lymphocytes. Also, hyprcoloc and stratified LD score regression analyses revealed that B cells at different developmental stages may be involved in mechanisms shared between two different diseases. Finally, two-sample MR analysis determined causal effects of asthma and rheumatoid arthritis on B-ALL.

CONCLUSIONS

Our research proved shared genetic architecture between autoimmune disorders and B-ALL and shed light on the potential mechanism that might involve in.

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.

The Relationship Between Systemic Lupus Erythematosus and Osteoporosis Based on Different Ethnic Groups: a Two-Sample Mendelian Randomization Analysis

The previous observational studies could not overcome the effects of confounding variables and reverse causality. We aimed to determine whether there is a causal relationship between systemic lupus erythematosus and osteoporosis in East Asian and European populations, respectively, by two-sample Mendelian Randomization analysis. We obtained and downloaded data from publicly available genome-wide association study databases and analyses for East Asian and European populations, including systemic lupus erythematosus (SLE), osteoporosis (OP), multisite bone mineral density (BMD), and OP with fracture. After screening for instrumental single-nucleotide polymorphisms (SNPs) significantly correlated to SLE, the inverse-variance weighted (IVW) method was used for calculating the ratio and 95% confidence interval, besides utilizing MR-Egger, weighted median, and weighted mode to assess the robustness of the primary outcome. Moreover, multiple analyses, including MR-PRESSO, MR-Egger intercept, Cochran's Q test, as well as "leave-one-out" sensitivity, were used for evaluating horizontal pleiotropy, heterogeneity, and stability. Finally, we exchanged exposure and outcome and performed a reverse MR analysis. IVW (OR = 1.05, 95% CI = 1.01-1.09, P = 0.009) indicated a significant positive correlation between genetically predicted SLE and OP in East Asians. Furthermore, neither heterogeneity nor horizontal pleiotropy was observed. In Europe, there was no significant genetically predicted causal relation between SLE and OP. Bi-directional MR analysis showed no reverse causality between SLE and OP. In the East Asian population, genetically predicted SLE may have had a positive causal relationship with OP. In Europe, there is insufficient evidence for a potential causal relation between SLE and OP or BMD and fracture, and the correlations currently observed may be attributed to a variety of confounder variables.

Genetic control of DNA methylation is largely shared across European and East Asian populations

DNA methylation is an ideal trait to study the extent of the shared genetic control across ancestries, effectively providing hundreds of thousands of model molecular traits with large QTL effect sizes. We investigate cis DNAm QTLs in three European (n = 3701) and two East Asian (n = 2099) cohorts to quantify the similarities and differences in the genetic architecture across populations. We observe 80,394 associated mQTLs (62.2% of DNAm probes with significant mQTL) to be significant in both ancestries, while 28,925 mQTLs (22.4%) are identified in only a single ancestry. mQTL effect sizes are highly conserved across populations, with differences in mQTL discovery likely due to differences in allele frequency of associated variants and differing linkage disequilibrium between causal variants and assayed SNPs. This study highlights the overall similarity of genetic control across ancestries and the value of ancestral diversity in increasing the power to detect associations and enhancing fine mapping resolution.

Association of genetic variation on X chromosome with systemic lupus erythematosus in both Thai and Chinese populations

OBJECTIVES

X chromosome has been considered as a risk factor for SLE, which is a prototype of autoimmune diseases with a significant sex difference (female:male ratio is around 9:1). Our study aimed at exploring the association of genetic variants in X chromosome and investigating the influence of trisomy X in the development of SLE.

METHODS

X chromosome-wide association studies were conducted using data from both Thai (835 patients with SLE and 2995 controls) and Chinese populations (1604 patients with SLE and 3324 controls). Association analyses were performed separately in females and males, followed by a meta-analysis of the sex-specific results. In addition, the dosage of X chromosome in females with SLE were also examined.

RESULTS

Our analyses replicated the association of TMEM187-IRAK1-MECP2, TLR7, PRPS2 and GPR173 loci with SLE. We also identified two loci suggestively associated with SLE. In addition, making use of the difference in linkage disequilibrium between Thai and Chinese populations, a synonymous variant in TMEM187 was prioritised as a likely causal variant. This variant located in an active enhancer of immune-related cells, with the risk allele associated with decreased expression level of TMEM187. More importantly, we identified trisomy X (47,XXX) in 5 of 2231 (0.22%) females with SLE. The frequency is significantly higher than that found in the female controls (0.08%; two-sided exact binomial test P=0.002).

CONCLUSION

Our study confirmed previous SLE associations in X chromosome, and identified two loci suggestively associated with SLE. More importantly, our study indicated a higher risk of SLE for females with trisomy X.

Genetic interrogation for sequence and copy number variants in systemic lupus erythematosus

Early-onset systemic lupus erythematosus presents with a more severe disease and is associated with a greater genetic burden, especially in patients from Black, Asian or Hispanic ancestries. Next-generation sequencing techniques, notably whole exome sequencing, have been extensively used in genomic interrogation studies to identify causal disease variants that are increasingly implicated in the development of autoimmunity. This Review discusses the known casual variants of polygenic and monogenic systemic lupus erythematosus and its implications under certain genetic disparities while suggesting an age-based sequencing strategy to aid in clinical diagnostics and patient management for improved patient care.

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.

Causal effects of endometriosis on SLE, RA and SS risk: evidence from meta-analysis and Mendelian randomization

BACKGROUND

Endometriosis is an underdiagnosed disorder that affects an estimated 6-10% of women of reproductive age. Endometriosis has been reported in epidemiological studies to be associated with autoimmune diseases. However, the relationship remains controversial.

METHODS

A meta-analysis of observational studies was undertaken to evaluate the risk of autoimmune diseases in patients with endometriosis. The relevant studies were retrieved via the databases Medline, Embase and Web of Science until July 20, 2023. Mendelian randomization (MR) was subsequently utilized to scrutinize the causal influence of genetic predisposition toward endometriosis on three autoimmune diseases.

RESULTS

The meta-analysis findings revealed a relationship between endometriosis and the onset of SLE (cohort studies: RR = 1.77, 95% confidence interval (CI): 1.47-2.13, I2 = 0%; Case-control and cross-sectional studies: OR = 5.23, 95% CI: 0.74-36.98, I2 = 98%), RA (cohort studies: RR = 2.18, 95% CI: 1.85-2.55, I2 = 92%; Case-control and cross-sectional studies: OR = 1.40, 95% CI: 1.19-1.64, I2 = 0%) and SS (cohort studies: RR = 1.49, 95% CI: 1.34-1.66, I2 = 0%). Similarly, in our MR study, the results of the inverse-variance-weighted (IVW) model suggested that genetic predisposition to endometriosis was causally associated with an increased risk for SLE (OR = 1.915, 95% CI: 1.204-3.045, p = 0.006) and RA (OR = 1.005, 95% CI: 1.001-1.009, p = 0.014).

CONCLUSIONS

Both our meta-analysis and MR study indicate that endometriosis increases the risk of autoimmune diseases. These findings not only broaden our understanding of the genetic mechanisms underlying the comorbidity of endometriosis and autoimmune diseases, but also offer a new strategy for autoimmune disease prevention.

Cell-type deconvolution of bulk-blood RNA-seq reveals biological insights into neuropsychiatric disorders

Genome-wide association studies (GWASs) have uncovered susceptibility loci associated with psychiatric disorders such as bipolar disorder (BP) and schizophrenia (SCZ). However, most of these loci are in non-coding regions of the genome, and the causal mechanisms of the link between genetic variation and disease risk is unknown. Expression quantitative trait locus (eQTL) analysis of bulk tissue is a common approach used for deciphering underlying mechanisms, although this can obscure cell-type-specific signals and thus mask trait-relevant mechanisms. Although single-cell sequencing can be prohibitively expensive in large cohorts, computationally inferred cell-type proportions and cell-type gene expression estimates have the potential to overcome these problems and advance mechanistic studies. Using bulk RNA-seq from 1,730 samples derived from whole blood in a cohort ascertained from individuals with BP and SCZ, this study estimated cell-type proportions and their relation with disease status and medication. For each cell type, we found between 2,875 and 4,629 eGenes (genes with an associated eQTL), including 1,211 that are not found on the basis of bulk expression alone. We performed a colocalization test between cell-type eQTLs and various traits and identified hundreds of associations that occur between cell-type eQTLs and GWASs but that are not detected in bulk eQTLs. Finally, we investigated the effects of lithium use on the regulation of cell-type expression loci and found examples of genes that are differentially regulated according to lithium use. Our study suggests that applying computational methods to large bulk RNA-seq datasets of non-brain tissue can identify disease-relevant, cell-type-specific biology of psychiatric disorders and psychiatric medication.

GWAS for systemic sclerosis identifies six novel susceptibility loci including one in the Fcγ receptor region

Here we report the largest Asian genome-wide association study (GWAS) for systemic sclerosis performed to date, based on data from Japanese subjects and comprising of 1428 cases and 112,599 controls. The lead SNP is in the FCGR/FCRL region, which shows a penetrating association in the Asian population, while a complete linkage disequilibrium SNP, rs10917688, is found in a cis-regulatory element for IRF8. IRF8 is also a significant locus in European GWAS for systemic sclerosis, but rs10917688 only shows an association in the presence of the risk allele of IRF8 in the Japanese population. Further analysis shows that rs10917688 is marked with H3K4me1 in primary B cells. A meta-analysis with a European GWAS detects 30 additional significant loci. Polygenic risk scores constructed with the effect sizes of the meta-analysis suggest the potential portability of genetic associations beyond populations. Prioritizing the top 5% of SNPs of IRF8 binding sites in B cells improves the fitting of the polygenic risk scores, underscoring the roles of B cells and IRF8 in the development of systemic sclerosis. The results also suggest that systemic sclerosis shares a common genetic architecture across populations.

Graves' disease and systemic lupus erythematosus: a Mendelian randomization study

Introduction

Previous observational studies have established a correlation between Graves' disease(GD) and systemic lupus erythematosus(SLE). However, whether a causal relationship exists between these two diseases remains unknown.We utilized Mendelian randomization to infer the causal association between GD and SLE.

Methods

This study employed GWAS summary statistics of GD and SLE in individuals of Asian descent. The random effect inverse variance weighted (IVW) method was utilized to aggregate the causal effect estimates of all SNPs. Cochran's Q values were computed to evaluate the heterogeneity among instrumental variables. Sensitivity analyses such as MR-Egger method, median weighting method, leave-one-out method, and MR-PRESSO method were used to test whether there was horizontal pleiotropy of instrumental variables.

Results

Our study found genetically predicted GD may increase risk of SLE (OR=1.17, 95% CI 0.99-1.40, p=0.069). Additionally, genetically predicted SLE elevated the risk of developing GD by 15% (OR=1.15, 95% CI 1.05-1.27, p= 0.004). After correcting for possible horizontal pleiotropy by excluding outlier SNPs, the results suggested that GD increased the risk of SLE (OR=1.27, 95% CI 1.09-1.48, p =0.018), while SLE also increased the risk of developing GD (OR=1.13, 95% CI 1.05-1.22, p =0.003).

Conclusion

The findings of the study indicate that there may be a correlation between GD and SLE, with each potentially increasing the risk of the other. These results have important implications for the screening and treatment of patients with co-morbidities in clinical settings, as well as for further research into the molecular mechanisms underlying the relationship between GD and SLE.

Autoimmune diseases and risk of non-Hodgkin lymphoma: A Mendelian randomisation study

Objective

To examine whether genetically predicted susceptibility to ten autoimmune diseases (Behçet's disease, coeliac disease, dermatitis herpetiformis, lupus, psoriasis, rheumatoid arthritis, sarcoidosis, Sjögren's syndrome, systemic sclerosis, and type 1 diabetes) is associated with risk of non-Hodgkin lymphoma (NHL).

Design

Two sample Mendelian randomization (MR) study.

Setting

Genome wide association studies (GWASs) of ten autoimmune diseases, NHL, and four NHL subtypes (i.e., follicular lymphoma, mature T/natural killer-cell lymphomas, non-follicular lymphoma, and other and unspecified types of NHL).

Analysis

We used data from the largest publicly available GWASs of European ancestry for each autoimmune disease, NHL, and NHL subtypes. For each autoimmune disease, we extracted single nucleotide polymorphisms (SNPs) strongly associated (P < 5×10-8) with that disease and that were independent of one another (R2 < 1×10-3) as genetic instruments. SNPs within the human leukocyte antigen region were not considered due to potential pleiotropy. Our primary MR analysis was the inverse-variance weighted analysis. Additionally, we conducted MR-Egger, weighted mode, and weighted median regression to address potential bias due to pleiotropy, and robust adjusted profile scores to address weak instrument bias. We carried out sensitivity analysis limited to the non-immune pathway for nominally significant findings. To account for multiple testing, we set the thresholds for statistical significance at P < 5×10-3.

Participants

The number of cases and controls identified in the relevant GWASs were 437 and 3,325 for Behçet's disease, 4,918 and 5,684 for coeliac disease, 435 and 341,188 for dermatitis herpetiformis, 4,576 and 8,039 for lupus, 11,988 and 275,335 for psoriasis, 22,350 and 74,823 for rheumatoid arthritis, 3,597 and 337,121 for sarcoidosis, 2,735 and 332,115 for Sjögren's syndrome, 9,095 and 17,584 for systemic sclerosis, 18,942 and 501,638 for type 1 diabetes, 2,400 and 410,350 for NHL; and 296 to 2,340 cases and 271,463 controls for NHL subtypes.

Exposures

Genetic variants predicting ten autoimmune diseases: Behçet's disease, coeliac disease, dermatitis herpetiformis, lupus, psoriasis, rheumatoid arthritis, sarcoidosis, Sjögren's syndrome, systemic sclerosis, and type 1 diabetes.

Main outcome measures

Estimated associations between genetically predicted susceptibility to ten autoimmune diseases and the risk of NHL.

Results

The variance of each autoimmune disease explained by the SNPs ranged from 0.3% to 3.1%. Negative associations between type 1 diabetes and sarcoidosis and the risk of NHL were observed (odds ratio [OR] 0.95, 95% confidence interval [CI]: 0.92 to 0.98, P = 5×10-3, and OR 0.92, 95% CI: 0.85 to 0.99, P = 2.8×10-2, respectively). These findings were supported by the sensitivity analyses accounting for potential pleiotropy and weak instrument bias. No significant associations were found between the other eight autoimmune diseases and NHL risk. Of the NHL subtypes, type 1 diabetes was most strongly associated with follicular lymphoma (OR 0.91, 95% CI: 0.86 to 0.96, P = 1×10-3), while sarcoidosis was most strongly associated with other and unspecified NHL (OR 0.86, 95% CI: 0.75 to 0.97, P = 1.8×10-2).

Conclusions

These findings suggest that genetically predicted susceptibility to type 1 diabetes, and to some extent sarcoidosis, might reduce the risk of NHL. However, future studies with different datasets, approaches, and populations are warranted to further examine the potential associations between these autoimmune diseases and the risk of NHL.

Causal relationship between systemic lupus erythematosus and primary liver cirrhosis based on two-sample bidirectional Mendelian randomization and transcriptome overlap analysis

BACKGROUND

Overlapping cases of systemic lupus erythematosus (SLE) and primary biliary cirrhosis (PBC) are rare and have not yet been fully proven to be accidental or have a common genetic basis.

METHODS

Two-sample bidirectional Mendelian randomization (MR) analysis was applied to explore the potential causal relationship between SLE and PBC. The heterogeneity and reliability of MR analysis were evaluated through Cochran's Q-test and sensitivity test, respectively. Next, transcriptome overlap analysis of SLE and PBC was performed using the Gene Expression Omnibus database to identify the potential mechanism of hub genes. Finally, based on MR analysis, the potential causal relationship between hub genes and SLE or PBC was validated again.

RESULTS

The MR analysis results indicated that SLE and PBC were both high-risk factors for the occurrence and development of the other party. On the one hand, MR analysis had heterogeneity, and on the other hand, it also had robustness. Nine hub genes were identified through transcriptome overlap analysis, and machine learning algorithms were used to verify their high recognition efficiency for SLE patients. Finally, based on MR analysis, it was verified that there was no potential causal relationship between the central gene SOCS3 and SLE, but it was a high-risk factor for the potential risk of PBC.

CONCLUSION

The two-sample bidirectional MR analysis revealed that SLE and PBC were high-risk factors for each other, indicating that they had similar genetic bases, which could to some extent overcome the limitation of insufficient overlap in case samples of SLE and PBC. The analysis of transcriptome overlapping hub genes provided a theoretical basis for the potential mechanisms and therapeutic targets of SLE with PBC overlapping cases.

The functional TNFAIP3 rs2230926T/G (Phe127Cys) variant confers risk to systemic lupus erythematosus in a Latin American population

TNFAIP3 is a classical systemic lupus erythematosus (SLE)-associated risk locus identified by genome-wide association studies (GWASs) and replicated by candidate gene association studies primarily in Caucasians and Asians. However, in Latin American populations, its role on SLE susceptibility is not known. We conducted a case-control study to evaluate whether the TNFAIP3 rs2230926T/G (Phe127Cys) variant is associated with risk of developing SLE in a cohort of Mexican patients. The TNFAIP3 rs2230926T/G variant was analyzed in 561 patients with SLE and 499 control subjects, using TaqMan probes. We found that the G allele was associated with susceptibility to SLE under the allelic (OR 2.09, p = 0.005) and genotypic (OR 2.14, p = 0.004) models. In conclusion, our results show that TNFAIP3 rs2230926T/G is a risk factor for the development of SLE in the Mexican population.

rs10924104 in the expression enhancer motif of CD58 confers susceptibility to human autoimmune diseases

CD58 plays roles in cell adhesion and co-stimulation with antigen presentation from major histocompatibility complex class II on antigen-presenting cells to T-cell antigen receptors on naïve T cells. CD58 reportedly contributes to the development of various human autoimmune diseases. Recently, genome-wide association studies (GWASs) identified CD58 as a susceptibility locus for autoimmune diseases such as systemic lupus erythematosus (SLE), multiple sclerosis (MS), and primary biliary cholangitis (PBC). However, the primary functional variant and molecular mechanisms of susceptibility to autoimmune diseases in the CD58 locus were not clarified. Here, rs10924104, located in the ZNF35-binding motif within the gene expression regulatory motif, was identified as the primary functional variant for SLE, MS, and PBC among genetic variants showing stronger linkage disequilibrium (LD) with GWAS-lead variants in the CD58 locus. Expression-quantitative trait locus (e-QTL) data for each distinct blood cell type and in vitro functional analysis using the CRISPR/Cas9 system corroborated the functional role of rs10924104 in the upregulation of CD58 transcription by the disease-risk allele. Additionally, the strength of disease susceptibility observed in the CD58 locus could be accounted for by the strength of LD between rs10924104 and each GWAS-lead variant. In conclusion, the present study demonstrated for the first time the existence of a shared autoimmune disease-related primary functional variant (i.e., rs10924104) that regulates the expression of CD58. Clarifying the molecular mechanism of disease susceptibility derived from such a shared genetic background is important for understanding human autoimmune diseases and human immunology.

Transcriptional profiling upon T cell stimulation reveals down-regulation of inflammatory pathways in T and B cells in SLE versus Sjögren's syndrome

Systemic lupus erythematosus (SLE) and primary Sjögren's syndrome (pSS) share clinical as well as pathogenic similarities. Although previous studies suggest various abnormalities in different immune cell compartments, dedicated cell-type specific transcriptomic signatures are often masked by patient heterogeneity. Here, we performed transcriptional profiling of isolated CD4, CD8, CD16 and CD19 lymphocytes from pSS and SLE patients upon T cell stimulation, in addition to a steady-state condition directly after blood drawing, in total comprising 581 sequencing samples. T cell stimulation, which induced a pronounced inflammatory response in all four cell types, gave rise to substantial re-modulation of lymphocyte subsets in the two autoimmune diseases compared to healthy controls, far exceeding the transcriptomic differences detected at steady-state. In particular, we detected cell-type and disease-specific down-regulation of a range of pro-inflammatory cytokine and chemokine pathways. Such differences between SLE and pSS patients are instrumental for selective immune targeting by future therapies.

Reframing health disparities in SLE: A critical reassessment of racial and ethnic differences in lupus disease outcomes

Health disparities in the prevalence and outcomes of systemic lupus erythematosus (SLE) are well documented across racial and ethnic groups. Similar to other chronic diseases, differences in disease severity among individuals with SLE are likely influenced by both genetic predisposition and multiple social determinants of health. However, research in SLE that jointly examines the genetic and environmental contributions to the disease course is limited, resulting in an incomplete understanding of the biologic and social mechanisms that underly health disparities. While research on health disparities can reveal inequalities and inform resource allocation to improve outcomes, research that relies on racial and ethnic categories to describe diverse groups of people can pose challenges. Additionally, results from research comparing outcomes across socially constructed groups without considering other contributing factors can be misleading. We herein comprehensively examine existing literature on health disparities in SLE, including both clinical studies that examine the relationship between self-reported race and ethnicity and disease outcomes and studies that explore the relationships between genomics and lupus outcomes. Having surveyed this body of research, we propose a framework for research examining health disparities in SLE, including ways to mitigate bias in future studies.

Dynamic regulatory elements in single-cell multimodal data implicate key immune cell states enriched for autoimmune disease heritability

In autoimmune diseases such as rheumatoid arthritis, the immune system attacks the body's own cells. Developing a precise understanding of the cell states where noncoding autoimmune risk variants impart causal mechanisms is critical to developing curative therapies. Here, to identify noncoding regions with accessible chromatin that associate with cell-state-defining gene expression patterns, we leveraged multimodal single-nucleus RNA and assay for transposase-accessible chromatin (ATAC) sequencing data across 28,674 cells from the inflamed synovial tissue of 12 donors. Specifically, we used a multivariate Poisson model to predict peak accessibility from single-nucleus RNA sequencing principal components. For 14 autoimmune diseases, we discovered that cell-state-dependent ('dynamic') chromatin accessibility peaks in immune cell types were enriched for heritability, compared with cell-state-invariant ('cs-invariant') peaks. These dynamic peaks marked regulatory elements associated with T peripheral helper, regulatory T, dendritic and STAT1+CXCL10+ myeloid cell states. We argue that dynamic regulatory elements can help identify precise cell states enriched for disease-critical genetic variation.

Translational immunology: Applying fundamental discoveries to human health and autoimmune diseases

Studying the human immune system is challenging. These challenges stem from the complexity of the immune system itself, the heterogeneity of the immune system between individuals, and the many factors that lead to this heterogeneity including the influence of genetics, environment, and immune experience. Studies of the human immune system in the context of disease are increased in complexity as multiple combinations and variations in immune pathways can lead to a single disease. Thus, although individuals with a disease may share clinical features, the underlying disease mechanisms and resulting pathophysiology can be diverse among individuals with the same disease diagnosis. This has consequences for the treatment of diseases, as no single therapy will work for everyone, therapeutic efficacy varies among patients, and targeting a single immune pathway is rarely 100% effective. This review discusses how to address these challenges by identifying and managing the sources of variation, improving access to high-quality, well-curated biological samples by building cohorts, applying new technologies such as single-cell omics and imaging technologies to interrogate samples, and bringing to bear computational expertise in conjunction with immunologists and clinicians to interpret those results. The review has a focus on autoimmune diseases, including rheumatoid arthritis, MS, systemic lupus erythematosus, and type 1 diabetes, but its recommendations are also applicable to studies of other immune-mediated diseases.

The Application of Genetic Risk Scores in Rheumatic Diseases: A Perspective

Modest effect sizes have limited the clinical applicability of genetic associations with rheumatic diseases. Genetic risk scores (GRSs) have emerged as a promising solution to translate genetics into useful tools. In this review, we provide an overview of the recent literature on GRSs in rheumatic diseases. We describe six categories for which GRSs are used: (a) disease (outcome) prediction, (b) genetic commonalities between diseases, (c) disease differentiation, (d) interplay between genetics and environmental factors, (e) heritability and transferability, and (f) detecting causal relationships between traits. In our review of the literature, we identified current lacunas and opportunities for future work. First, the shortage of non-European genetic data restricts the application of many GRSs to European populations. Next, many GRSs are tested in settings enriched for cases that limit the transferability to real life. If intended for clinical application, GRSs are ideally tested in the relevant setting. Finally, there is much to elucidate regarding the co-occurrence of clinical traits to identify shared causal paths and elucidate relationships between the diseases. GRSs are useful instruments for this. Overall, the ever-continuing research on GRSs gives a hopeful outlook into the future of GRSs and indicates significant progress in their potential applications.

CellGO: a novel deep learning-based framework and webserver for cell-type-specific gene function interpretation

Interpreting the function of genes and gene sets identified from omics experiments remains a challenge, as current pathway analysis tools often fail to consider the critical biological context, such as tissue or cell-type specificity. To address this limitation, we introduced CellGO. CellGO tackles this challenge by leveraging the visible neural network (VNN) and single-cell gene expressions to mimic cell-type-specific signaling propagation along the Gene Ontology tree within a cell. This design enables a novel scoring system to calculate the cell-type-specific gene-pathway paired active scores, based on which, CellGO is able to identify cell-type-specific active pathways associated with single genes. In addition, by aggregating the activities of single genes, CellGO extends its capability to identify cell-type-specific active pathways for a given gene set. To enhance biological interpretation, CellGO offers additional features, including the identification of significantly active cell types and driver genes and community analysis of pathways. To validate its performance, CellGO was assessed using a gene set comprising mixed cell-type markers, confirming its ability to discern active pathways across distinct cell types. Subsequent benchmarking analyses demonstrated CellGO's superiority in effectively identifying cell types and their corresponding cell-type-specific pathways affected by gene knockouts, using either single genes or sets of genes differentially expressed between knockout and control samples. Moreover, CellGO demonstrated its ability to infer cell-type-specific pathogenesis for disease risk genes. Accessible as a Python package, CellGO also provides a user-friendly web interface, making it a versatile and accessible tool for researchers in the field.

Trans-ethnic Mendelian randomization study of systemic lupus erythematosus and common female hormone-dependent malignancies

BACKGROUND

Observational research has reported that systemic lupus erythematosus (SLE) is related to common female hormone-dependent cancers, but the underlying causal effect remains undefined. This study aimed to explore the causal association of these conditions by Mendelian randomization (MR) analysis.

METHODS

We selected instrumental variables for SLE from genome-wide association studies (GWASs) conducted in European and East Asian populations. The genetic variants for female malignant neoplasms were obtained from corresponding ancestry GWASs. We utilized inverse variance weighted (IVW) as the primary analysis, followed by sensitivity analysis. Furthermore, we conducted multivariable MR (MVMR) to estimate direct effects by adjusting for the body mass index and estradiol. Finally, we implemented reverse direction MR analysis and gave a negative example to test the reliability of MR results.

RESULTS

We found SLE was significantly negatively associated with overall endometrial cancer risk (odds ratio [OR] = 0.961, 95% confidence interval [CI] = 0.935-0.987, P  = 3.57E-03) and moderately inversely related to endometrioid endometrial cancer (ENEC) (OR = 0.965, 95% CI = 0.936-0.995, P  = 0.024) risk in the European population by IVW. We replicated these results using other MR models and detected a direct effect by MVMR (overall endometrial cancer, OR = 0.962, 95% CI = 0.941-0.983, P  = 5.11E-04; ENEC, OR = 0.964, 95% CI = 0.940-0.989, P  = 0.005). Moreover, we revealed that SLE was correlated with decreased breast cancer risk (OR = 0.951, 95% CI = 0.918-0.986, P  = 0.006) in the East Asian population by IVW, and the effect was still significant in MVMR (OR = 0.934, 95% CI = 0.859-0.976, P  = 0.002). The statistical powers of positive MR results were all >0.9.

CONCLUSION

This finding suggests a possible causal effect of SLE on the risk of overall endometrial cancer and breast cancer in European and East Asian populations, respectively, by MR analysis, which compensates for inherent limitations of observational research.

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.

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.

Association of lipid lowering drugs and the risk of systemic lupus erythematosus: a drug target Mendelian randomization

Background and objective: An interaction between low-density lipoprotein level, lipid-lowering drugs, and systemic lupus erythematosus (SLE) was reported by previous studies. However, whether lipid-lowering drugs provided protective effect for reducing the risk of SLE was unclear. We aimed to clarify this causal relationship through a drug-target Mendelian randomization (MR) study. Methods: Genetic instruments-single nucleotide polymorphism (SNPs)-were utilized to proxy inhibition of the three genes-3-hydroxy-3-methylglutaryl-CoA reductase (HMGCR), proprotein convertase subtilisin/kexin type 9 (PCSK9), and Niemann-Pick C1-Like 1(NPC1L1), which was corresponded to three lipid-lowering drugs-statins, evolocumab, and ezetimibe. Low-density lipoprotein (LDL) cholesterol was selected as the biomarker for the measurement of the inhibitors of HMGCR, PCSK9, and NPC1L1, and the genetic data were acquired from the Global Lipids Genetics Consortium, which consisted of 1.3 million participants of European ancestry and 146.5 thousand participants of East Asian ancestry. The genetic dataset of SLE was acquired from two large-scale GWAS studies; one recruited 23,210 participants (7,219 SLE cases and 15,991 controls) of European ancestry and the other one recruited 12,653 participants (4,222 SLE cases and 8,431 controls) of Chinese ancestry. The primary analysis used the inverse variance weighted (IVW) method. Four additional sensitivity analyses, colocalization analysis, and stratification analysis were performed. Results: The primary analysis showed that inhibition of PCSK9 (evolocumab) was associated with a significantly lower risk of SLE [odds ratio (OR) 0.51, 95%CI 0.34 to 0.76, p = 0.001] in the European population. The secondary analyses had similar findings. Stratification analysis demonstrated that the preventive effect of PCSK9 inhibition for SLE was similar in both males and females. However, the results were not replicated in the East Asian population. The inhibition of HMGCR (statins) and NPC1L1 (ezetimibe) did not cause a lower risk of SLE. Conclusion: Evolocumab might provide a protective effect on the risk of SLE in the European population, but statins and ezetimibe might not have the protective effect. Further research is necessary to elucidate the specific mechanisms and potential therapeutic applications of PCSK9 inhibitors (evolocumab) in the context of SLE protection.

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.

Benchmarking of local genetic correlation estimation methods using summary statistics from genome-wide association studies

Local genetic correlation evaluates the correlation of additive genetic effects between different traits across the same genetic variants at a genomic locus. It has been proven informative for understanding the genetic similarities of complex traits beyond that captured by global genetic correlation calculated across the whole genome. Several summary-statistics-based approaches have been developed for estimating local genetic correlation, including $\rho$-hess, SUPERGNOVA and LAVA. However, there has not been a comprehensive evaluation of these methods to offer practical guidelines on the choices of these methods. In this study, we conduct benchmark comparisons of the performance of these three methods through extensive simulation and real data analyses. We focus on two technical difficulties in estimating local genetic correlation: sample overlaps across traits and local linkage disequilibrium (LD) estimates when only the external reference panels are available. Our simulations suggest the likelihood of incorrectly identifying correlated regions and local correlation estimation accuracy are highly dependent on the estimation of the local LD matrix. These observations are corroborated by real data analyses of 31 complex traits. Overall, our findings illuminate the distinct results yielded by different methods applied in post-genome-wide association studies (post-GWAS) local correlation studies. We underscore the sensitivity of local genetic correlation estimates and inferences to the precision of local LD estimation. These observations accentuate the vital need for ongoing refinement in methodologies.

Pervasive Sharing of Causal Genetic Risk Factors Contributes to Clinical and Molecular Overlap between Sjögren's Disease and Systemic Lupus Erythematosus

SjD (Sjögren's Disease) and SLE (Systemic Lupus Erythematosus) are similar diseases. There is extensive overlap between the two in terms of both clinical features and pathobiologic mechanisms. Shared genetic risk is a potential explanation of this overlap. In this study, we evaluated whether these diseases share causal genetic risk factors. We compared the causal genetic risk for SLE and SjD using three complementary approaches. First, we examined the published GWAS results for these two diseases by analyzing the predicted causal gene protein-protein interaction networks of both diseases. Since this method does not account for overlapping risk intervals, we examined whether such intervals also overlap. Third, we used two-sample Mendelian randomization (two sample MR) using GWAS summary statistics to determine whether risk variants for SLE are causal for SjD and vice versa. We found that both the putative causal genes and the genomic risk intervals for SLE and SjD overlap 28- and 130-times more than expected by chance (p < 1.1 × 10-24 and p < 1.1 × 10-41, respectively). Further, two sample MR analysis confirmed that alone or in aggregate, SLE is likely causal for SjD and vice versa. [SjD variants predicting SLE: OR = 2.56; 95% CI (1.98-3.30); p < 1.4 × 10-13, inverse-variance weighted; SLE variants predicting SjD: OR = 1.36; 95% CI (1.26-1.47); p < 1.6 × 10-11, inverse-variance weighted]. Notably, some variants have disparate impact in terms of effect size across disease states. Overlapping causal genetic risk factors were found for both diseases using complementary approaches. These observations support the hypothesis that shared genetic factors drive the clinical and pathobiologic overlap between these diseases. Our study has implications for both differential diagnosis and future genetic studies of these two conditions.

Associations Between Osteopontin Variants and Systemic Lupus Erythematosus: A Meta-Analysis

Objective: Osteopontin (OPN) increases T-cell proliferation, interferon production, and CD40 ligand expression, which leads to B-cell proliferation and antibody production. This study was designed to determine whether OPN variants are associated with susceptibility to systemic lupus erythematosus [SLE]. Methods: We searched the Medline, Embase, and KoreaMed databases for available articles. We performed a meta-analysis on the association of OPN 707 T/C (rs1126616) at exon 6, 1083 G/A (rs112772) at the 3'-untranslated region (3'-UTR), 1239 C/A (rs9138) at 3'-UTR, and 9250 T/C (rs11229919) variants in exon 7 with susceptibility to SLE. Results: Ten studies from 9 articles with 2175 SLE patients and 3233 controls were included. The meta-analysis showed a significant association between SLE and the 707 T allele of the OPN 707 T/C variant (odds ratio [OR] = 1.522, 95% confidence interval [CI] = 1.101-2.105, p = 0.044). Stratification by ethnicity indicated an association between the OPN 707 T/C variant and SLE in European and Arab populations. The meta-analysis also revealed a significant association between the OPN 9250 C allele and SLE in the Asian and Arab populations. A significant association was also identified between the +1239 C allele of the OPN 1239 C/A variant and SLE (OR = 1.192, 95% CI = 1.008-1.410, p = 0.040). The meta-analysis indicated no allelic association between SLE and OPN 1083 G/A and the OPN 1239 C/A variants. Conclusions: The OPN 707 T/C variant is associated with SLE susceptibility in European and Arab populations and the OPN 9250 T/C variant is associated with SLE susceptibility in Asian and Arab populations. In addition, associations were found between the OPN 1239 C/A variant and SLE.

Telomere length and its association with systemic lupus erythematosus in an Asian population: A Mendelian randomization study

OBJECTIVES

To investigate whether shorter telomere length is a causal risk factor for systemic lupus erythematosus (SLE) in the Asian population.

METHODS

We applied the two-sample Mendelian randomization (MR) method to the pooled statistics from a genome-wide association study (GWAS) of 6,707 SLE cases and 16,047 controls. We selected nine single-nucleotide polymorphisms (SNPs) with genome-wide significance as instrumental variables for telomere length. The main analysis was carried out by the random-effects inverse-variance weighted (IVW) method. Horizontal pleiotropy was evaluated by the intercept of MR-Egger regression.

RESULTS

A potentially causal relationship between longer genetically predicted telomere length and increased risk of systemic lupus erythematosus (OR = 1.72, 95%CI: 1.21, 2.46, p = 0.01) was observed. The MR-Egger regression demonstrated an intercept proximal to zero (intercept = 0.017, p = 0.69), which does not provide evidence of the presence of horizontal pleiotropy.

CONCLUSIONS

Our findings provided evidence supporting a potential causal relationship between longer telomere length and increased risk of systemic lupus erythematosus.

Functional evaluation of rare OASL variants by analysis of SLE patient-derived iPSCs

BACKGROUND

Systemic lupus erythematosus (SLE) is a chronic systemic autoimmune disease characterized by genetic heterogeneity and an interferon (IFN) signature. The overall landscapes of the heritability of SLE remains unclear.

OBJECTIVES

To identify and elucidate the biological functions of rare variants underlying SLE, we conducted analyses of patient-derived induced pluripotent stem cells (iPSCs) in combination with genetic analysis.

METHODS

Two familial SLE patient- and two healthy donor (HD)-derived iPSCs were established. Type 1 IFN-secreting dendritic cells (DCs) were differentiated from iPSCs. Genetic analyses of SLE-iPSCs, and 117 SLE patients and 107 HDs in the ImmuNexUT database were performed independently. Genome editing of the variants on iPSCs was performed with the CRISPR/Cas9 system.

RESULTS

Type 1 IFN secretion was significantly increased in DCs differentiated from SLE-iPSCs compared to HD-iPSCs. Genetic analyses revealed a rare variant in the 2'-5'-Oligoadenylate Synthetase Like (OASL) shared between SLE-iPSCs and another independent SLE patient, and significant accumulation of OASL variants among SLE patients (HD 0.93%, SLE 6.84%, OR 8.387) in the database. Genome editing of mutated OASL 202Q to wild-type 202 R or wild-type OASL 202 R to mutated 202Q resulted in reduced or enhanced Type 1 IFN secretion of DCs. Three other OASL variants (R60W, T261S and A447V) accumulated in SLE patients had also capacities to enhance Type 1 IFN secretion in response to dsRNA.

CONCLUSIONS

We established a patient-derived iPSC-based strategy to investigate the linkage of genotype and phenotype in autoimmune diseases. Detailed case-based investigations using patient-derived iPSCs provide information to unveil the heritability of the pathogenesis of autoimmune diseases.

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.

Translational implications of newly characterized pathogenic pathways in systemic lupus erythematosus

Improved characterization of relevant pathogenic pathways in systemic lupus erythematosus (SLE) has been further delineated over the last decades. This led to the development of targeted treatments including belimumab and anifrolumab, which recently became available in clinics. Therapeutic targets in SLE encompass interferon (IFN) signaling, B-T costimulation including immune checkpoints, and increasing modalities of B lineage targeting, such as chimeric antigen receptor (CAR) T cells directed against CD19 or sequential anti-B cell targeting. Patient profiling based on characterization of underlying molecular abnormalities, often performed through comprehensive omics analyses, has recently been shown to better predict patients' treatment responses and also holds promise to unravel key molecular mechanisms driving SLE. SLE carries two key signatures, namely the IFN and B lineage/plasma cell signatures. Recent advances in SLE treatments clearly indicate that targeting innate and adaptive immunity is successful in such a complex autoimmune disease. Although those signatures may interact at the molecular level and provide the basis for the first selective treatments in SLE, it remains to be clarified whether these distinct treatments show different treatment responses among certain patient subsets. In fact, notwithstanding the remarkable amount of novel clues for innovative SLE treatment, harmonization of big data within tailored treatment strategies will be instrumental to better understand and treat this challenging autoimmune disorder. This review will provide an overview of recent improvements in SLE pathogenesis, related insights by analyses of big data and machine learning as well as technical improvements in conducting clinical trials with the ultimate goal that translational research results in improved patient outcomes.

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.

Interferon-α regulates abnormally increased expression of RSAD2 in Th17 and Tfh cells in systemic lupus erythematosus patients

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease that often involves abnormal activation of regulatory IFN genes and regulation of B cells by CD4+ T cells. Radical S-adenosyl methionine domain containing 2 (RSAD2) is a viral suppressor protein regulated by type I IFN, and it has been proven to play an important regulatory role in SLE. However, the mechanism by which RSAD2 participates in the pathogenesis of SLE is unclear. In this study, we observed higher expression levels of RSAD2 in CD4+ T-cell subsets from the peripheral blood of SLE patients than in those from healthy controls by bioinformatics analysis and validation experiments. We analyzed the expression of RSAD2 in CD4+ T cells of patients with SLE and other autoimmune diseases. In addition, we found that the expression of RSAD2 in CD4+ T cells might be regulated by IFN-α, and RSAD2 significantly affected the differentiation of Th17 cells and T follicular helper (Tfh) cells. Our findings underlined that RSAD2 may promote B-cell activation by promoting the differentiation of Th17 and Tfh cells in SLE patients, a process that is regulated by IFN-α.

Pathogenesis of systemic lupus erythematosus: risks, mechanisms and therapeutic targets

Research elucidating the pathogenesis of systemic lupus erythematosus (SLE) has defined two critical families of mediators, type I interferon (IFN-I) and autoantibodies targeting nucleic acids and nucleic acid-binding proteins, as fundamental contributors to the disease. On the fertile background of significant genetic risk, a triggering stimulus, perhaps microbial, induces IFN-I, autoantibody production or most likely both. When innate and adaptive immune system cells are engaged and collaborate in the autoimmune response, clinical SLE can develop. This review describes recent data from genetic analyses of patients with SLE, along with current studies of innate and adaptive immune function that contribute to sustained IFN-I pathway activation, immune activation and autoantibody production, generation of inflammatory mediators and tissue damage. The goal of these studies is to understand disease mechanisms, identify therapeutic targets and stimulate development of therapeutics that can achieve improved outcomes for patients.

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.

Risk and Factors associated with disease manifestations in systemic lupus erythematosus - lupus nephritis (RIFLE-LN): a ten-year risk prediction strategy derived from a cohort of 1652 patients

Objectives

Lupus nephritis (LN) remains one of the most severe manifestations in patients with systemic lupus erythematosus (SLE). Onset and overall LN risk among SLE patients remains considerably difficult to predict. Utilizing a territory-wide longitudinal cohort of over 10 years serial follow-up data, we developed and validated a risk stratification strategy to predict LN risk among Chinese SLE patients - Risk and Factors associated with disease manifestations in systemic Lupus Erythematosus - Lupus Nephritis (RIFLE-LN).

Methods

Demographic and longitudinal data including autoantibody profiles, clinical manifestations, major organ involvement, LN biopsy results and outcomes were documented. Association analysis was performed to identify factors associated with LN. Regression modelling was used to develop a prediction model for 10-year risk of LN and thereafter validated.

Results

A total of 1652 patients were recruited: 1382 patients were assigned for training and validation of the RIFLE-LN model; while 270 were assigned for testing. The median follow-up duration was 21 years. In the training and validation cohort, 845 (61%) of SLE patients developed LN. Cox regression and log rank test showed significant positive association between male sex, age of SLE onset and anti-dsDNA positivity. These factors were thereafter used to develop RIFLE-LN. The algorithm was tested in 270 independent patients and showed good performance (AUC = 0·70).

Conclusion

By using male sex, anti-dsDNA positivity, age of SLE onset and SLE duration; RIFLE-LN can predict LN among Chinese SLE patients with good performance. We advocate its potential utility in guiding clinical management and disease monitoring. Further validation studies in independent cohorts are required.

Evaluating significance of European-associated index SNPs in the East Asian population for 31 complex phenotypes

BACKGROUND

Genome-wide association studies (GWASs) have identified many single-nucleotide polymorphisms (SNPs) associated with complex phenotypes in the European (EUR) population; however, the extent to which EUR-associated SNPs can be generalized to other populations such as East Asian (EAS) is not clear.

RESULTS

By leveraging summary statistics of 31 phenotypes in the EUR and EAS populations, we first evaluated the difference in heritability between the two populations and calculated the trans-ethnic genetic correlation. We observed the heritability estimates of some phenotypes varied substantially across populations and 53.3% of trans-ethnic genetic correlations were significantly smaller than one. Next, we examined whether EUR-associated SNPs of these phenotypes could be identified in EAS using the trans-ethnic false discovery rate method while accounting for winner's curse for SNP effect in EUR and difference of sample sizes in EAS. We found on average 54.5% of EUR-associated SNPs were also significant in EAS. Furthermore, we discovered non-significant SNPs had higher effect heterogeneity, and significant SNPs showed more consistent linkage disequilibrium and allele frequency patterns between the two populations. We also demonstrated non-significant SNPs were more likely to undergo natural selection.

CONCLUSIONS

Our study revealed the extent to which EUR-associated SNPs could be significant in the EAS population and offered deep insights into the similarity and diversity of genetic architectures underlying phenotypes in distinct ancestral groups.

Antibodies against Phosphorylcholine-Implications for Chronic Inflammatory Diseases

Atherosclerosis and its main consequence, cardiovascular disease (CVD) are nowadays regarded as chronic inflammatory disease conditions, and CVD is the main cause of death in the world. Other examples of chronic inflammation are rheumatic and other autoimmune conditions, but also diabetes, obesity, and even osteoarthritis among others. In addition, infectious diseases can have traits in common with these conditions. Systemic lupus erythematosus (SLE) is a prototypical autoimmune disease, where atherosclerosis is increased and the risk of CVD is very high. This is a clinical problem but could also shed light on the role of the immune system in atherosclerosis and CVD. Underlying mechanisms are of major interest and these are only partially known. Phosphorylcholine (PC) is a small lipid-related antigen, which is both a danger associated molecular pattern (DAMP), and a pathogen associated molecular pattern (PAMP). Antibodies against PC are ubiquitous and 5-10% of circulating IgM is IgM anti-PC. Anti-PC, especially IgM and IgG1 anti-PC, has been associated with protection in the chronic inflammatory conditions mentioned above, and develops during the first years of life, while being present at very low levels at birth. Animal experiments with immunization to raise anti-PC ameliorate atherosclerosis and other chronic inflammatory conditions. Potential mechanisms include anti-inflammatory, immune modulatory, clearance of dead cells and protection against infectious agents. An intriguing possibility is to raise anti-PC levels through immunization, to prevent and/or ameliorate chronic inflammation.