Spatial deconvolution of HER2-positive breast cancer delineates tumor-associated cell type interactions

Applications of single‑cell omics and spatial transcriptomics technologies in gastric cancer (Review)

Gastric cancer (GC) is a prominent contributor to global cancer-related mortalities, and a deeper understanding of its molecular characteristics and tumor heterogeneity is required. Single-cell omics and spatial transcriptomics (ST) technologies have revolutionized cancer research by enabling the exploration of cellular heterogeneity and molecular landscapes at the single-cell level. In the present review, an overview of the advancements in single-cell omics and ST technologies and their applications in GC research is provided. Firstly, multiple single-cell omics and ST methods are discussed, highlighting their ability to offer unique insights into gene expression, genetic alterations, epigenomic modifications, protein expression patterns and cellular location in tissues. Furthermore, a summary is provided of key findings from previous research on single-cell omics and ST methods used in GC, which have provided valuable insights into genetic alterations, tumor diagnosis and prognosis, tumor microenvironment analysis, and treatment response. In summary, the application of single-cell omics and ST technologies has revealed the levels of cellular heterogeneity and the molecular characteristics of GC, and holds promise for improving diagnostics, personalized treatments and patient outcomes in GC.

Deciphering the spatial landscape and plasticity of immunosuppressive fibroblasts in breast cancer

Although heterogeneity of FAP+ Cancer-Associated Fibroblasts (CAF) has been described in breast cancer, their plasticity and spatial distribution remain poorly understood. Here, we analyze trajectory inference, deconvolute spatial transcriptomics at single-cell level and perform functional assays to generate a high-resolution integrated map of breast cancer (BC), with a focus on inflammatory and myofibroblastic (iCAF/myCAF) FAP+ CAF clusters. We identify 10 spatially-organized FAP+ CAF-related cellular niches, called EcoCellTypes, which are differentially localized within tumors. Consistent with their spatial organization, cancer cells drive the transition of detoxification-associated iCAF (Detox-iCAF) towards immunosuppressive extracellular matrix (ECM)-producing myCAF (ECM-myCAF) via a DPP4- and YAP-dependent mechanism. In turn, ECM-myCAF polarize TREM2+ macrophages, regulatory NK and T cells to induce immunosuppressive EcoCellTypes, while Detox-iCAF are associated with FOLR2+ macrophages in an immuno-protective EcoCellType. FAP+ CAF subpopulations accumulate differently according to the invasive BC status and predict invasive recurrence of ductal carcinoma in situ (DCIS), which could help in identifying low-risk DCIS patients eligible for therapeutic de-escalation.

Dickkopf-1 drives perineural invasion via PI3K-AKT signaling pathway in head and neck squamous cancer

Perineural invasion (PNI) leads to the poor prognosis of head and neck squamous cancer (HNSCC) patients, but the mechanism of PNI remains unclear. Dickkopf-1 (DKK1), a secretory protein in the Wnt signaling pathway, was found indeed upregulated in HNSCC cells and tissues. Higher expression of DKK1 was statistically relevant to T stage, N stage, PNI, and poor prognosis of HNSCC. DKK1 overexpression enhanced the migration abilities of cancer cells. Moreover, DKK1-overexpressing cancer cells promoted cancer cells invasion of peripheral nerves in vitro and in vivo. Mechanistically, DKK1 could promote the PI3K-AKT signaling pathway. The migration abilities of neuroblastoma cells, which were enhanced by DKK1-overexpressing HNSCC cell lines, could be reversed by an inhibitor of Akt (MK2206). The association of DKK1 with PNI was also confirmed in HNSCC samples. Variables, including T stage, N stage, DKK1 expression, and PNI, were used to establish a nomogram to predict the survival probability and disease-free probability at 3 and 5 years. In summary, DKK1 can promote the PI3K-AKT signaling pathway in tumor cells and then could induce neuritogenesis and facilitate PNI. MK2206 may be a potential therapeutic target drug for HNSCC patients with PNI.

Benchmarking spatial clustering methods with spatially resolved transcriptomics data

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

Emerging measurements for tumor-infiltrating lymphocytes in breast cancer

Tumor-infiltrating lymphocytes are a general term for lymphocytes or immune cells infiltrating the tumor microenvironment. Numerous studies have demonstrated tumor-infiltrating lymphocytes to be robust prognostic and predictive biomarkers in breast cancer. Recently, immune checkpoint inhibitors, which directly target tumor-infiltrating lymphocytes, have become part of standard of care treatment for triple-negative breast cancer. Surprisingly, tumor-infiltrating lymphocytes quantified by conventional methods do not predict response to immune checkpoint inhibitors, which highlights the heterogeneity of tumor-infiltrating lymphocytes and the complexity of the immune network in the tumor microenvironment. Tumor-infiltrating lymphocytes are composed of diverse immune cell populations, including cytotoxic CD8-positive T lymphocytes, B cells and myeloid cells. Traditionally, tumor-infiltrating lymphocytes in tumor stroma have been evaluated by histology. However, the standardization of this approach is limited, necessitating the use of various novel technologies to elucidate the heterogeneity in the tumor microenvironment. This review outlines the evaluation methods for tumor-infiltrating lymphocytes from conventional pathological approaches that evaluate intratumoral and stromal tumor-infiltrating lymphocytes such as immunohistochemistry, to the more recent advancements in computer tissue imaging using artificial intelligence, flow cytometry sorting and multi-omics analyses using high-throughput assays to estimate tumor-infiltrating lymphocytes from bulk tumor using immune signatures or deconvolution tools. We also discuss higher resolution technologies that enable the analysis of tumor-infiltrating lymphocytes heterogeneity such as single-cell analysis and spatial transcriptomics. As we approach the era of personalized medicine, it is important for clinicians to understand these technologies.

Recent advances in differential expression analysis for single-cell RNA-seq and spatially resolved transcriptomic studies

Differential expression (DE) analysis is a necessary step in the analysis of single-cell RNA sequencing (scRNA-seq) and spatially resolved transcriptomics (SRT) data. Unlike traditional bulk RNA-seq, DE analysis for scRNA-seq or SRT data has unique characteristics that may contribute to the difficulty of detecting DE genes. However, the plethora of DE tools that work with various assumptions makes it difficult to choose an appropriate one. Furthermore, a comprehensive review on detecting DE genes for scRNA-seq data or SRT data from multi-condition, multi-sample experimental designs is lacking. To bridge such a gap, here, we first focus on the challenges of DE detection, then highlight potential opportunities that facilitate further progress in scRNA-seq or SRT analysis, and finally provide insights and guidance in selecting appropriate DE tools or developing new computational DE methods.

Adoptive T cell therapy for solid tumors: current landscape and future challenges

