Mapping systemic lupus erythematosus heterogeneity at the single-cell level

Single-cell analysis with childhood and adult systemic lupus erythematosus

Patients with systemic lupus erythematosus (SLE), a heterogeneous and chronic autoimmune disease, exhibit unique changes in the complex composition and transcriptional signatures of peripheral blood mononuclear cells (PBMCs). While the mechanism of pathogenesis for both childhood-onset SLE (cSLE) and adult-onset SLE (aSLE) remains unclear, cSLE patients are considered more unpredictable and dangerous than aSLE patients. In this study, we analysed single-cell RNA sequencing data (scRNA-seq) to profile the PBMC clusters of cSLE/aSLE patients and matched healthy donors and compared the PBMC composition and transcriptional variations between the two groups. Our analysis revealed that the PBMC composition and transcriptional variations in cSLE patients were similar to those in aSLE patients. Comparative single-cell transcriptome analysis between healthy donors and SLE patients revealed IFITM3, ISG15, IFI16 and LY6E as potential therapeutic targets for both aSLE and cSLE patients. Additionally, we observed that the percentage of pre-B cells (CD34-) was increased in cSLE patients, while the percentage of neutrophil cells was upregulated in aSLE patients. Notably, we found decreased expression of TPM2 in cSLE patients, and similarly, TMEM150B, IQSEC2, CHN2, LRP8 and USP46 were significantly downregulated in neutrophil cells from aSLE patients. Overall, our study highlights the differences in complex PBMC composition and transcriptional profiles between cSLE and aSLE patients, providing potential biomarkers that could aid in diagnosing SLE.

Single-Cell Analyses Reveal the Metabolic Heterogeneity and Plasticity of the Tumor Microenvironment during Head and Neck Squamous Cell Carcinoma Progression

Metabolic reprogramming is a hallmark of cancer. In addition to metabolic alterations in the tumor cells, multiple other metabolically active cell types in the tumor microenvironment (TME) contribute to the emergence of a tumor-specific metabolic milieu. Here, we defined the metabolic landscape of the TME during the progression of head and neck squamous cell carcinoma (HNSCC) by performing single-cell RNA sequencing on 26 human patient specimens, including normal tissue, precancerous lesions, early stage cancer, advanced-stage cancer, lymph node metastases, and recurrent tumors. The analysis revealed substantial heterogeneity at the transcriptional, developmental, metabolic, and functional levels in different cell types. SPP1+ macrophages were identified as a protumor and prometastatic macrophage subtype with high fructose and mannose metabolism, which was further substantiated by integrative analysis and validation experiments. An inhibitor of fructose metabolism reduced the proportion of SPP1+ macrophages, reshaped the immunosuppressive TME, and suppressed tumor growth. In conclusion, this work delineated the metabolic landscape of HNSCC at a single-cell resolution and identified fructose metabolism as a key metabolic feature of a protumor macrophage subpopulation. Significance: Fructose and mannose metabolism is a metabolic feature of a protumor and prometastasis macrophage subtype and can be targeted to reprogram macrophages and the microenvironment of head and neck squamous cell carcinoma.

Single-cell landscape of peripheral immune response in patients with anti-melanoma differentiation-associated gene 5 dermatomyositis

OBJECTIVE

Anti-melanoma differentiation-associated gene 5 (Anti-MDA5)-positive DM is a rare but life-threatening autoimmune disorder that is associated with a high risk of developing rapidly progressive interstitial lung disease. Current empirical therapies offer limited benefit in terms of patient survival, as little is known about the aetiology of anti-MDA5 DM. To best understand its immune landscape, we applied single-cell RNA sequencing (scRNA-seq) to peripheral blood samples from DM patients and healthy controls.

METHODS

Peripheral blood mononuclear cells (PBMCs) from eight DM patients (including three distinct subtypes of DM) and two healthy donors were sequenced using the 10X Genomics platform. Additional scRNA-seq data for four healthy donors were incorporated for further bioinformatic analysis.

RESULTS

Aberrantly increased proportions of CD14+ monocytes and plasma cells were observed in anti-MDA5 DM PBMC samples. Moreover, we found an overactivated type I IFN response and antiviral immunity in both innate and adaptive immune cells derived from anti-MDA5 DM patients that was positively correlated with disease severity. Importantly, a unique subset of CD14+ monocytes that highly expressed IFN alpha-inducible protein 27 (IFI27), a biomarker for viral infection, and IFN induced with helicase C domain 1 (IFIH1, which encodes MDA5) was specifically identified in anti-MDA5 DM samples for the first time.

CONCLUSION

Our study has illustrated the peripheral immune cell atlas of a number of DM subtypes, has provided compelling evidence for a viral infection-derived origin for anti-MDA5 DM, and has indicated potential targets for innovative therapeutic interventions.

Decoding the genetic landscape of juvenile dermatomyositis: insights from phosphorylation-associated single nucleotide polymorphisms

Genome-wide association studies (GWASs) have identified genetic susceptibility loci associated with juvenile dermatomyositis (JDM). Single nucleotide polymorphisms related to phosphorylation (phosSNPs) are critical nonsynonymous mutations exerting substantial influence on gene expression regulation. The aim of this study was to identify JDM susceptibility genes in the GWAS loci by the use of phosSNPs. We explored quantitative trait loci (QTLs) among the phosSNPs associated with JDM using data from eQTL (bulk tissues and single-cell) and pQTL studies. For gene expression and protein levels significantly influenced by JDM-associated phosSNPs, we assessed their associations with JDM through MR analyses. Additionally, we conducted differential expression gene analyses, incorporating single-cell transcriptomic profiling of 6 JDM cases and 11 juvenile controls (99,396 cells). We identified 31 phosSNPs situated in the 6p21 locus that were associated with JDM. Half of these phosSNPs showed effects on gene expression in various cells and circulating protein levels. In MR analyses, we established associations between the expression levels of pivotal JDM-associated genes, including MICB, C4A, HLA-DRB1, HLA-DRB5, and PSMB9, in skin, muscle, or blood cells and circulating levels of C4A, with JDM. Utilizing single-cell eQTL data, we identified a total of 276 association signals across 14 distinct immune cell types for 28 phosSNPs. Further insights were gained through single-cell differential expression analysis, revealing differential expression of PSMB9, HLA-A, HLA-B, HLA-C, HLA-DPB1, HLA-DQA1, HLA-DQB1, and HLA-DRB1 in immune cells. The present study pinpointed phosSNPs within susceptibility genes for JDM and unraveled the intricate relationships among these SNPs, gene expression levels, and JDM.

Single-cell omics: experimental workflow, data analyses and applications

Cells are the fundamental units of biological systems and exhibit unique development trajectories and molecular features. Our exploration of how the genomes orchestrate the formation and maintenance of each cell, and control the cellular phenotypes of various organismsis, is both captivating and intricate. Since the inception of the first single-cell RNA technology, technologies related to single-cell sequencing have experienced rapid advancements in recent years. These technologies have expanded horizontally to include single-cell genome, epigenome, proteome, and metabolome, while vertically, they have progressed to integrate multiple omics data and incorporate additional information such as spatial scRNA-seq and CRISPR screening. Single-cell omics represent a groundbreaking advancement in the biomedical field, offering profound insights into the understanding of complex diseases, including cancers. Here, we comprehensively summarize recent advances in single-cell omics technologies, with a specific focus on the methodology section. This overview aims to guide researchers in selecting appropriate methods for single-cell sequencing and related data analysis.

Molecular subtypes explain lupus epigenomic heterogeneity unveiling new regulatory genetic risk variants

The heterogeneity of systemic lupus erythematosus (SLE) can be explained by epigenetic alterations that disrupt transcriptional programs mediating environmental and genetic risk. This study evaluated the epigenetic contribution to SLE heterogeneity considering molecular and serological subtypes, genetics and transcriptional status, followed by drug target discovery. We performed a stratified epigenome-wide association studies of whole blood DNA methylation from 213 SLE patients and 221 controls. Methylation quantitative trait loci analyses, cytokine and transcription factor activity - epigenetic associations and methylation-expression correlations were conducted. New drug targets were searched for based on differentially methylated genes. In a stratified approach, a total of 974 differential methylation CpG sites with dependency on molecular subtypes and autoantibody profiles were found. Mediation analyses suggested that SLE-associated SNPs in the HLA region exert their risk through DNA methylation changes. Novel genetic variants regulating DNAm in disease or in specific molecular contexts were identified. The epigenetic landscapes showed strong association with transcription factor activity and cytokine levels, conditioned by the molecular context. Epigenetic signals were enriched in known and novel drug targets for SLE. This study reveals possible genetic drivers and consequences of epigenetic variability on SLE heterogeneity and disentangles the DNAm mediation role on SLE genetic risk and novel disease-specific meQTLs. Finally, novel targets for drug development were discovered.

