Voyager: exploratory single-cell genomics data analysis with geospatial statistics

Multimodal spatial proteomic profiling in acute myeloid leukemia

Acute myeloid leukemia (AML) resides in an immune-rich microenvironment, yet, immune-based therapies have faltered in eliciting durable responses. Bridging this paradox requires a comprehensive understanding of leukemic interactions within the bone marrow microenvironment. We optimized a high-throughput tissue-microarray-based pipeline for high-plex spatial immunofluorescence and mass cytometry imaging on a single slide, capturing immune, tumor, and structural components. Using unbiased clustering on the spatial K function, we unveiled the presence of tertiary lymphoid-like aggregates in bone marrow, which we validated using spatial transcriptomics and an independent proteomics approach. We then found validated TLS signatures predictive of outcomes in AML using an integrated public 480-patient transcriptomic dataset. By harnessing high-plex spatial proteomics, we open the possibility of discovering novel structures and interactions that underpin leukemic immune response. Further, our study's methodologies and resources can be adapted for other bone marrow diseases where decalcification and autofluorescence present challenges.

Cell-type deconvolution methods for spatial transcriptomics

Spatial transcriptomics is a powerful method for studying the spatial organization of cells, which is a critical feature in the development, function and evolution of multicellular life. However, sequencing-based spatial transcriptomics has not yet achieved cellular-level resolution, so advanced deconvolution methods are needed to infer cell-type contributions at each location in the data. Recent progress has led to diverse tools for cell-type deconvolution that are helping to describe tissue architectures in health and disease. In this Review, we describe the varied types of cell-type deconvolution methods for spatial transcriptomics, contrast their capabilities and summarize them in a web-based, interactive table to enable more efficient method selection.

Modeling the dynamics of EMT reveals genes associated with pan-cancer intermediate states and plasticity

Epithelial-mesenchymal transition (EMT) is a cell state transition co-opted by cancer that drives metastasis via stable intermediate states. Here we study EMT dynamics to identify marker genes of highly metastatic intermediate cells via mathematical modeling with single-cell RNA sequencing (scRNA-seq) data. Across multiple tumor types and stimuli, we identified genes consistently upregulated in EMT intermediate states, many previously unrecognized as EMT markers. Bayesian parameter inference of a simple EMT mathematical model revealed tumor-specific transition rates, providing a framework to quantify EMT progression. Consensus analysis of differential expression, RNA velocity, and model-derived dynamics highlighted SFN and NRG1 as key regulators of intermediate EMT. Independent validation confirmed SFN as an intermediate state marker. Our approach integrates modeling and inference to identify genes associated with EMT dynamics, offering biomarkers and therapeutic targets to modulate tumor-promoting cell state transitions driven by EMT.

Quantitative characterization of tissue states using multiomics and ecological spatial analysis

The spatial organization of cells in tissues underlies biological function, and recent advances in spatial profiling technologies have enhanced our ability to analyze such arrangements to study biological processes and disease progression. We propose MESA (multiomics and ecological spatial analysis), a framework drawing inspiration from ecological concepts to delineate functional and spatial shifts across tissue states. MESA introduces metrics to systematically quantify spatial diversity and identify hot spots, linking spatial patterns to phenotypic outcomes, including disease progression. Furthermore, MESA integrates spatial and single-cell multiomics data to facilitate an in-depth, molecular understanding of cellular neighborhoods and their spatial interactions within tissue microenvironments. Applying MESA to diverse datasets demonstrates additional insights it brings over prior methods, including newly identified spatial structures and key cell populations linked to disease states. Available as a Python package, MESA offers a versatile framework for quantitative decoding of tissue architectures in spatial omics across health and disease.

A factor-based analysis of individual human microglia uncovers regulators of an Alzheimer-related transcriptional signature

Human microglial heterogeneity is only beginning to be appreciated at the molecular level. Here, we present a large, single-cell atlas of expression signatures from 441,088 live microglia broadly sampled across a diverse set of brain regions and neurodegenerative and neuroinflammatory diseases obtained from 161 donors sampled at autopsy or during a neurosurgical procedure. Using single-cell hierarchical Poisson factorization (scHPF), we derived a 23-factor model for continuous gene expression signatures across microglia which capture specific biological processes (e.g., metabolism, phagocytosis, antigen presentation, inflammatory signaling, disease-associated states). Using external datasets, we evaluated the aspects of microglial phenotypes that are encapsulated in various in vitro and in vivo microglia models and identified and replicated the role of two factors in human postmortem tissue of Alzheimer's disease (AD). Further, we derived a complex network of transcriptional regulators for all factors, including regulators of an AD-related factor enriched for the mouse disease-associated microglia 2 (DAM2) signature: ARID5B, CEBPA, MITF, and PPARG. We replicated the role of these four regulators in the AD-related factor and then designed a multiplexed MERFISH panel to assess our microglial factors using spatial transcriptomics. We find that, unlike cells with high expression of the interferon-response factor, cells with high expression of the AD DAM2-like factor are widely distributed in neocortical tissue. We thus propose a novel analytic framework that provides a taxonomic approach for microglia that is more biologically interpretable and use it to uncover new therapeutic targets for AD.

