Three-dimensional intact-tissue sequencing of single-cell transcriptional states

Computational modeling for deciphering tissue microenvironment heterogeneity from spatially resolved transcriptomics

Spatial transcriptomics techniques, while measuring gene expression, retain spatial location information, aiding in situ studies of organismal tissue architecture and the progression of pathological processes. These techniques generate vast amounts of omics data, necessitating the development of computational methods to reveal the underlying tissue microenvironment heterogeneity. The main directions in spatial transcriptomics data analysis are spatial domain detection and spatial deconvolution, which can identify spatial functional regions and parse the distribution of cell types in spatial transcriptomics data by integrating single-cell transcriptomics data. In these two research directions, many computational methods have been successively proposed. This article will categorize them into three types: machine learning-based methods, probabilistic models-based methods, and deep learning-based methods. It will list and discuss the representative algorithms of each type along with their advantages and disadvantages and describe the datasets and evaluation metrics used to assess these computational methods, facilitating researchers in selecting suitable computational methods according to their research needs. Finally, combining the latest technological developments and the advantages and disadvantages of current algorithms, this article will look forward to the future directions of computational method development.

Computational methods for alignment and integration of spatially resolved transcriptomics data

Most of the complex biological regulatory activities occur in three dimensions (3D). To better analyze biological processes, it is essential not only to decipher the molecular information of numerous cells but also to understand how their spatial contexts influence their behavior. With the development of spatially resolved transcriptomics (SRT) technologies, SRT datasets are being generated to simultaneously characterize gene expression and spatial arrangement information within tissues, organs or organisms. To fully leverage spatial information, the focus extends beyond individual two-dimensional (2D) slices. Two tasks known as slices alignment and data integration have been introduced to establish correlations between multiple slices, enhancing the effectiveness of downstream tasks. Currently, numerous related methods have been developed. In this review, we first elucidate the details and principles behind several representative methods. Then we report the testing results of these methods on various SRT datasets, and assess their performance in representative downstream tasks. Insights into the strengths and weaknesses of each method and the reasons behind their performance are discussed. Finally, we provide an outlook on future developments. The codes and details of experiments are now publicly available at https://github.com/YangLabHKUST/SRT_alignment_and_integration.

Deciphering the spatiotemporal transcriptional landscape of intestinal diseases (Review)

The intestines are the largest barrier organ in the human body. The intestinal barrier plays a crucial role in maintaining the balance of the intestinal environment and protecting the intestines from harmful bacterial invasion. Single‑cell RNA sequencing technology allows the detection of the different cell types in the intestine in two dimensions and the exploration of cell types that have not been fully characterized. The intestinal mucosa is highly complex in structure, and its proper functioning is linked to multiple structures in the proximal‑distal intestinal and luminal‑mucosal axes. Spatial localization is at the core of the efforts to explore the interactions between the complex structures. Spatial transcriptomics (ST) is a method that allows for comprehensive tissue analysis and the acquisition of spatially separated genetic information from individual cells, while preserving their spatial location and interactions. This approach also prevents the loss of fragile cells during tissue disaggregation. The emergence of ST technology allows us to spatially dissect enzymatic processes and interactions between multiple cells, genes, proteins and signals in the intestine. This includes the exchange of oxygen and nutrients in the intestine, different gradients of microbial populations and the role of extracellular matrix proteins. This regionally precise approach to tissue studies is gaining more acceptance and is increasingly applied in the investigation of disease mechanisms related to the gastrointestinal tract. Therefore, this review summarized the application of ST in gastrointestinal diseases.

Spatial technologies to evaluate the HIV-1 reservoir and its microenvironment in the lymph node

The presence of the HIV-1 reservoir, a group of immune cells that contain intact, integrated, and replication-competent proviruses, is a major challenge to cure HIV-1. HIV-1 reservoir cells are largely unaffected by the cytopathic effects of viruses, antiviral immune responses, or antiretroviral therapy (ART). The HIV-1 reservoir is seeded early during HIV-1 infection and augmented during active viral replication. CD4+ T cells are the primary target for HIV-1 infection, and recent studies suggest that memory T follicular helper cells within the lymph node, more precisely in the B cell follicle, harbor integrated provirus, which contribute to viral rebound upon ART discontinuation. The B cell follicle, more specifically the germinal center, possesses a unique environment because of its distinct property of being partly immune privileged, potentially allowing HIV-1-infected cells within the lymph nodes to be protected from CD8+ T cells. This modified immune response in the germinal center of the follicle is potentially explained by the exclusion of CD8+ T cells and the presence of T regulatory cells at the junction of the follicle and extrafollicular region. The proviral makeup of HIV-1-infected cells is similar in lymph nodes and blood, suggesting trafficking between these compartments. Little is known about the cell-to-cell interactions, microenvironment of HIV-1-infected cells in the follicle, and trafficking between the lymph node follicle and other body compartments. Applying a spatiotemporal approach that integrates genomics, transcriptomics, and proteomics to investigate the HIV-1 reservoir and its neighboring cells in the lymph node has promising potential for informing HIV-1 cure efforts.

Single-cell omics: experimental workflow, data analyses and applications

Cells are the fundamental units of biological systems and exhibit unique development trajectories and molecular features. Our exploration of how the genomes orchestrate the formation and maintenance of each cell, and control the cellular phenotypes of various organismsis, is both captivating and intricate. Since the inception of the first single-cell RNA technology, technologies related to single-cell sequencing have experienced rapid advancements in recent years. These technologies have expanded horizontally to include single-cell genome, epigenome, proteome, and metabolome, while vertically, they have progressed to integrate multiple omics data and incorporate additional information such as spatial scRNA-seq and CRISPR screening. Single-cell omics represent a groundbreaking advancement in the biomedical field, offering profound insights into the understanding of complex diseases, including cancers. Here, we comprehensively summarize recent advances in single-cell omics technologies, with a specific focus on the methodology section. This overview aims to guide researchers in selecting appropriate methods for single-cell sequencing and related data analysis.

Unsupervised representation learning of chromatin images identifies changes in cell state and tissue organization in DCIS

Ductal carcinoma in situ (DCIS) is a pre-invasive tumor that can progress to invasive breast cancer, a leading cause of cancer death. We generate a large-scale tissue microarray dataset of chromatin images, from 560 samples from 122 female patients in 3 disease stages and 11 phenotypic categories. Using representation learning on chromatin images alone, without multiplexed staining or high-throughput sequencing, we identify eight morphological cell states and tissue features marking DCIS. All cell states are observed in all disease stages with different proportions, indicating that cell states enriched in invasive cancer exist in small fractions in normal breast tissue. Tissue-level analysis reveals significant changes in the spatial organization of cell states across disease stages, which is predictive of disease stage and phenotypic category. Taken together, we show that chromatin imaging represents a powerful measure of cell state and disease stage of DCIS, providing a simple and effective tumor biomarker.

