Fabian J. Theis

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

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.

Learning Discrete Cell and Niche Codes from Spatial Transcriptomics Using Dual Residual Vector Quantization

Creating bottom-up RNA transfer vehicles from synthetic protein assemblies.

Evolution guides biological systems to populate ecological niches, with viruses among the most successful examples of this principle. Viruses evolved over billions of years to efficiently transfer genetic information. Although viruses are highly diverse, most have converged towards remarkable similarity in the size and shape of their capsids1,2. By contrast, generative models for protein design enable the creation of protein architectures that are absent from nature3-5. Here we investigate whether protein assemblies designed by artificial intelligence can be functionalized to construct nucleic acid transport vehicles that are independent of evolutionary trajectories. By combining natural protein domains with synthetic protein assemblies, we create more than 100 bottom-up RNA transfer vehicles with unique sizes and shapes. These vehicles surpass the RNA transfer efficiency of widely used delivery vehicles by several orders of magnitude. In addition, we demonstrate that their tropism can be programmed by incorporation of computationally designed peptide binders and use them to deliver therapeutically relevant cargo RNAs into a wide range of cellular models. We show the in vivo biodistribution of one of these vehicles in a mouse at near-single-cell resolution, confirm its safety, and use it to perform a gene-editing treatment strategy for Duchenne muscular dystrophy in patient-derived cells and a pig. Our work demonstrates how proteins created by generative artificial intelligence can be harnessed for the rational engineering of RNA transport systems with the desired properties by overcoming the limitations of natural protein diversity.

Plasticity of human microglia and brain perivascular macrophages in aging and Alzheimer’s disease

Myeloid cells, including microglia and perivascular macrophages, are central to Alzheimer’s disease (AD) neurobiology, yet their role remains incompletely understood. We profiled 832,505 human myeloid cells from the prefrontal cortex of 1,607 donors spanning the lifespan and showing varying degrees of AD neuropathology. We delineated six subclasses comprising 13 transcriptionally distinct subtypes and identified adaptive changes associated with aging and AD progression. Here we show that a disease-associated microglial subtype, characterized by elevated GPNMB expression and enriched for polygenic AD risk, expands with AD pathology and shows increased phagocytic activity. We identify MITF as an upstream regulator required to maintain this microglial state. Cell–cell interaction analyses prioritize APOE–SORL1 and APOE–TREM2 signaling pairs associated with disease progression. Using human and mouse models, we demonstrate that the neuroprotective effects of this microglial subtype depend on TREM2. These findings provide mechanistic insights into myeloid cell function in aging and AD, aiding therapeutic discovery. Molecular profiling of myeloid cells from the prefrontal cortex of 1,607 donors with varying degrees of Alzheimer’s disease neuropathology delineates distinct myeloid subtypes and identifies changes associated with aging and disease progression.

A multimodal interrogation of Broca’s area in the pediatric and adult human brain

Language is a defining trait of our species, and disruptions in language acquisition can have profound consequences to the individuals affected. Uncovering the neurodevelopmental basis of this complex trait requires detailed molecular and cellular insights into the neocortical areas that support linguistic abilities. Here we performed joint gene expression and chromatin accessibility profiling at single-nucleus resolution (10x Genomics Single cell Multiome) and spatial transcriptomic profiling (Xenium high-plex in situ spatial transcriptomics) of Broca’s area alongside adjacent motor cortical areas. We profiled individuals from different ancestries (European and African) and developmental stages (infancy, childhood, adolescence, and adulthood). We provide a high-resolution dissection of the cellular and molecular architecture of Broca’s and motor cortical areas across early life stages and anchor the trajectories to the cellular states found in the adult human brain. We identify distinct area- and stage-specific cellular signatures, including a prominent role of glia populations and interneuron subtypes contributing to cytoarchitectonic specializations. Using longitudinal single cell spatial transcriptomic profiling, we orthogonally validate our consensus cell taxonomy and spatially resolve layer enrichment of neuronal and astrocyte subtypes that distinguish Broca’s area and motor cortex. We also uncover cell type-specific molecular signatures that distinguish cell developmental trajectories in these cortical areas, including an early molecular code established by differential expression of cadherin genes that might contribute to area-specific intercellular communication. We also identify cell type-specific vulnerabilities to language- related neurodevelopmental and neuropsychiatric disorders, with selective susceptibility of particular somatostatin-positive interneuron subtypes to ASD/ADHD. Finally, evolutionary analysis of differentially accessible regions between Broca’s area and motor cortex suggests that genetic mutations that might have contributed to the emergence of linguistic abilities accumulated over the course of million years following the divergence of human and chimpanzee lineages. Together, our study provides a comprehensive molecular, cellular and spatial definition of Broca’s area and motor cortex, laying the groundwork for investigations into unique aspects of human cognition and related neurodevelopmental and neuropsychiatric disorders.