Adoptive cell therapy (ACT) comprises different strategies to enhance the activity of T lymphocytes and other effector cells that orchestrate the antitumor immune response, including chimeric antigen receptor (CAR) T-cell therapy, T-cell receptor (TCR) gene-modified T cells, and therapy with tumor-infiltrating lymphocytes (TILs). The outstanding results of CAR-T cells in some hematologic malignancies have launched the investigation of ACT in patients with refractory solid malignancies. However, certain characteristics of solid tumors, such as their antigenic heterogeneity and immunosuppressive microenvironment, hamper the efficacy of antigen-targeted treatments. Other ACT modalities, such as TIL therapy, have emerged as promising new strategies. TIL therapy has shown safety and promising activity in certain immunogenic cancers, mainly advanced melanoma, with an exciting rationale for its combination with immune checkpoint inhibitors. However, the implementation of TIL therapy in clinical practice is hindered by several biological, logistic, and economic challenges. In this review, we aim to summarize the current knowledge, available clinical results, and potential areas of future research regarding the use of T cell therapy in patients with solid tumors.

Advanced insights on tumor-associated macrophages revealed by single-cell RNA sequencing: The intratumor heterogeneity, functional phenotypes, and cellular interactions

Single-cell RNA sequencing (scRNA-seq) is an emerging technology used for cellular transcriptome analysis. The application of scRNA-seq has led to profoundly advanced oncology research, continuously optimizing novel therapeutic strategies. Intratumor heterogeneity extensively consists of all tumor components, contributing to different tumor behaviors and treatment responses. Tumor-associated macrophages (TAMs), the core immune cells linking innate and adaptive immunity, play significant roles in tumor progression and resistance to therapies. Moreover, dynamic changes occur in TAM phenotypes and functions subject to the regulation of the tumor microenvironment. The heterogeneity of TAMs corresponding to the state of the tumor microenvironment has been comprehensively recognized using scRNA-seq. Herein, we reviewed recent research and summarized variations in TAM phenotypes and functions from a developmental perspective to better understand the significance of TAMs in the tumor microenvironment.

Spatial transcriptomics and the anatomical pathologist: Molecular meets morphology

In recent years anatomical pathology has been revolutionised by the incorporation of molecular findings into routine diagnostic practice, and in some diseases the presence of specific molecular alterations are now essential for diagnosis. Spatial transcriptomics describes a group of technologies that provide up to transcriptome-wide expression profiling while preserving the spatial origin of the data, with many of these technologies able to provide these data using a single tissue section. Spatial transcriptomics allows expression profiling of highly specific areas within a tissue section potentially to subcellular resolution, and allows correlation of expression data with morphology, tissue type and location relative to other structures. While largely still research laboratory-based, several spatial transcriptomics methods have now achieved compatibility with formalin-fixed paraffin-embedded tissue (FFPE), allowing their use in diagnostic tissue samples, and with further development potentially leading to their incorporation in routine anatomical pathology practice. This mini review provides an overview of spatial transcriptomics methods, with an emphasis on platforms compatible with FFPE tissue, approaches to assess the data and potential applications in anatomical pathology practice.

Spatial transcriptomics reveals prognosis-associated cellular heterogeneity in the papillary thyroid carcinoma microenvironment

BACKGROUND

Papillary thyroid carcinoma (PTC) is the most common malignant endocrine tumour, and its incidence and prevalence are increasing considerably. Cellular heterogeneity in the tumour microenvironment is important for PTC prognosis. Spatial transcriptomics is a powerful technique for cellular heterogeneity study.

METHODS

In conjunction with a clinical pathologist identification method, spatial transcriptomics was employed to characterise the spatial location and RNA profiles of PTC-associated cells within the tissue sections. The spatial RNA-clinical signature genes for each cell type were extracted and applied to outlining the distribution regions of specific cells on the entire section. The cellular heterogeneity of each cell type was further revealed by ContourPlot analysis, monocle analysis, trajectory analysis, ligand-receptor analysis and Gene Ontology enrichment analysis.

RESULTS

The spatial distribution region of tumour cells, typical and atypical follicular cells (FCs and AFCs) and immune cells were accurately and comprehensively identified in all five PTC tissue sections. AFCs were identified as a transitional state between FCs and tumour cells, exhibiting a higher resemblance to the latter. Three tumour foci were shared among all patients out of the 13 observed. Notably, tumour foci No. 2 displayed elevated expression levels of genes associated with lower relapse-free survival in PTC patients. We discovered key ligand-receptor interactions, including LAMB3-ITGA2, FN1-ITGA3 and FN1-SDC4, involved in the transition of PTC cells from FCs to AFCs and eventually to tumour cells. High expression of these patterns correlated with reduced relapse-free survival. In the tumour immune microenvironment, reduced interaction between myeloid-derived TGFB1 and TGFBR1 in tumour focus No. 2 contributed to tumourigenesis and increased heterogeneity. The spatial RNA-clinical analysis method developed here revealed prognosis-associated cellular heterogeneity in the PTC microenvironment.

CONCLUSIONS

The occurrence of tumour foci No. 2 and three enhanced ligand-receptor interactions in the AFC area/tumour foci reduced the relapse-free survival of PTC patients, potentially leading to improved prognostic strategies and targeted therapies for PTC patients.

Spatial Transcriptomic and Metabolomic Landscapes of Oral Submucous Fibrosis-Derived Oral Squamous Cell Carcinoma and its Tumor Microenvironment

In South and Southeast Asia, the habit of chewing betel nuts is prevalent, which leads to oral submucous fibrosis (OSF). OSF is a well-established precancerous lesion, and a portion of OSF cases eventually progress to oral squamous cell carcinoma (OSCC). However, the specific molecular mechanisms underlying the malignant transformation of OSCC from OSF are poorly understood. In this study, the leading-edge techniques of Spatial Transcriptomics (ST) and Spatial Metabolomics (SM) are integrated to obtain spatial location information of cancer cells, fibroblasts, and immune cells, as well as the transcriptomic and metabolomic landscapes in OSF-derived OSCC tissues. This work reveals for the first time that some OSF-derived OSCC cells undergo partial epithelial-mesenchymal transition (pEMT) within the in situ carcinoma (ISC) region, eventually acquiring fibroblast-like phenotypes and participating in collagen deposition. Complex interactions among epithelial cells, fibroblasts, and immune cells in the tumor microenvironment are demonstrated. Most importantly, significant metabolic reprogramming in OSF-derived OSCC, including abnormal polyamine metabolism, potentially playing a pivotal role in promoting tumorigenesis and immune evasion is discovered. The ST and SM data in this study shed new light on deciphering the mechanisms of OSF-derived OSCC. The work also offers invaluable clues for the prevention and treatment of OSCC.

Mapping Spatiotemporal Heterogeneity in Multifocal Breast Tumor Progression by Noninvasive Ultrasound Elastography-Guided Mass Spectrometry Imaging Strategy

