Altered expression of genes controlling metabolism characterizes the tissue response to immune injury in lupus.

To compare lupus pathogenesis in disparate tissues, we analyzed gene expression profiles of human discoid lupus erythematosus (DLE) and lupus nephritis (LN). We found common increases in myeloid cell-defining gene sets and decreases in genes controlling glucose and lipid metabolism in lupus-affected skin and kidney. Regression models in DLE indicated increased glycolysis was correlated with keratinocyte, endothelial, and inflammatory cell transcripts, and decreased tricarboxylic (TCA) cycle genes were correlated with the keratinocyte signature. In LN, regression models demonstrated decreased glycolysis and TCA cycle genes were correlated with increased endothelial or decreased kidney cell transcripts, respectively. Less severe glomerular LN exhibited similar alterations in metabolism and tissue cell transcripts before monocyte/myeloid cell infiltration in some patients. Additionally, changes to mitochondrial and peroxisomal transcripts were associated with specific cells rather than global signal changes. Examination of murine LN gene expression demonstrated metabolic changes were not driven by acute exposure to type I interferon and could be restored after immunosuppression. Finally, expression of HAVCR1, a tubule damage marker, was negatively correlated with the TCA cycle signature in LN models. These results indicate that altered metabolic dysfunction is a common, reversible change in lupus-affected tissues and appears to reflect damage downstream of immunologic processes

Identification and Characterization of Post-activated B Cells in Systemic Autoimmune Diseases

Autoimmune diseases (AID) such as systemic lupus erythematosus (SLE), primary Sjögren's syndrome (pSS), and rheumatoid arthritis (RA) are chronic inflammatory diseases in which abnormalities of B cell function play a central role. Although it is widely accepted that autoimmune B cells are hyperactive in vivo, a full understanding of their functional status in AID has not been delineated. Here, we present a detailed analysis of the functional capabilities of AID B cells and dissect the mechanisms underlying altered B cell function. Upon BCR activation, decreased spleen tyrosine kinase (Syk) and Bruton's tyrosine kinase (Btk) phosphorylation was noted in AID memory B cells combined with constitutive co-localization of CD22 and protein tyrosine phosphatase (PTP) non-receptor type 6 (SHP-1) along with hyporesponsiveness to TLR9 signaling, a Syk-dependent response. Similar BCR hyporesponsiveness was also noted specifically in SLE CD27- B cells together with increased PTP activities and increased transcripts for PTPN2, PTPN11, PTPN22, PTPRC, and PTPRO in SLE B cells. Additional studies revealed that repetitive BCR stimulation of normal B cells can induce BCR hyporesponsiveness and that tissue-resident memory B cells from AID patients also exhibited decreased responsiveness immediately ex vivo, suggesting that the hyporesponsive status can be acquired by repeated exposure to autoantigen(s) in vivo. Functional studies to overcome B cell hyporesponsiveness revealed that CD40 co-stimulation increased BCR signaling, induced proliferation, and downregulated PTP expression (PTPN2, PTPN22, and receptor-type PTPs). The data support the conclusion that hyporesponsiveness of AID and especially SLE B cells results from chronic in vivo stimulation through the BCR without T cell help mediated by CD40-CD154 interaction and is manifested by decreased phosphorylation of BCR-related proximal signaling molecules and increased PTPs. The hyporesponsiveness of AID B cells is similar to a form of functional anergy

Nucleic Acid-Sensing and Interferon-Inducible Pathways Show Differential Methylation in MZ Twins Discordant for Lupus and Overexpression in Independent Lupus Samples: Implications for Pathogenic Mechanism and Drug Targeting.

