Integrative spatial analysis reveals a multi-layered organization of glioblastoma.

Glioma contains malignant cells in diverse states. Here, we combine spatial transcriptomics, spatial proteomics, and computational approaches to define glioma cellular states and uncover their organization. We find three prominent modes of organization. First, gliomas are composed of small local environments, each typically enriched with one major cellular state. Second, specific pairs of states preferentially reside in proximity across multiple scales. This pairing of states is consistent across tumors. Third, these pairwise interactions collectively define a global architecture composed of five layers. Hypoxia appears to drive the layers, as it is associated with a long-range organization that includes all cancer cell states. Accordingly, tumor regions distant from any hypoxic/necrotic foci and tumors that lack hypoxia such as low-grade IDH-mutant glioma are less organized. In summary, we provide a conceptual framework for the organization of cellular states in glioma, highlighting hypoxia as a long-range tissue organizer.

Mapping the transcriptome: Realizing the full potential of spatial data analysis.

RNA sequencing in situ allows for whole-transcriptome characterization at high resolution, while retaining spatial information. These data present an analytical challenge for bioinformatics-how to leverage spatial information effectively? Properties of data with a spatial dimension require special handling, which necessitate a different set of statistical and inferential considerations when compared to non-spatial data. The geographical sciences primarily use spatial data and have developed methods to analye them. Here we discuss the challenges associated with spatial analysis and examine how we can take advantage of practice from the geographical sciences to realize the full potential of spatial information in transcriptomic datasets.

SnapShot: Spatial transcriptomics.

Spatially resolved transcriptomics methodologies using RNA sequencing principles have and will continue to contribute to decode the molecular landscape of tissues. Linking quantitative sequencing data with tissue morphology empowers profiling of cellular morphology and transcription over time and space in health and disease. To view this SnapShot, open or download the PDF.

Large-Scale Topological Changes Restrain Malignant Progression in Colorectal Cancer.

Widespread changes to DNA methylation and chromatin are well documented in cancer, but the fate of higher-order chromosomal structure remains obscure. Here we integrated topological maps for colon tumors and normal colons with epigenetic, transcriptional, and imaging data to characterize alterations to chromatin loops, topologically associated domains, and large-scale compartments. We found that spatial partitioning of the open and closed genome compartments is profoundly compromised in tumors. This reorganization is accompanied by compartment-specific hypomethylation and chromatin changes. Additionally, we identify a compartment at the interface between the canonical A and B compartments that is reorganized in tumors. Remarkably, similar shifts were evident in non-malignant cells that have accumulated excess divisions. Our analyses suggest that these topological changes repress stemness and invasion programs while inducing anti-tumor immunity genes and may therefore restrain malignant progression. Our findings call into question the conventional view that tumor-associated epigenomic alterations are primarily oncogenic.

A Structured Tumor-Immune Microenvironment in Triple Negative Breast Cancer Revealed by Multiplexed Ion Beam Imaging.

The immune system is critical in modulating cancer progression, but knowledge of immune composition, phenotype, and interactions with tumor is limited. We used multiplexed ion beam imaging by time-of-flight (MIBI-TOF) to simultaneously quantify in situ expression of 36 proteins covering identity, function, and immune regulation at sub-cellular resolution in 41 triple-negative breast cancer patients. Multi-step processing, including deep-learning-based segmentation, revealed variability in the composition of tumor-immune populations across individuals, reconciled by overall immune infiltration and enriched co-occurrence of immune subpopulations and checkpoint expression. Spatial enrichment analysis showed immune mixed and compartmentalized tumors, coinciding with expression of PD1, PD-L1, and IDO in a cell-type- and location-specific manner. Ordered immune structures along the tumor-immune border were associated with compartmentalization and linked to survival. These data demonstrate organization in the tumor-immune microenvironment that is structured in cellular composition, spatial arrangement, and regulatory-protein expression and provide a framework to apply multiplexed imaging to immune oncology.

A Mechanogenetic Toolkit for Interrogating Cell Signaling in Space and Time.