Spatiotemporal heterogeneity of tumors provides an escape mechanism for breast cancer cells, which can obstruct the investigation of tumor progression. While molecular profiling obtained from mass spectrometry imaging (MSI) is rich in biochemical information, it lacks the capacity for in vivo analysis. Ultrasound diagnosis has a high diagnostic accuracy but low chemical specificity. Here, we describe a noninvasive ultrasound elastography (UE)-guided MSI strategy (UEg-MSI) that integrates physical and biochemical characteristics of tumors acquired from both in vivo and in vitro imaging. Using UEg-MSI, both elasticity histopathology metabolism "fingerprints" and reciprocal crosstalk are revealed, indicating the intact, multifocal spatiotemporal heterogeneity of spontaneous tumorigenesis of the breast from early, middle, and late stages. Our results demonstrate a gradual increase in malignant degree of primary focus in cervical and thoracic mammary glands. This progression is characterized by increased stiffness according to elasticity scores, histopathological changes from hyperplasia to increased nests of neoplastic cells and necrotic areas, and regional metabolic heterogeneity and reprogramming at the spatiotemporal level. De novo fatty acid (FA) synthesis focused on independent (such as ω-9 FAs) and dependent (such as ω-6 FAs) dietary FA intake in the core cancerous nest areas in the middle and late stages of tumor or in the peripheral microareas in the early stage of the tumor. SM-Cer signaling pathway and GPs biosynthesis and degradation, as well as glycerophosphoinositol intensity, changed in multiple characteristic microareas. The UEg-MSI strategy holds the potential to expand MSI applications and enhance ultrasound-mediated cancer diagnosis. It offers new insight into early cancer discovery and the occurrence of metastasis.

Single-cell colocalization analysis using a deep generative model

Analyzing colocalization of single cells with heterogeneous molecular phenotypes is essential for understanding cell-cell interactions, and cellular responses to external stimuli and their biological functions in diseases and tissues. However, existing computational methodologies identified the colocalization patterns between predefined cell populations, which can obscure the molecular signatures arising from intercellular communication. Here, we introduce DeepCOLOR, a computational framework based on a deep generative model that recovers intercellular colocalization networks with single-cell resolution by the integration of single-cell and spatial transcriptomes. Along with colocalized population detection accuracy that is superior to existing methods in simulated dataset, DeepCOLOR identified plausible cell-cell interaction candidates between colocalized single cells and segregated cell populations defined by the colocalization relationships in mouse brain tissues, human squamous cell carcinoma samples, and human lung tissues infected with SARS-CoV-2. DeepCOLOR is applicable to studying cell-cell interactions behind various spatial niches. A record of this paper's transparent peer review process is included in the supplemental information.

Immune cell networking in solid tumors: focus on macrophages and neutrophils

The tumor microenvironment is composed of tumor cells, stromal cells and leukocytes, including innate and adaptive immune cells, and represents an ecological niche that regulates tumor development and progression. In general, inflammatory cells are considered to contribute to tumor progression through various mechanisms, including the formation of an immunosuppressive microenvironment. Macrophages and neutrophils are important components of the tumor microenvironment and can act as a double-edged sword, promoting or inhibiting the development of the tumor. Targeting of the immune system is emerging as an important therapeutic strategy for cancer patients. However, the efficacy of the various immunotherapies available is still limited. Given the crucial importance of the crosstalk between macrophages and neutrophils and other immune cells in the formation of the anti-tumor immune response, targeting these interactions may represent a promising therapeutic approach against cancer. Here we will review the current knowledge of the role played by macrophages and neutrophils in cancer, focusing on their interaction with other immune cells.

Breast Cancer Immunity: It is TIME for the Next Chapter

Our ability to interrogate the tumor immune microenvironment (TIME) at an ever-increasing granularity has uncovered critical determinants of disease progression. Not only do we now have a better understanding of the immune response in breast cancer, but it is becoming possible to leverage key mechanisms to effectively combat this disease. Almost every component of the immune system plays a role in enabling or inhibiting breast tumor growth. Building on early seminal work showing the involvement of T cells and macrophages in controlling breast cancer progression and metastasis, single-cell genomics and spatial proteomics approaches have recently expanded our view of the TIME. In this article, we provide a detailed description of the immune response against breast cancer and examine its heterogeneity in disease subtypes. We discuss preclinical models that enable dissecting the mechanisms responsible for tumor clearance or immune evasion and draw parallels and distinctions between human disease and murine counterparts. Last, as the cancer immunology field is moving toward the analysis of the TIME at the cellular and spatial levels, we highlight key studies that revealed previously unappreciated complexity in breast cancer using these technologies. Taken together, this article summarizes what is known in breast cancer immunology through the lens of translational research and identifies future directions to improve clinical outcomes.

Single-cell RNA-seq reveals the metabolic status of immune cells response to immunotherapy in triple-negative breast cancer

Immune checkpoint blockade (ICB) therapy offers promise in the treatment of triple-negative breast cancer (TNBC); however, its limited efficacy in certain TNBC patients poses a challenge. In this study, we elucidated the metabolic mechanism at 'sub-subtype' resolution underlying the non-response to ICB therapy in TNBC. Here, an analytic pipeline was developed to reveal the metabolic heterogeneity, which is correlated with the ICB outcomes, within each immune cell subtype. First, we identified metabolic 'sub-subtypes' within certain cell subtypes, predominantly T cell subsets, which are enriched in ICB non-responders and named as non-responder-enriched (NR-E) clusters. Notably, most of NR-E T metabolic cells exhibit globally higher metabolic activities compared to other cells within the same individual subtype. Further, we investigated the extra-cellular signals that trigger the metabolic status of NR-E T cells. In detail, the prediction of cell-to-cell communication indicated that NR-E T cells are regulated by plasmatic dendritic cells (pDCs) through TNFSF9, as well as by macrophages expressing SIGLEC9. In addition, we also validate the communication between TNFSF9+ pDCs and NR-E T cells utilizing deconvolution of spatial transcriptomics analysis. In summary, our research identified specific metabolic 'sub-subtypes' associated with ICB non-response and uncovered the mechanisms of their regulation in TNBC. And the proposed analytical pipeline can be used to examine metabolic heterogeneity within cell types that correlate with diverse phenotypes.

Mapping cancer biology in space: applications and perspectives on spatial omics for oncology

Technologies to decipher cellular biology, such as bulk sequencing technologies and single-cell sequencing technologies, have greatly assisted novel findings in tumor biology. Recent findings in tumor biology suggest that tumors construct architectures that influence the underlying cancerous mechanisms. Increasing research has reported novel techniques to map the tissue in a spatial context or targeted sampling-based characterization and has introduced such technologies to solve oncology regarding tumor heterogeneity, tumor microenvironment, and spatially located biomarkers. In this study, we address spatial technologies that can delineate the omics profile in a spatial context, novel findings discovered via spatial technologies in oncology, and suggest perspectives regarding therapeutic approaches and further technological developments.

Macrophage barrier in the tumor microenvironment and potential clinical applications