Systemic lupus erythematosus (SLE) is a chronic, multisystem, autoimmune inflammatory disease with genomic and non-genomic contributions to risk. We hypothesize that epigenetic factors are a significant contributor to SLE risk and may be informative for identifying pathogenic mechanisms and therapeutic targets. To test this hypothesis while controlling for genetic background, we performed an epigenome-wide analysis of DNA methylation in genomic DNA from whole blood in three pairs of female monozygotic (MZ) twins of European ancestry, discordant for SLE. Results were replicated on the same array in four cell types from a set of four Danish female MZ twin pairs discordant for SLE. Genes implicated by the epigenetic analyses were then evaluated in 10 independent SLE gene expression datasets from the Gene Expression Omnibus (GEO). There were 59 differentially methylated loci between unaffected and affected MZ twins in whole blood, including 11 novel loci. All but two of these loci were hypomethylated in the SLE twins relative to the unaffected twins. The genes harboring these hypomethylated loci exhibited increased expression in multiple independent datasets of SLE patients. This pattern was largely consistent regardless of disease activity, cell type, or renal tissue type. The genes proximal to CpGs exhibiting differential methylation (DM) in the SLE-discordant MZ twins and exhibiting differential expression (DE) in independent SLE GEO cohorts (DM-DE genes) clustered into two pathways: the nucleic acid-sensing pathway and the type I interferon pathway. The DM-DE genes were also informatically queried for potential gene-drug interactions, yielding a list of 41 drugs including a known SLE therapy. The DM-DE genes delineate two important biologic pathways that are not only reflective of the heterogeneity of SLE but may also correlate with distinct IFN responses that depend on the source, type, and location of nucleic acid molecules and the activated receptors in individual patients. Cell- and tissue-specific analyses will be critical to the understanding of genetic factors dysregulating the nucleic acid-sensing and IFN pathways and whether these factors could be appropriate targets for therapeutic intervention

Differential regulation of the interferon response in systemic lupus erythematosus distinguishes patients of Asian ancestry.

OBJECTIVES: Type I interferon (IFN) plays a role in the pathogenesis of systemic lupus erythematosus (SLE), but insufficient attention has been directed to the differences in IFN responses between ancestral populations. Here, we explored the expression of the interferon gene signatures (IGSs) in SLE patients of European ancestry (EA) and Asian ancestry (AsA). METHODS: We used gene set variation analysis with multiple IGS encompassing the response to both type 1 and type 2 IFN in isolated CD14+ monocytes, CD19+B cells, CD4+T cells and Natural Killer (NK) cells from patients with SLE stratified by self-identified ancestry. The expression of genes upstream of the IGS and influenced by lupus-associated risk alleles was also examined. Lastly, we employed machine learning (ML) models to assess the most important features classifying patients by disease activity. RESULTS: AsA patients with SLE exhibited greater enrichment in the IFN core and IFNA2 IGS compared with EA patients in all cell types examined and, in the presence and absence of autoantibodies. Overall, AsA patients with SLE demonstrated higher expression of genes upstream of the IGS than EA counterparts. ML with feature importance analysis indicated that IGS expression in NK cells, anti-dsDNA, complement levels and AsA status contributed to disease activity. CONCLUSIONS: AsA patients with SLE exhibited higher IGS than EA patients in all cell types regardless of autoantibody status, with enhanced expression of genetically associated genes upstream of the IGS potentially contributing. AsA, along with the IGS in NK cells, anti-dsDNA and complement, independently influenced SLE disease activity

Molecular endotypes of type 1 and type 2 SLE

Objective To character the molecular landscape of patients with type 1 and type 2 SLE by analysing gene expression profiles from peripheral blood.Methods Full transcriptomic RNA sequencing was carried out on whole blood samples from 18 subjects with SLE selected by the presence of manifestations typical of type 1 and type 2 SLE. The top 5000 row variance genes were analysed by Multiscale Embedded Gene Co-expression Network Analysis to generate gene co-expression modules that were functionally annotated and correlated with various demographic traits, clinical features and laboratory measures.Results Expression of specific gene co-expression modules correlated with individual features of type 1 and type 2 SLE and also effectively segregated samples from patients with type 1 SLE from those with type 2 SLE. Unique type 1 SLE enrichment included interferon, monocytes, T cells, cell cycle and neurotransmitter pathways, whereas unique type 2 SLE enrichment included B cells and metabolic and neuromuscular pathways. Gene co-expression modules of patients with type 2 SLE were identified in subsets of previously reported patients with inactive SLE and idiopathic fibromyalgia (FM) and also identified subsets of patients with active SLE with a greater frequency of severe fatigue.Conclusion Gene co-expression analysis successfully identified unique transcriptional patterns that segregate type 1 SLE from type 2 SLE and further identified type 2 molecular features in patients with inactive SLE or FM and with active SLE with severe fatigue

