Several Critical Cell Types, Tissues, and Pathways Are Implicated in Genome-Wide Association Studies for Systemic Lupus Erythematosus

Investigating the causal relationship and potential shared diagnostic genes between primary biliary cholangitis and systemic lupus erythematosus using bidirectional Mendelian randomization and transcriptomic analyses

Background

The co-occurrence of primary biliary cholangitis (PBC) and systemic lupus erythematosus (SLE) has been consistently reported in observational studies. Nevertheless, the underlying causal correlation between these two conditions still needs to be established.

Methods

We performed a bidirectional two-sample Mendelian randomization (MR) study to assess their causal association. Five MR analysis methods were utilized for causal inference, with inverse-variance weighted (IVW) selected as the primary method. The Mendelian Randomization Pleiotropy RESidual Sum and Outlier (MR-PRESSO) and the IVW Radial method were applied to exclude outlying SNPs. To assess the robustness of the MR results, five sensitivity analyses were carried out. Multivariable MR (MVMR) analysis was also employed to evaluate the effect of possible confounders. In addition, we integrated transcriptomic data from PBC and SLE, employing Weighted Gene Co-expression Network Analysis (WGCNA) to explore shared genes between the two diseases. Then, we used Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment methods to perform on the shared genes. The Least Absolute Shrinkage and Selection Operator (LASSO) regression algorithm was utilized to identify potential shared diagnostic genes. Finally, we verified the potential shared diagnostic genes in peripheral blood mononuclear cells (PBMCs)-specific cell populations of SLE patients by single-cell analysis.

Results

Our MR study provided evidence that PBC had a causal relationship with SLE (IVW, OR: 1.347, 95% CI: 1.276 - 1.422, P < 0.001) after removing outliers (MR-PRESSO, rs35464393, rs3771317; IVW Radial, rs11065987, rs12924729, rs3745516). Conversely, SLE also had a causal association with PBC (IVW, OR: 1.225, 95% CI: 1.141 - 1.315, P < 0.001) after outlier correction (MR-PRESSO, rs11065987, rs3763295, rs7774434; IVW Radial, rs2297067). Sensitivity analyses confirmed the robustness of the MR findings. MVMR analysis indicated that body mass index (BMI), smoking and drinking were not confounding factors. Moreover, bioinformatic analysis identified PARP9, ABCA1, CEACAM1, and DDX60L as promising diagnostic biomarkers for PBC and SLE. These four genes are highly expressed in CD14+ monocytes in PBMCs of SLE patients and potentially associated with innate immune responses and immune activation.

Conclusion

Our study confirmed the bidirectional causal relationship between PBC and SLE and identified PARP9, ABCA1, CEACAM1, and DDX60L genes as the most potentially shared diagnostic genes between the two diseases, providing insights for the exploration of the underlying mechanisms of these disorders.

Training vs. Tolerance: The Yin/Yang of the Innate Immune System

For almost nearly a century, memory functions have been attributed only to acquired immune cells. Lately, this paradigm has been challenged by an increasing number of studies revealing that innate immune cells are capable of exhibiting memory-like features resulting in increased responsiveness to subsequent challenges, a process known as trained immunity (known also as innate memory). In contrast, the refractory state of endotoxin tolerance has been defined as an immunosuppressive state of myeloid cells portrayed by a significant reduction in the inflammatory capacity. Both training as well tolerance as adaptive features are reported to be accompanied by epigenetic and metabolic alterations occurring in cells. While training conveys proper protection against secondary infections, the induction of endotoxin tolerance promotes repairing mechanisms in the cells. Consequently, the inappropriate induction of these adaptive cues may trigger maladaptive effects, promoting an increased susceptibility to secondary infections—tolerance, or contribute to the progression of the inflammatory disorder—trained immunity. This review aims at the discussion of these opposing manners of innate immune and non-immune cells, describing the molecular, metabolic and epigenetic mechanisms involved and interpreting the clinical implications in various inflammatory pathologies.