The tumor microenvironment (TME) constitutes a complex microenvironment comprising a diverse array of immune cells and stromal components. Within this intricate context, tumor-associated macrophages (TAMs) exhibit notable spatial heterogeneity. This heterogeneity contributes to various facets of tumor behavior, including immune response modulation, angiogenesis, tissue remodeling, and metastatic potential. This review summarizes the spatial distribution of macrophages in both the physiological environment and the TME. Moreover, this paper explores the intricate interactions between TAMs and diverse immune cell populations (T cells, dendritic cells, neutrophils, natural killer cells, and other immune cells) within the TME. These bidirectional exchanges form a complex network of immune interactions that influence tumor immune surveillance and evasion strategies. Investigating TAM heterogeneity and its intricate interactions with different immune cell populations offers potential avenues for therapeutic interventions. Additionally, this paper discusses therapeutic strategies targeting macrophages, aiming to uncover novel approaches for immunotherapy. Video Abstract.

Spatially contrastive variational autoencoder for deciphering tissue heterogeneity from spatially resolved transcriptomics

Recent advances in spatially resolved transcriptomics (SRT) have brought ever-increasing opportunities to characterize expression landscape in the context of tissue spatiality. Nevertheless, there still exist multiple challenges to accurately detect spatial functional regions in tissue. Here, we present a novel contrastive learning framework, SPAtially Contrastive variational AutoEncoder (SpaCAE), which contrasts transcriptomic signals of each spot and its spatial neighbors to achieve fine-grained tissue structures detection. By employing a graph embedding variational autoencoder and incorporating a deep contrastive strategy, SpaCAE achieves a balance between spatial local information and global information of expression, enabling effective learning of representations with spatial constraints. Particularly, SpaCAE provides a graph deconvolutional decoder to address the smoothing effect of local spatial structure on expression's self-supervised learning, an aspect often overlooked by current graph neural networks. We demonstrated that SpaCAE could achieve effective performance on SRT data generated from multiple technologies for spatial domains identification and data denoising, making it a remarkable tool to obtain novel insights from SRT studies.

Profiling the heterogeneity of colorectal cancer consensus molecular subtypes using spatial transcriptomics

The consensus molecular subtypes (CMS) of colorectal cancer (CRC) is the most widely-used gene expression-based classification and has contributed to a better understanding of disease heterogeneity and prognosis. Nevertheless, CMS intratumoral heterogeneity restricts its clinical application, stressing the necessity of further characterizing the composition and architecture of CRC. Here, we used Spatial Transcriptomics (ST) in combination with single-cell RNA sequencing (scRNA-seq) to decipher the spatially resolved cellular and molecular composition of CRC. In addition to mapping the intratumoral heterogeneity of CMS and their microenvironment, we identified cell communication events in the tumor-stroma interface of CMS2 carcinomas. This includes tumor growth-inhibiting as well as -activating signals, such as the potential regulation of the ETV4 transcriptional activity by DCN or the PLAU-PLAUR ligand-receptor interaction. Our study illustrates the potential of ST to resolve CRC molecular heterogeneity and thereby help advance personalized therapy.

STEM enables mapping of single-cell and spatial transcriptomics data with transfer learning

Profiling spatial variations of cellular composition and transcriptomic characteristics is important for understanding the physiology and pathology of tissues. Spatial transcriptomics (ST) data depict spatial gene expression but the currently dominating high-throughput technology is yet not at single-cell resolution. Single-cell RNA-sequencing (SC) data provide high-throughput transcriptomic information at the single-cell level but lack spatial information. Integrating these two types of data would be ideal for revealing transcriptomic landscapes at single-cell resolution. We develop the method STEM (SpaTially aware EMbedding) for this purpose. It uses deep transfer learning to encode both ST and SC data into a unified spatially aware embedding space, and then uses the embeddings to infer SC-ST mapping and predict pseudo-spatial adjacency between cells in SC data. Semi-simulation and real data experiments verify that the embeddings preserved spatial information and eliminated technical biases between SC and ST data. We apply STEM to human squamous cell carcinoma and hepatic lobule datasets to uncover the localization of rare cell types and reveal cell-type-specific gene expression variation along a spatial axis. STEM is powerful for mapping SC and ST data to build single-cell level spatial transcriptomic landscapes, and can provide mechanistic insights into the spatial heterogeneity and microenvironments of tissues.

A TCF4-dependent gene regulatory network confers resistance to immunotherapy in melanoma

To better understand intrinsic resistance to immune checkpoint blockade (ICB), we established a comprehensive view of the cellular architecture of the treatment-naive melanoma ecosystem and studied its evolution under ICB. Using single-cell, spatial multi-omics, we showed that the tumor microenvironment promotes the emergence of a complex melanoma transcriptomic landscape. Melanoma cells harboring a mesenchymal-like (MES) state, a population known to confer resistance to targeted therapy, were significantly enriched in early on-treatment biopsies from non-responders to ICB. TCF4 serves as the hub of this landscape by being a master regulator of the MES signature and a suppressor of the melanocytic and antigen presentation transcriptional programs. Targeting TCF4 genetically or pharmacologically, using a bromodomain inhibitor, increased immunogenicity and sensitivity of MES cells to ICB and targeted therapy. We thereby uncovered a TCF4-dependent regulatory network that orchestrates multiple transcriptional programs and contributes to resistance to both targeted therapy and ICB in melanoma.

Inferring super-resolution tissue architecture by integrating spatial transcriptomics with histology

Spatial transcriptomics (ST) has demonstrated enormous potential for generating intricate molecular maps of cells within tissues. Here we present iStar, a method based on hierarchical image feature extraction that integrates ST data and high-resolution histology images to predict spatial gene expression with super-resolution. Our method enhances gene expression resolution to near-single-cell levels in ST and enables gene expression prediction in tissue sections where only histology images are available.

Strategies for studying immune and non-immune human and canine mammary gland cancer tumour infiltrate

The tumour microenvironment (TME) is usually defined as a cell environment associated with tumours or cancerous stem cells where conditions are established affecting tumour development and progression through malignant cell interaction with non-malignant cells. The TME is made up of endothelial, immune and non-immune cells, extracellular matrix (ECM) components and signalling molecules acting specifically on tumour and non-tumour cells. Breast cancer (BC) is the commonest malignant neoplasm worldwide and the main cause of mortality in women globally; advances regarding BC study and understanding it are relevant for acquiring novel, personalised therapeutic tools. Studying canine mammary gland tumours (CMGT) is one of the most relevant options for understanding BC using animal models as they share common epidemiological, clinical, pathological, biological, environmental, genetic and molecular characteristics with human BC. In-depth, detailed investigation regarding knowledge of human BC-related TME and in its canine model is considered extremely relevant for understanding changes in TME composition during tumour development. This review addresses important aspects concerned with different methods used for studying BC- and CMGT-related TME that are important for developing new and more effective therapeutic strategies for attacking a tumour during specific evolutionary stages.

Transformer with convolution and graph-node co-embedding: An accurate and interpretable vision backbone for predicting gene expressions from local histopathological image