Pan-cancer single-cell dissection reveals phenotypically distinct B cell subtypes

Characterizing the compositional and phenotypic characteristics of tumor-infiltrating B cells (TIBs) is important for advancing our understanding of their role in cancer development. Here, we establish a comprehensive resource of human B cells by integrating single-cell RNA sequencing data of B cells from 649 patients across 19 major cancer types. We demonstrate substantial heterogeneity in their total abundance and subtype composition and observe immunoglobulin G (IgG)-skewness of antibody-secreting cell isotypes. Moreover, we identify stress-response memory B cells and tumor-associated atypical B cells (TAABs), two tumor-enriched subpopulations with prognostic potential, shared in a pan-cancer manner. In particular, TAABs, characterized by a high clonal expansion level and proliferative capacity as well as by close interactions with activated CD4 T cells in tumors, are predictive of immunotherapy response. Our integrative resource depicts distinct clinically relevant TIB subsets, laying a foundation for further exploration of functional commonality and diversity of B cells in cancer.

CD4+CD57+ senescent T cells as promoters of systemic lupus erythematosus pathogenesis and the therapeutic potential of senolytic BCL-2 inhibitor

Systemic lupus erythematosus (SLE) is a complex autoimmune disease characterized by persistent activation of immune cells and overproduction of autoantibodies. The accumulation of senescent T and B cells has been observed in SLE and other immune-mediated diseases. However, the exact mechanistic pathways contributing to this process in SLE remain incompletely understood. In this study, we found that in SLE patients: (1) the frequency of CD4+CD57+ senescent T cells was significantly elevated and positively correlated with disease activity; (2) the expression levels of B-lymphoma-2 (BCL-2) family and interferon-induced genes (ISGs) were significantly upregulated; and (3) in vitro, the cytokine IL-15 stimulation increased the frequency of senescent CD4+ T cells and upregulated the expression of BCL-2 family and ISGs. Further, treatment with ABT-263 (a senolytic BCL-2 inhibitor) in MRL/lpr mice resulted in decreased: (1) frequency of CD4+CD44hiCD62L-PD-1+CD153+ senescent CD4+ T cells; (2) frequency of CD19+CD11c+T-bet+ age-related B cells; (3) level of serum antinuclear antibody; (4) proteinuria; (5) frequency of Tfh cells; and (6) renal histopathological abnormalities. Collectively, these results indicated a dominant role for CD4+CD57+ senescent CD4+ T cells in the pathogenesis of SLE and senolytic BCL-2 inhibitor ABT-263 may be the potential treatment in ameliorating lupus phenotypes.

PredGCN: a Pruning-enabled Gene-Cell Net for automatic cell annotation of single cell transcriptome data

MOTIVATION

The annotation of cell types from single-cell transcriptomics is essential for understanding the biological identity and functionality of cellular populations. Although manual annotation remains the gold standard, the advent of automatic pipelines has become crucial for scalable, unbiased, and cost-effective annotations. Nonetheless, the effectiveness of these automatic methods, particularly those employing deep learning, significantly depends on the architecture of the classifier and the quality and diversity of the training datasets.

RESULTS

To address these limitations, we present a Pruning-enabled Gene-Cell Net (PredGCN) incorporating a Coupled Gene-Cell Net (CGCN) to enable representation learning and information storage. PredGCN integrates a Gene Splicing Net (GSN) and a Cell Stratification Net (CSN), employing a pruning operation (PrO) to dynamically tackle the complexity of heterogeneous cell identification. Among them, GSN leverages multiple statistical and hypothesis-driven feature extraction methods to selectively assemble genes with specificity for scRNA-seq data while CSN unifies elements based on diverse region demarcation principles, exploiting the representations from GSN and precise identification from different regional homogeneity perspectives. Furthermore, we develop a multi-objective Pareto pruning operation (Pareto PrO) to expand the dynamic capabilities of CGCN, optimizing the sub-network structure for accurate cell type annotation. Multiple comparison experiments on real scRNA-seq datasets from various species have demonstrated that PredGCN surpasses existing state-of-the-art methods, including its scalability to cross-species datasets. Moreover, PredGCN can uncover unknown cell types and provide functional genomic analysis by quantifying the influence of genes on cell clusters, bringing new insights into cell type identification and characterizing scRNA-seq data from different perspectives.

AVAILABILITY AND IMPLEMENTATION

The source code is available at https://github.com/IrisQi7/PredGCN and test data is available at https://figshare.com/articles/dataset/PredGCN/25251163.

Single-cell analysis in rheumatic and allergic diseases: insights for clinical practice

Since the advent of single-cell RNA sequencing (scRNA-seq) methodology, single-cell analysis has become a powerful tool for exploration of cellular networks and dysregulated immune responses in disease pathogenesis. Advanced bioinformatics tools have enabled the combined analysis of scRNA-seq data and information on various cell properties, such as cell surface molecular profiles, chromatin accessibility and spatial information, leading to a deeper understanding of pathology. This Review provides an overview of the achievements in single-cell analysis applied to clinical samples of rheumatic and allergic diseases, including rheumatoid arthritis, systemic lupus erythematosus, systemic sclerosis, allergic airway diseases and atopic dermatitis, with an expanded scope beyond peripheral blood cells to include local diseased tissues. Despite the valuable insights that single-cell analysis has provided into disease pathogenesis, challenges remain in translating single-cell findings into clinical practice and developing personalized treatment strategies. Beyond understanding the atlas of cellular diversity, we discuss the application of data obtained in each study to clinical practice, with a focus on identifying biomarkers and therapeutic targets.

Peripheral blood cells RNA-seq identifies differentially expressed gene network linked to lymphocyte subsets alterations and active lupus nephritis associated with declines in renal function

Background

The aim of this study was to investigate whether quantitative changes in lymphocyte subsets and gene expression in peripheral blood (PB) cells are related to the clinical manifestations and pathogenesis of lupus nephritis (LN).

Methods

We enrolled 95 pediatric-onset SLE patients with renal involvement who presented with 450 clinical episodes suspicious for LN flare. Percentages of lymphocyte subsets at each episode were determined. We stratified 55 of 95 patients as high or low subset group according to the median percentage of each lymphocyte subset and the association with changes in the eGFR (ΔeGFR) were analyzed. Peripheral blood bulk RNA-seq to identify differentially expressed genes (DEGs) in 9 active LN vs. 9 inactive LN patients and the DEG-derived network was constructed by Ingenuity Pathway Analysis (IPA).

Results

The mean ΔeGFR of low NK-low memory CD4+ T-high naive CD4+ T group (31.01 mL/min/1.73 m2) was significantly greater than that of high NK-high memory CD4+ T-low naive CD4+ T group (11.83 mL/min/1.73 m2; P = 0.0175). Kaplan-Meier analysis showed that the median time for ΔeGFR decline to mean ΔeGFR is approximately 10 years for high NK-high memory CD4+ T-low naive CD4+ T group and approximately 5 years for low NK-low memory CD4+ T-high naive CD4+ T group (log-rank test P = 0.0294).

Conclusions

Our study highlighted important connections between DEG-derived network, lymphocyte subset composition, and disease status of LN and GN. A novel scoring system based on lymphocyte subset proportions effectively stratified patients into groups with differential risks for declining renal function.

ETV5 promotes lupus pathogenesis and follicular helper T cell differentiation by inducing osteopontin expression

Follicular helper T (TFH) cells mediate germinal center reactions to generate high affinity antibodies against specific pathogens, and their excessive production is associated with the pathogenesis of systemic autoimmune diseases such as systemic lupus erythematosus (SLE). ETV5, a member of the ETS transcription factor family, promotes TFH cell differentiation in mice. In this study, we examined the role of ETV5 in the pathogenesis of lupus in mice and humans. T cell-specific deletion of Etv5 alleles ameliorated TFH cell differentiation and autoimmune phenotypes in lupus mouse models. Further, we identified SPP1 as an ETV5 target that promotes TFH cell differentiation in both mice and humans. Notably, extracellular osteopontin (OPN) encoded by SPP1 enhances TFH cell differentiation by activating the CD44-AKT signaling pathway. Furthermore, ETV5 and SPP1 levels were increased in CD4+ T cells from patients with SLE and were positively correlated with disease activity. Taken together, our findings demonstrate that ETV5 is a lupus-promoting transcription factor, and secreted OPN promotes TFH cell differentiation.