Tools capable of imaging and perturbing mechanical signaling pathways with fine spatiotemporal resolution have been elusive, despite their importance in diverse cellular processes. The challenge in developing a mechanogenetic toolkit (i.e., selective and quantitative activation of genetically encoded mechanoreceptors) stems from the fact that many mechanically activated processes are localized in space and time yet additionally require mechanical loading to become activated. To address this challenge, we synthesized magnetoplasmonic nanoparticles that can image, localize, and mechanically load targeted proteins with high spatiotemporal resolution. We demonstrate their utility by investigating the cell-surface activation of two mechanoreceptors: Notch and E-cadherin. By measuring cellular responses to a spectrum of spatial, chemical, temporal, and mechanical inputs at the single-molecule and single-cell levels, we reveal how spatial segregation and mechanical force cooperate to direct receptor activation dynamics. This generalizable technique can be used to control and understand diverse mechanosensitive processes in cell signaling. VIDEO ABSTRACT.

Spatially resolved clonal copy number alterations in benign and malignant tissue.

Defining the transition from benign to malignant tissue is fundamental to improving early diagnosis of cancer1. Here we use a systematic approach to study spatial genome integrity in situ and describe previously unidentified clonal relationships. We used spatially resolved transcriptomics2 to infer spatial copy number variations in >120,000 regions across multiple organs, in benign and malignant tissues. We demonstrate that genome-wide copy number variation reveals distinct clonal patterns within tumours and in nearby benign tissue using an organ-wide approach focused on the prostate. Our results suggest a model for how genomic instability arises in histologically benign tissue that may represent early events in cancer evolution. We highlight the power of capturing the molecular and spatial continuums in a tissue context and challenge the rationale for treatment paradigms, including focal therapy.

40†years of the homeobox: mechanisms of Hox spatial-temporal collinearity in vertebrates.

Animal body plans are established during embryonic development by the Hox genes. This patterning process relies on the differential expression of Hox genes along the head-to-tail axis. Hox spatial collinearity refers to the relationship between the organization of Hox genes in clusters and the differential Hox expression, whereby the relative order of the Hox genes within a cluster mirrors the spatial sequence of expression in the developing embryo. In vertebrates, the cluster organization is also associated with the timing of Hox activation, which harmonizes Hox expression with the progressive emergence of axial tissues. Thereby, in vertebrates, Hox temporal collinearity is intimately linked to Hox spatial collinearity. Understanding the mechanisms contributing to Hox temporal and spatial collinearity is thus key to the comprehension of vertebrate patterning. Here, we provide an overview of the main discoveries pertaining to the mechanisms of Hox spatial-temporal collinearity.

Benchmarking spatial clustering methods with spatially resolved transcriptomics data.

Spatial clustering, which shares an analogy with single-cell clustering, has expanded the scope of tissue physiology studies from cell-centroid to structure-centroid with spatially resolved transcriptomics (SRT) data. Computational methods have undergone remarkable development in recent years, but a comprehensive benchmark study is still lacking. Here we present a benchmark study of 13 computational methods on 34 SRT data (7 datasets). The performance was evaluated on the basis of accuracy, spatial continuity, marker genes detection, scalability, and robustness. We found existing methods were complementary in terms of their performance and functionality, and we provide guidance for selecting appropriate methods for given scenarios. On testing additional 22 challenging datasets, we identified challenges in identifying noncontinuous spatial domains and limitations of existing methods, highlighting their inadequacies in handling recent large-scale tasks. Furthermore, with 145 simulated data, we examined the robustness of these methods against four different factors, and assessed the impact of pre- and postprocessing approaches. Our study offers a comprehensive evaluation of existing spatial clustering methods with SRT data, paving the way for future advancements in this rapidly evolving field.

The universal visitation law of human mobility.