Integration tools for scRNA-seq data and spatial transcriptomics sequencing data

Numerous methods have been developed to integrate spatial transcriptomics sequencing data with single-cell RNA sequencing (scRNA-seq) data. Continuous development and improvement of these methods offer multiple options for integrating and analyzing scRNA-seq and spatial transcriptomics data based on diverse research inquiries. However, each method has its own advantages, limitations and scope of application. Researchers need to select the most suitable method for their research purposes based on the actual situation. This review article presents a compilation of 19 integration methods sourced from a wide range of available approaches, serving as a comprehensive reference for researchers to select the suitable integration method for their specific research inquiries. By understanding the principles of these methods, we can identify their similarities and differences, comprehend their applicability and potential complementarity, and lay the foundation for future method development and understanding. This review article presents 19 methods that aim to integrate scRNA-seq data and spatial transcriptomics data. The methods are classified into two main groups and described accordingly. The article also emphasizes the incorporation of High Variance Genes in annotating various technologies, aiming to obtain biologically relevant information aligned with the intended purpose.

TWCOM: an R package for inference of cell-cell communication on spatially resolved transcriptomics data

Summary

The inference of cell-cell communication is important, as it unveils the intricate cellular behaviors at the molecular level, providing crucial insights essential for understanding complex biological processes and informing targeted interventions in various pathological contexts. Here, we present TWCOM, an R package that implements a Tweedie distribution-based model for accurate cell-cell communication inference. Operating under a generalized additive model framework, TWCOM adeptly handles both single-cell resolution and spot-based spatially resolved transcriptomics data, providing a versatile tool for robust biological sample analysis.

Availability and implementation

The R package TWCOM is available at https://github.com/dongyuanwu/TWCOM. Comprehensive documentation is included with the package.

Categorization of 31 computational methods to detect spatially variable genes from spatially resolved transcriptomics data

In the analysis of spatially resolved transcriptomics data, detecting spatially variable genes (SVGs) is crucial. Numerous computational methods exist, but varying SVG definitions and methodologies lead to incomparable results. We review 31 state-of-the-art methods, categorizing SVGs into three types: overall, cell-type-specific, and spatial-domain-marker SVGs. Our review explains the intuitions underlying these methods, summarizes their applications, and categorizes the hypothesis tests they use in the trade-off between generality and specificity for SVG detection. We discuss challenges in SVG detection and propose future directions for improvement. Our review offers insights for method developers and users, advocating for category-specific benchmarking.

Automated single-cell proteomics providing sufficient proteome depth to study complex biology beyond cell type classifications

The recent technological and computational advances in mass spectrometry-based single-cell proteomics have pushed the boundaries of sensitivity and throughput. However, reproducible quantification of thousands of proteins within a single cell remains challenging. To address some of those limitations, we present a dedicated sample preparation chip, the proteoCHIP EVO 96 that directly interfaces with the Evosep One. This, in combination with the Bruker timsTOF demonstrates double the identifications without manual sample handling and the newest generation timsTOF Ultra identifies up to 4000 with an average of 3500 protein groups per single HEK-293T without a carrier or match-between runs. Our workflow spans 4 orders of magnitude, identifies over 50 E3 ubiquitin-protein ligases, and profiles key regulatory proteins upon small molecule stimulation. This study demonstrates that the proteoCHIP EVO 96-based sample preparation with the timsTOF Ultra provides sufficient proteome depth to study complex biology beyond cell-type classifications.

A point cloud segmentation framework for image-based spatial transcriptomics

Recent progress in image-based spatial RNA profiling enables to spatially resolve tens to hundreds of distinct RNA species with high spatial resolution. It presents new avenues for comprehending tissue organization. In this context, the ability to assign detected RNA transcripts to individual cells is crucial for downstream analyses, such as in-situ cell type calling. Yet, accurate cell segmentation can be challenging in tissue data, in particular in the absence of a high-quality membrane marker. To address this issue, we introduce ComSeg, a segmentation algorithm that operates directly on single RNA positions and that does not come with implicit or explicit priors on cell shape. ComSeg is applicable in complex tissues with arbitrary cell shapes. Through comprehensive evaluations on simulated and experimental datasets, we show that ComSeg outperforms existing state-of-the-art methods for in-situ single-cell RNA profiling and in-situ cell type calling. ComSeg is available as a documented and open source pip package at https://github.com/fish-quant/ComSeg .

Spatial Dissection of the Immune Landscape of Solid Tumors to Advance Precision Medicine

Chemotherapeutics, radiation, targeted therapeutics, and immunotherapeutics each demonstrate clinical benefits for a small subset of patients with solid malignancies. Immune cells infiltrating the tumor and the surrounding stroma play a critical role in shaping cancer progression and modulating therapy response. They do this by interacting with the other cellular and molecular components of the tumor microenvironment. Spatial multi-omics technologies are rapidly evolving. Currently, such technologies allow high-throughput RNA and protein profiling and retain geographical information about the tumor microenvironment cellular architecture and the functional phenotype of tumor, immune, and stromal cells. An in-depth spatial characterization of the heterogeneous tumor immune landscape can improve not only the prognosis but also the prediction of therapy response, directing cancer patients to more tailored and efficacious treatments. This review highlights recent advancements in spatial transcriptomics and proteomics profiling technologies and the ways these technologies are being applied for the dissection of the immune cell composition in solid malignancies in order to further both basic research in oncology and the implementation of precision treatments in the clinic.

An interpretable Bayesian clustering approach with feature selection for analyzing spatially resolved transcriptomics data

Recent breakthroughs in spatially resolved transcriptomics (SRT) technologies have enabled comprehensive molecular characterization at the spot or cellular level while preserving spatial information. Cells are the fundamental building blocks of tissues, organized into distinct yet connected components. Although many non-spatial and spatial clustering approaches have been used to partition the entire region into mutually exclusive spatial domains based on the SRT high-dimensional molecular profile, most require an ad hoc selection of less interpretable dimensional-reduction techniques. To overcome this challenge, we propose a zero-inflated negative binomial mixture model to cluster spots or cells based on their molecular profiles. To increase interpretability, we employ a feature selection mechanism to provide a low-dimensional summary of the SRT molecular profile in terms of discriminating genes that shed light on the clustering result. We further incorporate the SRT geospatial profile via a Markov random field prior. We demonstrate how this joint modeling strategy improves clustering accuracy, compared with alternative state-of-the-art approaches, through simulation studies and 3 real data applications.