Active juvenile systemic lupus erythematosus is associated with distinct NK cell transcriptional and phenotypic alterations

While adaptive immune responses have been studied extensively in SLE (systemic lupus erythematosus), there is limited and contradictory evidence regarding the contribution of natural killer (NK) cells to disease pathogenesis. There is even less evidence about the role of NK cells in the more severe phenotype with juvenile-onset (J)SLE. In this study, analysis of the phenotype and function of NK cells in a large cohort of JSLE patients demonstrated that total NK cells, as well as perforin and granzyme A expressing NK cell populations, were significantly diminished in JSLE patients compared to age- and sex-matched healthy controls. The reduction in NK cell frequency was associated with increased disease activity, and transcriptomic analysis of NK populations from active and low disease activity JSLE patients versus healthy controls confirmed that disease activity was the main driver of differential NK cell gene expression. Pathway analysis of differentially expressed genes revealed an upregulation of interferon-α responses and a downregulation of exocytosis in active disease compared to healthy controls. Further gene set enrichment analysis also demonstrated an overrepresentation of the apoptosis pathway in active disease. This points to increased propensity for apoptosis as a potential factor contributing to NK cell deficiency in JSLE.

Ragas: integration and enhanced visualization for single cell subcluster analysis

SUMMARY

Subcluster analysis is a powerful means to improve clustering and characterization of single cell RNA-Seq data. However, there are no existing tools to systematically integrate results from multiple subclusters, which creates hurdles for accurate data quantification, visualization, and interpretation in downstream analysis. To address this issue, we developed Ragas, an R package that integrates multi-level subclustering objects for streamlined analysis and visualization. A new data structure was implemented to seamlessly connect and assemble miscellaneous single cell analyses from different levels of subclustering, along with several new or enhanced visualization functions. Moreover, a re-projection algorithm was developed to integrate nearest-neighbor graphs from multiple subclusters in order to maximize their separability on the combined cell embeddings, which significantly improved the presentation of rare and homogeneous subpopulations.

AVAILABILITY AND IMPLEMENTATION

The Ragas package and its documentation can be accessed through https://github.com/jig4003/Ragas and its source code is also available at https://zenodo.org/records/11244921.

Uncoupling interferons and the interferon signature explains clinical and transcriptional subsets in SLE

Systemic lupus erythematosus (SLE) displays a hallmark interferon (IFN) signature. Yet, clinical trials targeting type I IFN (IFN-I) have shown variable efficacy, and blocking IFN-II failed to treat SLE. Here, we show that IFN type levels in SLE vary significantly across clinical and transcriptional endotypes. Whereas skin involvement correlated with IFN-I alone, systemic features like nephritis associated with co-elevation of IFN-I, IFN-II, and IFN-III, indicating additive IFN effects in severe SLE. Notably, while high IFN-II/-III levels without IFN-I had a limited effect on disease activity, IFN-II was linked to IFN-I-independent transcriptional profiles (e.g., OXPHOS and CD8+GZMH+ cells), and IFN-III enhanced IFN-induced gene expression when co-elevated with IFN-I. Moreover, dysregulated IFNs do not explain the IFN signature in 64% of patients or clinical manifestations including cytopenia, serositis, and anti-phospholipid syndrome, implying IFN-independent endotypes in SLE. This study sheds light on mechanisms underlying SLE heterogeneity and the variable response to IFN-targeted therapies in clinical trials.

Type I interferons induce an epigenetically distinct memory B cell subset in chronic viral infection

Memory B cells (MBCs) are key providers of long-lived immunity against infectious disease, yet in chronic viral infection, they do not produce effective protection. How chronic viral infection disrupts MBC development and whether such changes are reversible remain unknown. Through single-cell (sc)ATAC-seq and scRNA-seq during acute versus chronic lymphocytic choriomeningitis viral infection, we identified a memory subset enriched for interferon (IFN)-stimulated genes (ISGs) during chronic infection that was distinct from the T-bet+ subset normally associated with chronic infection. Blockade of IFNAR-1 early in infection transformed the chromatin landscape of chronic MBCs, decreasing accessibility at ISG-inducing transcription factor binding motifs and inducing phenotypic changes in the dominating MBC subset, with a decrease in the ISG subset and an increase in CD11c+CD80+ cells. However, timing was critical, with MBCs resistant to intervention at 4 weeks post-infection. Together, our research identifies a key mechanism to instruct MBC identity during viral infection.

Unraveling transcriptomic signatures and dysregulated pathways in systemic lupus erythematosus across disease states

OBJECTIVES

This study aims to elucidate the transcriptomic signatures and dysregulated pathways in patients with Systemic Lupus Erythematosus (SLE), with a particular focus on those persisting during disease remission.

METHODS

We conducted bulk RNA-sequencing of peripheral blood mononuclear cells (PBMCs) from a well-defined cohort comprising 26 remission patients meeting the Low Lupus Disease Activity State (LLDAS) criteria, 76 patients experiencing disease flares, and 15 healthy controls. To elucidate immune signature changes associated with varying disease states, we performed extensive analyses, including the identification of differentially expressed genes and pathways, as well as the construction of protein-protein interaction networks.

RESULTS

Several transcriptomic features recovered during remission compared to the active disease state, including down-regulation of plasma and cell cycle signatures, as well as up-regulation of lymphocytes. However, specific innate immune response signatures, such as the interferon (IFN) signature, and gene modules involved in chromatin structure modification, persisted across different disease states. Drug repurposing analysis revealed certain drug classes that can target these persistent signatures, potentially preventing disease relapse.

CONCLUSION

Our comprehensive transcriptomic study revealed gene expression signatures for SLE in both active and remission states. The discovery of gene expression modules persisting in the remission stage may shed light on the underlying mechanisms of vulnerability to relapse in these patients, providing valuable insights for their treatment.

Single-cell transcriptomic analysis of hematopoietic progenitor cells from patients with systemic lupus erythematosus reveals interferon-inducible reprogramming in early progenitors

Introduction

Immune cells that contribute to the pathogenesis of systemic lupus erythematosus (SLE) derive from adult hematopoietic stem and progenitor cells (HSPCs) within the bone marrow (BM). For this reason, we reasoned that fundamental abnormalities in SLE can be traced to a BM-derived HSPC inflammatory signature.

Methods

BM samples from four SLE patients, six healthy controls, and two umbilical cord blood (CB) samples were used. CD34+ cells were isolated from BM and CB samples, and single-cell RNA-sequencing was performed.

Results

A total of 426 cells and 24,473 genes were used in the analysis. Clustering analysis resulted in seven distinct clusters of cell types. Mutually exclusive markers, which were characteristic of each cell type, were identified. We identified three HSPC subpopulations, one of which consisted of proliferating cells (MKI67 expressing cells), one T-like, one B-like, and two myeloid-like progenitor subpopulations. Differential expression analysis revealed i) cell cycle-associated signatures, in healthy BM of HSPC clusters 3 and 4 when compared with CB, and ii) interferon (IFN) signatures in SLE BM of HSPC clusters 3 and 4 and myeloid-like progenitor cluster 5 when compared with healthy controls. The IFN signature in SLE appeared to be deregulated following TF regulatory network analysis and differential alternative splicing analysis between SLE and healthy controls in HSPC subpopulations.

Discussion

This study revealed both quantitative-as evidenced by decreased numbers of non-proliferating early progenitors-and qualitative differences-characterized by an IFN signature in SLE, which is known to drive loss of function and depletion of HSPCs. Chronic IFN exposure affects early hematopoietic progenitors in SLE, which may account for the immune aberrancies and the cytopenias in SLE.

Different Immunologic Profiles Are Associated With Distinct Clinical Phenotypes in Longitudinally Observed Patients With Systemic Lupus Erythematosus

OBJECTIVE

The aim of this study was to determine the immunologic profile associated with disease flares in patients with systemic lupus erythematosus (SLE) and to investigate the clinical significance of any differences observed between patients during and following a flare.

METHODS

Multiparameter flow cytometry was used to examine 47 immune populations within the peripheral blood of 16 healthy controls, 25 patients with clinically quiescent SLE, and 46 patients with SLE experiencing a flare at baseline and at 6- and 12-month follow-up visits. Unsupervised clustering was used to identify patients with similar immune profiles and to track changes over time. Parametric or nonparametric statistics were used when appropriate to assess the association of cellular phenotypes with clinical and laboratory parameters.

RESULTS