Human mobility impacts many aspects of a city, from its spatial structure1-3 to its response to an epidemic4-7. It is also ultimately key to social interactions8, innovation9,10 and productivity11. However, our quantitative understanding of the aggregate movements of individuals remains incomplete. Existing models-such as the gravity law12,13 or the radiation model14-concentrate on the purely spatial dependence of mobility flows and do not capture the varying frequencies of recurrent visits to the same locations. Here we reveal a simple and robust scaling law that captures the temporal and spatial spectrum of population movement on the basis of large-scale mobility data from diverse cities around the globe. According to this law, the number of visitors to any location decreases as the inverse square of the product of their visiting frequency and travel distance. We further show that the spatio-temporal flows to different locations give rise to prominent spatial clusters with an area distribution that follows Zipf's law15. Finally, we build an individual mobility model based on exploration and preferential return to provide a mechanistic explanation for the discovered scaling law and the emerging spatial structure. Our findings corroborate long-standing conjectures in human geography (such as central place theory16 and Weber's theory of emergent optimality10) and allow for predictions of recurrent flows, providing a basis for applications in urban planning, traffic engineering and the mitigation of epidemic diseases.

Latent feature extraction with a prior-based self-attention framework for spatial transcriptomics.

Rapid advances in spatial transcriptomics (ST) have revolutionized the interrogation of spatial heterogeneity and increase the demand for comprehensive methods to effectively characterize spatial domains. As a prerequisite for ST data analysis, spatial domain characterization is a crucial step for downstream analyses and biological implications. Here we propose a prior-based self-attention framework for spatial transcriptomics (PAST), a variational graph convolutional autoencoder for ST, which effectively integrates prior information via a Bayesian neural network, captures spatial patterns via a self-attention mechanism, and enables scalable application via a ripple walk sampler strategy. Through comprehensive experiments on data sets generated by different technologies, we show that PAST can effectively characterize spatial domains and facilitate various downstream analyses, including ST visualization, spatial trajectory inference and pseudotime analysis. Also, we highlight the advantages of PAST for multislice joint embedding and automatic annotation of spatial domains in newly sequenced ST data. Compared with existing methods, PAST is the first ST method that integrates reference data to analyze ST data. We anticipate that PAST will open up new avenues for researchers to decipher ST data with customized reference data, which expands the applicability of ST technology.

Bacterial coexistence driven by motility and spatial competition.

Elucidating elementary mechanisms that underlie bacterial diversity is central to ecology1,2 and microbiome research3. Bacteria are known to coexist by metabolic specialization4, cooperation5 and cyclic warfare6-8. Many species are also motile9, which is studied in terms of mechanism10,11, benefit12,13, strategy14,15, evolution16,17 and ecology18,19. Indeed, bacteria often compete for nutrient patches that become available periodically or by random disturbances2,20,21. However, the role of bacterial motility in coexistence remains unexplored experimentally. Here we show that-for mixed bacterial populations that colonize nutrient patches-either population outcompetes the other when low in relative abundance. This inversion of the competitive hierarchy is caused by active segregation and spatial exclusion within the patch: a small fast-moving population can outcompete a large fast-growing population by impeding its migration into the patch, while a small fast-growing population can outcompete a large fast-moving population by expelling it from the initial contact area. The resulting spatial segregation is lost for weak growth-migration trade-offs and a lack of virgin space, but is robust to population ratio, density and chemotactic ability, and is observed in both laboratory and wild strains. These findings show that motility differences and their trade-offs with growth are sufficient to promote diversity, and suggest previously undescribed roles for motility in niche formation and collective expulsion-containment strategies beyond individual search and survival.

Spatial Organization of Single mRNPs at Different Stages of the Gene Expression Pathway.

mRNAs form ribonucleoprotein complexes (mRNPs) by association with proteins that are crucial for mRNA metabolism. While the mRNP proteome has been well characterized, little is known about mRNP organization. Using a single-molecule approach, we show that mRNA conformation changes depending on its cellular localization and translational state. Compared to nuclear mRNPs and lncRNPs, association with ribosomes decompacts individual mRNAs, while pharmacologically dissociating ribosomes or sequestering them into stress granules leads to increased compaction. Moreover, translating mRNAs rarely show co-localized 5' and 3' ends, indicating either that mRNAs are not translated in a closed-loop configuration, or that mRNA circularization is transient, suggesting that a stable closed-loop conformation is not a universal state for all translating mRNAs.

