Back

Decoding Spatial Programs in Human Glioblastoma Through Surprisal Information-Theoretical Analysis

Bao, S.; Long, G.; Ghose, S.; Wang, N.; Li, H.; Zhong, M.; Matthews, L.; Lu, Y.; Sheng, J.; Tu, Z.; Escobar, W.; Gopal, P.; McGuone, D.; Erson Omay, Z. E.; Alok, K.; Zhang, D.; DiStasio, M.; Tesileanu, M.; Yang, M.; Li, K.; Moliterno, J.; Heath, J. R.; Raredon, M. S. B.; Remacle, F.; Levine, R. D.; Zhou, J.; Fan, R.

2026-07-08 systems biology
10.64898/2026.06.23.733938 bioRxiv
Show abstract

Glioblastoma (GBM) is a highly heterogeneous and invasive brain tumor in which complex interactions among tumor cells, immune cells, and neurons shape disease progression and therapeutic resistance. Resolving spatial patterns programs in glioblastoma (GBM) requires analytical approaches that go beyond variance-driven embeddings. Here, we applied thermodynamic surprisal analysis, an information-theoretic decomposition framework, to spatial transcriptomic sequencing from human GBM specimens to identify orthogonal constraint modes that capture dominant and previously hidden spatial programs. Surprisal analysis revealed structured patterns in the data that are not highlighted by standard approaches such as PCA or cell type deconvolution, including immune-activation-associated signatures as well as spatial programs consistent with synapse remodeling. Coupling surprisal decomposition with spatial ligand-receptor interaction analysis, NICHES, along with multiplexed protein imaging connected these spatial hidden modes to reveal potential communication networks. Together, these results position surprisal analysis as a powerful, complementary strategy for interrogating spatial tumor architecture, enabling discovery of non-obvious spatial programs and interactions that are obscured by variance-based methods.

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.