Spatial-Omics Methods and Applications

Traditional tissue profiling approaches have evolved from bulk studies to single-cell analysis over the last decade; however, the spatial context in tissues and microenvironments has always been lost. Over the last 5 years, spatial technologies have emerged that enabled researchers to investigate tissues in situ for proteins and transcripts without losing anatomy and histology. The field of spatial-omics enables highly multiplexed analysis of biomolecules like RNAs and proteins in their native spatial context-and has matured from initial proof-of-concept studies to a thriving field with widespread applications from basic research to translational and clinical studies. While there has been wide adoption of spatial technologies, there remain challenges with the standardization of methodologies, sample compatibility, throughput, resolution, and ease of use. In this chapter, we discuss the current state of the field and highlight technological advances and limitations.

Dendritic cells type 1 control the formation, maintenance, and function of tertiary lymphoid structures in cancer

Tertiary lymphoid structures (TLS) are organized immune cell aggregates that arise in chronic inflammatory conditions. In cancer, TLS are associated with better prognosis and enhanced response to immunotherapy, making these structures attractive therapeutic targets. However, the mechanisms regulating TLS formation and maintenance in cancer are incompletely understood. Using spatial transcriptomics and multiplex imaging across various human tumors, we found an enrichment of mature dendritic cells (DC) expressing high levels of CCR7 in TLS, prompting us to investigate the role of DC in the formation and maintenance of TLS in solid tumors. To address this, we developed a novel murine model of non-small cell lung cancer (NSCLC) that forms mature TLS, containing B cell follicles with germinal centers and T cell zones with T follicular helper cells (TFH) and TCF1+PD-1+ progenitor exhausted CD8+ T cells (Tpex). Here we show that, during the early stages of tumor development, TLS formation relies on IFNγ-driven maturation of the conventional DC type 1 (cDC1) subset, their migration to tumor-draining lymph nodes (tdLN), and recruitment of activated T cells to the tumor site. As tumors progress, TLS maintenance becomes independent of T cell egress from tdLN, coinciding with a significant reduction of cDC1 migration to tdLN. Instead, mature cDC1 accumulate within intratumoral CCR7 ligand-enriched stromal hubs. Notably, timed depletion of cDC1 or disruption of their migration to these stromal hubs after TLS are formed alters TLS maintenance. Importantly, we found that cDC1-mediated antigen presentation to both CD4+ and CD8+ T cells and intact CD40 signaling, is critical for the maintenance of TLS, the preservation of the TFH cell pool, the formation of germinal center and the production of tumor-specific IgG antibodies. These findings underscore the key role of mature cDC1 in establishing and maintaining functional TLS within tumor lesions and highlight the potential for cDC1-targeting therapies as a promising strategy to enhance TLS function and improve anti-tumor immunity in patients with cancer.

Spatial and single-nucleus transcriptomic analysis of genetic and sporadic forms of Alzheimer's disease

The pathogenesis of Alzheimer's disease (AD) depends on environmental and heritable factors, with its molecular etiology still unclear. Here we present a spatial transcriptomic (ST) and single-nucleus transcriptomic survey of late-onset sporadic AD and AD in Down syndrome (DSAD). Studying DSAD provides an opportunity to enhance our understanding of the AD transcriptome, potentially bridging the gap between genetic mouse models and sporadic AD. We identified transcriptomic changes that may underlie cortical layer-preferential pathology accumulation. Spatial co-expression network analyses revealed transient and regionally restricted disease processes, including a glial inflammatory program dysregulated in upper cortical layers and implicated in AD genetic risk and amyloid-associated processes. Cell-cell communication analysis further contextualized this gene program in dysregulated signaling networks. Finally, we generated ST data from an amyloid AD mouse model to identify cross-species amyloid-proximal transcriptomic changes with conformational context.

SpottedPy quantifies relationships between spatial transcriptomic hotspots and uncovers environmental cues of epithelial-mesenchymal plasticity in breast cancer

Spatial transcriptomics is revolutionizing the exploration of intratissue heterogeneity in cancer, yet capturing cellular niches and their spatial relationships remains challenging. We introduce SpottedPy, a Python package designed to identify tumor hotspots and map spatial interactions within the cancer ecosystem. Using SpottedPy, we examine epithelial-mesenchymal plasticity in breast cancer and highlight stable niches associated with angiogenic and hypoxic regions, shielded by CAFs and macrophages. Hybrid and mesenchymal hotspot distribution follows transformation gradients reflecting progressive immunosuppression. Our method offers flexibility to explore spatial relationships at different scales, from immediate neighbors to broader tissue modules, providing new insights into tumor microenvironment dynamics.

