New Computational Framework Tracks Cell Migration Across Space and Time

Scientists have unveiled a new computational framework that brings together two of the most powerful technologies in modern biology—single-cell RNA sequencing and spatial transcriptomics—to reconstruct how populations of cells redistribute and change state across time within living tissues. The work, published in BMC Bioinformatics, addresses one of the most persistent gaps in computational biology: the difficulty of watching cells move, settle, and transform inside complex tissue environments over the course of development or disease progression. By combining time-series gene expression profiles with spatially resolved molecular maps, the framework generates temporally resolved pictures of cellular organization that neither technology can provide on its own.

Single-cell RNA sequencing has transformed biology by allowing researchers to measure the complete transcriptional activity of thousands of individual cells in a single experiment. This technology reveals extraordinary diversity within tissues, distinguishing cell types and even subtle cell states that bulk measurements would average away. However, scRNA-seq carries a fundamental limitation: during the dissociation process required to prepare samples, cells are physically separated from their tissue context. The resulting data describe what each cell is doing, but not where it was located. Spatial transcriptomics was developed precisely to solve this problem, preserving the spatial organization of cells within tissue sections while still capturing gene expression information. Yet even the most advanced spatial platforms typically profile fewer genes per cell or provide lower molecular resolution than dissociative single-cell methods, and most experiments capture only a single moment in time.

The new framework, developed by Mahdi Pursalim, Kaveh Kavousi, and Parisa Shooshtari, tackles the integration problem directly. The approach operates in two complementary stages. In the first stage, the method performs reference-based spatial deconvolution, a computational procedure that uses scRNA-seq data as a reference atlas to estimate which cell types are present in each spot or region of a spatial transcriptomics measurement. Because spatial platforms often measure mixtures of cells within each capture location, deconvolution is essential for untangling these mixed signals and recovering the underlying cellular composition. The second stage applies single-cell spatial mapping, transferring individual cells from the scRNA-seq dataset onto their most likely spatial positions within the tissue, guided by the gene expression similarities between the single-cell profiles and the spatial measurements.

What distinguishes this framework from previous integration efforts is its explicit treatment of time. The method is designed to work with time-series data—sets of scRNA-seq and spatial transcriptomics measurements collected at multiple, consecutive time points. By applying the deconvolution and mapping procedures independently at each time point and then comparing the results across consecutive stages, the framework quantifies how spatial cell-type distributions shift over the course of a biological process. These temporal changes are summarized using complementary abundance-weighted spatial descriptors, quantitative measures that capture both the proportion of each cell type and its spatial arrangement within the tissue. In parallel, the framework characterizes inferred cell-state dynamics by tracking changes in dominant cell-type assignments across consecutive time points, providing a computational estimate of how the state composition of the tissue evolves.

Recognizing that any computational method is only as trustworthy as its validation, the researchers subjected the spatial abundance component of their framework to rigorous simulation-based benchmarking. They generated synthetic spatial datasets with known, predefined cell-type compositions, using these simulated ground truths to quantitatively compare their approach against alternative spatial deconvolution and mapping methods. This evaluation strategy is critical in a field where real experimental data rarely come with complete knowledge of the true cellular makeup of every spatial location. The benchmarking results support the framework’s ability to recover simulated spatial cell-type abundance patterns, offering reassurance that the estimates it produces reflect genuine biological structure rather than computational artifacts.

To demonstrate the framework’s biological utility and generalizability, the team applied it to two systems at opposite ends of the biological spectrum: human glioblastoma organoids and embryonic chicken heart development. Glioblastoma organoids are laboratory-grown, three-dimensional structures derived from human cells that recapitulate key features of aggressive brain tumors, making them valuable models for studying cancer progression in a controlled setting. The framework revealed temporal patterns of spatial cell-type organization and inferred cell-state composition within these organoids that could not be detected through static analyses alone, offering a dynamic view of how the tumor-like cellular ecosystem reorganizes over time.

The second application, embryonic chicken heart development, represents one of the most dramatic examples of coordinated cellular behavior in biology. The heart forms through precisely choreographed movements of diverse cell populations, and disruptions to these processes underlie many congenital heart defects. When the researchers applied their framework to time-series data from developing chicken hearts, it captured temporal patterns in spatial cell-type redistribution and inferred cell-state dynamics that illuminate how the cellular architecture of this vital organ emerges during embryogenesis. Together, the two case studies—one disease-focused, one developmental—demonstrate that the framework is not tailored to a single biological context but can be applied broadly across systems where time-series single-cell and spatial data are available.