Five clusters of patients were identified that variably contained patients with active and quiescent SLE, and that had distinct clinical phenotypes. Patients characterized by increased T peripheral helper, activated B, and age-associated B cells were the most likely to be flaring at baseline, as well as the most likely to remain active or flare over the subsequent year if they acquired or retained this phenotype at follow-up. In contrast, patients who had increased T helper (Th) cells in the absence of B cell changes, or who had increased Th1 cells and innate immune populations, mostly developed quiescent SLE on follow-up. A significant proportion of patients with SLE had depletion of many immune populations at flare and only showed increases in these populations post-flare.

CONCLUSION

Cellular phenotyping of patients with SLE reveals several distinct immunologic profiles that may help to stratify patients with regard to prognosis and treatment.

Immunological Characteristics of Hepatic Dendritic Cells in Patients and Mouse Model with Liver Echinococcus multilocularis Infection

The cestode Echinococcus multilocularis, which mainly dwells in the liver, leads to a serious parasitic liver disease called alveolar echinococcosis (AE). Despite the increased attention drawn to the immunosuppressive microenvironment formed by hepatic AE tissue, the immunological characteristics of hepatic dendritic cells (DCs) in the AE liver microenvironment have not been fully elucidated. Here, we profiled the immunophenotypic characteristics of hepatic DC subsets in both clinical AE patients and a mouse model. Single-cell RNA sequencing (scRNA-Seq) analysis of four AE patient specimens revealed that greater DC numbers were present within perilesional liver tissues and that the distributions of cDC and pDC subsets in the liver and periphery were different. cDCs highly expressed the costimulatory molecule CD86, the immune checkpoint molecule CD244, LAG3, CTLA4, and the checkpoint ligand CD48, while pDCs expressed these genes at low frequencies. Flow cytometric analysis of hepatic DC subsets in an E. multilocularis infection mouse model demonstrated that the number of cDCs significantly increased after parasite infection, and a tolerogenic phenotype characterized by a decrease in CD40 and CD80 expression levels was observed at an early stage, whereas an activated phenotype characterized by an increase in CD86 expression levels was observed at a late stage. Moreover, the expression profiles of major immune checkpoint molecules (CD244 and LAG3) and ligands (CD48) on hepatic DC subsets in a mouse model exhibited the same pattern as those in AE patients. Notably, the cDC and pDC subsets in the E. multilocularis infection group exhibited higher expression levels of PD-L1 and CD155 than those in the control group, suggesting the potential of these subsets to impair T cell function. These findings may provide valuable information for investigating the role of hepatic DC subsets in the AE microenvironment and guiding DC targeting treatments for AE.

Skin in the game: a review of single-cell and spatial transcriptomics in dermatological research

Single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics (ST) are two emerging research technologies that uniquely characterize gene expression microenvironments on a cellular or subcellular level. The skin, a clinically accessible tissue composed of diverse, essential cell populations, serves as an ideal target for these high-resolution investigative approaches. Using these tools, researchers are assembling a compendium of data and discoveries in healthy skin as well as a range of dermatologic pathophysiologies, including atopic dermatitis, psoriasis, and cutaneous malignancies. The ongoing advancement of single-cell approaches, coupled with anticipated decreases in cost with increased adoption, will reshape dermatologic research, profoundly influencing disease characterization, prognosis, and ultimately clinical practice.

B-cell immune dysregulation with low soluble CD22 levels in refractory seronegative myasthenia gravis

Myasthenia gravis (MG), primarily caused by acetylcholine receptor (AChR) autoantibodies, is a chronic autoimmune disorder causing severe muscle weakness and fatigability. In particular, seronegative MG constitutes 10%-15% of MG cases and presents diagnostic challenges especially in early-onset female patients who often show severe disease and resistance to immunosuppressive therapy. Furthermore, the immunopathology of seronegative MG remains unclear. Thus, in this study, we aimed to elucidate the pathogenic mechanism of seronegative MG using scRNA-seq analysis and plasma proteome analysis; in particular, we investigated the relationship between immune dysregulation status and disease severity in refractory seronegative MG. Employing single-cell RNA-sequencing and plasma proteome analyses, we analyzed peripheral blood samples from 30 women divided into three groups: 10 healthy controls, 10 early-onset AChR-positive MG, and 10 refractory early-onset seronegative MG patients, both before and after intravenous immunoglobulin treatment. The disease severity was evaluated using the MG-Activities of Daily Living (ADL), MG composite (MGC), and revised 15-item MG-Quality of Life (QOL) scales. We observed numerical abnormalities in multiple immune cells, particularly B cells, in patients with refractory seronegative MG, correlating with disease activity. Notably, severe MG cases had fewer regulatory T cells without functional abnormalities. Memory B cells were found to be enriched in peripheral blood cells compared with naïve B cells. Moreover, plasma proteome analysis indicated significantly lower plasma protein levels of soluble CD22, expressed in the lineage of B-cell maturation (including mature B cells and memory B cells), in refractory seronegative MG patients than in healthy donors or patients with AChR-positive MG. Soluble CD22 levels were correlated with disease severity, B-cell frequency, and RNA expression levels of CD22. In summary, this study elucidates the immunopathology of refractory seronegative MG, highlighting immune disorders centered on B cells and diminished soluble CD22 levels. These insights pave the way for novel MG treatment strategies focused on B-cell biology.

T cell expressions of aberrant gene signatures and Co-inhibitory receptors (Co-IRs) as predictors of renal damage and lupus disease activity

BACKGROUND

Systemic lupus erythematosus (SLE) is distinguished by an extensive range of clinical heterogeneity with unpredictable disease flares and organ damage. This research investigates the potential of aberrant signatures on T cell genes, soluble Co-IRs/ligands, and Co-IRs expression on T cells as biomarkers for lupus disease parameters.

METHODS

Comparative transcriptome profiling analysis of non-renal and end-stage renal disease (ESRD) phenotypes of SLE was performed using CD4 + and CD8 + cDNA microarrays of sorted T cells. Comparing the expression of Co-IRs on T cells and serum soluble mediators among healthy and SLE phenotypes.

RESULTS

SLE patients with ESRD were downregulated CD38, PLEK, interferon-γ, CX3CR1, FGFBP2, and SLCO4C1 transcripts on CD4 + and CD8 + T cells simultaneously and NKG7, FCRL6, GZMB/H, FcγRIII, ITGAM, Fas ligand, TBX21, LYN, granulysin, CCL4L1, CMKLR1, HLA-DRβ, KIR2DL3, and KLRD1 in CD8 T cells. Pathway enrichment and PPI network analyses revealed that the overwhelming majority of Differentially Expressed Genes (DEGs) have been affiliated with novel cytotoxic, antigen presentation, and chemokine-cell migration signature pathways. CD8 + GZMK + T cells that are varied in nature, including CD161 + Mucosal-associated invariant T (MAIT) cells and CD161- aged-associated T (Taa) cells and CD161-GZMK + GZMB + T cells might account for a higher level of GZMK in CD8 + T cells associated with ESRD. SLE patients have higher TIGIT + , PD1 + , and lower CD127 + cell percentages on CD4 + T cells, higher TIM3 + , TIGIT + , HLA-DR + cell frequency, and lower MFI expression of CD127, CD160 in CD8 T cells. Co-IRs expression in T cells was correlated with soluble PD-1, PDL-2, and TIM3 levels, as well as SLE disease activity, clinical phenotypes, and immune-therapy responses.

CONCLUSION

The signature of dysfunctional pathways defines a distinct immunity pattern in LN ESRD patients. Expression levels of Co-IRs in peripheral blood T cells and serum levels of soluble PD1/PDL-2/TIM3 can serve as biomarkers for evaluating clinical parameters and therapeutic responses.

Single-Cell Profiling of Bone Marrow B Cells and Early B Cell Developmental Disorders Associated With Systemic Lupus Erythematosus

OBJECTIVE

The peripheral B cell compartment is heavily disturbed in systemic lupus erythematosus (SLE), but whether B cells develop aberrantly in the bone marrow (BM) is largely unknown.

METHODS

We performed single-cell RNA/B cell receptor (BCR) sequencing and immune profiling of BM B cells and classified patients with SLE into two groups: early B cell (Pro-B and Pre-B) normal (EBnor) and EB defective/low (EBlo) groups.

RESULTS

The SLE-EBlo group exhibited more severe disease activity and proinflammatory status, overaction of type I interferon signaling and metabolic pathways within the B cell compartment, and aberrant BCR repertoires compared with the SLE-EBnor group. Moreover, in one patient with SLE who was initially classified in the SLE-EBlo group, early B cell deficiency and associated abnormalities were largely rectified in a second BM sample at the remission phase.

CONCLUSION

In summary, this study suggests that early B cell loss in BM defines a unique pathological state in a subset of patients with SLE that may play an active role in the dysregulated autoimmune responses.