Trained immunity as a possible newcomer in autoinflammatory and autoimmune diseases pathophysiology

Autoimmune disorders have been well characterized over the years and many pathways-but not all of them-have been found to explain their pathophysiology. Autoinflammatory disorders, on the other hand, are still hiding most of their molecular and cellular mechanisms. During the past few years, a newcomer has challenged the idea that only adaptive immunity could display memory response. Trained immunity is defined by innate immune responses that are faster and stronger to a second stimulus than to the first one, being the same or not. In response to the trained immunity inducer, and through metabolic and epigenetic changes of hematopoietic stem and progenitor cells in the bone marrow that are transmitted to their cellular progeny (peripheral trained immunity), or directly of tissue-resident cells (local innate immunity), innate cells responsiveness and functions upon stimulation are improved in the long-term. Innate immunity can be beneficial, but it could also be detrimental when maladaptive. Here, we discuss how trained immunity could contribute to the physiopathology of autoimmune and autoinflammatory diseases.

Biological insights into systemic lupus erythematosus through an immune cell-specific transcriptome-wide association study

OBJECTIVE

Genome-wide association studies (GWAS) have identified >100 risk loci for systemic lupus erythematosus (SLE), but the disease genes at most loci remain unclear, hampering translation of these genetic discoveries. We aimed to prioritise genes underlying the 110 SLE loci that were identified in the latest East Asian GWAS meta-analysis.

METHODS

We built gene expression predictive models in blood B cells, CD4+ and CD8+ T cells, monocytes, natural killer cells and peripheral blood cells of 105 Japanese individuals. We performed a transcriptome-wide association study (TWAS) using data from the latest genome-wide association meta-analysis of 208 370 East Asians and searched for candidate genes using TWAS and three data-driven computational approaches.

RESULTS

TWAS identified 171 genes for SLE (p<1.0×10-5); 114 (66.7%) showed significance only in a single cell type; 127 (74.3%) were in SLE GWAS loci. TWAS identified a strong association between CD83 and SLE (p<7.7×10-8). Meta-analysis of genetic associations in the existing 208 370 East Asian and additional 1498 cases and 3330 controls found a novel single-variant association at rs72836542 (OR=1.11, p=4.5×10-9) around CD83. For the 110 SLE loci, we identified 276 gene candidates, including 104 genes at recently-identified SLE novel loci. We demonstrated in vitro that putative causal variant rs61759532 exhibited an allele-specific regulatory effect on ACAP1, and that presence of the SLE risk allele decreased ACAP1 expression.

CONCLUSIONS

Cell-level TWAS in six types of immune cells complemented SLE gene discovery and guided the identification of novel genetic associations. The gene findings shed biological insights into SLE genetic associations.

IFN-γ Mediates the Development of Systemic Lupus Erythematosus

Objective

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease that can affect all organs in the body. It is characterized by overexpression of antibodies against autoantigen. Although previous bioinformatics analyses have identified several genetic factors underlying SLE, they did not discriminate between naive and individuals exposed to anti-SLE drugs. Here, we evaluated specific genes and pathways in active and recently diagnosed SLE population.

Methods

GSE46907 matrix downloaded from Gene Expression Omnibus (GEO) was analyzed using R, Metascape, STRING, and Cytoscape to identify differentially expressed genes (DEGs), enrichment pathways, protein-protein interaction (PPI), and hub genes between naive SLE individuals and healthy controls.

Results

A total of 134 DEGs were identified, in which 29 were downregulated, whereas 105 were upregulated in active and newly diagnosed SLE cases. GO term analysis revealed that transcriptional induction of the DEGs was particularly enhanced in response to secretion of interferon-γ and interferon-α and regulation of cytokine production innate immune responses among others. KEGG pathway analysis showed that the expression of DEGs was particularly enhanced in interferon signaling, IFN antiviral responses by activated genes, class I major histocompatibility complex (MHC-I) mediated antigen processing and presentation, and amyloid fiber formation. STAT1, IRF7, MX1, OASL, ISG15, IFIT3, IFIH1, IFIT1, OAS2, and GBP1 were the top 10 DEGs.

