Multi-scale modeling of human tissues from spatial transcriptomics with TERRA

Abstract

Spatial transcriptomics maps gene expression at cellular resolution, revealing how cells organize into multicellular niches. Yet computational analyses remain dataset-specific, without a transferable representation of tissue organization that generalizes across datasets, tasks and tissues or predicts how tissues behave under perturbation. We present TERRA, a foundation model pretrained on 112 million human cells profiled by spatial transcriptomics. From a single pretrained backbone, TERRA yields embeddings at the scale of cells, the genes they express and the neighborhoods in which they reside, and supports spatial in silico perturbation, all applied zero-shot to unseen tissues. At the cell level, in newly generated spatial data for developing pancreas, TERRA identified an islet-associated capillary state which we posit represents a developmental precursor of the mature islet microvasculature. At the gene level, in untreated kidney sections, in silico knockout of immune-checkpoint targets predicted a gene program of immune-checkpoint-blockade-associated nephrotoxicity, which we validated in treatment-exposed tissue and recovered in blood. At the neighborhood level, TERRA mapped macrophages across tissues to identify recurring cross-organ niches, which we term archetypes, including a tumor-boundary niche associated with poor prognosis in kidney cancer. Together, TERRA captures the spatial and multicellular logic of human tissue and predicts, in silico, its response to perturbation, providing a multi-scale framework for tissue biology, therapeutic development and clinical application.