Vascular damage in systemic lupus erythematosus

Vascular disease is a major cause of morbidity and mortality in patients with systemic autoimmune diseases, particularly systemic lupus erythematosus (SLE). Although comorbid cardiovascular risk factors are frequently present in patients with SLE, they do not explain the high burden of premature vascular disease. Profound innate and adaptive immune dysregulation seems to be the primary driver of accelerated vascular damage in SLE. In particular, evidence suggests that dysregulation of type 1 interferon (IFN-I) and aberrant neutrophils have key roles in the pathogenesis of vascular damage. IFN-I promotes endothelial dysfunction directly via effects on endothelial cells and indirectly via priming of immune cells that contribute to vascular damage. SLE neutrophils are vasculopathic in part because of their increased ability to form immunostimulatory neutrophil extracellular traps. Despite improvements in clinical care, cardiovascular disease remains the leading cause of mortality among patients with SLE, and treatments that improve vascular outcomes are urgently needed. Improved understanding of the mechanisms of vascular injury in inflammatory conditions such as SLE could also have implications for common cardiovascular diseases, such as atherosclerosis and hypertension, and may ultimately lead to personalized therapeutic approaches to the prevention and treatment of this potentially fatal complication.

Age-associated CD4+ T cells with B cell-promoting functions are regulated by ZEB2 in autoimmunity

Aging is a significant risk factor for autoimmunity, and many autoimmune diseases tend to onset during adulthood. We conducted an extensive analysis of CD4+ T cell subsets from 354 patients with autoimmune disease and healthy controls via flow cytometry and bulk RNA sequencing. As a result, we identified a distinct CXCR3midCD4+ effector memory T cell subset that expands with age, which we designated "age-associated T helper (THA) cells." THA cells exhibited both a cytotoxic phenotype and B cell helper functions, and these features were regulated by the transcription factor ZEB2. Consistent with the highly skewed T cell receptor usage of THA cells, gene expression in THA cells from patients with systemic lupus erythematosus reflected disease activity and was affected by treatment with a calcineurin inhibitor. Moreover, analysis of single-cell RNA sequencing data revealed that THA cells infiltrate damaged organs in patients with autoimmune diseases. Together, our characterization of THA cells may facilitate improved understanding of the relationship between aging and autoimmune diseases.

Enhanced GATA4 expression in senescent systemic lupus erythematosus monocytes promotes high levels of IFNα production

Enhanced interferon α (IFNα) production has been implicated in the pathogenesis of systemic lupus erythematosus (SLE). We previously reported IFNα production by monocytes upon activation of the stimulator of IFN genes (STING) pathway was enhanced in patients with SLE. We investigated the mechanism of enhanced IFNα production in SLE monocytes. Monocytes enriched from the peripheral blood of SLE patients and healthy controls (HC) were stimulated with 2'3'-cyclic GAMP (2'3'-cGAMP), a ligand of STING. IFNα positive/negative cells were FACS-sorted for RNA-sequencing analysis. Gene expression in untreated and 2'3'-cGAMP-stimulated SLE and HC monocytes was quantified by real-time PCR. The effect of GATA binding protein 4 (GATA4) on IFNα production was investigated by overexpressing GATA4 in monocytic U937 cells by vector transfection. Chromatin immunoprecipitation was performed to identify GATA4 binding target genes in U937 cells stimulated with 2'3'-cGAMP. Differentially expressed gene analysis of cGAS-STING stimulated SLE and HC monocytes revealed the enrichment of gene sets related to cellular senescence in SLE. CDKN2A, a marker gene of cellular senescence, was upregulated in SLE monocytes at steady state, and its expression was further enhanced upon STING stimulation. GATA4 expression was upregulated in IFNα-positive SLE monocytes. Overexpression of GATA4 enhanced IFNα production in U937 cells. GATA4 bound to the enhancer region of IFIT family genes and promoted the expressions of IFIT1, IFIT2, and IFIT3, which promote type I IFN induction. SLE monocytes with accelerated cellular senescence produced high levels of IFNα related to GATA4 expression upon activation of the cGAS-STING pathway.

Imbalance of programmed cell death patterns mediated by dendritic cell subsets in systemic lupus erythematosus and lupus nephritis

OBJECTIVES

Abnormal programmed cell death in immune cells is associated with autoimmune diseases, but the patterns of programmed cell death in systemic lupus erythematosus (SLE) and especially lupus nephritis (LN) remain unclear. This study aims to explore the association between SLE, LN, and immune cell death patterns.

METHODS

Bulk RNA sequencing (bulk RNA-seq) and single-cell RNA sequencing (scRNA-seq) data were downloaded from the Gene Expression Omnibus (GEO) database. Bioinformatic analysis was conducted to explore the expression levels of genes related to 3 cell death patterns in peripheral blood mononuclear cells of SLE patients. Key cell subsets involved in the imbalance of cell death patterns were identified through scRNA-seq. Immunofluorescence was used to detect the expression levels of receptor interacting serine/threonine kinase 3 (RIPK3), mixed-lineage kinase domain-like protein (MLKL), phosphorylated MLKL (pMLKL), caspase 1 (CASP1), CD1c molecule (CD1C), C-type lectin domain containing 9A (CLEC9A), and X-C motif chemokine receptor 1 (XCR1) in dendritic cells (DC). scRNA-seq was performed on kidney tissues collected from LN patients and healthy controls (HC) at the Third Xiangya Hospital of Central South University, followed by bioinformatic analysis to identify key cell subsets involved in the imbalance of cell death patterns. Pseudotime analysis and ligand-receptor analysis were used to explore the differentiation direction and cell communication of different DC subsets. Transient transfection was used to transfect RAW264.7 cells with empty plasmid, empty plasmid+dsDNA (HSV-DNA), empty plasmid+200 μmol/L tert-butyl hydroperoxide (TBHP), stimulator of interferon genes (STING) shRNA plasmid, STING shRNA plasmid+dsDNA (HSV-DNA), and STING shRNA plasmid+200 μmol/L TBHP. Annexin V-mCherry and SYTOX Green staining were used to detect cell death in each group. Western blotting was used to detect the activation of CASP1, gasdermin D (GSDMD), RIPK3, and MLKL in each group.

RESULTS

Bioinformatic analysis showed an imbalance in 3 cell death patterns in SLE and LN patients: Pro-inflammatory pyroptosis and necroptosis were activated, while anti-inflammatory apoptosis was inhibited. The key cell subsets involved were DC subsets, particularly focusing on CLEC9A+cDC1. Immunofluorescence results showed that the expression levels of RIPK3, MLKL, and CASP1 in DCs were higher in the SLE group compared to the HC group. pMLKL and CASP1 expression levels in renal cDC1 marked by CLEC9A and XCR1 were higher in the LN group than in the HC group. Pseudotime analysis and ligand-receptor analysis suggested that the CLEC9A+cDC1 subset in LN kidney tissues originated from peripheral circulation. Annexin V-mCherry and SYTOX Green staining results showed that the number of dead cells decreased in the STING shRNA transfection group compared to the empty plasmid group in RAW264.7 cells. Western blotting results showed that the activation of CASP1, GSDMD, RIPK3, and MLKL was decreased in the STING shRNA transfection group compared to the empty plasmid group.

CONCLUSIONS

This study provides novel insights into the role of CLEC9A+cDC1 in the imbalance of cell death patterns in SLE and LN.

CD169+ classical monocyte as an important participant in Graves' ophthalmopathy through CXCL12-CXCR4 axis

Patients with Graves' disease (GD) can develop Graves' ophthalmopathy (GO), but the underlying pathological mechanisms driving this development remain unclear. In our study, which included patients with GD and GO, we utilized single-cell RNA sequencing (scRNA-seq) and multiplatform analyses to investigate CD169+ classical monocytes, which secrete proinflammatory cytokines and are expanded through activated interferon signaling. We found that CD169+ clas_mono was clinically significant in predicting GO progression and prognosis, and differentiated into CD169+ macrophages that promote inflammation, adipogenesis, and fibrosis. Our murine model of early-stage GO showed that CD169+ classical monocytes accumulated in orbital tissue via the Cxcl12-Cxcr4 axis. Further studies are needed to investigate whether targeting circulating monocytes and the Cxcl12-Cxcr4 axis could alleviate GO progression.

Blood immunophenotyping identifies distinct kidney histopathology and outcomes in patients with lupus nephritis