Conclusions

Our findings suggest that interferon-related gene expression and pathways are common features for SLE pathogenesis, and IFN-γ and IFN-γ-inducible GBP1 gene in naive SLE were emphasized. Together, the identified genes and cellular pathways have expanded our understanding on the mechanism underlying development of SLE. They have also opened a new frontier on potential biomarkers for diagnosis, biotherapy, and prognosis for SLE.

The epigenetic face of lupus: Focus on antigen-presenting cells

In recent years, epigenetic mechanisms became widely known due to their ability to regulate and maintain physiological processes such as cell growth, development, differentiation and genomic stability. When dysregulated, epigenetic mechanisms, may introduce gene expression changes and disturbance in immune homeostasis leading to autoimmune diseases. Systemic lupus erythematosus (SLE), the most extensively studied autoimmune disorder, has already been correlated with epigenetic modifications, especially in T cells. Since these cell rely on antigen presentation, it may be assumed that erroneous activity of antigen-presenting cells (APCs), culminates in T cell abnormalities. In this review we summarize and discuss the epigenetic modifications in SLE affected APCs, with the focus on dendritic cells (DCs), B cells and monocytes. Unravelling this aspect of SLE pathogenesis, might result in identification of new disease biomarkers and putative therapeutic approaches.

Overall Downregulation of mRNAs and Enrichment of H3K4me3 Change Near Genome-Wide Association Study Signals in Systemic Lupus Erythematosus: Cell-Specific Effects

This study was designed to define gene expression and H3K4me3 histone modifications in T cells, B cells, and monocytes in systemic lupus erythematosus (SLE). Array studies of total peripheral blood mononuclear cells have demonstrated gene expression signatures related to neutrophils, interferon, and other inflammatory pathways. It is not clear how consistent these effects are across different cell types. In this study, RNA-seq and chromatin immunoprecipitation-seq were utilized to identify gene expression patterns and H3K4me3 histone modifications related to promoter activation in SLE. Across the three cell types, there was 55% concordance for gene expression changes related to SLE. Key conserved pathways were ribosome biogenesis among upregulated genes and heat shock response among downregulated genes. ETS family transcription factors (TFs) and STAT1 were revealed as common regulators by position weight matrices. When epigenetic changes were leveraged with gene expression, the pivotal TFs ATF3 and FOS were defined with ATF3 also cross-referencing with gene expression-identified TFs. Genome-wide association study (GWAS) single nucleotide polymorphisms associated with SLE were cross-referenced with both mRNA and H3K4me3 changes in SLE. Baseline mRNA expression and H3K4me3 peak height was higher at sites that cross-referenced with GWAS signals, however, all three cell types exhibited an overall decrease in expression of GWAS-associated RNAs differentially expressed in SLE. H3K4me3 changes in SLE were also enriched in GWAS-associated sites. In summary, the SLE disease process is associated with both shared and cell-specific changes in gene expression and epigenetics. Surprisingly, GWAS-associated RNAs were overall markedly decreased across all three cell types. TF analysis identified ATF3, FOS, STAT1, and ETS family members as critical, all pathways with a recognized relationship to the SLE disease process. GWAS signals clearly mark both cell-type specific changes in SLE as well as concordant changes across all three cell types. Interpretation of single nucleotide polymorphism effects in SLE will require tissue-specific mechanistic studies and therapeutics will require mechanistic studies in multiple cell types.

The Potential Role of Trained Immunity in Autoimmune and Autoinflammatory Disorders

During induction of trained immunity, monocytes and macrophages undergo a functional and transcriptional reprogramming toward increased activation. Important rewiring of cellular metabolism of the myeloid cells takes place during induction of trained immunity, including a shift toward glycolysis induced through the mTOR pathway, as well as glutaminolysis and cholesterol synthesis. Subsequently, this leads to modulation of the function of epigenetic enzymes, resulting in important changes in chromatin architecture that enables increased gene transcription. However, in addition to the beneficial effects of trained immunity as a host defense mechanism, we hypothesize that trained immunity also plays a deleterious role in the induction and/or maintenance of autoimmune and autoinflammatory diseases if inappropriately activated.