The impact of package selection and versioning on single-cell RNA-seq analysis

Standard single-cell RNA-sequencing analysis (scRNA-seq) workflows consist of converting raw read data into cell-gene count matrices through sequence alignment, followed by analyses including filtering, highly variable gene selection, dimensionality reduction, clustering, and differential expression analysis. Seurat and Scanpy are the most widely-used packages implementing such workflows, and are generally thought to implement individual steps similarly. We investigate in detail the algorithms and methods underlying Seurat and Scanpy and find that there are, in fact, considerable differences in the outputs of Seurat and Scanpy. The extent of differences between the programs is approximately equivalent to the variability that would be introduced in benchmarking scRNA-seq datasets by sequencing less than 5% of the reads or analyzing less than 20% of the cell population. Additionally, distinct versions of Seurat and Scanpy can produce very different results, especially during parts of differential expression analysis. Our analysis highlights the need for users of scRNA-seq to carefully assess the tools on which they rely, and the importance of developers of scientific software to prioritize transparency, consistency, and reproducibility for their tools.

A perspective on FAIR quality control in multiplexed imaging data processing

Multiplexed imaging approaches are getting increasingly adopted for imaging of large tissue areas, yielding big imaging datasets both in terms of the number of samples and the size of image data per sample. The processing and analysis of these datasets is complex owing to frequent technical artifacts and heterogeneous profiles from a high number of stained targets To streamline the analysis of multiplexed images, automated pipelines making use of state-of-the-art algorithms have been developed. In these pipelines, the output quality of one processing step is typically dependent on the output of the previous step and errors from each step, even when they appear minor, can propagate and confound the results. Thus, rigorous quality control (QC) at each of these different steps of the image processing pipeline is of paramount importance both for the proper analysis and interpretation of the analysis results and for ensuring the reusability of the data. Ideally, QC should become an integral and easily retrievable part of the imaging datasets and the analysis process. Yet, limitations of the currently available frameworks make integration of interactive QC difficult for large multiplexed imaging data. Given the increasing size and complexity of multiplexed imaging datasets, we present the different challenges for integrating QC in image analysis pipelines as well as suggest possible solutions that build on top of recent advances in bioimage analysis.

Giotto Suite: a multi-scale and technology-agnostic spatial multi-omics analysis ecosystem

Emerging spatial omics technologies continue to advance the molecular mapping of tissue architecture and the investigation of gene regulation and cellular crosstalk, which in turn provide new mechanistic insights into a wide range of biological processes and diseases. Such technologies provide an increasingly large amount of information content at multiple spatial scales. However, representing and harmonizing diverse spatial datasets efficiently, including combining multiple modalities or spatial scales in a scalable and flexible manner, remains a substantial challenge. Here, we present Giotto Suite, a suite of open-source software packages that underlies a fully modular and integrated spatial data analysis toolbox. At its core, Giotto Suite is centered around an innovative and technology-agnostic data framework embedded in the R software environment, which allows the representation and integration of virtually any type of spatial omics data at any spatial resolution. In addition, Giotto Suite provides both scalable and extensible end-to-end solutions for data analysis, integration, and visualization. Giotto Suite integrates molecular, morphology, spatial, and annotated feature information to create a responsive and flexible workflow for multi-scale, multi-omic data analyses, as demonstrated here by applications to several state-of-the-art spatial technologies. Furthermore, Giotto Suite builds upon interoperable interfaces and data structures that bridge the established fields of genomics and spatial data science, thereby enabling independent developers to create custom-engineered pipelines. As such, Giotto Suite creates an immersive ecosystem for spatial multi-omic data analysis.

Spatial and single-nucleus transcriptomic analysis of genetic and sporadic forms of Alzheimer's Disease

The pathogenesis of Alzheimer's disease (AD) depends on environmental and heritable factors, with remarkable differences evident between individuals at the molecular level. Here we present a transcriptomic survey of AD using spatial transcriptomics (ST) and single-nucleus RNA-seq in cortical samples from early-stage AD, late-stage AD, and AD in Down Syndrome (AD in DS) donors. Studying AD in DS provides an opportunity to enhance our understanding of the AD transcriptome, potentially bridging the gap between genetic mouse models and sporadic AD. Our analysis revealed spatial and cell-type specific changes in disease, with broad similarities in these changes between sAD and AD in DS. We performed additional ST experiments in a disease timecourse of 5xFAD and wildtype mice to facilitate cross-species comparisons. Finally, amyloid plaque and fibril imaging in the same tissue samples used for ST enabled us to directly link changes in gene expression with accumulation and spread of pathology.