Multiplexed in situ RNA imaging by combFISH

Multiplexed in situ RNA imaging offers new opportunities for gene expression profiling by providing high-throughput spatial information. In this work, we present a cyclic combinatorial fluorescent in situ hybridization (combFISH) assay to achieve multiplexed detection of RNA in cell cultures and tissues. Specifically, multiplexing is achieved through cyclic interrogation of barcode sequences on the rolling circle amplicons generated from the padlock probe assay by using sets of combinatorial detection probes. Theoretically, combFISH can detect 64 genes in three hybridization cycles by combinatorial barcoding using 12 fluorescently labeled detection probes. Our method eliminates sequencing-by-ligation (SBL) chemistry in the in situ sequencing protocol and directly uses RNA as targets for ligation, making it more straightforward. We showed that our method works in fresh-frozen and formalin-fixed paraffin-embedded tissue sections. With its straightforward protocols, we expect our method to be adopted by the scientific community and extended to clinical settings.

Accurate and efficient integrative reference-informed spatial domain detection for spatial transcriptomics

Spatially resolved transcriptomics (SRT) studies are becoming increasingly common and large, offering unprecedented opportunities in mapping complex tissue structures and functions. Here we present integrative and reference-informed tissue segmentation (IRIS), a computational method designed to characterize tissue spatial organization in SRT studies through accurately and efficiently detecting spatial domains. IRIS uniquely leverages single-cell RNA sequencing data for reference-informed detection of biologically interpretable spatial domains, integrating multiple SRT slices while explicitly considering correlations both within and across slices. We demonstrate the advantages of IRIS through in-depth analysis of six SRT datasets encompassing diverse technologies, tissues, species and resolutions. In these applications, IRIS achieves substantial accuracy gains (39-1,083%) and speed improvements (4.6-666.0) in moderate-sized datasets, while representing the only method applicable for large datasets including Stereo-seq and 10x Xenium. As a result, IRIS reveals intricate brain structures, uncovers tumor microenvironment heterogeneity and detects structural changes in diabetes-affected testis, all with exceptional speed and accuracy.

New tools to study renal fibrogenesis

PURPOSE OF REVIEW

Kidney fibrosis is a key pathological aspect and outcome of chronic kidney disease (CKD). The advent of multiomic analyses using human kidney tissue, enabled by technological advances, marks a new chapter of discovery in fibrosis research of the kidney. This review highlights the rapid advancements of single-cell and spatial multiomic techniques that offer new avenues for exploring research questions related to human kidney fibrosis development.

RECENT FINDINGS

We recently focused on understanding the origin and transition of myofibroblasts in kidney fibrosis using single-cell RNA sequencing (scRNA-seq) [1] . We analysed cells from healthy human kidneys and compared them to patient samples with CKD. We identified PDGFRα+/PDGFRβ+ mesenchymal cells as the primary cellular source of extracellular matrix (ECM) in human kidney fibrosis. We found several commonly shared cell states of fibroblasts and myofibroblasts and provided insights into molecular regulators. Novel single-cell and spatial multiomics tools are now available to shed light on cell lineages, the plasticity of kidney cells and cell-cell communication in fibrosis.

SUMMARY

As further single-cell and spatial multiomic approaches are being developed, opportunities to apply these methods to human kidney tissues expand similarly. Careful design and optimisation of the multiomic experiments are needed to answer questions related to cell lineages, plasticity and cell-cell communication in kidney fibrosis.

Foundation models in molecular biology

Determining correlations between molecules at various levels is an important topic in molecular biology. Large language models have demonstrated a remarkable ability to capture correlations from large amounts of data in the field of natural language processing as well as image generation, and correlations captured from data using large language models can also be applicable to solving a wide range of specific tasks, hence large language models are also referred to as foundation models. The massive amount of data that exists in the field of molecular biology provides an excellent basis for the development of foundation models, and the recent emergence of foundation models in the field of molecular biology has really pushed the entire field forward. We summarize the foundation models developed based on RNA sequence data, DNA sequence data, protein sequence data, single-cell transcriptome data, and spatial transcriptome data respectively, and further discuss the research directions for the development of foundation models in molecular biology.

SPRITE: improving spatial gene expression imputation with gene and cell networks

MOTIVATION

Spatially resolved single-cell transcriptomics have provided unprecedented insights into gene expression in situ, particularly in the context of cell interactions or organization of tissues. However, current technologies for profiling spatial gene expression at single-cell resolution are generally limited to the measurement of a small number of genes. To address this limitation, several algorithms have been developed to impute or predict the expression of additional genes that were not present in the measured gene panel. Current algorithms do not leverage the rich spatial and gene relational information in spatial transcriptomics. To improve spatial gene expression predictions, we introduce Spatial Propagation and Reinforcement of Imputed Transcript Expression (SPRITE) as a meta-algorithm that processes predictions obtained from existing methods by propagating information across gene correlation networks and spatial neighborhood graphs.

RESULTS

SPRITE improves spatial gene expression predictions across multiple spatial transcriptomics datasets. Furthermore, SPRITE predicted spatial gene expression leads to improved clustering, visualization, and classification of cells. SPRITE can be used in spatial transcriptomics data analysis to improve inferences based on predicted gene expression.

AVAILABILITY AND IMPLEMENTATION

The SPRITE software package is available at https://github.com/sunericd/SPRITE. Code for generating experiments and analyses in the manuscript is available at https://github.com/sunericd/sprite-figures-and-analyses.

Single-cell analysis in rheumatic and allergic diseases: insights for clinical practice

Since the advent of single-cell RNA sequencing (scRNA-seq) methodology, single-cell analysis has become a powerful tool for exploration of cellular networks and dysregulated immune responses in disease pathogenesis. Advanced bioinformatics tools have enabled the combined analysis of scRNA-seq data and information on various cell properties, such as cell surface molecular profiles, chromatin accessibility and spatial information, leading to a deeper understanding of pathology. This Review provides an overview of the achievements in single-cell analysis applied to clinical samples of rheumatic and allergic diseases, including rheumatoid arthritis, systemic lupus erythematosus, systemic sclerosis, allergic airway diseases and atopic dermatitis, with an expanded scope beyond peripheral blood cells to include local diseased tissues. Despite the valuable insights that single-cell analysis has provided into disease pathogenesis, challenges remain in translating single-cell findings into clinical practice and developing personalized treatment strategies. Beyond understanding the atlas of cellular diversity, we discuss the application of data obtained in each study to clinical practice, with a focus on identifying biomarkers and therapeutic targets.

