# Biologically informed deep learning to query gene programs in single-cell atlases

**Authors:** Mohammad Lotfollahi, Sergei Rybakov, Karin Hrovatin, Soroor Hediyeh-zadeh, Carlos Talavera-López, Alexander V. Misharin, Fabian J. Theis

PMC · DOI: 10.1038/s41556-022-01072-x · 2023-02-02

## TL;DR

The paper introduces ExpiMap, a deep-learning model that improves the interpretability of single-cell data analysis by mapping cells to biologically meaningful gene programs.

## Contribution

ExpiMap introduces a novel deep-learning architecture that maps single-cell data to interpretable gene programs, enhancing biological understanding.

## Key findings

- ExpiMap performs comparably to existing methods but provides additional interpretability through gene programs.
- The model successfully analyzes perturbation responses in various tissues and species.
- ExpiMap resolves patient responses to coronavirus disease 2019 treatments across cell types.

## Abstract

The increasing availability of large-scale single-cell atlases has enabled the detailed description of cell states. In parallel, advances in deep learning allow rapid analysis of newly generated query datasets by mapping them into reference atlases. However, existing data transformations learned to map query data are not easily explainable using biologically known concepts such as genes or pathways. Here we propose expiMap, a biologically informed deep-learning architecture that enables single-cell reference mapping. ExpiMap learns to map cells into biologically understandable components representing known ‘gene programs’. The activity of each cell for a gene program is learned while simultaneously refining them and learning de novo programs. We show that expiMap compares favourably to existing methods while bringing an additional layer of interpretability to integrative single-cell analysis. Furthermore, we demonstrate its applicability to analyse single-cell perturbation responses in different tissues and species and resolve responses of patients who have coronavirus disease 2019 to different treatments across cell types.

Lotfollahi et al. present ExpiMap, a deep-learning model enabling interpretable reference mapping of RNA sequencing data using biologically defined entities, offering end-to-end analysis from dataset integration to functional interpretation.

## Linked entities

- **Diseases:** coronavirus disease 2019 (MONDO:0100096)

## Full-text entities

- **Genes:** Gcg (glucagon) [NCBI Gene 14526] {aka GLP-1, Glu, PPG}, Ifnb1 (interferon beta 1, fibroblast) [NCBI Gene 15977] {aka IFN-beta, IFNB, If1da1, Ifb}, CD19 (CD19 molecule) [NCBI Gene 930] {aka B4, CVID3}, FPR1 (formyl peptide receptor 1) [NCBI Gene 2357] {aka FMLP, FPR}, IFNB1 (interferon beta 1) [NCBI Gene 3456] {aka IFB, IFF, IFN-beta, IFNB}, JCHAIN (joining chain of multimeric IgA and IgM) [NCBI Gene 3512] {aka IGCJ, IGJ, JCH}, CD14 (CD14 molecule) [NCBI Gene 929], GPR166P (G protein-coupled receptor 166, pseudogene) [NCBI Gene 442206] {aka GPCR, PGR9}, Cdkn2a (cyclin dependent kinase inhibitor 2A) [NCBI Gene 12578] {aka ARF-INK4a, Arf, INK4a-ARF, Ink4a/Arf, MTS1, Pctr1}, ANXA1 (annexin A1) [NCBI Gene 301] {aka ANX1, LPC1}, MS4A1 (membrane spanning 4-domains A1) [NCBI Gene 931] {aka B1, Bp35, CD20, CVID5, FMC7, LEU-16}, Cfap126 (cilia and flagella associated protein 126) [NCBI Gene 75472] {aka 1700009P17Rik, Flattop, Fltp}, TMSB4X (thymosin beta 4 X-linked) [NCBI Gene 7114] {aka FX, PTMB4, TB4X, TMSB4}, RNF130 (ring finger protein 130) [NCBI Gene 55819] {aka G1RP, G1RZFP, GOLIATH, GP}, IFIH1 (interferon induced with helicase C domain 1) [NCBI Gene 64135] {aka AGS7, Hlcd, IDDM19, IMD95, MDA-5, MDA5}, CD74 (CD74 molecule) [NCBI Gene 972] {aka CLIP, DHLAG, HLADG, II, Ia-GAMMA, p33}, IL6R (interleukin 6 receptor) [NCBI Gene 3570] {aka CD126, HIES5, IL-1Ra, IL-6R, IL-6R-1, IL-6RA}, CD4 (CD4 molecule) [NCBI Gene 920] {aka CD4mut, IMD79, Leu-3, OKT4D, T4}, EBF1 (EBF transcription factor 1) [NCBI Gene 1879] {aka COE1, EBF, O/E-1, OLF1}, MME (membrane metalloendopeptidase) [NCBI Gene 4311] {aka CALLA, CD10, CMT2T, NEP, SCA43, SFE}, IFNG (interferon gamma) [NCBI Gene 3458] {aka IFG, IFI, IMD69}, IFNA1 (interferon alpha 1) [NCBI Gene 3439] {aka IFL, IFN, IFN-ALPHA, IFN-alphaD, IFNA13, IFNA@}, TNFRSF8 (TNF receptor superfamily member 8) [NCBI Gene 943] {aka CD30, D1S166E, Ki-1}, IFNA8 (interferon alpha 8) [NCBI Gene 3445] {aka IFN-alphaB}, FCGR3A (Fc gamma receptor IIIa) [NCBI Gene 2214] {aka CD16-II, CD16A, FCG3, FCGR3, FCRIIIA, FcGRIIIA}, RIGI (RNA sensor RIG-I) [NCBI Gene 23586] {aka DDX58, RIG-I, RIG1, RLR-1, SGMRT2}, CXCL10 (C-X-C motif chemokine ligand 10) [NCBI Gene 3627] {aka C7, IFI10, INP10, IP-10, SCYB10, crg-2}, CD8A (CD8 subunit alpha) [NCBI Gene 925] {aka CD8, CD8alpha, IMD116, Leu2, p32}, CD52 (CD52 molecule) [NCBI Gene 1043] {aka CDW52, EDDM5, HE5}, CD79A (CD79a molecule) [NCBI Gene 973] {aka IGA, IGAlpha, MB-1, MB1}
- **Diseases:** T2D (MESH:D003924), inflammatory (MESH:D007249), GP (OMIM:614201), Pancreas (MESH:D010190), INTERFERON_SIGNALING (MESH:C535530), UPR (MESH:D011488), COVID (MESH:D000086382), viral infection (MESH:D014777), cancer (MESH:D009369), diabetes (MESH:D003920), scVI (MESH:C567857), severe COVID (MESH:D045169), infection (MESH:D007239), Lupus (MESH:D008180), ANTIGEN PRESENTATION (MESH:C535887), cytokine storm syndrome (MESH:D000080424)
- **Species:** Homo sapiens (human, species) [taxon 9606], Mus musculus (house mouse, species) [taxon 10090]

## Figures

16 figures with captions in the complete paper: https://tomesphere.com/paper/PMC9928587/full.md

---
Source: https://tomesphere.com/paper/PMC9928587