Back

ST-Assign: a probabilistic model for joint cell type identification in spatial transcriptomics and single-cell RNA sequencing data

Geras, A.; Szczurek, E.

2023-05-30 bioinformatics
10.1101/2023.05.29.542559 bioRxiv
Show abstract

Understanding the intricate composition of tissues in complex living organisms is crucial for unraveling the mechanisms underlying health and disease. This study addresses the challenge of dissecting cell types within tissues by integrating information from two powerful experimental techniques: single-cell RNA-sequencing (scRNA-seq) and spatial transcriptomics (ST). While scRNA-seq offers insights into transcriptional heterogeneity at the cellular level, ST provides spatial information within tissues. Current methods for cell-type annotation in scRNA-seq and mixture decomposition in ST data are often conducted independently, resulting in reduced statistical power and accuracy. To bridge this gap, we propose ST-Assign, a novel hierarchical Bayesian probabilistic model that jointly performs cell-type annotation in scRNA-seq data and cell-type mixture decomposition in ST data. ST-Assign accounts for shared variables such as gene expression profiles and leverages prior knowledge about marker genes, amplifying statistical strength and mitigating experimental noise. The models excellent performance is demonstrated on simulated and real mouse brain data, showcasing accurate cell-type mixture decomposition and cell-type assignment. In comparison to existing tools, ST-Assign demonstrates superior capabilities, particularly in the task of assigning cell types to individual cells. ST-Assign enables exploring the spatial composition of cell types and holds the potential for enhancing our comprehension of diverse biological systems.

Matching journals

The top 5 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.