Protocol for enhancing visualization clarity for categorical spatial datasets using Spaco

In categorical data visualization, appropriate color arrangements can avoid perceptual ambiguity and help perceive underlying data patterns. We introduce a protocol to assign contrastive colors to neighboring categories using both Python and R packages. We describe steps for calculating the interlacement between clusters and generating a proper color palette and calculating color contrast. We then detail procedures for aligning cluster interlacement and color contrast to get an optimized cluster-color assignment, achieving clear categorical visualization. For complete details on the use and execution of this protocol, please refer to Jing et al.1.

STimage-1K4M: A histopathology image-gene expression dataset for spatial transcriptomics

Recent advances in multi-modal algorithms have driven and been driven by the increasing availability of large image-text datasets, leading to significant strides in various fields, including computational pathology. However, in most existing medical image-text datasets, the text typically provides high-level summaries that may not sufficiently describe sub-tile regions within a large pathology image. For example, an image might cover an extensive tissue area containing cancerous and healthy regions, but the accompanying text might only specify that this image is a cancer slide, lacking the nuanced details needed for in-depth analysis. In this study, we introduce STimage-1K4M, a novel dataset designed to bridge this gap by providing genomic features for sub-tile images. STimage-1K4M contains 1,149 images derived from spatial transcriptomics data, which captures gene expression information at the level of individual spatial spots within a pathology image. Specifically, each image in the dataset is broken down into smaller sub-image tiles, with each tile paired with 15,000 - 30,000 dimensional gene expressions. With 4,293,195 pairs of sub-tile images and gene expressions, STimage-1K4M offers unprecedented granularity, paving the way for a wide range of advanced research in multi-modal data analysis an innovative applications in computational pathology, and beyond.

Advances in spatial transcriptomics and its applications in cancer research

Malignant tumors have increasing morbidity and high mortality, and their occurrence and development is a complicate process. The development of sequencing technologies enabled us to gain a better understanding of the underlying genetic and molecular mechanisms in tumors. In recent years, the spatial transcriptomics sequencing technologies have been developed rapidly and allow the quantification and illustration of gene expression in the spatial context of tissues. Compared with the traditional transcriptomics technologies, spatial transcriptomics technologies not only detect gene expression levels in cells, but also inform the spatial location of genes within tissues, cell composition of biological tissues, and interaction between cells. Here we summarize the development of spatial transcriptomics technologies, spatial transcriptomics tools and its application in cancer research. We also discuss the limitations and challenges of current spatial transcriptomics approaches, as well as future development and prospects.

Tissue module discovery in single-cell-resolution spatial transcriptomics data via cell-cell interaction-aware cell embedding

Computational methods are desired for single-cell-resolution spatial transcriptomics (ST) data analysis to uncover spatial organization principles for how individual cells exert tissue-specific functions. Here, we present ST data analysis via interaction-aware cell embedding (SPACE), a deep-learning method for cell-type identification and tissue module discovery from single-cell-resolution ST data by learning a cell representation that captures its gene expression profile and interactions with its spatial neighbors. SPACE identified spatially informed cell subtypes defined by their special spatial distribution patterns and distinct proximal-interacting cell types. SPACE also automatically discovered "cell communities"-tissue modules with discernible boundaries and a uniform spatial distribution of constituent cell types. For each cell community, SPACE outputs a characteristic proximal cell-cell interaction network associated with physiological processes, which can be used to refine ligand-receptor-based intercellular signaling analyses. We envision that SPACE can be used in large-scale ST projects to understand how proximal cell-cell interactions contribute to emergent biological functions within cell communities. A record of this paper's transparent peer review process is included in the supplemental information.

Computational analysis of super-resolved in situ sequencing data reveals genes modified by immune-tumor contact events

Cancer cells can manipulate immune cells and escape from the immune system response. Quantifying the molecular changes that occur when an immune cell touches a tumor cell can increase our understanding of the underlying mechanisms. Recently, it became possible to perform such measurements in situ-for example, using expansion sequencing, which enabled in situ sequencing of genes with super-resolution. We systematically examined whether individual immune cells from specific cell types express genes differently when in physical proximity to individual tumor cells. First, we demonstrated that a dense mapping of genes in situ can be used for the segmentation of cell bodies in 3D, thus improving our ability to detect likely touching cells. Next, we used three different computational approaches to detect the molecular changes that are triggered by proximity: differential expression analysis, tree-based machine learning classifiers, and matrix factorization analysis. This systematic analysis revealed tens of genes, in specific cell types, whose expression separates immune cells that are proximal to tumor cells from those that are not proximal, with a significant overlap between the different detection methods. Remarkably, an order of magnitude more genes are triggered by proximity to tumor cells in CD8 T cells compared to CD4 T cells, in line with the ability of CD8 T cells to directly bind major histocompatibility complex (MHC) class I on tumor cells. Thus, in situ sequencing of an individual biopsy can be used to detect genes likely involved in immune-tumor cell-cell interactions. The data used in this manuscript and the code of the InSituSeg, machine learning, cNMF, and Moran's I methods are publicly available at doi:10.5281/zenodo.7497981.

Dissecting tumor microenvironment from spatially resolved transcriptomics data by heterogeneous graph learning

Spatially resolved transcriptomics (SRT) has enabled precise dissection of tumor-microenvironment (TME) by analyzing its intracellular molecular networks and intercellular cell-cell communication (CCC). However, lacking computational exploration of complicated relations between cells, genes, and histological regions, severely limits the ability to interpret the complex structure of TME. Here, we introduce stKeep, a heterogeneous graph (HG) learning method that integrates multimodality and gene-gene interactions, in unraveling TME from SRT data. stKeep leverages HG to learn both cell-modules and gene-modules by incorporating features of diverse nodes including genes, cells, and histological regions, allows for identifying finer cell-states within TME and cell-state-specific gene-gene relations, respectively. Furthermore, stKeep employs HG to infer CCC for each cell, while ensuring that learned CCC patterns are comparable across different cell-states through contrastive learning. In various cancer samples, stKeep outperforms other tools in dissecting TME such as detecting bi-potent basal populations, neoplastic myoepithelial cells, and metastatic cells distributed within the tumor or leading-edge regions. Notably, stKeep identifies key transcription factors, ligands, and receptors relevant to disease progression, which are further validated by the functional and survival analysis of independent clinical data, thereby highlighting its clinical prognostic and immunotherapy applications.