Inferring gene expressions from histopathological images has long been a fascinating yet challenging task, primarily due to the substantial disparities between the two modality. Existing strategies using local or global features of histological images are suffering model complexity, GPU consumption, low interpretability, insufficient encoding of local features, and over-smooth prediction of gene expressions among neighboring sites. In this paper, we develop TCGN (Transformer with Convolution and Graph-Node co-embedding method) for gene expression estimation from H&E-stained pathological slide images. TCGN comprises a combination of convolutional layers, transformer encoders, and graph neural networks, and is the first to integrate these blocks in a general and interpretable computer vision backbone. Notably, TCGN uniquely operates with just a single spot image as input for histopathological image analysis, simplifying the process while maintaining interpretability. We validate TCGN on three publicly available spatial transcriptomic datasets. TCGN consistently exhibited the best performance (with median PCC 0.232). TCGN offers superior accuracy while keeping parameters to a minimum (just 86.241 million), and it consumes minimal memory, allowing it to run smoothly even on personal computers. Moreover, TCGN can be extended to handle bulk RNA-seq data while providing the interpretability. Enhancing the accuracy of omics information prediction from pathological images not only establishes a connection between genotype and phenotype, enabling the prediction of costly-to-measure biomarkers from affordable histopathological images, but also lays the groundwork for future multi-modal data modeling. Our results confirm that TCGN is a powerful tool for inferring gene expressions from histopathological images in precision health applications.

Spatial transcriptomics analysis identifies a unique tumor-promoting function of the meningeal stroma in melanoma leptomeningeal disease

Leptomeningeal disease (LMD) remains a rapidly lethal complication for late-stage melanoma patients. The inaccessible nature of the disease site and lack of understanding of the biology of this unique metastatic site are major barriers to developing efficacious therapies for patients with melanoma LMD. Here, we characterize the tumor microenvironment of the leptomeningeal tissues and patient-matched extra-cranial metastatic sites using spatial transcriptomic analyses with in vitro and in vivo validation. We show the spatial landscape of melanoma LMD to be characterized by a lack of immune infiltration and instead exhibit a higher level of stromal involvement. We show that the tumor-stroma interactions at the leptomeninges activate pathways implicated in tumor-promoting signaling, mediated through upregulation of SERPINA3 at the tumor-stroma interface. Our functional experiments establish that the meningeal stroma is required for melanoma cells to survive in the CSF environment and that these interactions lead to a lack of MAPK inhibitor sensitivity in the tumor. We show that knocking down SERPINA3 or inhibiting the downstream IGR1R/PI3K/AKT axis results in re-sensitization of the tumor to MAPK-targeting therapy and tumor cell death in the leptomeningeal environment. Our data provides a spatial atlas of melanoma LMD, identifies the tumor-promoting role of meningeal stroma, and demonstrates a mechanism for overcoming microenvironment-mediated drug resistance unique to this metastatic site.

Spatial transcriptomics of B cell and T cell receptors reveals lymphocyte clonal dynamics

The spatial distribution of lymphocyte clones within tissues is critical to their development, selection, and expansion. We have developed spatial transcriptomics of variable, diversity, and joining (VDJ) sequences (Spatial VDJ), a method that maps B cell and T cell receptor sequences in human tissue sections. Spatial VDJ captures lymphocyte clones that match canonical B and T cell distributions and amplifies clonal sequences confirmed by orthogonal methods. We found spatial congruency between paired receptor chains, developed a computational framework to predict receptor pairs, and linked the expansion of distinct B cell clones to different tumor-associated gene expression programs. Spatial VDJ delineates B cell clonal diversity and lineage trajectories within their anatomical niche. Thus, Spatial VDJ captures lymphocyte spatial clonal architecture across tissues, providing a platform to harness clonal sequences for therapy.

Benchmarking computational methods to identify spatially variable genes and peaks

Spatially resolved transcriptomics offers unprecedented insight by enabling the profiling of gene expression within the intact spatial context of cells, effectively adding a new and essential dimension to data interpretation. To efficiently detect spatial structure of interest, an essential step in analyzing such data involves identifying spatially variable genes. Despite researchers having developed several computational methods to accomplish this task, the lack of a comprehensive benchmark evaluating their performance remains a considerable gap in the field. Here, we present a systematic evaluation of 14 methods using 60 simulated datasets generated by four different simulation strategies, 12 real-world transcriptomics, and three spatial ATAC-seq datasets. We find that spatialDE2 consistently outperforms the other benchmarked methods, and Moran's I achieves competitive performance in different experimental settings. Moreover, our results reveal that more specialized algorithms are needed to identify spatially variable peaks.

Multimodal Imaging of Pancreatic Cancer Microenvironment in Response to an Antiglycolytic Drug

Lipid metabolism and glycolysis play crucial roles in the progression and metastasis of cancer, and the use of 3-bromopyruvate (3-BP) as an antiglycolytic agent has shown promise in killing pancreatic cancer cells. However, developing an effective strategy to avoid chemoresistance requires the ability to probe the interaction of cancer drugs with complex tumor-associated microenvironments (TAMs). Unfortunately, no robust and multiplexed molecular imaging technology is currently available to analyze TAMs. In this study, the simultaneous profiling of three protein biomarkers using SERS nanotags and antibody-functionalized nanoparticles in a syngeneic mouse model of pancreatic cancer (PC) is demonstrated. This allows for comprehensive information about biomarkers and TAM alterations before and after treatment. These multimodal imaging techniques include surface-enhanced Raman spectroscopy (SERS), immunohistochemistry (IHC), polarized light microscopy, second harmonic generation (SHG) microscopy, fluorescence lifetime imaging microscopy (FLIM), and untargeted liquid chromatography and mass spectrometry (LC-MS) analysis. The study reveals the efficacy of 3-BP in treating pancreatic cancer and identifies drug treatment-induced lipid species remodeling and associated pathways through bioinformatics analysis.

Charting the tumor microenvironment with spatial profiling technologies

In recent years technologies that can achieve readouts at cellular resolution such as single-cell RNA sequencing (scRNA-seq) have provided a comprehensive characterization of the cellular proportions and phenotypes that populate the tumor microenvironment (TME). However, because of the sample dissociation steps required by these protocols, they fail to capture information related to the intricate spatial context in which cells operate as well as their dense networks of interactions. Spatial profiling technologies have recently emerged as a valuable way to investigate the physical organization of cells crowding the TME in intact tissues. In this review we first discuss how spatial profiling technologies have propelled TME characterization, and then explore their potential to improve both diagnosis and prognosis for cancer patients in the clinic.

THItoGene: a deep learning method for predicting spatial transcriptomics from histological images

Spatial transcriptomics unveils the complex dynamics of cell regulation and transcriptomes, but it is typically cost-prohibitive. Predicting spatial gene expression from histological images via artificial intelligence offers a more affordable option, yet existing methods fall short in extracting deep-level information from pathological images. In this paper, we present THItoGene, a hybrid neural network that utilizes dynamic convolutional and capsule networks to adaptively sense potential molecular signals in histological images for exploring the relationship between high-resolution pathology image phenotypes and regulation of gene expression. A comprehensive benchmark evaluation using datasets from human breast cancer and cutaneous squamous cell carcinoma has demonstrated the superior performance of THItoGene in spatial gene expression prediction. Moreover, THItoGene has demonstrated its capacity to decipher both the spatial context and enrichment signals within specific tissue regions. THItoGene can be freely accessed at https://github.com/yrjia1015/THItoGene.

