AUTHOR=Li Yixi , Ma Chiyu , Liao Shengyou , Qi Suwen , Meng Shuhui , Cai Wanxia , Dai Weier , Cao Rui , Dong Xiangnan , Krämer Bernhard K. , Yun Chen , Hocher Berthold , Hong Xiaoping , Liu Dongzhou , Tang Donge , He Jingquan , Yin Lianghong , Dai Yong TITLE=Combined proteomics and single cell RNA-sequencing analysis to identify biomarkers of disease diagnosis and disease exacerbation for systemic lupus erythematosus JOURNAL=Frontiers in Immunology VOLUME=13 YEAR=2022 URL=https://www.frontiersin.org/journals/immunology/articles/10.3389/fimmu.2022.969509 DOI=10.3389/fimmu.2022.969509 ISSN=1664-3224 ABSTRACT=Introduction

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease for which there is no cure. Effective diagnosis and precise assessment of disease exacerbation remains a major challenge.

Methods

We performed peripheral blood mononuclear cell (PBMC) proteomics of a discovery cohort, including patients with active SLE and inactive SLE, patients with rheumatoid arthritis (RA), and healthy controls (HC). Then, we performed a machine learning pipeline to identify biomarker combinations. The biomarker combinations were further validated using enzyme-linked immunosorbent assays (ELISAs) in another cohort. Single-cell RNA sequencing (scRNA-seq) data from active SLE, inactive SLE, and HC PBMC samples further elucidated the potential immune cellular sources of each of these PBMC biomarkers.

Results

Screening of the PBMC proteome identified 1023, 168, and 124 proteins that were significantly different between SLE vs. HC, SLE vs. RA, and active SLE vs. inactive SLE, respectively. The machine learning pipeline identified two biomarker combinations that accurately distinguished patients with SLE from controls and discriminated between active and inactive SLE. The validated results of ELISAs for two biomarker combinations were in line with the discovery cohort results. Among them, the six-protein combination (IFIT3, MX1, TOMM40, STAT1, STAT2, and OAS3) exhibited good performance for SLE disease diagnosis, with AUC of 0.723 and 0.815 for distinguishing SLE from HC and RA, respectively. Nine-protein combination (PHACTR2, GOT2, L-selectin, CMC4, MAP2K1, CMPK2, ECPAS, SRA1, and STAT2) showed a robust performance in assessing disease exacerbation (AUC=0.990). Further, the potential immune cellular sources of nine PBMC biomarkers, which had the consistent changes with the proteomics data, were elucidated by PBMC scRNAseq.

Discussion

Unbiased proteomic quantification and experimental validation of PBMC samples from two cohorts of patients with SLE were identified as biomarker combinations for diagnosis and activity monitoring. Furthermore, the immune cell subtype origin of the biomarkers in the transcript expression level was determined using PBMC scRNAseq. These findings present valuable PBMC biomarkers associated with SLE and may reveal potential therapeutic targets.