Back

A transcriptional patient map of systemic lupus erythematosus reveals disease-related multicellular immune programs conserved between blood and kidney

Linares-Blanco, J.; Schäfer, P. S. L.; Zimmermann, L.; Melo Ferreira, R.; Toro-Dominguez, D.; Carmona-Saez, P.; Tanevski, J.; Alarcon Riquelme, M. E.; Eadon, M. T.; Ramirez Flores, R. O.; Saez-Rodriguez, J.

2026-05-01 systems biology
10.64898/2026.04.28.721379 bioRxiv
Show abstract

Systemic lupus erythematosus (SLE) shows marked clinical and molecular heterogeneity, yet patient stratification often relies on gene expression signatures lacking multicellular context. Here we construct a transcriptional patient map of SLE by analyzing 1,167 total samples (783 SLE, 384 healthy controls) across different resolutions, including single-cell and bulk blood as well as spatially resolved kidney tissue transcriptomes. Using an unsupervised approach we inferred patient-level transcriptomic immune programs from two independent single-cell RNA sequencing cohorts of peripheral blood mononuclear cells (PBMCs), capturing both differences between SLE and health as well as within-SLE heterogeneity. Specifically, we identified four conserved programs comprising two multicellular inflammatory programs driven by interferon and TNF/NFkB activity across immune cells, and two cell type-specific programs reflecting CD8 T cell cytotoxicity and a CD4 T cell naive-to-effector state. Functional analysis of these programs revealed a rewiring of both cell-to-cell interactions and task allocation across cell types during disease activation. In addition, mapping these programs onto an external longitudinal blood transcriptomic cohort predicted flare risk and identified candidate blood protein biomarkers detectable by proteomics. Finally, we showed that these blood programs were enriched in immune-infiltrated glomerular regions from kidney biopsies of individuals with lupus nephritis using spatially resolved transcriptomic data, thereby linking systemic immune programs to local tissue pathology. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=115 SRC="FIGDIR/small/721379v1_ufig1.gif" ALT="Figure 1"> View larger version (27K): org.highwire.dtl.DTLVardef@1e75f00org.highwire.dtl.DTLVardef@10e0a63org.highwire.dtl.DTLVardef@cc194forg.highwire.dtl.DTLVardef@191b4d2_HPS_FORMAT_FIGEXP M_FIG C_FIG

Matching journals

The top 6 journals account for 50% of the predicted probability mass.

50% of probability mass above

"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.