Lupus nephritis (LN) is a frequent manifestation of systemic lupus erythematosus, and fewer than half of patients achieve complete renal response with standard immunosuppressants. Identifying non-invasive, blood-based pathologic immune alterations associated with renal injury could aid therapeutic decisions. Here, we used mass cytometry immunophenotyping of peripheral blood mononuclear cells in 145 patients with biopsy-proven LN and 40 healthy controls to evaluate the heterogeneity of immune activation in patients with LN and to identify correlates of renal parameters and treatment response. Unbiased analysis identified 3 immunologically distinct groups of patients with LN that were associated with different patterns of histopathology, renal cell infiltrates, urine proteomic profiles, and treatment response at one year. Patients with enriched circulating granzyme B+ T cells at baseline showed more severe disease and increased numbers of activated CD8 T cells in the kidney, yet they had the highest likelihood of treatment response. A second group characterized primarily by a high type I interferon signature had a lower likelihood of response to therapy, while a third group appeared immunologically inactive by immunophenotyping at enrollment but with chronic renal injuries. Main immune profiles could be distilled down to 5 simple cytometric parameters that recapitulate several of the associations, highlighting the potential for blood immune profiling to translate to clinically useful non-invasive metrics to assess immune-mediated disease in LN.

UNC93B1 variants underlie TLR7-dependent autoimmunity

UNC93B1 is critical for trafficking and function of nucleic acid-sensing Toll-like receptors (TLRs) TLR3, TLR7, TLR8, and TLR9, which are essential for antiviral immunity. Overactive TLR7 signaling induced by recognition of self-nucleic acids has been implicated in systemic lupus erythematosus (SLE). Here, we report UNC93B1 variants (E92G and R336L) in four patients with early-onset SLE. Patient cells or mouse macrophages carrying the UNC93B1 variants produced high amounts of TNF-α and IL-6 and upon stimulation with TLR7/TLR8 agonist, but not with TLR3 or TLR9 agonists. E92G causes UNC93B1 protein instability and reduced interaction with TLR7, leading to selective TLR7 hyperactivation with constitutive type I IFN signaling. Thus, UNC93B1 regulates TLR subtype-specific mechanisms of ligand recognition. Our findings establish a pivotal role for UNC93B1 in TLR7-dependent autoimmunity and highlight the therapeutic potential of targeting TLR7 in SLE.

Comprehensive analysis of lactate-related gene profiles and immune characteristics in lupus nephritis

Objectives

The most frequent cause of kidney damage in systemic lupus erythematosus (SLE) is lupus nephritis (LN), which is also a significant risk factor for morbidity and mortality. Lactate metabolism and protein lactylation might be related to the development of LN. However, there is still a lack of relative research to prove the hypothesis. Hence, this study was conducted to screen the lactate-related biomarkers for LN and analyze the underlying mechanism.

Methods

To identify differentially expressed genes (DEGs) in the training set (GSE32591, GSE127797), we conducted a differential expression analysis (LN samples versus normal samples). Then, module genes were mined using WGCNA concerning LN. The overlapping of DEGs, critical module genes, and lactate-related genes (LRGs) was used to create the lactate-related differentially expressed genes (LR-DEGs). By using a machine-learning algorithm, ROC, and expression levels, biomarkers were discovered. We also carried out an immune infiltration study based on biomarkers and GSEA.

Results

A sum of 1259 DEGs was obtained between LN and normal groups. Then, 3800 module genes in reference to LN were procured. 19 LR-DEGs were screened out by the intersection of DEGs, key module genes, and LRGs. Moreover, 8 pivotal genes were acquired via two machine-learning algorithms. Subsequently, 3 biomarkers related to lactate metabolism were obtained, including COQ2, COQ4, and NDUFV1. And these three biomarkers were enriched in pathways 'antigen processing and presentation' and 'NOD-like receptor signaling pathway'. We found that Macrophages M0 and T cells regulatory (Tregs) were associated with these three biomarkers as well.

Conclusion

Overall, the results indicated that lactate-related biomarkers COQ2, COQ4, and NDUFV1 were associated with LN, which laid a theoretical foundation for the diagnosis and treatment of LN.

Type I IFN in Glomerular Disease: Scarring beyond the STING

The field of nephrology has recently directed a considerable amount of attention towards the stimulator of interferon genes (STING) molecule since it appears to be a potent driver of chronic kidney disease (CKD). STING and its activator, the cyclic GMP-AMP synthase (cGAS), along with intracellular RIG-like receptors (RLRs) and toll-like receptors (TLRs), are potent inducers of type I interferon (IFN-I) expression. These cytokines have been long recognized as part of the mechanism used by the innate immune system to battle viral infections; however, their involvement in sterile inflammation remains unclear. Mounting evidence pointing to the involvement of the IFN-I pathway in sterile kidney inflammation provides potential insights into the complex interplay between the innate immune system and damage to the most sensitive segment of the nephron, the glomerulus. The STING pathway is often cited as one cause of renal disease not attributed to viral infections. Instead, this pathway can recognize and signal in response to host-derived nucleic acids, which are also recognized by RLRs and TLRs. It is still unclear, however, whether the development of renal diseases depends on subsequent IFN-I induction or other processes involved. This review aims to explore the main endogenous inducers of IFN-I in glomerular cells, to discuss what effects autocrine and paracrine signaling have on IFN-I induction, and to identify the pathways that are implicated in the development of glomerular damage.

Single-cell resolution of longitudinal blood transcriptome profiles in rheumatoid arthritis, systemic lupus erythematosus and healthy control pregnancies

OBJECTIVES

Comparative longitudinal analyses of cellular composition and peripheral blood gene expression in Rheumatoid arthritis (RA), systemic lupus erythematosus (SLE) and healthy pregnancies.

METHODS

In total, 335 whole blood samples from 84 RA, SLE and healthy controls before pregnancy, at each trimester, 6 weeks, 6 months and 12 months post partum were analysed. We combined bulk and single cell RNA analyses for cell-type estimation, validated by flow cytometry, before combining this in a cell-type adjusted analysis for an improved resolution of unrecognised gene expression changes associated with RA and SLE pregnancies.

RESULTS

Patients were well regulated throughout pregnancy, and few had pregnancy complications. In SLE, the interferon signature was augmented during pregnancy, and the pregnancy signature was continued post partum. An altered cell type composition strongly influences the profile. In the pregnancy signature, transcripts involved in galactosylation potentially altering the effector functions of autoantibodies became more evident. Several genes in the adjusted RA signature are expressed in mucosal associated invariant T cells.

CONCLUSION

We found distinct RA, SLE and pregnancy signatures, and no expression patterns could be attributed to medication or disease activity. Our results support the need for close postpartum follow-up of patients with SLE. Gene expression patterns in RA were closer to healthy controls than to SLE, and primarily became evident after cell-type adjustment. Adjusting for cell abundance unravelled gene expression signatures less associated with variation in cell-composition and highlighted genes with expression profiles associated with changes in specialised cell populations.

Single-cell transcriptome landscape of circulating CD4+ T cell populations in autoimmune diseases

CD4+ T cells are key mediators of various autoimmune diseases; however, their role in disease progression remains unclear due to cellular heterogeneity. Here, we evaluated CD4+ T cell subpopulations using decomposition-based transcriptome characterization and canonical clustering strategies. This approach identified 12 independent gene programs governing whole CD4+ T cell heterogeneity, which can explain the ambiguity of canonical clustering. In addition, we performed a meta-analysis using public single-cell datasets of over 1.8 million peripheral CD4+ T cells from 953 individuals by projecting cells onto the reference and cataloging cell frequency and qualitative alterations of the populations in 20 diseases. The analyses revealed that the 12 transcriptional programs were useful in characterizing each autoimmune disease and predicting its clinical status. Moreover, genetic variants associated with autoimmune diseases showed disease-specific enrichment within the 12 gene programs. The results collectively provide a landscape of single-cell transcriptomes of CD4+ T cell subpopulations involved in autoimmune disease.

ProBDNF contributed to patrolling monocyte infiltration and renal damage in systemic lupus erythematosus

Monocyte aberrations have been increasingly recognized as contributors to renal damage in systemic lupus erythematosus (SLE), however, recognition of the underlying mechanisms and modulating strategies is at an early stage. Our studies have demonstrated that brain-derived neurotrophic factor precursor (proBDNF) drives the progress of SLE by perturbing antibody-secreting B cells, and proBDNF facilitates pro-inflammatory responses in monocytes. By utilizing peripheral blood from patients with SLE, GEO database and spontaneous MRL/lpr lupus mice, we demonstrated in the present study that CX3CR1+ patrolling monocytes (PMo) numbers were decreased in SLE. ProBDNF was specifically expressed in CX3CR1+ PMo and was closely correlated with disease activity and the degree of renal injury in SLE patients. In MRL/lpr mice, elevated proBDNF was found in circulating PMo and the kidney, and blockade of proBDNF restored the balance of circulating and kidney-infiltrating PMo. This blockade also led to the reversal of pro-inflammatory responses in monocytes and a noticeable improvement in renal damage in lupus mice. Overall, the results indicate that the upregulation of proBDNF in PMo plays a crucial role in their infiltration into the kidney, thereby contributing to nephritis in SLE. Targeting of proBDNF offers a potential therapeutic role in modulating monocyte-driven renal damage in SLE.