Simulating multiple variability in spatially resolved transcriptomics with scCube

A pressing challenge in spatially resolved transcriptomics (SRT) is to benchmark the computational methods. A widely-used approach involves utilizing simulated data. However, biases exist in terms of the currently available simulated SRT data, which seriously affects the accuracy of method evaluation and validation. Herein, we present scCube ( https://github.com/ZJUFanLab/scCube ), a Python package for independent, reproducible, and technology-diverse simulation of SRT data. scCube not only enables the preservation of spatial expression patterns of genes in reference-based simulations, but also generates simulated data with different spatial variability (covering the spatial pattern type, the resolution, the spot arrangement, the targeted gene type, and the tissue slice dimension, etc.) in reference-free simulations. We comprehensively benchmark scCube with existing single-cell or SRT simulators, and demonstrate the utility of scCube in benchmarking spot deconvolution, gene imputation, and resolution enhancement methods in detail through three applications.

Complete spatially resolved gene expression is not necessary for identifying spatial domains

Spatially resolved transcriptomics (SRT) technologies have revolutionized the study of tissue organization. We introduce a graph convolutional network with an attention and positive emphasis mechanism, termed BINARY, relying exclusively on binarized SRT data to accurately delineate spatial domains. BINARY outperforms existing methods across various SRT data types while using significantly less input information. Our study suggests that precise gene expression quantification may not always be essential, inspiring further exploration of the broader applications of spatially resolved binarized gene expression data.

Gene count normalization in single-cell imaging-based spatially resolved transcriptomics

BACKGROUND

Recent advances in imaging-based spatially resolved transcriptomics (im-SRT) technologies now enable high-throughput profiling of targeted genes and their locations in fixed tissues. Normalization of gene expression data is often needed to account for technical factors that may confound underlying biological signals.

RESULTS

Here, we investigate the potential impact of different gene count normalization methods with different targeted gene panels in the analysis and interpretation of im-SRT data. Using different simulated gene panels that overrepresent genes expressed in specific tissue regions or cell types, we demonstrate how normalization methods based on detected gene counts per cell differentially impact normalized gene expression magnitudes in a region- or cell type-specific manner. We show that these normalization-induced effects may reduce the reliability of downstream analyses including differential gene expression, gene fold change, and spatially variable gene analysis, introducing false positive and false negative results when compared to results obtained from gene panels that are more representative of the gene expression of the tissue's component cell types. These effects are not observed with normalization approaches that do not use detected gene counts for gene expression magnitude adjustment, such as with cell volume or cell area normalization.

CONCLUSIONS

We recommend using non-gene count-based normalization approaches when feasible and evaluating gene panel representativeness before using gene count-based normalization methods if necessary. Overall, we caution that the choice of normalization method and gene panel may impact the biological interpretation of the im-SRT data.

Accurate Identification of Spatial Domain by Incorporating Global Spatial Proximity and Local Expression Proximity

Accurate identification of spatial domains is essential in the analysis of spatial transcriptomics data in order to elucidate tissue microenvironments and biological functions. However, existing methods only perform domain segmentation based on local or global spatial relationships between spots, resulting in an underutilization of spatial information. To this end, we propose SECE, a deep learning-based method that captures both local and global relationships among spots and aggregates their information using expression similarity and spatial similarity. We benchmarked SECE against eight state-of-the-art methods on six real spatial transcriptomics datasets spanning four different platforms. SECE consistently outperformed other methods in spatial domain identification accuracy. Moreover, SECE produced spatial embeddings that exhibited clearer patterns in low-dimensional visualizations and facilitated a more accurate trajectory inference.

Developmental origin of oligodendrocytes determines their function in the adult brain

In the mouse embryonic forebrain, developmentally distinct oligodendrocyte progenitor cell populations and their progeny, oligodendrocytes, emerge from three distinct regions in a spatiotemporal gradient from ventral to dorsal. However, the functional importance of this oligodendrocyte developmental heterogeneity is unknown. Using a genetic strategy to ablate dorsally derived oligodendrocyte lineage cells (OLCs), we show here that the areas in which dorsally derived OLCs normally reside in the adult central nervous system become populated and myelinated by OLCs of ventral origin. These ectopic oligodendrocytes (eOLs) have a distinctive gene expression profile as well as subtle myelination abnormalities. The failure of eOLs to fully assume the role of the original dorsally derived cells results in locomotor and cognitive deficits in the adult animal. This study reveals the importance of developmental heterogeneity within the oligodendrocyte lineage and its importance for homeostatic brain function.

iIMPACT: integrating image and molecular profiles for spatial transcriptomics analysis

Current clustering analysis of spatial transcriptomics data primarily relies on molecular information and fails to fully exploit the morphological features present in histology images, leading to compromised accuracy and interpretability. To overcome these limitations, we have developed a multi-stage statistical method called iIMPACT. It identifies and defines histology-based spatial domains based on AI-reconstructed histology images and spatial context of gene expression measurements, and detects domain-specific differentially expressed genes. Through multiple case studies, we demonstrate iIMPACT outperforms existing methods in accuracy and interpretability and provides insights into the cellular spatial organization and landscape of functional genes within spatial transcriptomics data.

Large-scale neurophysiology and single-cell profiling in human neuroscience

Advances in large-scale single-unit human neurophysiology, single-cell RNA sequencing, spatial transcriptomics and long-term ex vivo tissue culture of surgically resected human brain tissue have provided an unprecedented opportunity to study human neuroscience. In this Perspective, we describe the development of these paradigms, including Neuropixels and recent brain-cell atlas efforts, and discuss how their convergence will further investigations into the cellular underpinnings of network-level activity in the human brain. Specifically, we introduce a workflow in which functionally mapped samples of human brain tissue resected during awake brain surgery can be cultured ex vivo for multi-modal cellular and functional profiling. We then explore how advances in human neuroscience will affect clinical practice, and conclude by discussing societal and ethical implications to consider. Potential findings from the field of human neuroscience will be vast, ranging from insights into human neurodiversity and evolution to providing cell-type-specific access to study and manipulate diseased circuits in pathology. This Perspective aims to provide a unifying framework for the field of human neuroscience as we welcome an exciting era for understanding the functional cytoarchitecture of the human brain.

ST-CellSeg: Cell segmentation for imaging-based spatial transcriptomics using multi-scale manifold learning

Spatial transcriptomics has gained popularity over the past decade due to its ability to evaluate transcriptome data while preserving spatial information. Cell segmentation is a crucial step in spatial transcriptomic analysis, as it enables the avoidance of unpredictable tissue disentanglement steps. Although high-quality cell segmentation algorithms can aid in the extraction of valuable data, traditional methods are frequently non-spatial, do not account for spatial information efficiently, and perform poorly when confronted with the problem of spatial transcriptome cell segmentation with varying shapes. In this study, we propose ST-CellSeg, an image-based machine learning method for spatial transcriptomics that uses manifold for cell segmentation and is novel in its consideration of multi-scale information. We first construct a fully connected graph which acts as a spatial transcriptomic manifold. Using multi-scale data, we then determine the low-dimensional spatial probability distribution representation for cell segmentation. Using the adjusted Rand index (ARI), normalized mutual information (NMI), and Silhouette coefficient (SC) as model performance measures, the proposed algorithm significantly outperforms baseline models in selected datasets and is efficient in computational complexity.

Gene representation bias in spatial transcriptomics

For sequencing-based spatial transcriptomics data, the gene-spot count matrix is highly sparse. This feature is similar to scRNA-seq. The goal of this paper is to identify whether there exist genes that are frequently under-detected in Visium compared to bulk RNA-seq, and the underlying potential mechanism of under-detection in Visium. We collected paired Visium and bulk RNA-seq data for 28 human samples and 19 mouse samples, which covered diverse tissue sources. We compared the two data types and observed that there indeed exists a collection of genes frequently under-detected in Visium compared to bulk RNA-seq. We performed a motif search to examine the last 350 bp of the frequently under-detected genes, and we observed that the poly (T) motif was significantly enriched in genes identified from both human and mouse data, which matches with our previous finding about frequently under-detected genes in scRNA-seq. We hypothesized that the poly (T) motif may be able to form a hairpin structure with the poly (A) tails of their mRNA transcripts, making it difficult for their mRNA transcripts to be captured during Visium library preparation.

A practical guide to spatial transcriptomics

Spatial transcriptomics is revolutionizing modern biology, offering researchers an unprecedented ability to unravel intricate gene expression patterns within tissues. From pioneering techniques to newly commercialized platforms, the field of spatial transcriptomics has evolved rapidly, ushering in a new era of understanding across various disciplines, from developmental biology to disease research. This dynamic expansion is reflected in the rapidly growing number of technologies and data analysis techniques developed and introduced. However, the expanding landscape presents a considerable challenge for researchers, especially newcomers to the field, as staying informed about these advancements becomes increasingly complex. To address this challenge, we have prepared an updated review with a particular focus on technologies that have reached commercialization and are, therefore, accessible to a broad spectrum of potential new users. In this review, we present the fundamental principles of spatial transcriptomic methods, discuss the challenges in data analysis, provide insights into experimental considerations, offer information about available resources for spatial transcriptomics, and conclude with a guide for method selection and a forward-looking perspective. Our aim is to serve as a guiding resource for both experienced users and newcomers navigating the complex realm of spatial transcriptomics in this era of rapid development. We intend to equip researchers with the necessary knowledge to make informed decisions and contribute to the cutting-edge research that spatial transcriptomics offers.

Dissecting Spatiotemporal Structures in Spatial Transcriptomics via Diffusion-Based Adversarial Learning

Recent advancements in spatial transcriptomics (ST) technologies offer unprecedented opportunities to unveil the spatial heterogeneity of gene expression and cell states within tissues. Despite these capabilities of the ST data, accurately dissecting spatiotemporal structures (e.g., spatial domains, temporal trajectories, and functional interactions) remains challenging. Here, we introduce a computational framework, PearlST (partial differential equation [PDE]-enhanced adversarial graph autoencoder of ST), for accurate inference of spatiotemporal structures from the ST data using PDE-enhanced adversarial graph autoencoder. PearlST employs contrastive learning to extract histological image features, integrates a PDE-based diffusion model to enhance characterization of spatial features at domain boundaries, and learns the latent low-dimensional embeddings via Wasserstein adversarial regularized graph autoencoders. Comparative analyses across multiple ST datasets with varying resolutions demonstrate that PearlST outperforms existing methods in spatial clustering, trajectory inference, and pseudotime analysis. Furthermore, PearlST elucidates functional regulations of the latent features by linking intercellular ligand-receptor interactions to most contributing genes of the low-dimensional embeddings, as illustrated in a human breast cancer dataset. Overall, PearlST proves to be a powerful tool for extracting interpretable latent features and dissecting intricate spatiotemporal structures in ST data across various biological contexts.

SpatialPrompt: spatially aware scalable and accurate tool for spot deconvolution and domain identification in spatial transcriptomics

Efficiently mapping of cell types in situ remains a major challenge in spatial transcriptomics. Most spot deconvolution tools ignore spatial coordinate information and perform extremely slow on large datasets. Here, we introduce SpatialPrompt, a spatially aware and scalable tool for spot deconvolution and domain identification. SpatialPrompt integrates gene expression, spatial location, and single-cell RNA sequencing (scRNA-seq) dataset as reference to accurately infer cell-type proportions of spatial spots. SpatialPrompt uses non-negative ridge regression and graph neural network to efficiently capture local microenvironment information. Our extensive benchmarking analysis on Visium, Slide-seq, and MERFISH datasets demonstrated superior performance of SpatialPrompt over 15 existing tools. On mouse hippocampus dataset, SpatialPrompt achieves spot deconvolution and domain identification within 2 minutes for 50,000 spots. Overall, domain identification using SpatialPrompt was 44 to 150 times faster than existing methods. We build a database housing 40 plus curated scRNA-seq datasets for seamless integration with SpatialPrompt for spot deconvolution.

Spotless, a reproducible pipeline for benchmarking cell type deconvolution in spatial transcriptomics

Spatial transcriptomics (ST) technologies allow the profiling of the transcriptome of cells while keeping their spatial context. Since most commercial untargeted ST technologies do not yet operate at single-cell resolution, computational methods such as deconvolution are often used to infer the cell type composition of each sequenced spot. We benchmarked 11 deconvolution methods using 63 silver standards, 3 gold standards, and 2 case studies on liver and melanoma tissues. We developed a simulation engine called synthspot to generate silver standards from single-cell RNA-sequencing data, while gold standards are generated by pooling single cells from targeted ST data. We evaluated methods based on their performance, stability across different reference datasets, and scalability. We found that cell2location and RCTD are the top-performing methods, but surprisingly, a simple regression model outperforms almost half of the dedicated spatial deconvolution methods. Furthermore, we observe that the performance of all methods significantly decreased in datasets with highly abundant or rare cell types. Our results are reproducible in a Nextflow pipeline, which also allows users to generate synthetic data, run deconvolution methods and optionally benchmark them on their dataset (https://github.com/saeyslab/spotless-benchmark).

A complementary approach for neocortical cytoarchitecture inspection with cellular resolution imaging at whole brain scale

Cytoarchitecture, the organization of cells within organs and tissues, serves as a crucial anatomical foundation for the delineation of various regions. It enables the segmentation of the cortex into distinct areas with unique structural and functional characteristics. While traditional 2D atlases have focused on cytoarchitectonic mapping of cortical regions through individual sections, the intricate cortical gyri and sulci demands a 3D perspective for unambiguous interpretation. In this study, we employed fluorescent micro-optical sectioning tomography to acquire architectural datasets of the entire macaque brain at a resolution of 0.65 μm × 0.65 μm × 3 μm. With these volumetric data, the cortical laminar textures were remarkably presented in appropriate view planes. Additionally, we established a stereo coordinate system to represent the cytoarchitectonic information as surface-based tomograms. Utilizing these cytoarchitectonic features, we were able to three-dimensionally parcel the macaque cortex into multiple regions exhibiting contrasting architectural patterns. The whole-brain analysis was also conducted on mice that clearly revealed the presence of barrel cortex and reflected biological reasonability of this method. Leveraging these high-resolution continuous datasets, our method offers a robust tool for exploring the organizational logic and pathological mechanisms of the brain's 3D anatomical structure.

Benchmarking mapping algorithms for cell-type annotating in mouse brain by integrating single-nucleus RNA-seq and Stereo-seq data

Limited gene capture efficiency and spot size of spatial transcriptome (ST) data pose significant challenges in cell-type characterization. The heterogeneity and complexity of cell composition in the mammalian brain make it more challenging to accurately annotate ST data from brain. Many algorithms attempt to characterize subtypes of neuron by integrating ST data with single-nucleus RNA sequencing (snRNA-seq) or single-cell RNA sequencing. However, assessing the accuracy of these algorithms on Stereo-seq ST data remains unresolved. Here, we benchmarked 9 mapping algorithms using 10 ST datasets from four mouse brain regions in two different resolutions and 24 pseudo-ST datasets from snRNA-seq. Both actual ST data and pseudo-ST data were mapped using snRNA-seq datasets from the corresponding brain regions as reference data. After comparing the performance across different areas and resolutions of the mouse brain, we have reached the conclusion that both robust cell-type decomposition and SpatialDWLS demonstrated superior robustness and accuracy in cell-type annotation. Testing with publicly available snRNA-seq data from another sequencing platform in the cortex region further validated our conclusions. Altogether, we developed a workflow for assessing suitability of mapping algorithm that fits for ST datasets, which can improve the efficiency and accuracy of spatial data annotation.

MNMST: topology of cell networks leverages identification of spatial domains from spatial transcriptomics data

Advances in spatial transcriptomics provide an unprecedented opportunity to reveal the structure and function of biology systems. However, current algorithms fail to address the heterogeneity and interpretability of spatial transcriptomics data. Here, we present a multi-layer network model for identifying spatial domains in spatial transcriptomics data with joint learning. We demonstrate that spatial domains can be precisely characterized and discriminated by the topological structure of cell networks, facilitating identification and interpretability of spatial domains, which outperforms state-of-the-art baselines. Furthermore, we prove that network model offers an effective and efficient strategy for integrative analysis of spatial transcriptomics data from various platforms.

Multi-modal domain adaptation for revealing spatial functional landscape from spatially resolved transcriptomics

Spatially resolved transcriptomics (SRT) has emerged as a powerful tool for investigating gene expression in spatial contexts, providing insights into the molecular mechanisms underlying organ development and disease pathology. However, the expression sparsity poses a computational challenge to integrate other modalities (e.g. histological images and spatial locations) that are simultaneously captured in SRT datasets for spatial clustering and variation analyses. In this study, to meet such a challenge, we propose multi-modal domain adaption for spatial transcriptomics (stMDA), a novel multi-modal unsupervised domain adaptation method, which integrates gene expression and other modalities to reveal the spatial functional landscape. Specifically, stMDA first learns the modality-specific representations from spatial multi-modal data using multiple neural network architectures and then aligns the spatial distributions across modal representations to integrate these multi-modal representations, thus facilitating the integration of global and spatially local information and improving the consistency of clustering assignments. Our results demonstrate that stMDA outperforms existing methods in identifying spatial domains across diverse platforms and species. Furthermore, stMDA excels in identifying spatially variable genes with high prognostic potential in cancer tissues. In conclusion, stMDA as a new tool of multi-modal data integration provides a powerful and flexible framework for analyzing SRT datasets, thereby advancing our understanding of intricate biological systems.

Spatial Multi-Omics in Alzheimer's Disease: A Multi-Dimensional Approach to Understanding Pathology and Progression

Alzheimer's Disease (AD) presents a complex neuropathological landscape characterized by hallmark amyloid plaques and neurofibrillary tangles, leading to progressive cognitive decline. Despite extensive research, the molecular intricacies contributing to AD pathogenesis are inadequately understood. While single-cell omics technology holds great promise for application in AD, particularly in deciphering the understanding of different cell types and analyzing rare cell types and transcriptomic expression changes, it is unable to provide spatial distribution information, which is crucial for understanding the pathological processes of AD. In contrast, spatial multi-omics research emerges as a promising and comprehensive approach to analyzing tissue cells, potentially better suited for addressing these issues in AD. This article focuses on the latest advancements in spatial multi-omics technology and compares various techniques. Additionally, we provide an overview of current spatial omics-based research results in AD. These technologies play a crucial role in facilitating new discoveries and advancing translational AD research in the future. Despite challenges such as balancing resolution, increasing throughput, and data analysis, the application of spatial multi-omics holds immense potential in revolutionizing our understanding of human disease processes and identifying new biomarkers and therapeutic targets, thereby potentially contributing to the advancement of AD research.

Rapid, biochemical tagging of cellular activity history in vivo

Intracellular calcium (Ca2+) is ubiquitous to cell signaling across all biology. While existing fluorescent sensors and reporters can detect activated cells with elevated Ca2+ levels, these approaches require implants to deliver light to deep tissue, precluding their noninvasive use in freely-behaving animals. Here we engineered an enzyme-catalyzed approach that rapidly and biochemically tags cells with elevated Ca2+ in vivo. Ca2+-activated Split-TurboID (CaST) labels activated cells within 10 minutes with an exogenously-delivered biotin molecule. The enzymatic signal increases with Ca2+ concentration and biotin labeling time, demonstrating that CaST is a time-gated integrator of total Ca2+ activity. Furthermore, the CaST read-out can be performed immediately after activity labeling, in contrast to transcriptional reporters that require hours to produce signal. These capabilities allowed us to apply CaST to tag prefrontal cortex neurons activated by psilocybin, and to correlate the CaST signal with psilocybin-induced head-twitch responses in untethered mice.

Spatial multi-omics at subcellular resolution via high-throughput in situ pairwise sequencing

Technology for spatial multi-omics aids the discovery of new insights into cellular functions and disease mechanisms. Here we report the development and applicability of multi-omics in situ pairwise sequencing (MiP-seq), a method for the simultaneous detection of DNAs, RNAs, proteins and biomolecules at subcellular resolution. Compared with other in situ sequencing methods, MiP-seq enhances decoding capacity and reduces sequencing and imaging costs while maintaining the efficacy of detection of gene mutations, allele-specific expression and RNA modifications. MiP-seq can be integrated with in vivo calcium imaging and Raman imaging, which enabled us to generate a spatial multi-omics atlas of mouse brain tissues and to correlate gene expression with neuronal activity and cellular biochemical fingerprints. We also report a sequential dilution strategy for resolving optically crowded signals during in situ sequencing. High-throughput in situ pairwise sequencing may facilitate the multidimensional analysis of molecular and functional maps of tissues.

Cross-domain information fusion for enhanced cell population delineation in single-cell spatial-omics data

Cell population delineation and identification is an essential step in single-cell and spatial-omics studies. Spatial-omics technologies can simultaneously measure information from three complementary domains related to this task: expression levels of a panel of molecular biomarkers at single-cell resolution, relative positions of cells, and images of tissue sections, but existing computational methods for performing this task on single-cell spatial-omics datasets often relinquish information from one or more domains. The additional reliance on the availability of "atlas" training or reference datasets limits cell type discovery to well-defined but limited cell population labels, thus posing major challenges for using these methods in practice. Successful integration of all three domains presents an opportunity for uncovering cell populations that are functionally stratified by their spatial contexts at cellular and tissue levels: the key motivation for employing spatial-omics technologies in the first place. In this work, we introduce Cell Spatio- and Neighborhood-informed Annotation and Patterning (CellSNAP), a self-supervised computational method that learns a representation vector for each cell in tissue samples measured by spatial-omics technologies at the single-cell or finer resolution. The learned representation vector fuses information about the corresponding cell across all three aforementioned domains. By applying CellSNAP to datasets spanning both spatial proteomic and spatial transcriptomic modalities, and across different tissue types and disease settings, we show that CellSNAP markedly enhances de novo discovery of biologically relevant cell populations at fine granularity, beyond current approaches, by fully integrating cells' molecular profiles with cellular neighborhood and tissue image information.

Spatial Transcriptomics: A Powerful Tool in Disease Understanding and Drug Discovery

Recent advancements in modern science have provided robust tools for drug discovery. The rapid development of transcriptome sequencing technologies has given rise to single-cell transcriptomics and single-nucleus transcriptomics, increasing the accuracy of sequencing and accelerating the drug discovery process. With the evolution of single-cell transcriptomics, spatial transcriptomics (ST) technology has emerged as a derivative approach. Spatial transcriptomics has emerged as a hot topic in the field of omics research in recent years; it not only provides information on gene expression levels but also offers spatial information on gene expression. This technology has shown tremendous potential in research on disease understanding and drug discovery. In this article, we introduce the analytical strategies of spatial transcriptomics and review its applications in novel target discovery and drug mechanism unravelling. Moreover, we discuss the current challenges and issues in this research field that need to be addressed. In conclusion, spatial transcriptomics offers a new perspective for drug discovery.

Spatial transcriptomics in health and disease

The ability to localize hundreds of macromolecules to discrete locations, structures and cell types in a tissue is a powerful approach to understand the cellular and spatial organization of an organ. Spatially resolved transcriptomic technologies enable mapping of transcripts at single-cell or near single-cell resolution in a multiplex manner. The rapid development of spatial transcriptomic technologies has accelerated the pace of discovery in several fields, including nephrology. Its application to preclinical models and human samples has provided spatial information about new cell types discovered by single-cell sequencing and new insights into the cell-cell interactions within neighbourhoods, and has improved our understanding of the changes that occur in response to injury. Integration of spatial transcriptomic technologies with other omics methods, such as proteomics and spatial epigenetics, will further facilitate the generation of comprehensive molecular atlases, and provide insights into the dynamic relationships of molecular components in homeostasis and disease. This Review provides an overview of current and emerging spatial transcriptomic methods, their applications and remaining challenges for the field.

The broken Alzheimer's disease genome

The complex pathobiology of late-onset Alzheimer's disease (AD) poses significant challenges to therapeutic and preventative interventions. Despite these difficulties, genomics and related disciplines are allowing fundamental mechanistic insights to emerge with clarity, particularly with the introduction of high-resolution sequencing technologies. After all, the disrupted processes at the interface between DNA and gene expression, which we call the broken AD genome, offer detailed quantitative evidence unrestrained by preconceived notions about the disease. In addition to highlighting biological pathways beyond the classical pathology hallmarks, these advances have revitalized drug discovery efforts and are driving improvements in clinical tools. We review genetic, epigenomic, and gene expression findings related to AD pathogenesis and explore how their integration enables a better understanding of the multicellular imbalances contributing to this heterogeneous condition. The frontiers opening on the back of these research milestones promise a future of AD care that is both more personalized and predictive.

scBSP: A fast and accurate tool for identifying spatially variable genes from spatial transcriptomic data

Spatially resolved transcriptomics have enabled the inference of gene expression patterns within two and three-dimensional space, while introducing computational challenges due to growing spatial resolutions and sparse expressions. Here, we introduce scBSP, an open-source, versatile, and user-friendly package designed for identifying spatially variable genes in large-scale spatial transcriptomics. scBSP implements sparse matrix operation to significantly increase the computational efficiency in both computational time and memory usage, processing the high-definition spatial transcriptomics data for 19,950 genes on 181,367 spots within 10 seconds. Applied to diverse sequencing data and simulations, scBSP efficiently identifies spatially variable genes, demonstrating fast computational speed and consistency across various sequencing techniques and spatial resolutions for both two and three-dimensional data with up to millions of cells. On a sample with hundreds of thousands of sports, scBSP identifies SVGs accurately in seconds to on a typical desktop computer.