Spatially informed cell-type deconvolution for spatial transcriptomics

STModule: identifying tissue modules to uncover spatial components and characteristics of transcriptomic landscapes

Here we present STModule, a Bayesian method developed to identify tissue modules from spatially resolved transcriptomics that reveal spatial components and essential characteristics of tissues. STModule uncovers diverse expression signals in transcriptomic landscapes such as cancer, intraepithelial neoplasia, immune infiltration, outcome-related molecular features and various cell types, which facilitate downstream analysis and provide insights into tumor microenvironments, disease mechanisms, treatment development, and histological organization of tissues. STModule captures a broader spectrum of biological signals compared to other methods and detects novel spatial components. The tissue modules characterized by gene sets demonstrate greater robustness and transferability across different biopsies. STModule: https://github.com/rwang-z/STModule.git .

Profiling immune cell tissue niches in the spatial -omics era

Immune responses require complex, spatially coordinated interactions between immune cells and their tissue environment. For decades, we have imaged tissue sections to visualize a limited number of immune-related macromolecules in situ, functioning as surrogates for cell types or processes of interest. However, this inevitably provides a limited snapshot of the tissue's immune landscape. Recent developments in high-throughput spatial -omics technologies, particularly spatial transcriptomics, and its application to human samples has facilitated a more comprehensive understanding of tissue immunity by mapping fine-grained immune cell states to their precise tissue location while providing contextual information about their immediate cellular and tissue environment. These data provide opportunities to investigate mechanisms underlying the spatial distribution of immune cells and its functional implications, including the identification of immune niches, although the criteria used to define this term have been inconsistent. Here, we review recent technological and analytic advances in multiparameter spatial profiling, focusing on how these methods have generated new insights in translational immunology. We propose a 3-step framework for the definition and characterization of immune niches, which is powerfully facilitated by new spatial profiling methodologies. Finally, we summarize current approaches to analyze adaptive immune repertoires and lymphocyte clonal expansion in a spatially resolved manner.

Techniques and analytic workflow for spatial transcriptomics and its application to allergy and inflammation

Spatial profiling, through single-cell gene-level expression data paired with cell localization, offers unprecedented biologic insights within the intact spatial context of cells in healthy and diseased tissue, adding a novel dimension to data interpretation. This review summarizes recent developments in this field, its application to allergy and inflammation, and recent single-cell resolution platforms designed for spatial transcriptomics with a focus on data processing and analyses for efficient biologic interpretation of data. By preserving spatial context, these technologies provide critical insights into tissue architecture and cellular interactions that are unattainable with traditional transcriptomics methods, such as revealing localized inflammatory cell network in atopic dermatitis and T-cell interactions in the lung in chronic obstructive pulmonary disease. Spatial profiling offers opportunities for discovering novel biomarkers, defining compartmentalization of immune responses within tissues and individual diseases, and accelerating novel discoveries toward a greater understanding of fundamental disease mechanisms and, eventually, toward the development of future targeted therapies.

LETSmix: a spatially informed and learning-based domain adaptation method for cell-type deconvolution in spatial transcriptomics