Machine learning application in autoimmune diseases: State of art and future prospectives

Autoimmune diseases are a group of disorders resulting from an alteration of immune tolerance, characterized by the formation of autoantibodies and the consequent development of heterogeneous clinical manifestations. Diagnosing autoimmune diseases is often complicated, and the available prognostic tools are limited. Machine learning allows us to analyze large amounts of data and carry out complex calculations quickly and with minimal effort. In this work, we examine the literature focusing on the use of machine learning in the field of the main systemic (systemic lupus erythematosus and rheumatoid arthritis) and organ-specific autoimmune diseases (type 1 diabetes mellitus, autoimmune thyroid, gastrointestinal, and skin diseases). From our analysis, interesting applications of machine learning emerged for developing algorithms useful in the early diagnosis of disease or prognostic models (risk of complications, therapeutic response). Subsequent studies and the creation of increasingly rich databases to be supplied to the algorithms will eventually guide the clinician in the diagnosis, allowing intervention when the pathology is still in an early stage and immediately directing towards a correct therapeutic approach.

Proximity extension assay proteomics and renal single cell transcriptomics uncover novel urinary biomarkers for active lupus nephritis

OBJECTIVE

To identify urinary biomarkers that can distinguish active renal involvement in Lupus Nephritis (LN), a severe manifestation of systemic lupus erythematosus (SLE).

METHODS

Urine from 117 subjects, comprised of inactive SLE, active non-renal lupus, active LN, and healthy controls, were subjected to Proximity Extension Assay (PEA) based comprehensive proteomics followed by ELISA validation in an independent, ethnically diverse cohort. Proteomic data is also cross-referenced to renal transcriptomic data to elucidate cellular origins of biomarkers.

RESULTS

Systems biology analyses revealed progressive activation of cytokine signaling, chemokine activity and coagulation pathways, with worsening renal disease. In addition to validating 30 previously reported biomarkers, this study uncovers several novel candidates. Following ELISA validation in an independent cohort of different ethnicity, the six most discriminatory biomarkers for active LN were urinary ICAM-2, FABP4, FASLG, IGFBP-2, SELE and TNFSF13B/BAFF, with ROC AUC ≥80%, with most correlating strongly with clinical disease activity. Transcriptomic analyses of LN kidneys mapped the likely origin of these proteins to intra-renal myeloid cells (CXCL16, IL-1RT2, TNFSF13B/BAFF), T/NK cells (FASLG), leukocytes (ICAM2) and endothelial cells (SELE).

CONCLUSION

In addition to confirming the diagnostic potential of urine ALCAM, CD163, MCP1, SELL, ICAM1, VCAM1, NGAL and TWEAK for active LN, this study adds urine ICAM-2, FABP4, FASLG, IGFBP-2, SELE, and TNFSF13B/BAFF as additional markers that warrant systematic validation in larger cross-sectional and longitudinal cohorts.

Multimodal Single-Cell Sequencing of B Cells in Primary Sjögren's Syndrome

OBJECTIVE

B cells are important in the pathogenesis of primary Sjögren's syndrome (pSS). Patients positive for Sjögren's syndrome antigen A/Sjögren syndrome antigen B (SSA/SSB) autoantibodies are more prone to systemic disease manifestations and adverse outcomes. We aimed to determine the role of B cell composition, gene expression, and B cell receptor usage in pSS subgroups stratified for SSA/SSB antibodies.

METHODS

Over 230,000 B cells were isolated from peripheral blood of patients with pSS (n = 6 SSA-, n = 8 SSA+ single positive and n = 10 SSA/SSB+ double positive) and four healthy controls and processed for single-cell RNA sequencing (scRNA-seq) and single-cell variable, diversity, and joining (VDJ) gene sequencing (scVDJ-seq).

RESULTS

We show that SSA/SSB+ patients present the highest and lowest proportion of naïve and memory B cells, respectively, and the highest up-regulation of interferon-induced genes across all B cell subtypes. Differential usage of IGHV showed that IGHV1-69 and IGHV4-30-4 were more often used in all pSS subgroups compared with controls. Memory B cells from SSA/SSB+ patients displayed a higher proportion of cells with unmutated VDJ transcripts compared with other pSS patient groups and controls, indicating altered somatic hypermutation processes. Comparison with previous studies revealed heterogeneous clonotype pools, with little overlap in CDR3 sequences. Joint analysis using scRNA-seq and scVDJ-seq data allowed unsupervised stratification of patients with pSS and identified novel parameters that correlated to disease manifestations and antibody status.

CONCLUSION

We describe heterogeneity and molecular characteristics in B cells from patients with pSS, providing clues to intrinsic differences in B cells that affect the phenotype and outcome and allowing stratification of patients with pSS at improved resolution.

Adverse effects of gestational diabetes mellitus on fetal monocytes revealed by single-cell RNA sequencing

Gestational diabetes mellitus (GDM), the most prevalent metabolic disorder during pregnancy, has long-term risks of metabolic diseases in offspring. However, the underlying mechanisms remain unclear. Here, we analyzed single-cell transcriptional data of cord blood mononuclear cells (CBMCs) from fetuses of healthy and GDM mothers, peripheral blood mononuclear cells from children and adolescents, and coronary plaques myeloid cells from atherosclerosis. Our results demonstrated that monocytes in cord blood were characterized with down-regulated proinflammatory-related pathways and up-regulated proliferation-related pathways. And monocytes in cord blood from GDM mothers were featured with expanded CXCL8+IL1B+ subclusters, enhanced crosstalk with neutrophil granulocytes and augmented adhesive and phagocytic abilities. Interestingly, CXCL8+IL1B+ monocytes influenced by GDM had transcriptome similarity with those of coronary plaques myeloid cells from individuals with atherosclerotic cardiovascular disease. Collectively, our data reveal adverse impact of maternal GDM environment on fetal monocytes and propose potential mechanisms between maternal GDM and offspring atherosclerosis.

Immune and molecular landscape behind non-response to Mycophenolate Mofetil and Azathioprine in lupus nephritis therapy

Lupus nephritis (LN) represents one of the most severe complications of systemic lupus erythematosus, leading to end-stage kidney disease in worst cases. Current first-line therapies for LN, including mycophenolate mofetil (MMF) and azathioprine (AZA), fail to induce long-term remission in 60-70% of the patients, evidencing the urgent need to delve into the molecular knowledge-gap behind the non-response to these therapies. A longitudinal cohort of treated LN patients including clinical, cellular and transcriptomic data, was analyzed. Gene-expression signatures behind non-response to different drugs were revealed by differential expression analysis. Drug-specific non-response mechanisms and cell proportion differences were identified. Blood cell subsets mediating non-response were described using single-cell RNASeq data. We show that AZA and MMF non-response implicates different cells and regulatory functions. Mechanistic models were used to suggest add-on therapies to improve their current performance. Our results provide new insights into the molecular mechanisms associated with treatment failures in LN.

Effect of BAFF blockade on the B cell receptor repertoire and transcriptome in a mouse model of systemic lupus erythematosus

Introduction

Systemic lupus erythematosus (SLE) is a heterogeneous autoimmune disease. Anti-B-cell-activating factor (BAFF) therapy effectively depletes B cells and reduces SLE disease activity. This research aimed to evaluate the effect of BAFF blockade on B cell receptor (BCR) repertoire and gene expression.

Methods

Through next-generation sequencing, we analyzed gene expression and BCR repertoire in MRL/lpr mice that received long-term anti-BAFF therapy. Based on gene expression profiles, we predicted the relative proportion of immune cells using ImmuCellAI-mouse, validating our predictions via flow cytometry and FluoroSpot.

Results

The loss of BCR repertoire diversity and richness, along with increased clonality and differential frequency distribution of the immunoglobulin heavy chain variable (IGHV) segment gene usage, were observed in BAFF-blockade mice. Meanwhile, the distribution of complementarity-determining region 3 (CDR3) length and CDR3 amino acid usage remained unaffected. BAFF blockade resulted in extensive changes in gene expression, particularly that of genes related to B cells and immunoglobulins. Besides, the tumor necrosis factor (TNF)-α responses and interferon (IFN)-α/γ were downregulated, consistent with the decrease in IFN-γ and TNF-α serum levels following anti-BAFF therapy. In addition, BAFF blockade significantly reduced B cell subpopulations and plasmacytoid dendritic cells, and caused the depletion of antibody-secreting cells.

