Single-cell multi-omics analysis reveals IFN-driven alterations in T lymphocytes and natural killer cells in systemic lupus erythematosus

An abnormal increase in CD26(-)CD28(-) cytotoxic effector CD4 and CD8 T cell populations in patients with systemic lupus erythematosus

CD26 is a human T cell costimulatory molecule as well as a T cell subset marker, and the increase of CD26+ T cells in inflamed tissues and peripheral blood has been reported in diverse autoimmune diseases. In contrast, our group has previously shown that levels of circulating CD26+ T cells are decreased in patients with systemic lupus erythematosus (SLE), although the role of reduced CD26 T cell surface expression in SLE pathology remains to be elucidated. In the present study, we conducted CD26-based T cell subset analyses utilizing peripheral blood mononuclear cells from 57 SLE patients and 31 healthy adult volunteers. We show that the increase in the CD26(-) T cell population reflects the abnormal expansion of CD26(-)CD28(-) cytotoxic subsets of both CD8 T cells and CD4 T cells in SLE patients. Single-cell RNA sequencing analysis of the CD26(-)CD28(-) CD4 and CD8 T cell populations reveals unique characteristics with similarities to natural killer T cells. In addition, the level of CD26(-)CD28(-) T cells is increased in some active-stage SLE patients with renal manifestation. Meanwhile, the effect of prednisolone treatment on these populations varies from patient to patient, with levels of these cytotoxic effector populations still being elevated in some inactive-stage SLE patients. Taken together, our data suggest that analysis of these populations in SLE may be a useful tool to classify this markedly heterogeneous condition.

Multimodal hierarchical classification of CITE-seq data delineates immune cell states across lineages and tissues

Single-cell RNA sequencing (scRNA-seq) is invaluable for profiling cellular heterogeneity and transcriptional states, but transcriptomic profiles do not always delineate subsets defined by surface proteins. Cellular indexing of transcriptomes and epitopes (CITE-seq) enables simultaneous profiling of single-cell transcriptomes and surface proteomes; however, accurate cell-type annotation requires a classifier that integrates multimodal data. Here, we describe multimodal classifier hierarchy (MMoCHi), a marker-based approach for accurate cell-type classification across multiple single-cell modalities that does not rely on reference atlases. We benchmark MMoCHi using sorted T lymphocyte subsets and annotate a cross-tissue human immune cell dataset. MMoCHi outperforms leading transcriptome-based classifiers and multimodal unsupervised clustering in its ability to identify immune cell subsets that are not readily resolved and to reveal subset markers. MMoCHi is designed for adaptability and can integrate annotation of cell types and developmental states across diverse lineages, samples, or modalities.

Advanced methods and novel biomarkers in autoimmune diseases ‑ a review of the recent years progress in systemic lupus erythematosus

There are several autoimmune and rheumatic diseases affecting different organs of the human body. Multiple sclerosis (MS) mainly affects brain, rheumatoid arthritis (RA) mainly affects joints, Type 1 diabetes (T1D) mainly affects pancreas, Sjogren's syndrome (SS) mainly affects salivary glands, while systemic lupus erythematosus (SLE) affects almost every organ of the body. Autoimmune diseases are characterized by production of autoantibodies, activation of immune cells, increased expression of pro-inflammatory cytokines, and activation of type I interferons. Despite improvements in treatments and diagnostic tools, the time it takes for the patients to be diagnosed is too long, and the main treatment for these diseases is still non-specific anti-inflammatory drugs. Thus, there is an urgent need for better biomarkers, as well as tailored, personalized treatment. This review focus on SLE and the organs affected in this disease. We have used the results from various rheumatic and autoimmune diseases and the organs involved with an aim to identify advanced methods and possible biomarkers to be utilized in the diagnosis of SLE, disease monitoring, and response to treatment.

Targeted multiomics in childhood-onset SLE reveal distinct biological phenotypes associated with disease activity: results from an explorative study

OBJECTIVE

To combine targeted transcriptomic and proteomic data in an unsupervised hierarchical clustering method to stratify patients with childhood-onset SLE (cSLE) into similar biological phenotypes, and study the immunological cellular landscape that characterises the clusters.

