STIE: Single-cell level deconvolution, convolution, and clustering in spatial transcriptomics by aligning spot level transcriptome to nuclear morphology

In spot-based spatial transcriptomics, spots that are of the same size and printed at the fixed location cannot precisely capture the actual randomly located single cells, therefore failing to profile the transcriptome at the single-cell level. The current studies primarily focused on enhancing the spot resolution in size via computational imputation or technical improvement, however, they largely overlooked that single-cell resolution, i.e., resolution in cellular or even smaller size, does not equal single-cell level. Using both real and simulated spatial transcriptomics data, we demonstrated that even the high-resolution spatial transcriptomics still has a large number of spots partially covering multiple cells simultaneously, revealing the intrinsic non-single-cell level of spot-based spatial transcriptomics regardless of spot size. To this end, we present STIE, an EM algorithm that aligns the spatial transcriptome to its matched histology image-based nuclear morphology and recovers missing cells from up to ∼70% gap area between spots via the nuclear morphological similarity and neighborhood information, thereby achieving the real single-cell level and whole-slide scale deconvolution/convolution and clustering for both low- and high-resolution spots. On both real and simulation spatial transcriptomics data, STIE characterizes the cell-type specific gene expression variation and demonstrates the outperforming concordance with the single-cell RNAseq-derived cell type transcriptomic signatures compared to the other spot- and subspot-level methods. Furthermore, STIE enabled us to gain novel insights that failed to be revealed by the existing methods due to the lack of single-cell level, for instance, lower actual spot resolution than its reported spot size, the additional contribution of cellular morphology to cell typing beyond transcriptome, unbiased evaluation of cell type colocalization, superior power of high-resolution spot in distinguishing nuanced cell types, and spatially resolved cell-cell interactions at the single-cell level other than spot level. The STIE code is publicly available as an R package at https://github.com/zhushijia/STIE.

[1]  C. Danko,et al.  Cell type and gene expression deconvolution with BayesPrism enables Bayesian integrative analysis across bulk and single-cell RNA sequencing in oncology , 2022, Nature Cancer.

[2]  Lani F. Wu,et al.  Integrative spatial analysis of cell morphologies and transcriptional states with MUSE , 2022, Nature Biotechnology.

[3]  Evan Z. Macosko,et al.  Deep learning and alignment of spatially resolved single-cell transcriptomes with Tangram , 2021, Nature Methods.

[4]  Mingyao Li,et al.  SpaGCN: Integrating gene expression, spatial location and histology to identify spatial domains and spatially variable genes by graph convolutional network , 2021, Nature Methods.

[5]  Å. Borg,et al.  Spatial deconvolution of HER2-positive breast cancer delineates tumor-associated cell type interactions , 2021, Nature Communications.

[6]  M. Unser,et al.  DeepImageJ: A user-friendly environment to run deep learning models in ImageJ , 2021, Nature Methods.

[7]  Stephen R. Williams,et al.  A single-cell and spatially resolved atlas of human breast cancers , 2021, Nature Genetics.

[8]  Raphael Gottardo,et al.  Spatial transcriptomics at subspot resolution with BayesSpace , 2021, Nature Biotechnology.

[9]  Romain F. Laine,et al.  Democratising deep learning for microscopy with ZeroCostDL4Mic , 2021, Nature Communications.

[10]  Guocheng Yuan,et al.  SpatialDWLS: accurate deconvolution of spatial transcriptomic data , 2021, Genome Biology.

[11]  J. Lundeberg,et al.  Spatial Transcriptomics to define transcriptional patterns of zonation and structural components in the mouse liver , 2020, Nature Communications.

[12]  Cindy C. Guo,et al.  High-Spatial-Resolution Multi-Omics Sequencing via Deterministic Barcoding in Tissue , 2020, Cell.

[13]  Joseph Bergenstråhle,et al.  Single-cell and spatial transcriptomics enables probabilistic inference of cell type topography , 2020, Communications Biology.

[14]  Lihua Zhang,et al.  Inference and analysis of cell-cell communication using CellChat , 2020, Nature Communications.

[15]  Holger Heyn,et al.  Seeded NMF regression to Deconvolute Spatial Transcriptomics Spots with Single-Cell Transcriptomes , 2020 .

[16]  Q. Nguyen,et al.  stLearn: integrating spatial location, tissue morphology and gene expression to find cell types, cell-cell interactions and spatial trajectories within undissociated tissues , 2020, bioRxiv.

[17]  Rafael A. Irizarry,et al.  Robust decomposition of cell type mixtures in spatial transcriptomics , 2020, Nature Biotechnology.