VistoSeg: Processing utilities for high-resolution images for spatially resolved transcriptomics data

Spatially resolved transcriptomics (SRT) is a growing field that links gene expression to anatomical context. SRT approaches that use next-generation sequencing (NGS) combine RNA sequencing with histological or fluorescent imaging to generate spatial maps of gene expression in intact tissue sections. These technologies directly couple gene expression measurements with high-resolution histological or immunofluorescent images that contain rich morphological information about the tissue under study. While broad access to NGS-based spatial transcriptomic technology is now commercially available through the Visium platform from the vendor 10× Genomics, computational tools for extracting image-derived metrics for integration with gene expression data remain limited. We developed VistoSeg as a MATLAB pipeline to process, analyze and interactively visualize the high-resolution images generated in the Visium platform. VistoSeg outputs can be easily integrated with accompanying transcriptomic data to facilitate downstream analyses in common programing languages including R and Python. VistoSeg provides user-friendly tools for integrating image-derived metrics from histological and immunofluorescent images with spatially resolved gene expression data. Integration of this data enhances the ability to understand the transcriptional landscape within tissue architecture. VistoSeg is freely available at http://research.libd.org/VistoSeg/.

Spatial Transcriptomics: Emerging Technologies in Tissue Gene Expression Profiling

In this Perspective, we discuss the current status and advances in spatial transcriptomics technologies, which allow high-resolution mapping of gene expression in intact cell and tissue samples. Spatial transcriptomics enables the creation of high-resolution maps of gene expression patterns within their native spatial context, adding an extra layer of information to the bulk sequencing data. Spatial transcriptomics has expanded significantly in recent years and is making a notable impact on a range of fields, including tissue architecture, developmental biology, cancer, and neurodegenerative and infectious diseases. The latest advancements in spatial transcriptomics have resulted in the development of highly multiplexed methods, transcriptomic-wide analysis, and single-cell resolution utilizing diverse technological approaches. In this Perspective, we provide a detailed analysis of the molecular foundations behind the main spatial transcriptomics technologies, including methods based on microdissection, in situ sequencing, single-molecule FISH, spatial capturing, selection of regions of interest, and single-cell or nuclei dissociation. We contextualize the detection and capturing efficiency, strengths, limitations, tissue compatibility, and applications of these techniques as well as provide information on data analysis. In addition, this Perspective discusses future directions and potential applications of spatial transcriptomics, highlighting the importance of the continued development to promote widespread adoption of these techniques within the research community.

NeST: nested hierarchical structure identification in spatial transcriptomic data

Spatial gene expression in tissue is characterized by regions in which particular genes are enriched or depleted. Frequently, these regions contain nested inside them subregions with distinct expression patterns. Segmentation methods in spatial transcriptomic (ST) data extract disjoint regions maximizing similarity over the greatest number of genes, typically on a particular spatial scale, thus lacking the ability to find region-within-region structure. We present NeST, which extracts spatial structure through coexpression hotspots-regions exhibiting localized spatial coexpression of some set of genes. Coexpression hotspots identify structure on any spatial scale, over any possible subset of genes, and are highly explainable. NeST also performs spatial analysis of cell-cell interactions via ligand-receptor, identifying active areas de novo without restriction of cell type or other groupings, in both two and three dimensions. Through application on ST datasets of varying type and resolution, we demonstrate the ability of NeST to reveal a new level of biological structure.

Harnessing computational spatial omics to explore the spatial biology intricacies

Spatially resolved transcriptomics (SRT) has unlocked new dimensions in our understanding of intricate tissue architectures. However, this rapidly expanding field produces a wealth of diverse and voluminous data, necessitating the evolution of sophisticated computational strategies to unravel inherent patterns. Two distinct methodologies, gene spatial pattern recognition (GSPR) and tissue spatial pattern recognition (TSPR), have emerged as vital tools in this process. GSPR methodologies are designed to identify and classify genes exhibiting noteworthy spatial patterns, while TSPR strategies aim to understand intercellular interactions and recognize tissue domains with molecular and spatial coherence. In this review, we provide a comprehensive exploration of SRT, highlighting crucial data modalities and resources that are instrumental for the development of methods and biological insights. We address the complexities and challenges posed by the use of heterogeneous data in developing GSPR and TSPR methodologies and propose an optimal workflow for both. We delve into the latest advancements in GSPR and TSPR, examining their interrelationships. Lastly, we peer into the future, envisaging the potential directions and perspectives in this dynamic field.

Spatially Resolved Transcriptomics Technology Facilitates Cancer Research

Single cell RNA sequencing (scRNA-seq) provides a great convenience for studying tumor occurrence and development for its ability to study gene expression at the individual cell level. However, patient-derived tumor tissues are composed of multiple types of cells including tumor cells and adjacent non-malignant cells such as stromal cells and immune cells. The spatial locations of various cells in situ tissues plays a pivotal role in the occurrence and development of tumors, which cannot be elucidated by scRNA-seq alone. Spatially resolved transcriptomics (SRT) technology emerges timely to explore the unrecognized relationship between the spatial background of a particular cell and its functions, and is increasingly used in cancer research. This review provides a systematic overview of the SRT technologies that are developed, in particular the more widely used cutting-edge SRT technologies based on next-generation sequencing (NGS). In addition, the main achievements by SRT technologies in precisely unveiling the underappreciated spatial locations on gene expression and cell function with unprecedented high-resolution in cancer research are emphasized, with the aim of developing more effective clinical therapeutics oriented to a deeper understanding of the interaction between tumor cells and surrounding non-malignant cells.

Machine learning reveals genetic modifiers of the immune microenvironment of cancer

Heritability in the immune tumor microenvironment (iTME) has been widely observed yet remains largely uncharacterized. Here, we developed a machine learning approach to map iTME modifiers within loci from genome-wide association studies (GWASs) for breast cancer (BrCa) incidence. A random forest model was trained on a positive set of immune-oncology (I-O) targets, and then used to assign I-O target probability scores to 1,362 candidate genes in linkage disequilibrium with 155 BrCa GWAS loci. Cluster analysis of the most probable candidates revealed two subfamilies of genes related to effector functions and adaptive immune responses, suggesting that iTME modifiers impact multiple aspects of anticancer immunity. Two of the top ranking BrCa candidates, LSP1 and TLR1, were orthogonally validated as iTME modifiers using BrCa patient biopsies and comparative mapping studies, respectively. Collectively, these data demonstrate a robust and flexible framework for functionally fine-mapping GWAS risk loci to identify translatable therapeutic targets.

Novel insights from spatial transcriptome analysis in solid tumors