Differential regulation of the interferon response in systemic lupus erythematosus distinguishes patients of Asian ancestry

Objectives Type I interferon (IFN) plays a role in the pathogenesis of systemic lupus erythematosus (SLE), but insufficient attention has been directed to the differences in IFN responses between ancestral populations. Here, we explored the expression of the interferon gene signatures (IGSs) in SLE patients of European ancestry (EA) and Asian ancestry (AsA).Methods We used gene set variation analysis with multiple IGS encompassing the response to both type 1 and type 2 IFN in isolated CD14+ monocytes, CD19+B cells, CD4+T cells and Natural Killer (NK) cells from patients with SLE stratified by self-identified ancestry. The expression of genes upstream of the IGS and influenced by lupus-associated risk alleles was also examined. Lastly, we employed machine learning (ML) models to assess the most important features classifying patients by disease activity.Results AsA patients with SLE exhibited greater enrichment in the IFN core and IFNA2 IGS compared with EA patients in all cell types examined and, in the presence and absence of autoantibodies. Overall, AsA patients with SLE demonstrated higher expression of genes upstream of the IGS than EA counterparts. ML with feature importance analysis indicated that IGS expression in NK cells, anti-dsDNA, complement levels and AsA status contributed to disease activity.Conclusions AsA patients with SLE exhibited higher IGS than EA patients in all cell types regardless of autoantibody status, with enhanced expression of genetically associated genes upstream of the IGS potentially contributing. AsA, along with the IGS in NK cells, anti-dsDNA and complement, independently influenced SLE disease activity

Identification of alterations in macrophage activation associated with disease activity in systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is characterized by abnormalities in B cell and T cell function, but the role of disturbances in the activation status of macrophages (Mϕ) has not been well described in human patients. To address this, gene expression profiles from isolated lymphoid and myeloid populations were analyzed to identify differentially expressed (DE) genes between healthy controls and patients with either inactive or active SLE. While hundreds of DE genes were identified in B and T cells of active SLE patients, there were no DE genes found in B or T cells from patients with inactive SLE compared to healthy controls. In contrast, large numbers of DE genes were found in myeloid cells (MC) from both active and inactive SLE patients. Among the DE genes were several known to play roles in Mϕ activation and polarization, including the M1 genes STAT1 and SOCS3 and the M2 genes STAT3, STAT6, and CD163. M1-associated genes were far more frequent in data sets from active versus inactive SLE patients. To characterize the relationship between Mϕ activation and disease activity in greater detail, weighted gene co-expression network analysis (WGCNA) was used to identify modules of genes associated with clinical activity in SLE patients. Among these were disease activity-correlated modules containing activation signatures of predominantly M1-associated genes. No disease activity-correlated modules were enriched in M2-associated genes. Pathway and upstream regulator analysis of DE genes from both active and inactive SLE MC were cross-referenced with high-scoring hits from the drug discovery Library of Integrated Network-based Cellular Signatures (LINCS) to identify new strategies to treat both stages of SLE. A machine learning approach employing MC gene modules and a generalized linear model was able to predict the disease activity status in unrelated gene expression data sets. In summary, altered MC gene expression is characteristic of both active and inactive SLE. However, disease activity is associated with an alteration in the activation of MC, with a bias toward the M1 proinflammatory phenotype. These data suggest that while hyperactivity of B cells and T cells is associated with active SLE, MC potentially direct flare-ups and remission by altering their activation status toward the M1 state