The authors are careful to frame their results with appropriate scientific caution, a nuance that matters for how the work should be interpreted. The framework captures what the researchers describe as spatiotemporal cell-type redistribution and inferred cell-state dynamics—computational measures derived from molecular data. These inferred quantities should be understood as descriptive computational measures rather than direct evidence of physical cell migration or lineage-validated state transitions. In other words, observing that the abundance of a cell type increases in a particular spatial region between two time points does not prove that cells physically traveled there; the change could also reflect differential proliferation, cell death, or changes in the sampling of the tissue. Similarly, apparent shifts in cell-state composition represent computational inferences rather than experimentally validated transitions of individual cells between states.

This caution reflects a broader truth about computational modeling of biological dynamics. Single-cell and spatial transcriptomics are destructive measurements—each sample is consumed in the process of profiling it—so reconstructing temporal behavior necessarily involves stitching together observations from different specimens assumed to represent the same underlying process. Despite this inherent limitation, the value of such models is considerable. They generate testable hypotheses about which cell populations move where, when specific state transitions occur, and how the spatial organization of tissues changes during development or disease. These hypotheses can then guide targeted follow-up experiments, such as live imaging or lineage tracing, that directly observe the behaviors the models predict.

The publication arrives at a moment when the experimental technologies underpinning the framework are advancing rapidly. Spatial transcriptomics platforms continue to increase in resolution, throughput, and affordability, and time-course studies combining scRNA-seq with spatial profiling are becoming increasingly common in developmental biology, cancer research, and regenerative medicine. What has lagged behind is the computational toolkit for integrating these data modalities across time in a principled, quantitative, and benchmarked way. By providing a broadly applicable framework with validated spatial abundance estimation, temporally resolved descriptors of cellular organization, and demonstrated performance across biologically distinct systems, the researchers offer the community a practical resource for systematic characterization of dynamic cellular behaviors. The work points toward a future in which the choreography of cells—their movements, positions, and transformations—can be reconstructed computationally from molecular snapshots, deepening understanding of how tissues build themselves, how tumors evolve, and how these processes might ultimately be modulated for therapeutic benefit.

The methodological choices behind the framework reflect practical constraints inherent to current spatial transcriptomics platforms. Reference-based deconvolution is necessary because many widely used spatial assays capture expression from groups of cells within each measured location, meaning the raw signal represents a mixture rather than a single cell. By anchoring this unmixing step to scRNA-seq references, the framework leverages the superior molecular resolution of dissociative profiling while retaining the positional information that spatial measurements provide. The complementary use of single-cell spatial mapping then adds a second layer of inference, assigning individual cells from the reference to plausible tissue positions based on expression similarity.

The benchmarking strategy deserves particular attention from readers evaluating the method. Simulation-based validation, in which synthetic spatial datasets are generated with known cell-type compositions serving as ground truth, allows quantitative comparison against alternative deconvolution and mapping approaches in a way that real tissue data cannot. Because no experimental measurement of a real tissue comes with a complete, verified inventory of every cell type at every location, simulations provide the only setting where estimation accuracy can be measured directly. The reported performance of the spatial abundance component against competing methods therefore offers a meaningful, if partial, assessment of reliability.

The two biological applications also illustrate the diversity of questions the framework can address. Glioblastoma organoids and the embryonic chicken heart differ not only in species and biological context but in the character of their cellular dynamics: one models the progressive reorganization of a tumor-like cellular ecosystem, while the other captures the choreographed emergence of organ architecture during embryogenesis. That a single computational pipeline can extract temporally resolved patterns of spatial cell-type organization and inferred cell-state composition from both systems suggests the approach is not dependent on tissue-specific assumptions, an important property for a tool intended for broad community use.

The authors’ framing of their results as descriptive computational measures rather than direct evidence of physical cell migration or lineage-validated state transitions is a notable feature of the work. Such interpretive restraint is uncommon and valuable in a field where language about cell movement can easily outpace what molecular snapshot data can actually demonstrate. The framework is published as open access under a Creative Commons license, and the article appeared as a citable, peer-reviewed accepted manuscript carrying a permanent DOI, with a final Version of Record to follow. The research was supported in part by the Children’s Health Research Institute and the Ontario Institute for Cancer Research.

Subject of Research: Integrative computational modeling of cell migration and tissue organization using time-series single-cell RNA sequencing and spatial transcriptomics data

Article Title: Systems-level modeling of cell migration using spatially-resolved single-cell data

Article References: Pursalim, M., Kavousi, K., & Shooshtari, P. (2026). Systems-level modeling of cell migration using spatially-resolved single-cell data. BMC Bioinformatics. https://doi.org/10.1186/s12859-026-06649-z

Image Credits: AI Generated

DOI: 10.1186/s12859-026-06649-z

Keywords: cell migration, spatial transcriptomics, single-cell RNA sequencing, spatiotemporal modeling, deconvolution, cell-state dynamics, glioblastoma organoids, chicken heart development, time-series analysis, computational biology, tissue organization, cell types