Neuroimaging meta regression for coordinate based meta analysis data with a spatial model.

Coordinate-based meta-analysis combines evidence from a collection of neuroimaging studies to estimate brain activation. In such analyses, a key practical challenge is to find a computationally efficient approach with good statistical interpretability to model the locations of activation foci. In this article, we propose a generative coordinate-based meta-regression (CBMR) framework to approximate a smooth activation intensity function and investigate the effect of study-level covariates (e.g. year of publication, sample size). We employ a spline parameterization to model the spatial structure of brain activation and consider four stochastic models for modeling the random variation in foci. To examine the validity of CBMR, we estimate brain activation on 20 meta-analytic datasets, conduct spatial homogeneity tests at the voxel level, and compare the results to those generated by existing kernel-based and model-based approaches. Keywords: generalized linear models; meta-analysis; spatial statistics; statistical modeling.

Spatial distribution of poultry farms using point pattern modelling: A method to address livestock environmental impacts and disease transmission risks.

The distribution of farm locations and sizes is paramount to characterize patterns of disease spread. With some regions undergoing rapid intensification of livestock production, resulting in increased clustering of farms in peri-urban areas, measuring changes in the spatial distribution of farms is crucial to design effective interventions. However, those data are not available in many countries, their generation being resource-intensive. Here, we develop a farm distribution model (FDM), which allows the prediction of locations and sizes of poultry farms in countries with scarce data. The model combines (i) a Log-Gaussian Cox process model to simulate the farm distribution as a spatial Poisson point process, and (ii) a random forest model to simulate farm sizes (i.e. the number of animals per farm). Spatial predictors were used to calibrate the FDM on intensive broiler and layer farm distributions in Bangladesh, Gujarat (Indian state) and Thailand. The FDM yielded realistic farm distributions in terms of spatial clustering, farm locations and sizes, while providing insights on the factors influencing these distributions. Finally, we illustrate the relevance of modelling realistic farm distributions in the context of epidemic spread by simulating pathogen transmission on an array of spatial distributions of farms. We found that farm distributions generated from the FDM yielded spreading patterns consistent with simulations using observed data, while random point patterns underestimated the probability of large outbreaks. Indeed, spatial clustering increases vulnerability to epidemics, highlighting the need to account for it in epidemiological modelling studies. As the FDM maintains a realistic distribution of farm location and sizes, its use to inform mathematical models of disease transmission is particularly relevant for regions where these data are not available.

A quantitative and spatial analysis of cell cycle regulators during the fission yeast cycle.

We have carried out a systems-level analysis of the spatial and temporal dynamics of cell cycle regulators in the fission yeast Schizosaccharomyces pombe. In a comprehensive single-cell analysis, we have precisely quantified the levels of 38 proteins previously identified as regulators of the G2 to mitosis transition and of 7 proteins acting at the G1- to S-phase transition. Only 2 of the 38 mitotic regulators exhibit changes in concentration at the whole-cell level: the mitotic B-type cyclin Cdc13, which accumulates continually throughout the cell cycle, and the regulatory phosphatase Cdc25, which exhibits a complex cell cycle pattern. Both proteins show similar patterns of change within the nucleus as in the whole cell but at higher concentrations. In addition, the concentrations of the major fission yeast cyclin-dependent kinase (CDK) Cdc2, the CDK regulator Suc1, and the inhibitory kinase Wee1 also increase in the nucleus, peaking at mitotic onset, but are constant in the whole cell. The significant increase in concentration with size for Cdc13 supports the view that mitotic B-type cyclin accumulation could act as a cell size sensor. We propose a two-step process for the control of mitosis. First, Cdc13 accumulates in a size-dependent manner, which drives increasing CDK activity. Second, from mid-G2, the increasing nuclear accumulation of Cdc25 and the counteracting Wee1 introduce a bistability switch that results in a rapid rise of CDK activity at the end of G2 and thus, brings about an orderly progression into mitosis.

