Track papers across labs, topics, and trends.
Search papers, labs, and topics across Lattice.
To enable cohort-scale comparison of single-cell distributions while preserving intrinsic manifold geometry, the authors develop Pheno-GS, an unbalanced geodesic optimal transport framework tailored for sparse, noisy, and massive biological datasets. Traditional geodesic OT approaches fail to scale to all-pairs patient comparisons and break down when cell manifolds are unevenly sampled or disconnected. By combining graph connectivity regularization, KL-penalized unbalanced transport, and a batched solver that evaluates all pairs in a single heat diffusion pass, Pheno-GS achieves over a 200× speedup compared to standard Geodesic Sinkhorn across 500 distributions.
All-pairs geodesic optimal transport across patient cell manifolds can now run over 200× faster by solving all distributions simultaneously in a single batched heat diffusion pass.
High-throughput single-cell data is now collected across large patient cohorts. Understanding patient-level heterogeneity from cellular-level data motivates phenoscaping: embedding each single-cell distribution as a"datapoint,"with distances given by optimal transport (OT). Computing geometry-aware OT at this scale, between all pairs of patient datasets, remains an open challenge, since existing methods either rely on Euclidean ground metrics that distort manifold structure or fail under sparse, unevenly sampled, or large-scale data. We present \textbf{Pheno-GS} (Phenoscape-scale Geodesic Sinkhorn), which computes accurate, scalable geodesic transport distances under noisy, unbalanced, large-scale settings via three components: ($1$) graph connectivity regularization for well-defined geodesics on sparse/disconnected manifolds; ($2$) an unbalanced OT formulation via KL marginal penalties; and ($3$) a batched matrix algorithm computing all pairwise distances in one heat diffusion (over $200 \times$ faster than Geodesic Sinkhorn for $500$ distributions). We validate Pheno-GS on synthetic benchmarks and a CyTOF perturbation dataset.