Spatial transcriptomics (ST) enables the study of gene expression in spatial context, but many ST technologies face challenges due to limited resolution, leading to cell mixtures at each spot. We present LETSmix to deconvolve cell types by integrating spatial correlations through a tailored LETS filter, which leverages layer annotations, expression similarities, image texture features, and spatial coordinates to refine ST data. Additionally, LETSmix employs a mixup-augmented domain adaptation strategy to address discrepancies between ST and reference single-cell RNA sequencing data. Comprehensive evaluations across diverse ST platforms and tissue types demonstrate its high accuracy in estimating cell-type proportions and spatial patterns, surpassing existing methods (URL: https://github.com/ZhanYangen/LETSmix ).

Spatial transcriptomic analysis of tumor microenvironment in esophageal squamous cell carcinoma with HIV infection

BACKGROUND

Human Immunodeficiency Virus (HIV) is one of the most prevalent viruses, causing significant immune depletion in affected individuals. Current treatments can control HIV and prolong patients' lives, but new challenges have emerged. Increasing incidence of cancers occur in HIV patients. Esophageal squamous cell carcinoma (ESCC) is one of the most common cancers observed in HIV patients. However, the spatial cellular characteristics of HIV-related ESCC have not been explored, and the differences between HIV-ESCC and typical ESCC remain unclear.

METHODS

We performed spatial transcriptome sequencing on HIV-ESCC samples to depict the microenvironment and employed cell communication analysis and multiplex immunofluorescence to investigate the molecular mechanism in HIV-ESCC.

RESULTS

We found that HIV-ESCC exhibited a unique cellular composition, with fibroblasts and epithelial cells intermixed throughout the tumor tissue, lacking obvious spatial separation, while other cell types were sparse. Besides, HIV-ESCC exhibited an immune desert phenotype, characterized by a low degree of immune cell infiltration, with only a few SPP1+ macrophages showing immune resistance functions. Cell communication analysis and multiplex immunofluorescence staining revealed that tumor fibroblasts in HIV-ESCC interact with CD44+ epithelial cells via COL1A2, promoting the expression of PIK3R1 in epithelial cells. This interaction activates the PI3K-AKT signaling pathway, which contributes to the progression of HIV-ESCC.

CONCLUSIONS

Our findings depict the spatial microenvironment of HIV-ESCC and elucidate a molecular mechanism in the progression of HIV-ESCC. This will provide us insights into the molecular basis of HIV-ESCC and potential treatment strategies.

Systematic inference of super-resolution cell spatial profiles from histology images

Inferring cell spatial profiles from histology images is critical for cancer diagnosis and treatment in clinical settings. In this study, we report a weakly-supervised deep-learning method, HistoCell, to directly infer super-resolution cell spatial profiles consisting of cell types, cell states and their spatial network from histology images at the single-nucleus-level. Benchmark analysis demonstrates that HistoCell robustly achieves state-of-the-art performance in terms of cell type/states prediction solely from histology images across multiple cancer tissues. HistoCell can significantly enhance the deconvolution accuracy for the spatial transcriptomics data and enable accurate annotation of subtle cancer tissue architectures. Moreover, HistoCell is applied to de novo discovery of clinically relevant spatial organization indicators, including prognosis and drug response biomarkers, across diverse cancer types. HistoCell also enable image-based screening of cell populations that drives phenotype of interest, and is applied to discover the cell population and corresponding spatial organization indicators associated with gastric malignant transformation risk. Overall, HistoCell emerges as a powerful and versatile tool for cancer studies in histology image-only cohorts.

Spatial analysis of a complete DIPG-infiltrated brainstem reveals novel ligand-receptor mediators of tumour-to-TME crosstalk

Previous studies have highlighted the capacity of brain cancer cells to functionally interact with the tumour microenvironment (TME). This TME-cancer crosstalk crucially contributes to tumour cell invasion and disease progression. In this study, we performed spatial transcriptomic sequencing analysis of a complete annotated tumour-infiltrated brainstem from a single diffuse intrinsic pontine glioma (DIPG) patient. Gene signatures from ten sequential tumour regions were analysed to assess mechanisms of disease progression and oncogenic interactions with the TME. We identified four distinct tumour subpopulations and assessed respective ligand-receptor pairs that actively promote DIPG tumour progression via crosstalk with endothelial, neuronal and immune cell communities. Our analysis found potential targetable mediators of tumour-to-TME communication, including members of the complement component system and the neuropeptide/GPCR ligand-receptor pair ADCYAP1-ADCYAP1R1. These interactions could influence DIPG tumour progression and represent novel therapeutic targets.

Deconvolution and inference of spatial communication through optimization algorithm for spatial transcriptomics

Spatial transcriptomics technologies can capture gene expression at spatial loci. However, at certain resolutions, the obtained gene expression reflects the sum of either a heterogeneous or homogeneous set of cells, rather than individual cell. This limitation gives rise to the deconvolution algorithm to make cell-type inferences at each location. Yet, the vast majority of deconvolution methods that have been developed ignore the spatial information of the tissue and the communications between the cells or spots. To overcome these afflictions, we proposed a deconvolution method, non-negative least squares-based and optimization search-based deconvolution (NODE), that combines cell-type-specific information from single-cell RNA sequencing (scRNA-seq) and intercellular communications in tissue. NODE deconvolution algorithm, incorporating the spatial information of the tissue, allows us to quantify intercellular communications at the same instant. NODE can not only utilize optimization method to infer the deconvolution results of spatial transcriptomics data and reduce the probability of overfitting situations, but also make reasonable inferences for spatial communications. Subsequently, we applied NODE to four datasets to validate the correctness of the NODE deconvolution results and compare them with existing deconvolution algorithms. NODE also inferred spatial communications and validated them in tissue development of human heart.

Precise gene expression deconvolution in spatial transcriptomics with STged

Spatially resolved transcriptomics (SRT) has transformed tissue biology by linking gene expression profiles with spatial information. However, sequencing-based SRT methods aggregate signals from multiple cell types within capture locations ("spots"), masking cell-type-specific gene expression patterns. Traditional cell-type deconvolution methods estimate cell compositions within spots but fail to resolve cell-type-specific gene expression, limiting their ability to uncover critical biological processes such as cellular interactions and microenvironmental dynamics. Here, we present STged (spatial transcriptomic gene expression deconvolution), a novel computational framework that goes beyond traditional deconvolution by reconstructing cell-type-specific gene expression profiles from mixed spots. STged integrates graph-based spatial correlations and reference-derived gene signatures using a non-negative least-squares regression framework, achieving precise and biologically meaningful deconvolution. Comprehensive simulations show that STged consistently outperforms existing methods in accuracy and robustness. Applications to human pancreatic ductal adenocarcinoma and human squamous cell carcinoma datasets reveal its capacity to identify microenvironment-specific highly variable genes, reconstruct spatial cell-cell communication networks, and resolve tissue architecture at near-single-cell resolution. In mouse kidney tissues, STged uncovers dynamic spatial gene expression patterns and distinct gene programs, advancing our understanding of tissue heterogeneity and cellular dynamics.

Searching for the cellular underpinnings of the selective vulnerability to tauopathic insults in Alzheimer's disease

Neurodegenerative diseases such as Alzheimer's disease exhibit pathological changes in the brain that proceed in a stereotyped and regionally specific fashion. However, the cellular underpinnings of regional vulnerability are poorly understood, in part because whole-brain maps of a comprehensive collection of cell types have been inaccessible. Here, we deployed a recent cell-type mapping pipeline, Matrix Inversion and Subset Selection (MISS), to determine the brain-wide distributions of pan-hippocampal and neocortical cells in the mouse, and then used these maps to identify general principles of cell-type-based selective vulnerability in PS19 mouse models. We found that hippocampal glutamatergic neurons as a whole were significantly positively associated with regional tau deposition, suggesting vulnerability, while cortical glutamatergic and GABAergic neurons were negatively associated. We also identified oligodendrocytes as the single-most strongly negatively associated cell type. Further, cell-type distributions were more predictive of end-time-point tau pathology than AD-risk-gene expression. Using gene ontology analysis, we found that the genes that are directly correlated to tau pathology are functionally distinct from those that constitutively embody the vulnerable cells. In short, we have elucidated cell-type correlates of tau deposition across mouse models of tauopathy, advancing our understanding of selective cellular vulnerability at a whole-brain level.

Tumor-associated stroma shapes the spatial tumor immune microenvironment of primary Ewing sarcomas

To date, few studies have detailed the tumor microenvironment (TME) of Ewing sarcoma (EwS). The TME has a vital role in cancer survival and progression with implications in drug resistance and immune escape. By performing spatially resolved transcriptomic analysis of primary treatment-naïve EwS samples, we discovered greater stromal enrichment in localized EwS tumors compared to metastatic EwS tumors. Through spatial ligand-receptor analysis, we show that the stromal enriched regions harbor unique extracellular matrix related cytokines, immune recruitment and proinflammatory microenvironmental signals, implying EwS stroma may play an anti-tumorigenic role by acting as an immune recruitment center. All EwS tumors expressed pro-tumorigenic MIF-CD74 immune signaling, suggesting a potential immune-evasive mechanism and immunotherapy target. Our findings provide insight into tumor cell/stromal cell interactions in EwS and serve as a valuable resource for further investigations in the tumor immune microenvironment of EwS.

Single-cell spatial immune profiling for precision immunotherapy in Lynch syndrome

Lynch syndrome (LS) is the most common hereditary colorectal cancer (CRC) predisposition syndrome, characterized by a high mutational burden and microsatellite instability-high (MSI-H) tumors. Immunology of LS-associated CRC (LS-CRC) is unique, with significant implications for treatment. Despite well-established knowledge of LS immunology, immunotherapy dose and treatment response can vary significantly based on local tumor immunity and specific germline pathogenic variant of LS genes. This variability necessitates tailored surveillance strategies and new personalised immunotherapy approaches for LS patients. LS-CRC often benefits from immunotherapy due to the distinct tumor microenvironment (TME) and the variety of tumor infiltrating lymphocytes (TILs). This perspective discusses a novel approach of analysing spatial TILs at a single-cell level using tumor whole slide images (WSIs) that accounts for the distinct TME of LS-CRC. By emphasizing the necessity of personalized medicine in hereditary cancer syndromes, the future research and clinical practices that enhance patient outcomes through precision oncology is inspired.

Starfysh integrates spatial transcriptomic and histologic data to reveal heterogeneous tumor-immune hubs

Spatially resolved gene expression profiling provides insight into tissue organization and cell-cell crosstalk; however, sequencing-based spatial transcriptomics (ST) lacks single-cell resolution. Current ST analysis methods require single-cell RNA sequencing data as a reference for rigorous interpretation of cell states, mostly do not use associated histology images and are not capable of inferring shared neighborhoods across multiple tissues. Here we present Starfysh, a computational toolbox using a deep generative model that incorporates archetypal analysis and any known cell type markers to characterize known or new tissue-specific cell states without a single-cell reference. Starfysh improves the characterization of spatial dynamics in complex tissues using histology images and enables the comparison of niches as spatial hubs across tissues. Integrative analysis of primary estrogen receptor (ER)-positive breast cancer, triple-negative breast cancer (TNBC) and metaplastic breast cancer (MBC) tissues led to the identification of spatial hubs with patient- and disease-specific cell type compositions and revealed metabolic reprogramming shaping immunosuppressive hubs in aggressive MBC.

Spatial integration of multi-omics single-cell data with SIMO

Technical limitations in spatial and single-cell omics sequencing pose challenges for capturing and describing multimodal information at the spatial scale. To address this, we develop SIMO, a computational method designed for the Spatial Integration of Multi-Omics datasets through probabilistic alignment. Unlike previous tools, SIMO not only integrates spatial transcriptomics with single-cell RNA-seq but expands beyond, enabling integration across multiple single-cell modalities, such as chromatin accessibility and DNA methylation, which have not been co-profiled spatially before. We benchmark SIMO on simulated datasets, demonstrating its high accuracy and robustness. Further application on biological datasets reveals SIMO's ability to detect topological patterns of cells and their regulatory modes across multiple omics layers. Through comprehensive analysis of real-world data, SIMO uncovers multimodal spatial heterogeneity, offering deeper insights into the spatial organization and regulation of biological molecules. These findings position SIMO as a powerful tool for advancing spatial biology by revealing previously inaccessible multimodal insights.

Flexible analysis of spatial transcriptomics data (FAST): a deconvolution approach

MOTIVATION

Spatial transcriptomics is a state-of-art technique that allows researchers to study gene expression patterns in tissues over the spatial domain. As a result of technical limitations, the majority of spatial transcriptomics techniques provide bulk data for each sequencing spot. Consequently, in order to obtain high-resolution spatial transcriptomics data, performing deconvolution becomes essential. Most existing deconvolution methods rely on reference data (e.g., single-cell data), which may not be available in real applications. Current reference-free methods encounter limitations due to their dependence on distribution assumptions, reliance on marker genes, or the absence of leveraging histology and spatial information. Consequently, there is a critical need for the development of highly flexible, robust, and user-friendly reference-free deconvolution methods capable of unifying or leveraging case-specific information in the analysis of spatial transcriptomics data.

RESULTS

We propose a novel reference-free method based on regularized non-negative matrix factorization (NMF), named Flexible Analysis of Spatial Transcriptomics (FAST), that can effectively incorporate gene expression data, spatial, and histology information into a unified deconvolution framework. Compared to existing methods, FAST imposes fewer distribution assumptions, utilizes the spatial structure information of tissues, and encourages interpretable factorization results. These features enable greater flexibility and accuracy, making FAST an effective tool for deciphering the complex cell-type composition of tissues and advancing our understanding of various biological processes and diseases. Extensive simulation studies have shown that FAST outperforms other existing reference-free methods. In real data applications, FAST is able to uncover the underlying tissue structures and identify the corresponding marker genes.

Statistical identification of cell type-specific spatially variable genes in spatial transcriptomics

An essential task in spatial transcriptomics is identifying spatially variable genes (SVGs). Here, we present Celina, a statistical method for systematically detecting cell type-specific SVGs (ct-SVGs)-a subset of SVGs exhibiting distinct spatial expression patterns within specific cell types. Celina utilizes a spatially varying coefficient model to accurately capture each gene's spatial expression pattern in relation to the distribution of cell types across tissue locations, ensuring effective type I error control and high power. Celina proves powerful compared to existing methods in single-cell resolution spatial transcriptomics and stands as the only effective solution for spot-resolution spatial transcriptomics. Applied to five real datasets, Celina uncovers ct-SVGs associated with tumor progression and patient survival in lung cancer, identifies metagenes with unique spatial patterns linked to cell proliferation and immune response in kidney cancer, and detects genes preferentially expressed near amyloid-β plaques in an Alzheimer's model.

KanCell: dissecting cellular heterogeneity in biological tissues through integrated single-cell and spatial transcriptomics

KanCell is a deep learning model based on Kolmogorov-Arnold networks (KAN) designed to enhance cellular heterogeneity analysis by integrating single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics (ST) data. ST technologies provide insights into gene expression within tissue context, revealing cellular interactions and microenvironments. To fully leverage this potential, effective computational models are crucial. We evaluate KanCell on both simulated and real datasets from technologies such as STARmap, Slide-seq, Visium, and Spatial Transcriptomics. Our results demonstrate that KanCell outperforms existing methods across metrics like PCC, SSIM, COSSIM, RMSE, JSD, ARS, and ROC, with robust performance under varying cell numbers and background noise. Real-world applications on human lymph nodes, hearts, melanoma, breast cancer, dorsolateral prefrontal cortex, and mouse embryo brains confirmed its reliability. Compared with traditional approaches, KanCell effectively captures non-linear relationships and optimizes computational efficiency through KAN, providing an accurate and efficient tool for ST. By improving data accuracy and resolving cell type composition, KanCell reveals cellular heterogeneity, clarifies disease microenvironments, and identifies therapeutic targets, addressing complex biological challenges.

Comprehensive single-cell atlas of colorectal neuroendocrine tumors with liver metastases: unraveling tumor microenvironment heterogeneity between primary lesions and metastases

BACKGROUND

Colorectal neuroendocrine tumors with liver metastases (CRNELM) are associated with a poorer prognosis compared to their nonmetastatic counterparts. A comprehensive understanding of the tumor microenvironment (TME) heterogeneity between primary lesions (PL) and liver metastases (LM) could provide crucial insights for enhancing clinical management strategies for these patients.

METHODS

We utilized single-cell RNA sequencing to analyze fresh tissue samples from CRNELM patients, aiming to elucidate the variations in TME between PL and LM. Complementary multidimensional validation was achieved through spatial transcriptomics, bulk RNA sequencing, and multiplex immunohistochemistry/immunofluorescence.

RESULTS

Our single-cell RNA sequencing analysis revealed that LM harboured a higher proportion of CD8 + T cells, CD4 + T cells, NK cells, NKT cells, and B cells exhibiting a stress-like phenotype compared to PL. RGS5 + pericytes may play a role in the stress-like phenotype observed in immune cells within LM. MCs in PL (PL_MCs) and LM (LM_MCs) exhibit distinct activation of tumor-associated signaling pathways. Notably, COLEC11 + matrix cancer-associated fibroblasts (COLEC11_mCAFs) were found to be significantly associated with LM_MCs. Cell communication analysis unveiled potential targetable receptor-ligand interactions between COLEC11_mCAFs and LM_MCs. Multidimensional validation confirmed the prominence of the characteristic stress-like phenotypes, including HSPA6_CD8_Tstr, HSPA6_NK, and COLEC11_mCAFs in LM. Moreover, a higher abundance of COLEC11_mCAFs correlated with poorer survival rates in the neuroendocrine tumor patient cohort.

CONCLUSION

Overall, our study provides the first single-cell analysis of the cellular and molecular differences between PL and LM in CRNELM patients. We identified distinct cell subsets and receptor-ligand interactions that may drive TME discrepancies and support metastatic tumor growth. These insights highlight potential therapeutic targets and inform strategies for better managing CRNELM patients.

Pairpot: a database with real-time lasso-based analysis tailored for paired single-cell and spatial transcriptomics

Paired single-cell and spatially resolved transcriptomics (SRT) data supplement each other, providing in-depth insights into biological processes and disease mechanisms. Previous SRT databases have limitations in curating sufficient single-cell and SRT pairs (SC-SP pairs) and providing real-time heuristic analysis, which hinder the effort to uncover potential biological insights. Here, we developed Pairpot (http://pairpot.bioxai.cn), a database tailored for paired single-cell and SRT data with real-time heuristic analysis. Pairpot curates 99 high-quality pairs including 1,425,656 spots from 299 datasets, and creates the association networks. It constructs the curated pairs by integrating multiple slices and establishing potential associations between single-cell and SRT data. On this basis, Pairpot adopts semi-supervised learning that enables real-time heuristic analysis for SC-SP pairs where Lasso-View refines the user-selected SRT domains within milliseconds, Pair-View infers cell proportions of spots based on user-selected cell types in real-time and Layer-View displays SRT slices using a 3D hierarchical layout. Experiments demonstrated Pairpot's efficiency in identifying heterogeneous domains and cell proportions.

Spatial Transcriptomics: Biotechnologies, Computational Tools, and Neuroscience Applications

Spatial transcriptomics (ST) represents a revolutionary approach in molecular biology, providing unprecedented insights into the spatial organization of gene expression within tissues. This review aims to elucidate advancements in ST technologies, their computational tools, and their pivotal applications in neuroscience. It is begun with a historical overview, tracing the evolution from early image-based techniques to contemporary sequence-based methods. Subsequently, the computational methods essential for ST data analysis, including preprocessing, cell type annotation, spatial clustering, detection of spatially variable genes, cell-cell interaction analysis, and 3D multi-slices integration are discussed. The central focus of this review is the application of ST in neuroscience, where it has significantly contributed to understanding the brain's complexity. Through ST, researchers advance brain atlas projects, gain insights into brain development, and explore neuroimmune dysfunctions, particularly in brain tumors. Additionally, ST enhances understanding of neuronal vulnerability in neurodegenerative diseases like Alzheimer's and neuropsychiatric disorders such as schizophrenia. In conclusion, while ST has already profoundly impacted neuroscience, challenges remain issues such as enhancing sequencing technologies and developing robust computational tools. This review underscores the transformative potential of ST in neuroscience, paving the way for new therapeutic insights and advancements in brain research.

Epitranscriptomic rRNA fingerprinting reveals tissue-of-origin and tumor-specific signatures

Mammalian ribosomal RNA (rRNA) molecules are highly abundant RNAs, decorated with over 220 rRNA modifications. Previous works have shown that some rRNA modification types can be dynamically regulated; however, how and when the mammalian rRNA modification landscape is remodeled remains largely unexplored. Here, we employ direct RNA sequencing to chart the human and mouse rRNA epitranscriptome across tissues, developmental stages, cell types, and disease. Our analyses reveal multiple rRNA sites that are differentially modified in a tissue- and/or developmental stage-specific manner, including previously unannotated modified sites. We demonstrate that rRNA modification patterns can be used for tissue and cell-type identification, which we hereby term "epitranscriptomic fingerprinting." We then explore rRNA modification patterns in normal-tumor matched samples from lung cancer patients, finding that epitranscriptomic fingerprinting accurately classifies clinical samples into normal and tumor groups from only 250 reads per sample, demonstrating the potential of rRNA modifications as diagnostic biomarkers.

SpatialDeX Is a Reference-Free Method for Cell-Type Deconvolution of Spatial Transcriptomics Data in Solid Tumors

The rapid development of spatial transcriptomics (ST) technologies has enabled transcriptome-wide profiling of gene expression in tissue sections. Despite the emergence of single-cell resolution platforms, most ST sequencing studies still operate at a multicell resolution. Consequently, deconvolution of cell identities within the spatial spots has become imperative for characterizing cell-type-specific spatial organization. To this end, we developed Spatial Deconvolution Explorer (SpatialDeX), a regression model-based method for estimating cell-type proportions in tumor ST spots. SpatialDeX exhibited comparable performance to reference-based methods and outperformed other reference-free methods with simulated ST data. Using experimental ST data, SpatialDeX demonstrated superior performance compared with both reference-based and reference-free approaches. Additionally, a pan-cancer clustering analysis on tumor spots identified by SpatialDeX unveiled distinct tumor progression mechanisms both within and across diverse cancer types. Overall, SpatialDeX is a valuable tool for unraveling the spatial cellular organization of tissues from ST data without requiring single-cell RNA-seq references. Significance: The development of a reference-free method for deconvolving the identity of cells in spatial transcriptomics datasets enables exploration of tumor architecture to gain deeper insights into the dynamics of the tumor microenvironment.

Single-cell RNA sequencing and spatial transcriptome analysis in bladder cancer: Current status and future perspectives

Background

Bladder cancer is one of the most prevalent malignancies, and the mechanisms underlying its progression and the role of the tumor microenvironment (TME) are unclear. Recent advancements in single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics (ST) enable detailed analysis of the cellular heterogeneity, gene expression, and cell-cell interactions in bladder diseases.

Methodology

We conducted a comprehensive search for recent articles that have investigated bladder diseases using scRNA-seq and ST.

Results

scRNA-seq and ST have led to significant discoveries in bladder disease research. These technologies have enabled the identification of multiple molecular subtypes within individual tumors and of the mechanisms of treatment resistance. Additionally, molecular differences based on gender have been explored, explaining the heterogeneity of the incidence and progression of bladder cancer. These findings deepen our understanding of the pathology of bladder diseases and highlight the transformative potential of scRNA-seq and ST in identifying novel biomarkers and therapeutic targets.

Conclusions

Integrating scRNA-seq and ST has considerably enhanced our understanding of tumor heterogeneity and the tumor microenvironment within tissues. These insights may lead to the development of personalized therapies and the improvement of patient outcomes. Several challenges, such as technical limitations and access difficulties, need to be addressed for the future clinical application of these technologies.

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.

ImSpiRE: image feature-aided spatial resolution enhancement method

The resolution of most spatially resolved transcriptomic technologies usually cannot attain the single-cell level, limiting their applications in biological discoveries. Here, we introduce ImSpiRE, an image feature-aided spatial resolution enhancement method for in situ capturing spatial transcriptome. Taking the information stored in histological images, ImSpiRE solves an optimal transport problem to redistribute the expression profiles of spots to construct new transcriptional profiles with enhanced resolution, together with extending the gene expression profiles into unmeasured regions. Applications to multiple datasets confirm that ImSpiRE can enhance spatial resolution to the subspot level while contributing to the discovery of tissue domains, signaling communication patterns, and spatiotemporal characterization.

Granuloma Annulare Exhibits Mixed Immune and Macrophage Polarization Profiles with Spatial Transcriptomics

Granuloma annulare (GA) is an idiopathic condition characterized by granulomatous inflammation in the skin. Prior studies have suggested that GA develops from various triggers, leading to a complex interplay involving innate and adaptive immunity, tissue remodeling, and fibrosis. Macrophages are the major immune cells comprising GA granulomas; however, the molecular drivers and inflammatory signaling cascade behind macrophage activation are poorly understood. Histologically, GA exhibits both palisaded and interstitial patterns on histology; however, the molecular composition of GA at the spatial level remains unexplored. GA is a condition without Food and Drug Administration-approved therapies despite the significant impact of GA on QOL. Spatial transcriptomics is a valuable tool for profiling localized, genome-wide gene expression changes across tissues, with emerging applications in clinical medicine. To improve our understanding of the spatially localized gene expression patterns underlying GA, we profiled the spatial gene expression landscape from 6 patients with GA. Our findings revealed mixed T helper 1 and T helper 2 signals comprising the GA microenvironment and spatially distinct M1 and M2 macrophage polarization characteristics. IFN-γ and TNF signals emerged as important regulators of GA granulomatous inflammation, and IL-32 emerged as a key driver of granulomatous inflammation. Overall, our spatial transcriptomics data indicate that GA exhibits mixed immune and macrophage polarization.

Spatial analyses revealed CXCL5 and SLC6A14 as the markers of microvascular invasion in intrahepatic cholangiocarcinoma

BACKGROUND

Microvascular invasion (MVI) is a critical prognostic factor in intrahepatic cholangiocarcinoma (ICC), strongly associated with postoperative recurrence. However, the phenotypic features and spatial organization of MVI remain inadequately understood.

METHODS

We performed a spatial transcriptomic analysis on 29,632 spots from six ICC samples, manually delineating MVI clusters using the cloupe software. Key biomarkers were identified and validated in an independent cohort of 135 ICC patients. Functional and survival analyses were conducted to assess clinical relevance, and cell-cell communication pathways were investigated.

RESULTS

MVI regions exhibited heightened proliferation, angiogenesis, and epithelial-mesenchymal transition, driven by increased expression of transcription factors SOX10, ZEB1, and SNAI2. CXCL5 and SLC6A14 were identified as potential MVI biomarkers and showed high expression in tumor-invasive areas. Serum CXCL5 demonstrated strong predictive power for vascular invasion (AUC = 0.92) and intrahepatic metastasis (AUC = 0.96). High expression of both CXCL5 and SLC6A14 was associated with the worst survival outcomes. MVI regions were enriched with immunosuppressive MRC1+ macrophages and exhibited elevated immune checkpoint expression, including HAVCR2 and TIGHT, indicative of immune resistance. Cell-cell communication analysis revealed CXCL5-CXCR2 and LGALS9-HAVCR2 as key ligand-receptor pairs contributing to the immunosuppressive microenvironment.

CONCLUSIONS

This study identifies CXCL5 and SLC6A14 as key biomarkers of MVI, highlighting their roles in tumor proliferation, immune resistance, and poor clinical outcomes. These findings provide valuable insights into the spatial organization of MVI and its contribution to ICC progression, offering potential therapeutic targets.

Repurposing large-format microarrays for scalable spatial transcriptomics

Spatiomolecular analyses are key to study tissue functions and malfunctions. However, we lack profiling tools for spatial transcriptomics that are easy to adopt, low cost and scalable in terms of sample size and number. Here, we describe a method, Array-seq, to repurpose classical oligonucleotide microarrays for spatial transcriptomics profiling. We generate Array-seq slides from microarrays carrying custom-design probes that contain common sequences flanking unique barcodes at known coordinates. Then we perform a simple, two-step reaction that produces mRNA capture probes across all spots on the microarray. We demonstrate that Array-seq yields spatial transcriptomes with high detection sensitivity and localization specificity using histological sections from mouse tissues as test systems. Moreover, we show that the large surface area of Array-seq slides yields spatial transcriptomes (i) at high throughput by profiling multi-organ sections, (ii) in three dimensions by processing serial sections from one sample, and (iii) across whole human organs. Thus, by combining classical DNA microarrays and next-generation sequencing, we have created a simple and flexible platform for spatiomolecular studies of small-to-large specimens at scale.

NLSDeconv: an efficient cell-type deconvolution method for spatial transcriptomics data

SUMMARY

Spatial transcriptomics (ST) allows gene expression profiling within intact tissue samples but lacks single-cell resolution. This necessitates computational deconvolution methods to estimate the contributions of distinct cell types. This article introduces NLSDeconv, a novel cell-type deconvolution method based on non-negative least squares, along with an accompanying Python package. Benchmarking against 18 existing deconvolution methods on various ST datasets demonstrates NLSDeconv's competitive statistical performance and superior computational efficiency.

AVAILABILITY AND IMPLEMENTATION

NLSDeconv is freely available at https://github.com/tinachentc/NLSDeconv as a Python package.

Advances and applications in single-cell and spatial genomics

The applications of single-cell and spatial technologies in recent times have revolutionized the present understanding of cellular states and the cellular heterogeneity inherent in complex biological systems. These advancements offer unprecedented resolution in the examination of the functional genomics of individual cells and their spatial context within tissues. In this review, we have comprehensively discussed the historical development and recent progress in the field of single-cell and spatial genomics. We have reviewed the breakthroughs in single-cell multi-omics technologies, spatial genomics methods, and the computational strategies employed toward the analyses of single-cell atlas data. Furthermore, we have highlighted the advances made in constructing cellular atlases and their clinical applications, particularly in the context of disease. Finally, we have discussed the emerging trends, challenges, and opportunities in this rapidly evolving field.

Spall: accurate and robust unveiling cellular landscapes from spatially resolved transcriptomics data using a decomposition network

Recent developments in spatially resolved transcriptomics (SRT) enable the characterization of spatial structures for different tissues. Many decomposition methods have been proposed to depict the cellular distribution within tissues. However, existing computational methods struggle to balance spatial continuity in cell distribution with the preservation of cell-specific characteristics. To address this, we propose Spall, a novel decomposition network that integrates scRNA-seq data with SRT data to accurately infer cell type proportions. Spall introduced the GATv2 module, featuring a flexible dynamic attention mechanism to capture relationships between spots. This improves the identification of cellular distribution patterns in spatial analysis. Additionally, Spall incorporates skip connections to address the loss of cell-specific information, thereby enhancing the prediction capability for rare cell types. Experimental results show that Spall outperforms the state-of-the-art methods in reconstructing cell distribution patterns on multiple datasets. Notably, Spall reveals tumor heterogeneity in human pancreatic ductal adenocarcinoma samples and delineates complex tissue structures, such as the laminar organization of the mouse cerebral cortex and the mouse cerebellum. These findings highlight the ability of Spall to provide reliable low-dimensional embeddings for downstream analyses, offering new opportunities for deciphering tissue structures.

Integration of Imaging-based and Sequencing-based Spatial Omics Mapping on the Same Tissue Section via DBiTplus

Spatially mapping the transcriptome and proteome in the same tissue section can significantly advance our understanding of heterogeneous cellular processes and connect cell type to function. Here, we present Deterministic Barcoding in Tissue sequencing plus (DBiTplus), an integrative multi-modality spatial omics approach that combines sequencing-based spatial transcriptomics and image-based spatial protein profiling on the same tissue section to enable both single-cell resolution cell typing and genome-scale interrogation of biological pathways. DBiTplus begins with in situ reverse transcription for cDNA synthesis, microfluidic delivery of DNA oligos for spatial barcoding, retrieval of barcoded cDNA using RNaseH, an enzyme that selectively degrades RNA in an RNA-DNA hybrid, preserving the intact tissue section for high-plex protein imaging with CODEX. We developed computational pipelines to register data from two distinct modalities. Performing both DBiT-seq and CODEX on the same tissue slide enables accurate cell typing in each spatial transcriptome spot and subsequently image-guided decomposition to generate single-cell resolved spatial transcriptome atlases. DBiTplus was applied to mouse embryos with limited protein markers but still demonstrated excellent integration for single-cell transcriptome decomposition, to normal human lymph nodes with high-plex protein profiling to yield a single-cell spatial transcriptome map, and to human lymphoma FFPE tissue to explore the mechanisms of lymphomagenesis and progression. DBiTplusCODEX is a unified workflow including integrative experimental procedure and computational innovation for spatially resolved single-cell atlasing and exploration of biological pathways cell-by-cell at genome-scale.

Single-cell RNA sequencing in stroke and traumatic brain injury: Current achievements, challenges, and future perspectives on transcriptomic profiling

Single-cell RNA sequencing (scRNA-seq) is a high-throughput transcriptomic approach with the power to identify rare cells, discover new cellular subclusters, and describe novel genes. scRNA-seq can simultaneously reveal dynamic shifts in cellular phenotypes and heterogeneities in cellular subtypes. Since the publication of the first protocol on scRNA-seq in 2009, this evolving technology has continued to improve, through the use of cell-specific barcodes, adoption of droplet-based systems, and development of advanced computational methods. Despite induction of the cellular stress response during the tissue dissociation process, scRNA-seq remains a popular technology, and commercially available scRNA-seq methods have been applied to the brain. Recent advances in spatial transcriptomics now allow the researcher to capture the positional context of transcriptional activity, strengthening our knowledge of cellular organization and cell-cell interactions in spatially intact tissues. A combination of spatial transcriptomic data with proteomic, metabolomic, or chromatin accessibility data is a promising direction for future research. Herein, we provide an overview of the workflow, data analyses methods, and pros and cons of scRNA-seq technology. We also summarize the latest achievements of scRNA-seq in stroke and acute traumatic brain injury, and describe future applications of scRNA-seq and spatial transcriptomics.

SMART: spatial transcriptomics deconvolution using marker-gene-assisted topic model

While spatial transcriptomics offer valuable insights into gene expression patterns within the spatial context of tissue, many technologies do not have a single-cell resolution. Here, we present SMART, a marker gene-assisted deconvolution method that simultaneously infers the cell type-specific gene expression profile and the cellular composition at each spot. Using multiple datasets, we show that SMART outperforms the existing methods in realistic settings. It also provides a two-stage approach to enhance its performance on cell subtypes. The covariate model of SMART enables the identification of cell type-specific differentially expressed genes across conditions, elucidating biological changes at a single-cell-type resolution.

Inferring allele-specific copy number aberrations and tumor phylogeography from spatially resolved transcriptomics

Analyzing somatic evolution within a tumor over time and across space is a key challenge in cancer research. Spatially resolved transcriptomics (SRT) measures gene expression at thousands of spatial locations in a tumor, but does not directly reveal genomic aberrations. We introduce CalicoST, an algorithm to simultaneously infer allele-specific copy number aberrations (CNAs) and reconstruct spatial tumor evolution, or phylogeography, from SRT data. CalicoST identifies important classes of CNAs-including copy-neutral loss of heterozygosity and mirrored subclonal CNAs-that are invisible to total copy number analysis. Using nine patients' data from the Human Tumor Atlas Network, CalicoST achieves an average accuracy of 86%, approximately 21% higher than existing methods. CalicoST reconstructs a tumor phylogeography in three-dimensional space for two patients with multiple adjacent slices. CalicoST analysis of multiple SRT slices from a cancerous prostate organ reveals mirrored subclonal CNAs on the two sides of the prostate, forming a bifurcating phylogeography in both genetic and physical space.

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.

Deciphering spatial domains from spatially resolved transcriptomics through spatially regularized deep graph networks

BACKGROUND

Recent advancements in spatially resolved transcriptomics (SRT) have opened up unprecedented opportunities to explore gene expression patterns within spatial contexts. Deciphering spatial domains is a critical task in spatial transcriptomic data analysis, aiding in the elucidation of tissue structural heterogeneity and biological functions. However, existing spatial domain detection methods ignore the consistency of expression patterns and spatial arrangements between spots, as well as the severe gene dropout phenomenon present in SRT data, resulting in suboptimal performance in identifying tissue spatial heterogeneity.

RESULTS

In this paper, we introduce a novel framework, spatially regularized deep graph networks (SR-DGN), which integrates gene expression profiles with spatial information to learn spatially-consistent and informative spot representations. Specifically, SR-DGN employs graph attention networks (GAT) to adaptively aggregate gene expression information from neighboring spots, considering local expression patterns between spots. In addition, the spatial regularization constraint ensures the consistency of neighborhood relationships between physical and embedded spaces in an end-to-end manner. SR-DGN also employs cross-entropy (CE) loss to model gene expression states, effectively mitigating the impact of noisy gene dropouts.

CONCLUSIONS

Experimental results demonstrate that SR-DGN outperforms state-of-the-art methods in spatial domain identification across SRT data from different sequencing platforms. Moreover, SR-DGN is capable of recovering known microanatomical structures, yielding clearer low-dimensional visualizations and more accurate spatial trajectory inferences.

Spatiotemporal analysis of lung immune dynamics in lethal Coccidioides posadasii infection

Coccidioidomycosis, or Valley fever, is a lung disease caused by inhalation of Coccidioides fungi, prevalent in the Southwestern United States, Mexico, and parts of Central and South America. Annually, the United States reports 10,000-20,000 cases, although those numbers are expected to increase as climate change expands the fungal geographic range. While 60% of infections are asymptomatic, 40% symptomatic infections are often misdiagnosed due to similarities with bronchitis or pneumonia. A small subset of infection progress to severe illness, necessitating a better understanding of immune responses during lethal infection. Using single-cell RNA sequencing and spatial transcriptomics, we characterized lung responses during Coccidioides infection. We identified monocyte-derived Spp1-expressing macrophages as potential mediators of tissue remodeling and fibrosis, marked by high expression of profibrotic and proinflammatory transcripts. These macrophages showed elevated TGF-β and IL-6 signaling, pathways involved in fibrosis pathogenesis. Additionally, we observed significant neutrophil infiltration and defective lymphocyte responses, indicating severe adaptive immunity dysregulation in lethal, acute infection. These findings enhance our understanding of Coccidioides infection and suggest new therapeutic targets.IMPORTANCECoccidioidomycosis, commonly known as Valley fever, is a lung disease caused by the inhalation of Coccidioides fungi, which is prevalent in the Southwestern United States, Mexico, and parts of Central and South America. With climate change potentially expanding the geographic range of this fungus, understanding the immune responses during severe infections is crucial. Our study used advanced techniques to analyze lung responses during Coccidioides infection, identifying specific immune cells that may contribute to tissue damage and fibrosis. These findings provide new insights into the disease mechanisms and suggest potential targets for therapeutic intervention, which could improve outcomes for patients suffering from severe Valley fever.

Spatial dissection of tumour microenvironments in gastric cancers reveals the immunosuppressive crosstalk between CCL2+ fibroblasts and STAT3-activated macrophages

BACKGROUND

A spatially resolved, niche-level analysis of tumour microenvironments (TME) can provide insights into cellular interactions and their functional impacts in gastric cancers (GC).

OBJECTIVE

Our goal was to translate the spatial organisation of GC ecosystems into a functional landscape of cellular interactions involving malignant, stromal and immune cells.

DESIGN

We performed spatial transcriptomics on nine primary GC samples using the Visium platform to delineate the transcriptional landscape and dynamics of malignant, stromal and immune cells within the GC tissue architecture, highlighting cellular crosstalks and their functional consequences in the TME.

RESULTS

GC spatial transcriptomes with substantial cellular heterogeneity were delineated into six regional compartments. Specifically, the fibroblast-enriched TME upregulates epithelial-to-mesenchymal transformation and immunosuppressive response in malignant and TME cells, respectively. Cell type-specific transcriptional dynamics revealed that malignant and endothelial cells promote the cellular proliferations of TME cells, whereas the fibroblasts and immune cells are associated with procancer and anticancer immunity, respectively. Ligand-receptor analysis revealed that CCL2-expressing fibroblasts promote the tumour progression via JAK-STAT3 signalling and inflammatory response in tumour-infiltrated macrophages. CCL2+ fibroblasts and STAT3-activated macrophages are co-localised and their co-abundance was associated with unfavourable prognosis. We experimentally validated that CCL2+ fibroblasts recruit myeloid cells and stimulate STAT3 activation in recruited macrophages. The development of immunosuppressive TME by CCL2+ fibroblasts were also validated in syngeneic mouse models.

CONCLUSION

GC spatial transcriptomes revealed functional cellular crosstalk involving multiple cell types among which the interaction between CCL2+ fibroblasts and STAT3-activated macrophages plays roles in establishing immune-suppressive GC TME with potential clinical relevance.

Deep learning in integrating spatial transcriptomics with other modalities

Spatial transcriptomics technologies have been extensively applied in biological research, enabling the study of transcriptome while preserving the spatial context of tissues. Paired with spatial transcriptomics data, platforms often provide histology and (or) chromatin images, which capture cellular morphology and chromatin organization. Additionally, single-cell RNA sequencing (scRNA-seq) data from matching tissues often accompany spatial data, offering a transcriptome-wide gene expression profile of individual cells. Integrating such additional data from other modalities can effectively enhance spatial transcriptomics data, and, conversely, spatial transcriptomics data can supplement scRNA-seq with spatial information. Moreover, the rapid development of spatial multi-omics technology has spurred the demand for the integration of spatial multi-omics data to present a more detailed molecular landscape within tissues. Numerous deep learning (DL) methods have been developed for integrating spatial transcriptomics with other modalities. However, a comprehensive review of DL approaches for integrating spatial transcriptomics data with other modalities remains absent. In this study, we systematically review the applications of DL in integrating spatial transcriptomics data with other modalities. We first delineate the DL techniques applied in this integration and the key tasks involved. Next, we detail these methods and categorize them based on integrated modality and key task. Furthermore, we summarize the integration strategies of these integration methods. Finally, we discuss the challenges and future directions in integrating spatial transcriptomics with other modalities, aiming to facilitate the development of robust computational methods that more comprehensively exploit multimodal information.

Clustering-independent estimation of cell abundances in bulk tissues using single-cell RNA-seq data

Single-cell RNA sequencing has transformed the study of biological tissues by enabling transcriptomic characterizations of their constituent cell states. Computational methods for gene expression deconvolution use this information to infer the cell composition of related tissues profiled at the bulk level. However, current deconvolution methods are restricted to discrete cell types and have limited power to make inferences about continuous cellular processes such as cell differentiation or immune cell activation. We present ConDecon, a clustering-independent method for inferring the likelihood for each cell in a single-cell dataset to be present in a bulk tissue. ConDecon represents an improvement in phenotypic resolution and functionality with respect to regression-based methods. Using ConDecon, we discover the implication of neurodegenerative microglia inflammatory pathways in the mesenchymal transformation of pediatric ependymoma and characterize their spatial trajectories of activation. The generality of this approach enables the deconvolution of other data modalities, such as bulk ATAC-seq data.

Non-coding RNAs in oral cancer: Emerging biomarkers and therapeutic frontier

Around the world, oral cancer (OC) is a major public health problem, resulting in a significant number of deaths each year. Early detection and treatment are crucial for improving patient outcomes. Recent progress in DNA sequencing and transcriptome profiling has revealed extensive non-coding RNAs (ncRNAs) transcription, underscoring their regulatory importance. NcRNAs influence genomic transcription and translation and molecular signaling pathways, making them valuable for various clinical applications. Combining spatial transcriptomics (ST) and spatial metabolomics (SM) with single-cell RNA sequencing provides deeper insights into tumor microenvironments, enhancing diagnostic and therapeutic precision for OC. Additionally, the exploration of salivary biomarkers offers a non-invasive diagnostic avenue. This article explores the potential of ncRNAs as diagnostic and therapeutic tools for OC.

Integration of Single-Cell and Spatial Transcriptomic Data Reveals Spatial Architecture and Potential Biomarkers in Alzheimer's Disease

Alzheimer's disease (AD) is a complex neurodegenerative disorder characterized by the gradual loss of neurons and the accumulation of amyloid plaques and neurofibrillary tangles. Despite advancements in the understanding of AD's pathophysiology, the cellular organization and interactions in the prefrontal cortex (PFC) remain elusive. Eight single-cell RNA sequencing (scRNA-seq) datasets from both normal controls and individuals with AD were harmonized. Stringent preprocessing protocols were implemented to uphold dataset integrity. Unsupervised clustering and annotation revealed 22 distinct cell clusters corresponding to 19 unique cell types. The spatial architecture of the PFC region was constructed using the CARD tool. Further analyses encompassed trajectory examination of Oligodendrocyte subtypes, evaluation of regulon activity scores, and spot clustering within white matter regions (WM). Differential expression analysis and functional enrichment assays unveiled molecular signatures linked to AD progression and were validated using microarray data sourced from neurodegenerative disorder patients. Our investigation employs scRNA-seq and spatial transcriptomics to uncover the cellular atlas and spatial architecture of the human PFC in AD. Moreover, our results indicate that Oligodendrocytes are more prevalent in AD patients, showcasing diverse subtypes and spatial organization within WM regions. Each subtype appears to be associated with distinct biological processes and transcriptional regulators, shedding light on their involvement in AD pathology. Notably, the Oligodendrocyte_C6 subtype is linked to neurological damage in AD patients, characterized by heightened expression of genes involved in cell-cell connections, cell membrane stability, and myelination. Additionally, 12 target genes regulated by NFIA were identified, which are upregulated in AD patients and associated with disease progression. Elevated PLXDC2 expression in peripheral blood was also identified, suggesting its potential as a non-invasive biomarker for early AD detection. Our study provides novel insights into the role of Oligodendrocytes in AD and highlights the potential of PLXDC2 as a blood biomarker for non-invasive diagnosis and monitoring of AD patients.

Rapid response of lichen planus to baricitinib associated with suppression of cytotoxic CXCL13+CD8+ T cells

BACKGROUNDCutaneous lichen planus (LP) is a recalcitrant, difficult-to-treat, inflammatory skin disease characterized by pruritic, flat-topped, violaceous papules on the skin. Baricitinib is an oral Janus kinase (JAK) 1/2 inhibitor that interrupts the signaling pathway of IFN-γ, a cytokine implicated in the pathogenesis of LP.METHODSIn this phase II trial, 12 patients with cutaneous LP received 2 mg daily baricitinib for 16 weeks, accompanied by in-depth spatial, single-cell, and bulk transcriptomic profiling of pre- and posttreatment samples.RESULTSAn early and sustained clinical response was seen, with 83.3% of patients responsive at week 16. Our molecular data identified a unique, oligoclonal IFN-γ, CD8+, and CXCL13+ cytotoxic T cell population in LP skin and demonstrated a rapid decrease in IFN signature within 2 weeks of treatment, most prominently in the basal layer of the epidermis.CONCLUSIONThis study demonstrates the efficacy and molecular mechanisms of JAK inhibition in LP.TRIAL REGISTRATIONNCT05188521FUNDINGEli Lilly, Appignani Benefactor Funds, 5P30AR075043, Mayo Clinic Clinical Trials Stimulus Funds.

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.

Dual decoding of cell types and gene expression in spatial transcriptomics with PANDA

Sequencing-based spatial transcriptomics technologies have revolutionized our understanding of complex biological systems by enabling transcriptome profiling while preserving spatial context. However, spot-level expression measurements often amalgamate signals from diverse cells, obscuring potential heterogeneity. Existing methods aim to deconvolute spatial transcriptomics data into cell type proportions for each spot using single-cell RNA sequencing references but overlook cell-type-specific gene expression, essential for uncovering intra-type heterogeneity. We present PANDA (ProbAbilistic-based decoNvolution with spot-aDaptive cell type signAtures), a novel method that concurrently deciphers spot-level gene expression into both cell type proportions and cell-type-specific gene expression. PANDA integrates archetypal analysis to capture within-cell-type heterogeneity and dynamically learns cell type signatures for each spot during deconvolution. Simulations demonstrate PANDA's superior performance. Applied to real spatial transcriptomics data from diverse tissues, including tumor, brain, and developing heart, PANDA reconstructs spatial structures and reveals subtle transcriptional variations within specific cell types, offering a comprehensive understanding of tissue dynamics.

scDOT: optimal transport for mapping senescent cells in spatial transcriptomics

The low resolution of spatial transcriptomics data necessitates additional information for optimal use. We developed scDOT, which combines spatial transcriptomics and single cell RNA sequencing to improve the ability to reconstruct single cell resolved spatial maps and identify senescent cells. scDOT integrates optimal transport and expression deconvolution to learn non-linear couplings between cells and spots and to infer cell placements. Application of scDOT to lung spatial transcriptomics data improves on prior methods and allows the identification of the spatial organization of senescent cells, their neighboring cells and novel genes involved in cell-cell interactions that may be driving senescence.

eMCI: An Explainable Multimodal Correlation Integration Model for Unveiling Spatial Transcriptomics and Intercellular Signaling

Current integration methods for single-cell RNA sequencing (scRNA-seq) data and spatial transcriptomics (ST) data are typically designed for specific tasks, such as deconvolution of cell types or spatial distribution prediction of RNA transcripts. These methods usually only offer a partial analysis of ST data, neglecting the complex relationship between spatial expression patterns underlying cell-type specificity and intercellular cross-talk. Here, we present eMCI, an explainable multimodal correlation integration model based on deep neural network framework. eMCI leverages the fusion of scRNA-seq and ST data using different spot-cell correlations to integrate multiple synthetic analysis tasks of ST data at cellular level. First, eMCI can achieve better or comparable accuracy in cell-type classification and deconvolution according to wide evaluations and comparisons with state-of-the-art methods on both simulated and real ST datasets. Second, eMCI can identify key components across spatial domains responsible for different cell types and elucidate the spatial expression patterns underlying cell-type specificity and intercellular communication, by employing an attribution algorithm to dissect the visual input. Especially, eMCI has been applied to 3 cross-species datasets, including zebrafish melanomas, soybean nodule maturation, and human embryonic lung, which accurately and efficiently estimate per-spot cell composition and infer proximal and distal cellular interactions within the spatial and temporal context. In summary, eMCI serves as an integrative analytical framework to better resolve the spatial transcriptome based on existing single-cell datasets and elucidate proximal and distal intercellular signal transduction mechanisms over spatial domains without requirement of biological prior reference. This approach is expected to facilitate the discovery of spatial expression patterns of potential biomolecules with cell type and cell-cell communication specificity.

Age, sex, and cell type-resolved hypothalamic gene expression across the pubertal transition in mice

BACKGROUND

The hypothalamus plays a central role in regulating puberty. However, our knowledge of the postnatal gene regulatory networks that control the pubertal transition in males and females is incomplete. Here, we investigate the age-, sex- and cell-type-specific gene regulation in the hypothalamus across the pubertal transition.

METHODS

We used RNA-seq to profile hypothalamic gene expression in male and female mice at five time points spanning the onset of puberty (postnatal days (PD) 12, 22, 27, 32, and 37). By combining this data with hypothalamic single nuclei RNA-seq data from pre- and postpubertal mice, we assigned gene expression changes to their most likely cell types of origin. In our colony, pubertal onset occurs earlier in male mice, allowing us to focus on genes whose expression is dynamic across ages and offset between sexes, and to explore the bases of sex effects.

RESULTS

Our age-by-sex pattern of expression enriched for biological pathways involved hormone production, neuronal activation, and glial maturation. Additionally, we inferred a robust expansion of oligodendrocytes precursor cells into mature oligodendrocytes spanning the prepubertal (PD12) to peri-pubertal (PD27) timepoints. Using spatial transcriptomic data from postpubertal mice, we observed the lateral hypothalamic area and zona incerta were the most oligodendrocyte-rich regions and that these cells expressed genes known to be involved in pubertal regulation.

CONCLUSION

Together, by incorporating multiple biological timepoints and using sex as a variable, we identified gene and cell-type changes that may participate in orchestrating the pubertal transition and provided a resource for future studies of postnatal hypothalamic gene regulation.