Another Round of "Clue" to Uncover the Mystery of Complex Traits

A plethora of genetic association analyses have identified several genetic risk loci. Technological and statistical advancements have now led to the identification of not only common genetic variants, but also low-frequency variants, structural variants, and environmental factors, as well as multi-omics variations that affect the phenotypic variance of complex traits in a population, thus referred to as complex trait architecture. The concept of heritability, or the proportion of phenotypic variance due to genetic inheritance, has been studied for several decades, but its application is mainly in addressing the narrow sense heritability (or additive genetic component) from Genome-Wide Association Studies (GWAS). In this commentary, we reflect on our perspective on the complexity of understanding heritability for human traits in comparison to model organisms, highlighting another round of clues beyond GWAS and an alternative approach, investigating these clues comprehensively to help in elucidating the genetic architecture of complex traits.

Confirmation of five novel susceptibility loci for systemic lupus erythematosus (SLE) and integrated network analysis of 82 SLE susceptibility loci

We recently identified ten novel SLE susceptibility loci in Asians and uncovered several additional suggestive loci requiring further validation. This study aimed to replicate five of these suggestive loci in a Han Chinese cohort from Hong Kong, followed by meta-analysis (11,656 cases and 23,968 controls) on previously reported Asian and European populations, and to perform bioinformatic analyses on all 82 reported SLE loci to identify shared regulatory signatures. We performed a battery of analyses for these five loci, as well as joint analyses on all 82 SLE loci. All five loci passed genome-wide significance: MYNN (rs10936599, Pmeta = 1.92 × 10-13, OR = 1.14), ATG16L2 (rs11235604, Pmeta = 8.87 × 10 -12, OR = 0.78), CCL22 (rs223881, Pmeta = 5.87 × 10-16, OR = 0.87), ANKS1A (rs2762340, Pmeta = 4.93 × 10-15, OR = 0.87) and RNASEH2C (rs1308020, Pmeta = 2.96 × 10-19, OR = 0.84) and co-located with annotated gene regulatory elements. The novel loci share genetic signatures with other reported SLE loci, including effects on gene expression, transcription factor binding, and epigenetic characteristics. Most (56%) of the correlated (r2 > 0.8) SNPs from the 82 SLE loci were implicated in differential expression (9.81 × 10-198 < P < 5 × 10-3) of cis-genes. Transcription factor binding sites for p53, MEF2A and E2F1 were significantly (P < 0.05) over-represented in SLE loci, consistent with apoptosis playing a critical role in SLE. Enrichment analysis revealed common pathways, gene ontology, protein domains, and cell type-specific expression. In summary, we provide evidence of five novel SLE susceptibility loci. Integrated bioinformatics using all 82 loci revealed that SLE susceptibility loci share many gene regulatory features, suggestive of conserved mechanisms of SLE etiopathogenesis.

DNA-reactive B cells in lupus

IgG anti-DNA antibodies are both diagnostic and pathogenic for systemic lupus erythematosus (SLE). They contribute to tissue inflammation through direct tissue binding and to systemic inflammation through activation of Toll-like receptors by nucleic acid-containing immune complexes. IgG DNA-reactive antibodies originate when B cell tolerance mechanisms are impaired. The heterogeneous immune perturbations in SLE lead to the survival and activation of DNA-reactive B cells in various B cell subsets at distinct stages of B cell maturation and differentiation. We propose that the spectrum of B cell alterations and failed tolerance mechanisms for DNA-reactive B cells in lupus patients is best understood by studying genetic risk alleles. This implies that the B cells producing IgG anti-DNA antibodies and the failed tolerance mechanisms(s) will differ across patients. A better understanding of these differences should lead to better patient stratification, improved outcomes of clinical trials, and the identification of novel therapeutic targets.