Discussion

Our comparative BCR repertoire and transcriptome analyses of MRL/lpr mice subjected to BAFF blockade provide innovative insights into the molecular pathophysiology of SLE.

Lupus-associated innate receptors drive extrafollicular evolution of autoreactive B cells

In systemic lupus erythematosus, recent findings highlight the extrafollicular (EF) pathway as prominent origin of autoantibody-secreting cells (ASCs). CD21loCD11c+ B cells, associated with aging, infection, and autoimmunity, are contributors to autoreactive EF ASCs but have an obscure developmental trajectory. To study EF kinetics of autoreactive B cell in tissue, we adoptively transferred WT and gene knockout B cell populations into the 564Igi mice - an autoreactive host enriched with autoantigens and T cell help. Time-stamped analyses revealed TLR7 dependence in early escape of peripheral B cell tolerance and establishment of a pre-ASC division program. We propose CD21lo cells as precursors to EF ASCs due to their elevated TLR7 sensitivity and proliferative nature. Blocking receptor function reversed CD21 loss and reduced effector cell generation, portraying CD21 as a differentiation initiator and a possible target for autoreactive B cell suppression. Repertoire analysis further delineated proto-autoreactive B cell selection and receptor evolution toward self-reactivity. This work elucidates receptor and clonal dynamics in EF development of autoreactive B cells, and establishes modular, native systems to probe mechanisms of autoreactivity.

Upregulation of CX3CR1 expression in circulating T cells of systemic lupus erythematosus patients as a reflection of autoimmune status through characterization of cytotoxic capacity

OBJECTIVE

This study investigated CX3CR1 expression in human peripheral blood T lymphocytes and their subsets, exploring changes in SLE patients and its diagnostic potential.

METHODS

Peripheral blood samples from 31 healthy controls and 50 SLE patients were collected. RNA-Seq data from SLE patient PBMCs were used to analyze CX3CR1 expression in T cells. Flow cytometry determined CX3CR1-expressing T lymphocyte subset proportions in SLE patients and healthy controls. Subset composition and presence of GZMB, GPR56, and perforin in CX3CR1+ T lymphocytes were analyzed. T cell-clinical indicator correlations were assessed. ROC curves explored CX3CR1's diagnostic potential for SLE.

RESULTS

CX3CR1+CD8+ T cells exhibited higher GPR56, perforin, and GZMB expression than other T cell subsets. The proportion of CX3CR1+ was higher in TEMRA and lower in Tn and TCM. PMA activation reduced CX3CR1+ T cell proportions. Both RNA-Seq and flow cytometry revealed elevated CX3CR1+ T cell proportions in SLE patients. Significantly lower perforin+ and GPR56+ proportions were observed in CX3CR1+CD8+ T cells in SLE patients. CX3CR1+ T cells correlated with clinical indicators.

CONCLUSION

CX3CR1+ T cells display cytotoxic features, with heightened expression in CD8+ T cells, particularly in adult SLE patients. Increased CX3CR1 expression in SLE patient T cells suggests its potential as an adjunctive diagnostic marker for SLE.

Single-cell transcriptional profiling reveals aberrant gene expression patterns and cell states in autoimmune diseases

Multiple sclerosis(MS), primary Sjögren syndrome (pSS), and systemic lupus erythematosus (SLE) share numerous clinical symptoms and serological characteristics. We analyzed 153550 cells of scRNA-seq data of 17 treatment-naive patients (5 MS, 5 pSS, and 7 SLE) and 10 healthy controls, and we examined the enrichment of biological processes, differentially expressed genes (DEGs), immune cell types, and their subpopulations, and cell-cell communication in peripheral blood mononuclear cells (PBMCs). The percentage of B cells, megakaryocytes, monocytes, and proliferating T cells presented significant changes in autoimmune diseases. The enrichment of cell types based on gene expression revealed an elevated monocyte. MIF, MK, and GALECTIN signaling networks were obvious differences in autoimmune diseases. Taken together, our analysis provides a comprehensive map of the cell types and states of ADs patients at the single-cell level to understand better the pathogenesis and treatment of these ADs.

A comprehensive exploration of the heterogeneity of immune cells in Han and Zang systemic lupus erythematosus patients via single-cell RNA sequencing

Systemic Lupus Erythematosus (SLE) is an autoimmune sickness with unclear pathogenesis. The goal of this research was to reveal the heterogeneity of immune cells in SLE patients of Han and Zang nationality by single-cell RNA sequencing (scRNA-seq) and bioinformatics profiling.

METHODS

A total of 94,102 peripheral blood mononuclear cells (PBMCs) from six volunteers with SLE (3 Zang, 3 Han) and six healthy controls were first conducted through scRNA-seq analysis. The immune cell subsets in the pathogenesis of SLE were analyzed as well. Real-time quantitative PCR (RT-qPCR) was applied to confirm the results of sc-RNA seq analysis.

RESULTS

For the Tibetan samples, the ratios of Naïve CD4 RPS4Y1 cells, Naïve CD4 cells, Memory BC CD24 and Memory BC differed significantly between the SLE and control samples, while that of CD8 CTL MAL cells was significantly different between the two groups in Han nationality samples. Variable differentiation states of CD8 CTL MAL cells, CD8 CTL GZMK cells, and Naïve CD4 cells were detected through pseudotime analysis. Moreover, T-cell receptor (TCR) abundance was notably higher in Tibetan SLE specimens than that in controls, while B-cell receptor (BCR) abundance in Tibetan and Han samples was higher than in control groups.

CONCLUSIONS

In summary, the immune cellular heterogeneity of SLE patients both Han and Zang nationality was explored based on various bioinformatics approaches, providing new perspectives for immunological characteristics of SLE among different ethnic groups.

Immunohistochemical expression of perforin in adult systemic lupus erythematosus associated macrophage activation syndrome: Clinicohematological correlation and literature review

OBJECTIVE

To study the bone marrow (BM) immunohistomorphological characteristics in adult systemic lupus erythematosus (SLE) associated macrophage activation syndrome (SLE-MAS).

MATERIALS AND METHODS

Immunohistochemical (IHC) expression of CD3, CD8, perforin (PFN), and CD163 was studied on BM trephine biopsies from 30 cytopenic adult SLE cases (male: female = 1:5, age; 24 years, range; 19-32) and compared them with ten age matched controls. Clinicopathological parameters were compared among the cases likely (L) or unlikely (U) to have MAS using probability scoring criteria. The best cut off laboratory parameters to discriminate between the two were obtained through receiver operator curve (ROC) analysis.

RESULTS

MAS occurred in 12/30 (40%) cases and was more commonly associated with prior immunosuppressive therapy (p = .07), ≥ 3 system involvement (p = .09), lower fibrinogen (p < .01), increased triglyceride (p = .002), increased BM hemophagocytosis (p = .002), and higher MAS score [185 (176-203) vs. 105 (77-119), p < .01] than MAS-U subgroup. Although PFN+CD8+ T lymphocytes significantly decreased among cases than controls (p < .05), it was comparable between MAS-L and MAS-U subgroups. Fibrinogen (< 2.4 g/L, AUC; 0.93, p < .01), hemophagocytosis score (> 1.5, AUC; 0.71, p = .03), and an MAS probability score of ≥ 164 (AUC; 1, p < .01) discriminated MAS from those without MAS.

CONCLUSION

We noted a decrease in perforin mediated CD8 + T cell cytotoxicity in SLE. Immunohistochemical demonstration of the same along with histiocytic hemophagocytosis on BM biopsy may be useful adjunct in early diagnosis and management of MAS in SLE.

Atypical B cells consist of subsets with distinct functional profiles

Atypical B cells are a population of activated B cells that are commonly enriched in individuals with chronic immune activation but are also part of a normal immune response to infection or vaccination. To better define the role of atypical B cells in the human adaptive immune response, we performed single-cell sequencing of transcriptomes, cell surface markers, and B cell receptors in individuals with chronic exposure to the malaria parasite Plasmodium falciparum, a condition known to lead to accumulation of circulating atypical B cells. We identified three previously uncharacterized populations of atypical B cells with distinct transcriptional and functional profiles and observed marked differences among these three subsets in their ability to produce immunoglobulin G upon T-cell-dependent activation. Our findings help explain the conflicting observations in prior studies regarding the function of atypical B cells and highlight their different roles in the adaptive immune response in chronic inflammatory conditions.

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.