Since its first application in 2016, spatial transcriptomics has become a rapidly evolving technology in recent years. Spatial transcriptomics enables transcriptomic data to be acquired from intact tissue sections and provides spatial distribution information and remedies the disadvantage of single-cell RNA sequencing (scRNA-seq), whose data lack spatially resolved information. Presently, spatial transcriptomics has been widely applied to various tissue types, especially for the study of tumor heterogeneity. In this review, we provide a summary of the research progress in utilizing spatial transcriptomics to investigate tumor heterogeneity and the microenvironment with a focus on solid tumors. We summarize the research breakthroughs in various fields and perspectives due to the application of spatial transcriptomics, including cell clustering and interaction, cellular metabolism, gene expression, immune cell programs and combination with other techniques. As a combination of multiple transcriptomics, single-cell multiomics shows its superiority and validity in single-cell analysis. We also discuss the application prospect of single-cell multiomics, and we believe that with the progress of data integration from various transcriptomics, a multilayered subcellular landscape will be revealed.

Deep learning exploration of single-cell and spatially resolved cancer transcriptomics to unravel tumour heterogeneity

Tumour heterogeneity is one of the critical confounding aspects in decoding tumour growth. Malignant cells display variations in their gene transcription profiles and mutation spectra even when originating from a single progenitor cell. Single-cell and spatial transcriptomics sequencing have recently emerged as key technologies for unravelling tumour heterogeneity. Single-cell sequencing promotes individual cell-type identification through transcriptome-wide gene expression measurements of each cell. Spatial transcriptomics facilitates identification of cell-cell interactions and the structural organization of heterogeneous cells within a tumour tissue through associating spatial RNA abundance of cells at distinct spots in the tissue section. However, extracting features and analyzing single-cell and spatial transcriptomics data poses challenges. Single-cell transcriptome data is extremely noisy and its sparse nature and dropouts can lead to misinterpretation of gene expression and the misclassification of cell types. Deep learning predictive power can overcome data challenges, provide high-resolution analysis and enhance precision oncology applications that involve early cancer prognosis, diagnosis, patient survival estimation and anti-cancer therapy planning. In this paper, we provide a background to and review of the recent progress of deep learning frameworks to investigate tumour heterogeneity using both single-cell and spatial transcriptomics data types.

A hybrid machine learning and regression method for cell type deconvolution of spatial barcoding-based transcriptomic data

Spatial barcoding-based transcriptomic (ST) data require cell type deconvolution for cellular-level downstream analysis. Here we present SDePER, a hybrid machine learning and regression method, to deconvolve ST data using reference single-cell RNA sequencing (scRNA-seq) data. SDePER uses a machine learning approach to remove the systematic difference between ST and scRNA-seq data (platform effects) explicitly and efficiently to ensure the linear relationship between ST data and cell type-specific expression profile. It also considers sparsity of cell types per capture spot and across-spots spatial correlation in cell type compositions. Based on the estimated cell type proportions, SDePER imputes cell type compositions and gene expression at unmeasured locations in a tissue map with enhanced resolution. Applications to coarse-grained simulated data and four real datasets showed that SDePER achieved more accurate and robust results than existing methods, suggesting the importance of considering platform effects, sparsity and spatial correlation in cell type deconvolution.

Omics Technologies Improving Breast Cancer Research and Diagnostics

Breast cancer (BC) has yielded approximately 2.26 million new cases and has caused nearly 685,000 deaths worldwide in the last two years, making it the most common diagnosed cancer type in the world. BC is an intricate ecosystem formed by both the tumor microenvironment and malignant cells, and its heterogeneity impacts the response to treatment. Biomedical research has entered the era of massive omics data thanks to the high-throughput sequencing revolution, quick progress and widespread adoption. These technologies-liquid biopsy, transcriptomics, epigenomics, proteomics, metabolomics, pharmaco-omics and artificial intelligence imaging-could help researchers and clinicians to better understand the formation and evolution of BC. This review focuses on the findings of recent multi-omics-based research that has been applied to BC research, with an introduction to every omics technique and their applications for the different BC phenotypes, biomarkers, target therapies, diagnosis, treatment and prognosis, to provide a comprehensive overview of the possibilities of BC research.

Applications of multi-omics analysis in human diseases

Multi-omics usually refers to the crossover application of multiple high-throughput screening technologies represented by genomics, transcriptomics, single-cell transcriptomics, proteomics and metabolomics, spatial transcriptomics, and so on, which play a great role in promoting the study of human diseases. Most of the current reviews focus on describing the development of multi-omics technologies, data integration, and application to a particular disease; however, few of them provide a comprehensive and systematic introduction of multi-omics. This review outlines the existing technical categories of multi-omics, cautions for experimental design, focuses on the integrated analysis methods of multi-omics, especially the approach of machine learning and deep learning in multi-omics data integration and the corresponding tools, and the application of multi-omics in medical researches (e.g., cancer, neurodegenerative diseases, aging, and drug target discovery) as well as the corresponding open-source analysis tools and databases, and finally, discusses the challenges and future directions of multi-omics integration and application in precision medicine. With the development of high-throughput technologies and data integration algorithms, as important directions of multi-omics for future disease research, single-cell multi-omics and spatial multi-omics also provided a detailed introduction. This review will provide important guidance for researchers, especially who are just entering into multi-omics medical research.

Implications of tumour heterogeneity on cancer evolution and therapy resistance: lessons from breast cancer

Tumour heterogeneity is pervasive amongst many cancers and leads to disease progression, and therapy resistance. In this review, using breast cancer as an exemplar, we focus on the recent advances in understanding the interplay between tumour cells and their microenvironment using single cell sequencing and digital spatial profiling technologies. Further, we discuss the utility of lineage tracing methodologies in pre-clinical models of breast cancer, and how these are being used to unravel new therapeutic vulnerabilities and reveal biomarkers of breast cancer progression. © 2023 The Pathological Society of Great Britain and Ireland.

Primary Aldosteronism: Spatial Multiomics Mapping of Genotype-Dependent Heterogeneity and Tumor Expansion of Aldosterone-Producing Adenomas

BACKGROUND

Primary aldosteronism is frequently caused by an adrenocortical aldosterone-producing adenoma (APA) carrying a somatic mutation that drives aldosterone overproduction. APAs with a mutation in KCNJ5 (APA-KCNJ5MUT) are characterized by heterogeneous CYP11B2 (aldosterone synthase) expression, a particular cellular composition and larger tumor diameter than those with wild-type KCNJ5 (APA-KCNJ5WT). We exploited these differences to decipher the roles of transcriptome and metabolome reprogramming in tumor pathogenesis.

METHODS

Consecutive adrenal cryosections (7 APAs and 7 paired adjacent adrenal cortex) were analyzed by spatial transcriptomics (10x Genomics platform) and metabolomics (in situ matrix-assisted laser desorption/ionization mass spectrometry imaging) co-integrated with CYP11B2 immunohistochemistry.

RESULTS