Spatial Analysis of Drug-Susceptible and Multidrug-Resistant Cases of Tuberculosis, Ho Chi Minh City, Vietnam, 2020-2023.

We characterized the spatial distribution of drug-susceptible (DS) and multidrug-resistant (MDR) tuberculosis (TB) cases in Ho Chi Minh City, Vietnam, a major metropolis in southeastern Asia, and explored demographic and socioeconomic factors associated with local TB burden. Hot spots of DS and MDR TB incidence were observed in the central parts of Ho Chi Minh City, and substantial heterogeneity was observed across wards. Positive spatial autocorrelation was observed for both DS TB and MDR TB. Ward-level TB incidence was associated with HIV prevalence and the male proportion of the population. No ward-level demographic and socioeconomic indicators were associated with MDR TB case count relative to total TB case count. Our findings might inform spatially targeted TB control strategies and provide insights for generating hypotheses about the nature of the relationship between DS and MDR TB in Ho Chi Minh City and the wider southeastern region of Asia.

An introduction to bayesian spatial smoothing methods for disease mapping: modeling county firearm suicide mortality rates.

This article introduces bayesian spatial smoothing models for disease mapping-a specific application of small area estimation where the full universe of data is known-to a wider audience of public health professionals using firearm suicide as a motivating example. Besag, York, and Mollié (BYM) Poisson spatial and space-time smoothing models were fitted to firearm suicide counts for the years 2014-2018. County raw death rates in 2018 ranged from 0 to 24.81 deaths per 10 000 people. However, the highest mortality rate was highly unstable, based on only 2 deaths in a population of approximately 800, and 80.5% of contiguous US counties experienced fewer than 10 firearm suicide deaths and were thus suppressed. Spatially smoothed county firearm suicide mortality estimates ranged from 0.06 to 4.05 deaths per 10 000 people and could be reported for all counties. The space-time smoothing model produced similar estimates with narrower credible intervals as it allowed counties to gain precision from adjacent neighbors and their own counts in adjacent years. bayesian spatial smoothing methods are a useful tool for evaluating spatial health disparities in small geographies where small numbers can result in highly variable rate estimates, and new estimation techniques in R software have made fitting these models more accessible to researchers.

Advances in spatial transcriptomic data analysis.

Spatial transcriptomics is a rapidly growing field that promises to comprehensively characterize tissue organization and architecture at the single-cell or subcellular resolution. Such information provides a solid foundation for mechanistic understanding of many biological processes in both health and disease that cannot be obtained by using traditional technologies. The development of computational methods plays important roles in extracting biological signals from raw data. Various approaches have been developed to overcome technology-specific limitations such as spatial resolution, gene coverage, sensitivity, and technical biases. Downstream analysis tools formulate spatial organization and cell-cell communications as quantifiable properties, and provide algorithms to derive such properties. Integrative pipelines further assemble multiple tools in one package, allowing biologists to conveniently analyze data from beginning to end. In this review, we summarize the state of the art of spatial transcriptomic data analysis methods and pipelines, and discuss how they operate on different technological platforms.

SpaGCN: Integrating gene expression, spatial location and histology to identify spatial domains and spatially variable genes by graph convolutional network.

Recent advances in spatially resolved transcriptomics (SRT) technologies have enabled comprehensive characterization of gene expression patterns in the context of tissue microenvironment. To elucidate spatial gene expression variation, we present SpaGCN, a graph convolutional network approach that integrates gene expression, spatial location and histology in SRT data analysis. Through graph convolution, SpaGCN aggregates gene expression of each spot from its neighboring spots, which enables the identification of spatial domains with coherent expression and histology. The subsequent domain guided differential expression (DE) analysis then detects genes with enriched expression patterns in the identified domains. Analyzing seven SRT datasets using SpaGCN, we show it can detect genes with much more enriched spatial expression patterns than competing methods. Furthermore, genes detected by SpaGCN are transferrable and can be utilized to study spatial variation of gene expression in other datasets. SpaGCN is computationally fast, platform independent, making it a desirable tool for diverse SRT studies.