METHODS

Targeted whole blood gene expression and serum cytokines were determined in patients with cSLE, preselected on disease activity state (at diagnosis, Low Lupus Disease Activity State (LLDAS), flare). Unsupervised hierarchical clustering, agnostic to disease characteristics, was used to identify clusters with distinct biological phenotypes. Disease activity was scored by clinical SELENA-SLEDAI (Safety of Estrogens in Systemic Lupus Erythematosus National Assessment-Systemic Lupus Erythematosus Disease Activity Index). High-dimensional 40-colour flow cytometry was used to identify immune cell subsets.

RESULTS

Three unique clusters were identified, each characterised by a set of differentially expressed genes and cytokines, and by disease activity state: cluster 1 contained primarily patients in LLDAS, cluster 2 contained mainly treatment-naïve patients at diagnosis and cluster 3 contained a mixed group of patients, namely in LLDAS, at diagnosis and disease flare. The biological phenotypes did not reflect previous organ system involvement and over time, patients could move from one cluster to another. Healthy controls clustered together in cluster 1. Specific immune cell subsets, including CD11c+ B cells, conventional dendritic cells, plasmablasts and early effector CD4+ T cells, differed between the clusters.

CONCLUSION

Using a targeted multiomic approach, we clustered patients into distinct biological phenotypes that are related to disease activity state but not to organ system involvement. This supports a new concept where choice of treatment and tapering strategies are not solely based on clinical phenotype but includes measuring novel biological parameters.

Low-dose IL-2 reduces IL-21+ T cell frequency and induces anti-inflammatory gene expression in type 1 diabetes

Despite early clinical successes, the mechanisms of action of low-dose interleukin-2 (LD-IL-2) immunotherapy remain only partly understood. Here we examine the effects of interval administration of low-dose recombinant IL-2 (iLD-IL-2) in type 1 diabetes using high-resolution single-cell multiomics and flow cytometry on longitudinally-collected peripheral blood samples. Our results confirm that iLD-IL-2 selectively expands thymic-derived FOXP3+HELIOS+ regulatory T cells and CD56bright NK cells, and show that the treatment reduces the frequency of IL-21-producing CD4+ T cells and of two innate-like mucosal-associated invariant T and Vγ9Vδ2 CD8+ T cell subsets. The cellular changes induced by iLD-IL-2 associate with an anti-inflammatory gene expression signature, which remains detectable in all T and NK cell subsets analysed one month after treatment. These findings warrant investigations into the potential longer-term clinical benefits of iLD-IL-2 in immunotherapy.

Accurate Machine Learning Model to Diagnose Chronic Autoimmune Diseases Utilizing Information From B Cells and Monocytes

Heterogeneity and limited comprehension of chronic autoimmune disease pathophysiology cause accurate diagnosis a challenging process. With the increasing resources of single-cell sequencing data, a reasonable way could be found to address this issue. In our study, with the use of large-scale public single-cell RNA sequencing (scRNA-seq) data, analysis of dataset integration (3.1 × 10⁵ PBMCs from fifteen SLE patients and eight healthy donors) and cellular cross talking (3.8 × 10⁵ PBMCs from twenty-eight SLE patients and eight healthy donors) were performed to identify the most crucial information characterizing SLE. Our findings revealed that the interactions among the PBMC subpopulations of SLE patients may be weakened under the inflammatory microenvironment, which could result in abnormal emergences or variations in signaling patterns within PBMCs. In particular, the alterations of B cells and monocytes may be the most significant findings. Utilizing this powerful information, an efficient mathematical model of unbiased random forest machine learning was established to distinguish SLE patients from healthy donors via not only scRNA-seq data but also bulk RNA-seq data. Surprisingly, our mathematical model could also accurately identify patients with rheumatoid arthritis and multiple sclerosis, not just SLE, via bulk RNA-seq data (derived from 688 samples). Since the variations in PBMCs should predate the clinical manifestations of these diseases, our machine learning model may be feasible to develop into an efficient tool for accurate diagnosis of chronic autoimmune diseases.