We identified intratumoral transcriptional heterogeneity that delineated functionally distinct biological pathways. Common transcriptomic signatures were established across all APA specimens which encompassed 2 distinct transcriptional profiles in CYP11B2-immunopositive regions (CYP11B2-type 1 or 2). The CYP11B2-type 1 signature was characterized by zona glomerulosa gene markers and was detected in both APA-KCNJ5MUT and APA-KCNJ5WT. The CYP11B2-type 2 signature displayed markers of the zona fasciculata or reticularis and predominated in APA-KCNJ5MUT. Metabolites that promote oxidative stress and cell death accumulated in APA-KCNJ5WT. In contrast, antioxidant metabolites were abundant in APA-KCNJ5MUT. Finally, APA-like cell subpopulations-negative for CYP11B2 gene expression-were identified in adrenocortical tissue adjacent to APAs suggesting the existence of tumor precursor states.

CONCLUSIONS

Our findings provide insight into intra- and intertumoral transcriptional heterogeneity and support a role for prooxidant versus antioxidant systems in APA pathogenesis highlighting genotype-dependent capacities for tumor expansion.

Pan-cancer T cell atlas links a cellular stress response state to immunotherapy resistance

Tumor-infiltrating T cells offer a promising avenue for cancer treatment, yet their states remain to be fully characterized. Here we present a single-cell atlas of T cells from 308,048 transcriptomes across 16 cancer types, uncovering previously undescribed T cell states and heterogeneous subpopulations of follicular helper, regulatory and proliferative T cells. We identified a unique stress response state, TSTR, characterized by heat shock gene expression. TSTR cells are detectable in situ in the tumor microenvironment across various cancer types, mostly within lymphocyte aggregates or potential tertiary lymphoid structures in tumor beds or surrounding tumor edges. T cell states/compositions correlated with genomic, pathological and clinical features in 375 patients from 23 cohorts, including 171 patients who received immune checkpoint blockade therapy. We also found significantly upregulated heat shock gene expression in intratumoral CD4/CD8+ cells following immune checkpoint blockade treatment, particularly in nonresponsive tumors, suggesting a potential role of TSTR cells in immunotherapy resistance. Our well-annotated T cell reference maps, web portal and automatic alignment/annotation tool could provide valuable resources for T cell therapy optimization and biomarker discovery.

Advances in spatial transcriptomics and related data analysis strategies

Spatial transcriptomics technologies developed in recent years can provide various information including tissue heterogeneity, which is fundamental in biological and medical research, and have been making significant breakthroughs. Single-cell RNA sequencing (scRNA-seq) cannot provide spatial information, while spatial transcriptomics technologies allow gene expression information to be obtained from intact tissue sections in the original physiological context at a spatial resolution. Various biological insights can be generated into tissue architecture and further the elucidation of the interaction between cells and the microenvironment. Thus, we can gain a general understanding of histogenesis processes and disease pathogenesis, etc. Furthermore, in silico methods involving the widely distributed R and Python packages for data analysis play essential roles in deriving indispensable bioinformation and eliminating technological limitations. In this review, we summarize available technologies of spatial transcriptomics, probe into several applications, discuss the computational strategies and raise future perspectives, highlighting the developmental potential.

Spatially resolved multi-omics highlights cell-specific metabolic remodeling and interactions in gastric cancer

Mapping tumor metabolic remodeling and their spatial crosstalk with surrounding non-tumor cells can fundamentally improve our understanding of tumor biology, facilitates the designing of advanced therapeutic strategies. Here, we present an integration of mass spectrometry imaging-based spatial metabolomics and lipidomics with microarray-based spatial transcriptomics to hierarchically visualize the intratumor metabolic heterogeneity and cell metabolic interactions in same gastric cancer sample. Tumor-associated metabolic reprogramming is imaged at metabolic-transcriptional levels, and maker metabolites, lipids, genes are connected in metabolic pathways and colocalized in the heterogeneous cancer tissues. Integrated data from spatial multi-omics approaches coherently identify cell types and distributions within the complex tumor microenvironment, and an immune cell-dominated "tumor-normal interface" region where tumor cells contact adjacent tissues are characterized with distinct transcriptional signatures and significant immunometabolic alterations. Our approach for mapping tissue molecular architecture provides highly integrated picture of intratumor heterogeneity, and transform the understanding of cancer metabolism at systemic level.

Deciphering tumor ecosystems at super resolution from spatial transcriptomics with TESLA

Cell populations in the tumor microenvironment (TME), including their abundance, composition, and spatial location, are critical determinants of patient response to therapy. Recent advances in spatial transcriptomics (ST) have enabled the comprehensive characterization of gene expression in the TME. However, popular ST platforms, such as Visium, only measure expression in low-resolution spots and have large tissue areas that are not covered by any spots, which limits their usefulness in studying the detailed structure of TME. Here, we present TESLA, a machine learning framework for tissue annotation with pixel-level resolution in ST. TESLA integrates histological information with gene expression to annotate heterogeneous immune and tumor cells directly on the histology image. TESLA further detects unique TME features such as tertiary lymphoid structures, which represents a promising avenue for understanding the spatial architecture of the TME. Although we mainly illustrated the applications in cancer, TESLA can also be applied to other diseases.

Inferring ligand-receptor cellular networks from bulk and spatial transcriptomic datasets with BulkSignalR

The study of cellular networks mediated by ligand-receptor interactions has attracted much attention recently owing to single-cell omics. However, rich collections of bulk data accompanied with clinical information exists and continue to be generated with no equivalent in single-cell so far. In parallel, spatial transcriptomic (ST) analyses represent a revolutionary tool in biology. A large number of ST projects rely on multicellular resolution, for instance the Visium™ platform, where several cells are analyzed at each location, thus producing localized bulk data. Here, we describe BulkSignalR, a R package to infer ligand-receptor networks from bulk data. BulkSignalR integrates ligand-receptor interactions with downstream pathways to estimate statistical significance. A range of visualization methods complement the statistics, including functions dedicated to spatial data. We demonstrate BulkSignalR relevance using different datasets, including new Visium liver metastasis ST data, with experimental validation of protein colocalization. A comparison with other ST packages shows the significantly higher quality of BulkSignalR inferences. BulkSignalR can be applied to any species thanks to its built-in generic ortholog mapping functionality.

Reconstruction of the cell pseudo-space from single-cell RNA sequencing data with scSpace

Tissues are highly complicated with spatial heterogeneity in gene expression. However, the cutting-edge single-cell RNA-seq technology eliminates the spatial information of individual cells, which contributes to the characterization of cell identities. Herein, we propose single-cell spatial position associated co-embeddings (scSpace), an integrative method to identify spatially variable cell subpopulations by reconstructing cells onto a pseudo-space with spatial transcriptome references (Visium, STARmap, Slide-seq, etc.). We benchmark scSpace with both simulated and biological datasets, and demonstrate that scSpace can accurately and robustly identify spatially variated cell subpopulations. When employed to reconstruct the spatial architectures of complex tissue such as the brain cortex, the small intestinal villus, the liver lobule, the kidney, the embryonic heart, and others, scSpace shows promising performance on revealing the pairwise cellular spatial association within single-cell data. The application of scSpace in melanoma and COVID-19 exhibits a broad prospect in the discovery of spatial therapeutic markers.