[18]  Joseph Bergenstråhle,et al.  Spatially Resolved Transcriptomes—Next Generation Tools for Tissue Exploration , 2020, BioEssays : news and reviews in molecular, cellular and developmental biology.

[19]  Hao Chen,et al.  A Multi-Organ Nucleus Segmentation Challenge , 2020, IEEE Transactions on Medical Imaging.

[20]  D. Ribatti,et al.  Epithelial-Mesenchymal Transition in Cancer: A Historical Overview , 2020, Translational oncology.

[21]  J. Kleinman,et al.  Transcriptome-scale spatial gene expression in the human dorsolateral prefrontal cortex , 2020, Nature Neuroscience.

[22]  Joseph Bergenstråhle,et al.  Super-resolved spatial transcriptomics by deep data fusion , 2020, Nature Biotechnology.

[23]  Zev J. Gartner,et al.  ZipSeq: barcoding for real-time mapping of single cell transcriptomes , 2020, Nature Methods.

[24]  Richard Bonneau,et al.  High-definition spatial transcriptomics for in situ tissue profiling , 2019, Nature Methods.

[25]  Guo-Cheng Yuan,et al.  Accurate estimation of cell-type composition from gene expression data , 2019, Nature Communications.

[26]  Wei Guo,et al.  SCINA: A Semi-Supervised Subtyping Algorithm of Single Cells and Bulk Samples , 2019, Genes.

[27]  Guo-Cheng Yuan,et al.  Transcriptome-scale super-resolved imaging in tissues by RNA seqFISH+ , 2019, Nature.

[28]  Evan Z. Macosko,et al.  Slide-seq: A scalable technology for measuring genome-wide expression at high spatial resolution , 2019, Science.

[29]  D. Larsimont,et al.  Transcriptional output, cell-type densities, and normalization in spatial transcriptomics , 2018, bioRxiv.

[30]  Damian Szklarczyk,et al.  STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets , 2018, Nucleic Acids Res..

[31]  M. Gezmen-Karadag,et al.  The multiple functions and mechanisms of osteopontin. , 2018, Clinical biochemistry.

[32]  William E. Allen,et al.  Three-dimensional intact-tissue sequencing of single-cell transcriptional states , 2018, Science.

[33]  S. Orkin,et al.  Mapping the Mouse Cell Atlas by Microwell-Seq , 2018, Cell.

[34]  V. Weaver,et al.  α5β1-Integrin promotes tension-dependent mammary epithelial cell invasion by engaging the fibronectin synergy site , 2017, Molecular biology of the cell.

[35]  D. Schreiner,et al.  An alternative splicing switch shapes neurexin repertoires in principal neurons versus interneurons in the mouse hippocampus , 2016, eLife.

[36]  Jeffrey R Moffitt,et al.  High-performance multiplexed fluorescence in situ hybridization in culture and tissue with matrix imprinting and clearing , 2016, Proceedings of the National Academy of Sciences.

[37]  Cynthia C. Hession,et al.  Div-Seq: Single-nucleus RNA-Seq reveals dynamics of rare adult newborn neurons , 2016, Science.

[38]  Patrik L. Ståhl,et al.  Visualization and analysis of gene expression in tissue sections by spatial transcriptomics , 2016, Science.

[39]  X. Zhuang,et al.  Spatially resolved, highly multiplexed RNA profiling in single cells , 2015, Science.

[40]  Peng Li,et al.  Controlling cell–cell interactions using surface acoustic waves , 2014, Proceedings of the National Academy of Sciences.

[41]  L. Cai,et al.  Single-cell in situ RNA profiling by sequential hybridization , 2014, Nature Methods.

[42]  George M. Church,et al.  Highly Multiplexed Subcellular RNA Sequencing in Situ , 2014, Science.

[43]  J. Schwarzbauer,et al.  Mammary epithelial cell interactions with fibronectin stimulate epithelial-mesenchymal transition , 2013, Oncogene.

[44]  Dmitry I. Strokotov,et al.  Is there a difference between T- and B-lymphocyte morphology? , 2009, Journal of biomedical optics.

[45]  E. Dahl,et al.  Dual role of macrophage migration inhibitory factor (MIF) in human breast cancer , 2009, BMC Cancer.

[46]  S. Sleijfer,et al.  Association of an Extracellular Matrix Gene Cluster with Breast Cancer Prognosis and Endocrine Therapy Response , 2008, Clinical Cancer Research.

[47]  Allan R. Jones,et al.  Genome-wide atlas of gene expression in the adult mouse brain , 2007, Nature.