- 1West China School of Pharmacy, Sichuan University, Chengdu, China
- 2State Key Laboratory of Southwestern Chinese Medicine Resources, Chengdu University of Traditional Chinese Medicine, Chengdu, China
- 3Cardiovascular Department, Hospital of Chengdu University of Traditional Chinese Medicine, Chengdu, China
- 4Department of Clinical Medicine, Zhengzhou University, Zhengzhou, China
- 5School of Chinese Medicine, The University of Hong Kong, Hong Kong, Hong Kong SAR, China
- 6Key Laboratory of Drug-Targeting and Drug Delivery System of the Education Ministry, Sichuan Engineering Laboratory for Plant-Sourced Drug and Sichuan Research Center for Drug Precision Industrial Technology, Sichuan University, Chengdu, China
A comprehensive understanding of tumor heterogeneity, tumor microenvironment and the mechanisms of drug resistance is fundamental to advancing breast cancer research. While single-cell RNA sequencing has resolved the issue of “temporal dynamic expression” of genes at the single-cell level, the lack of spatial information still prevents us from gaining a comprehensive understanding of breast cancer. The introduction and application of spatial transcriptomics addresses this limitation. As the annual technical method of 2020, spatial transcriptomics preserves the spatial location of tissues and resolves RNA-seq data to help localize and differentiate the active expression of functional genes within a specific tissue region, enabling the study of spatial location attributes of gene locations and cellular tissue environments. In the context of breast cancer, spatial transcriptomics can assist in the identification of novel breast cancer subtypes and spatially discriminative features that show promise for individualized precise treatment. This article summarized the key technical approaches, recent advances in spatial transcriptomics and its applications in breast cancer, and discusses the limitations of current spatial transcriptomics methods and the prospects for future development, with a view to advancing the application of this technology in clinical practice.
1 Introduction
High-throughput sequencing (HTS), also known as next-generation sequencing (NGS), is a class of technologies that can sequence a large number of DNA or RNA samples simultaneously. Since its emergence, it has gradually become the most important tool for human genetic research due to its high sequencing efficiency, short time, low cost and ability to provide more accurate and richer genome sequence information (1). In the early development stage of HTS, it was mainly focused on DNA sequencing, especially in the field of genomics (2). With the continuous development of technology, RNA sequencing (RNA-seq) has gradually become a popular research direction, which promotes the research of transcriptomics and epigenetics by determining the sequence of RNA and analyzing the information of gene expression, transcript variants and variable splicing (3). With the rapid development of bioinformatics, materials science and computer science, RNA-seq technology has been applied at the single-cell level to analyze the transcriptome expression of a single cell, which is called single-cell RNA sequencing (scRNA-seq) (4). scRNA-seq technology can enable researchers to elucidate the ‘time-dynamic expression’ of genes at the single-cell level through the selection of samples collected at different time points, thus facilitating in-depth analysis of cellular heterogeneity and identification of different cell types (5, 6). However, the spatial information of the tissue samples is inevitably lost in the process. To address this issue, researchers have begun to focus on the spatial analysis of RNA-seq.
In 2016, the concept of spatial transcriptomics (ST) was first proposed (7). ST represents an emerging frontier technology that follows the development of scRNA-seq. It can preserve the spatial location of tissues and simultaneously analyze RNA-seq data of tissue sections to locate and distinguish the active expression of functional genes in specific tissue regions (8). As a current research hotspot of single-cellomics, ST has been rated as the annual technical method of 2020, allowing the location of genes to be assessed and the tissue environment and spatial properties of cells to be tracked (9). Various sequencing technologies including RNA-seq, scRNA-seq and ST are important tools in today’s biology and genomics research, bringing more choices to researchers. Each has its unique advantages and limitations for different research needs, as shown in Table 1.
Breast cancer is one of the most common cancers in the world, and its tumor heterogeneity and biological function require further analysis. The spatial location of cells in tissues is closely related to cell function. The tumor microenvironment (TME) of breast cancer is a highly structured ecosystem, including a variety of immune cells, cancer-associated fibroblasts (CAFs), endothelial cells and extracellular matrix (ECM) and other components (10). In the occurrence and development of breast cancer, in addition to changes in cell gene expression and function, it is often accompanied by disturbances in tissue structure and spatial location of cells (11). ST combines complex phenotypes and spatial information to quantify tumor spatial heterogeneity at the single-cell level, thus achieving new breast cancer subtype identification and spatially resolved feature analysis, which is of great significance for understanding the occurrence and development of breast cancer and providing new strategies for the treatment of breast cancer (12).
In this article, we reviewed the current methods and recent advances in ST. We also discussed the application of ST in breast cancer research, focusing on how ST can be used to study tumor heterogeneity and TME, and provide a reference for solving breast cancer drug resistance.
2 Technical approaches to ST
The concept of ST was originally derived fromfluorescently-labeled RNA sequences (FISH), which uses fluorescent probes in combination with targets to achieve in-situ detection of genes for qualitative and quantitative analysis and relative position analysis. In a broad sense, technologies that can simultaneously obtain tissue or cell transcriptomics information and spatial information are classified as ST, including the initial in situ imaging technology and the subsequent development of laser microdissection technology and spatial indexing technology (13, 14). After the combination of microarray technology and next-generation sequencing (NGS) technology, ST goes to the detection of a large number of spatial spot information, which greatly improves the detection efficiency. The current ST sequencing generally refers to the high-throughput ST technology based on NGS technology and microarray technology.
2.1 Laser capture microdissection-based approaches
To obtain spatial gene expression information, the target region can first be isolated from the sample, and then the dissected single region is placed separately in a test tube for RNA extraction and sequencing. The separation of target regions can be achieved by frozen section, RNA tomography (tomo-seq), laser capture microdissection (LCM) and other methods (15–17). Among them, LCM is more favored by researchers due to its advantages of high speed, high precision and versatility.
LCM is a technique for cutting, separating and collecting selected samples under a microscope using a micromanipulation system. It works by combining a laser with an inverted microscope and connecting it to a computer for direct visual analysis (18, 19). The combination of LCM and scRNA-seq allows RNA sequencing of individual cells or tissue regions (20). On this basis, Chen et al. constructed a geographical position sequencing (GEO-seq) method by integrating and optimizing LCM and scRNA-seq, which can obtain a small amount of cell transcriptomics information while retaining the original location information of cells, so as to achieve high-throughput gene expression (21). Geo-seq enables high-resolution ST by using a technique that tags mRNA in situ with spatially barcoded probes to map gene expression across tissue sections.
Additionally, transcriptomics in vivo analysis (TIVA) enables the examination of ST profiles in living cells. However, its analytical throughput is limited, with the capacity to analyze only a few single cells at a time (22). NICHE-seq employs photoactivation to conduct transcriptomics analysis on individual cells within a specific niche, but its applicability is limited to model organisms due to the necessity of genetic engineering (23). ProximID is another approach to examine the niche within the cellular microenvironment. It innovatively identifies protein-protein interactions within cells by using proximity labeling and high-throughput sequencing, allowing for spatially resolved mapping of protein interactions. Nevertheless, it remains constrained by the limitation of low flux (24).
In conclusion, as shown in Figure 1A, this technology can be used to obtain ST information at the cellular level; however, its analytical throughput is typically low, and the operation is complex and time-consuming.
Figure 1. Advantages and disadvantages of different spatial transcriptomics technologies. (The green color shows the advantages of the technology, while the yellow color shows the disadvantages.) (A) Laser capture microdissection-based approaches. (B) In situ hybridization-based approaches. (C) In situ sequencing-based approaches. (D) Next-generation sequencing-based approaches.
2.2 In situ imaging-based approaches
Imaging-based ST technology employs microscopy to observe mRNA following the combination of a chromogenic group with a probe for in situ imaging. In situ hybridization (ISH) and in situ sequencing (ISS) are the two most commonly utilized methods.
ISH is continuously hybridized in the tissue by a marker probe that is complementary to the target sequence, and then imaged and quantitatively positioned under a microscope (25). Single-molecule RNA fluorescence in situ hybridization (smFISH) employs multiple shorter probes to target discrete regions of the transcriptomics, thereby enabling quantitative measurement of the transcriptomics (26). Subsequently, the advent of RNAscope technology enabled the simultaneous detection of up to 12 RNA targets, facilitating the identification of a smaller number of mRNA and protein targets (27). Continuous fluorescence in situ hybridization increases the signal intensity of FISH by performing multiple rounds of hybridization of a series of probes with transcripts to achieve high-throughput analysis (28–30). Furthermore, a novel technology designated as DNA Microscopy was put forth in 2019 (31). This method is based on thermodynamic entropy and does not necessitate the use of sophisticated, costly optical instrumentation. The sample itself, in conjunction with the requisite reagents, is capable of providing the requisite spatial information as part of the chemical reaction, thus assisting in the clarification of the spatial structure information present in cells and tissues (32).
In summary, ISH offers the benefits of high resolution and extensive flux detection; however, there are still constraints, including prolonged processing times, significant errors, and elevated costs (Figure 1B).
In 2013, the inaugural publication of ISS technology was released (33). In complete tissue sections, mRNA is initially transcribed into cDNA, and subsequently, a single-stranded DNA molecule, referred to as a lock-probe in situ sequencing (Lock-Seq) method, is employed to target identified genes. In comparison to ISH, ISS exhibits superior subcellular resolution; however, the number of targets that it can examine is restricted. Subsequently, FISSEQ and ExSeq enhanced the capabilities of ISS by employing two query bases to sequence circular and RCA-amplified cDNAs, as opposed to utilizing a single query base per probe to sequence gene barcodes (34, 35). The advantage of STARmap is that a second primer is added to the site next to Lock-Seq, thereby avoiding the reverse transcription step and improving efficiency and reducing noise (36). Furthermore, the advent of technologies such as BaristaSeq, INSTA-Seq, and HybISS have revolutionized ISS by improving spatial resolution, throughput, and sensitivity. BaristaSeq enables high-throughput sequencing directly from tissue sections with enhanced precision, INSTA-Seq provides rapid, single-molecule sequencing for faster analysis, and HybISS combines hybridization and ISS for precise ST (37–39). These advancements allow for more detailed, scalable, and efficient mapping of gene expression across tissues, opening new possibilities for studying complex tissue architectures and their molecular dynamics in health and disease (Figure 1C).
2.3 NGS-based approaches
The ST techniques described above are based on the separation of known tissue regions of interest, or in situ visualization of RNA molecules by hybridization or sequencing, while NGS technology is based on in situ capture of transcripts, followed by ectopic sequencing. This method was developed from the conceptual innovation of the scRNA method.
In 2016, Ståhl PL et al. successfully obtained spatial-resolved full-transcriptomics information using ST technology. This technique captures polyadenylated RNA on the slide of the spatial barcode microarray before reverse transcription, thereby ensuring that each transcript can be mapped back to its original point through a unique positional molecular barcode (40). Subsequently, this innovative approach was acquired by 10× Genomics and commercialized as a Visium space gene expression solution (41).
In 2019, high-throughput ST sequencing emerged, with the most representative method being high definition spatial transcriptomics (HDST) (42–44). HDST involves embedding silica gel magnetic beads with a diameter of 2 μm on a slide, with a specific spatial barcode connected to the surface of the magnetic beads to capture the mRNA of the cells at the corresponding position for reverse transcription and transcriptomics sequencing.
As shown in Figure 1D, the limitation of ST is that the analysis location is a small area of multiple cells, which is not a genuine analysis of a single cell. Consequently, subsequent new technology has focused on maintaining the resolution of a single cell or subcellular level. Seq-Scope, based on the Illumina sequencing platform, allows for the visualization of transcriptomic heterogeneity at the cellular and subcellular levels in a range of tissues, with a sub-micron resolution (45). Sci-Space combines single-cell gene expression differences with spatial backgrounds, preserving single-cell resolution (46). Furthermore, Stereo-Seq integrates DNA nanospheres (DNB) and in situ RNA capture technology to facilitate high-throughput transcriptomics analysis of tissue sections at unparalleled nanoscale resolutions. The visual field can be extended to the centimeter level, exhibiting high sensitivity and uniform capture rate (47). Stereo-Seq is anticipated to become a fundamental tool for analyzing transcriptomics heterogeneity in complex tissues and organisms.
2.4 Computer reconstruction of spatial data
The high cost of ST precludes its extensive use in clinical applications, particularly when evaluated against the backdrop of alternative experimental methods. Currently, the advent of information technology has given rise to a novel approach for calculating spatial analysis gene expression datasets. This method employs computer algorithms to simulate the spatial morphology of reconstructed tissues based on single-cell transcriptomics data. To illustrate, Seurat, an R package that is frequently employed in single-cell data analysis (48). By integrating spatial information and gene expression data, Seurat is able to reconstruct spatial transcriptional profiles of tissues, revealing the spatial distribution and functional status of cells in tissues. It can process in situ sequencing data of spatially tagged genes for spatial clustering, differential expression analysis and cell type annotation. In addition, Seurat supports multimodal data integration to help researchers understand the dynamics of gene expression at the spatial level and promote in-depth studies of disease mechanisms and tissue structure (49, 50).
Furthermore, the optimization of deep learning calculation methods can be conducted repeatedly in the context of continuous modeling and verification. Additionally, the prediction of ST expression can be achieved through the utilization of existing histological images (51, 52). The toolboxes Starfysh, DeepST, MISTy, standR, and others are based on deep generative models. Such methods can not only analyze the spatial dynamic characteristics of complex tissues, but also identify spatial hubs composed of different subtypes of cells. They have become invaluable for investigating the structural and functional interactions of disparate spatial backgrounds in ST.
Starfysh is designed for processing spatial transcriptome data from tissue sections. It is able to automatically identify and classify spatial regions and reveal the spatial distribution of different cell populations in tissues by combining image data and gene expression information. Using deep neural networks, Starfysh accurately combines spatial information from tissue sections with gene expression patterns to achieve a comprehensive analysis of spatial structure and function, which is widely used to study spatial heterogeneity of tissues and dynamic changes in the microenvironment (53). DeepST aims to predict the spatial location of cells in tissue sections by learning complex patterns in spatial transcriptomic data. The tool is able to automatically learn spatial associations between cells from the data by combining gene expression data with tissue section image information and using a convolutional neural network (CNN) model for spatial reconstruction and prediction. DeepST is mainly used to study complex tissue structures such as cell type distribution, spatial dynamics features, and the TME, and helps to reveal the interactions between cell function and tissue morphology (54). MISTy accurately predicts the spatial location of cells and their transcriptional status in different spatial regions by analyzing high-resolution images of tissue sections in combination with spatial transcriptome data. It can reveal the multilevel structure and dynamic changes of spatial gene expression, and is particularly suitable for spatial gene expression analysis of the tumor immune microenvironment and complex tissues (55). StandR is mainly used to solve the problem of correlation between spatial structure and gene expression of complex tissues. It is capable of identifying cell subtypes and their distribution patterns in different spatial locations through deep learning modeling of spatial transcriptome data combined with tissue section image information. standR not only supports large-scale data analysis, but also effectively identifies spatial interactions among cells, and is widely used in the study of tissue differentiation, tumor heterogeneity and disease-related microenvironments (56).
The innovative calculation method can also address the issue of ST being challenging to differentiate between similar cell populations, thereby enhancing the efficiency of sample capture. Artificial intelligence (AI) can identify different cell types and tissue regions in HE staining images and assign cell type labels to each spatial location. Combined with spatial transcriptome data, this annotation information helps to resolve the spatial layout of different cell populations in tissues and their association with gene expression (57–59). Zubair A et al. employed cell type-specific immunofluorescence markers obtained from mouse brain tissue sections and utilized artificial intelligence to annotate the output of HE tissue images, which markedly enhanced the recognition rate of clinically pertinent immune cell infiltration in tissues and broadened the scope of ST applications (60).
In short, laser capture microdissection-based approaches and in situ imaging-based approaches usually require expensive equipment (e.g. laser microdissectors, high-resolution imaging microscopes) and are therefore more costly. Both approaches also require more time for data acquisition and are less efficient, especially for complex tissue sections. NGS-based approaches and computer reconstruction of spatial data, on the other hand, are relatively less costly as they rely on high-throughput equipment and computational power, but do not require expensive hardware support. They are more efficient when dealing with large numbers of samples and can quickly process and analyze large spatial data sets. In terms of accuracy of results, laser capture microdissection-based approaches and in situ imaging-based approaches provide high spatial resolution and accuracy, and are suitable for fine analysis at the cellular level. NGS-based approaches and computer reconstruction, although capable of handling large data sets, may not provide the same precise spatial localization as the former two and are therefore slightly less accurate in terms of result accuracy. Table 2 shows a comparison of the advantages and disadvantages of the different technological tools.
3 The workflow of ST
Despite the advent of numerous technical approaches to ST analysis, the majority of these methods are hindered by significant limitations, including intricate operational complexity, suboptimal detection accuracy, and constrained detection throughput. Consequently, they are not yet suitable for widespread implementation. At present, 10× Genomics Visium is the only commercially available ST sequencing scheme that offers high precision and high throughput. The technology has a straightforward workflow, a relatively short processing time, and a wide range of potential applications. The key steps in the Visium experiment are illustrated in Figure 2.
Figure 2. Experimentation workflow of the spatial transcriptomics. (Created with BioRender.com).
3.1 Sample and slice preparation
Initially, fresh tissue samples were fixed in isopentane and frozen in liquid nitrogen. This process was employed to safeguard the quality of RNA and the integrity of tissue morphology, thereby ensuring an accurate reflection of cellular and genetic distribution in space. Optimal cutting temperature compound (OCT) is a fixative agent. Following rapid freezing, the tissue samples must be embedded with OCT embedding agent. This has the dual benefit of reducing wrinkles and fragmentation during slicing and increasing the continuity of tissues. Prior to slicing with a freezing microtome, the tissue and 10 × slide were placed in a −20°C freezing microtome box for approximately 30 minutes. The thickness of the slice was typically 10 μm, with a maximum size limitation (61).
3.2 RNA capture and cDNA synthesis
After the preparation of the slices, HE staining or IF staining was performed and photographed to detect the quality of the slices and patches. Subsequently, tissue permeabilization was performed and fluorescent cDNA was synthesized. The optimal permeabilization time was then identified according to the fluorescence intensity. The slices were then processed according to the optimal permeabilization time. The addition of permeabilization reagents to the slices resulted in the mRNA of the cells being released vertically down to the slide and bound to the capture sequence of the spatial microarray. cDNA was then synthesized by reverse transcription of the captured mRNA.
3.3 Library construction, high-throughput sequencing and data analysis
Following the acquisition of cDNA, the library was constructed and sequenced using an illumina high-throughput sequencing instrument, thereby obtaining the requisite sequencing data. The fundamental aspect of ST data analysis is the implementation of dimensionality reduction and clustering techniques, which are based on the gene expression variations observed across different spatial units on each chip. This is then correlated with the corresponding tissue image. Concurrently, the spatial position of each gene expression within the tissue can be determined.
In general, Seurat is employed for the analysis, visualization and integration of spatial data sets. The field of spatial analysis encompasses a range of techniques, including normalization, dimensional reduction and clustering, spatially-variable features, interactive visualization, and working with multiple slices. The steps of ST data preprocessing are analogous to those of typical scRNA-seq. To address the discrepancy in sequencing depth between data points, the data must be normalized. In contrast to scRNA-seq, the molecular plot of a spatial dataset can exhibit considerable variation, particularly in instances where there are notable differences in cell density across the tissue, resulting in a substantial degree of heterogeneity. Sctransform facilitates the normalization of data, enabling the detection of high variation features and the storage of data in SCT format (62).
3.4 Integration of ST and scRNA-seq data
10× Genomics Visium technology is a powerful ST tool that enables efficient spatial localization of gene expression in tissue sections (63, 64). However, there are some potential limitations of this technology, especially when dealing with different types of breast cancer samples. Firstly, the limited spatial resolution of the Visium technique, typically 55 µm, means that it may not be able to capture fine-scale differences between cellular subpopulations in the TME, particularly interactions between tumor cells and immune or stromal cells. Secondly, Visium is dependent on the quality and integrity of the tissue sections, which may affect the accuracy of the data for some breast cancer samples, especially those with high tissue stiffness or heterogeneity. In addition, the technique may be affected by background noise and technical bias when processing complex tumor samples, requiring multiple optimizations and algorithmic adjustments to improve data availability and resolution. Therefore, the Visium technology still needs to be combined with other technologies (e.g., laser microdissection or scRNA-seq) to obtain more accurate spatial transcriptomic data in studies of different breast cancer subtypes (65, 66).
By anchoring and integrating ST and scRNA-seq data, a three-dimensional transcriptomics map of the target tissue can be obtained. In this process, since the noise patterns of the representation space and the single-cell data set are very different, it is generally recommended to use the integral method rather than the deconvolution method (67).
The emerging ST technology, combined with microscopic imaging and scRNA-seq, can reveal the gene expression and arrangement distribution of individual cells at subcellular spatial resolution. This allows for the maximization of cellular heterogeneity, the discovery of new cell populations and cell subsets, and the clarification of the regulatory mechanisms that govern cellular status. This, in turn, provides technical support for a range of applications (68, 69).
4 Application of ST in the treatment of breast cancer
Conventional scRNA-seq fails to provide spatial information about cells in tumor tissue, which is particularly critical for breast cancer treatment (70). Tumor heterogeneity and microenvironment have a significant impact on treatment response, and single-cell data cannot reveal gene expression differences and cellular interactions between different regions (71). The application of spatial transcriptome sequencing in breast cancer provides an important tool to study the TME and its interaction with tumor cells. By combining gene expression and spatial information, this technology can reveal the spatial distribution and functional status of different cell populations within a tumor, helping researchers to understand tumor heterogeneity, immune escape mechanisms, and the process of tumor progression. Spatial transcriptomic data in breast cancer samples can help identify changes in immune cells, stromal cells, and blood vessels around the tumor, providing targets for precision therapy. For example, it has been found that the spatial relationship between immune cells and tumor cells may influence the tumor’s ability to escape immunity, and thus the response to treatment (72, 73). In addition, the spatial transcriptome could be used to monitor molecular signatures in different subtypes of breast cancer, which in turn could facilitate the development of personalized treatment strategies.
4.1 Tumor heterogeneity
Tumor heterogeneity refers to the changes in molecular biology or genes during tumor growth, resulting in differences in growth rate, invasion ability, drug sensitivity and prognosis of different tumor cells (74, 75). Despite the prevailing view that breast cancer has a monoclonal origin, the occurrence and development of this disease is the result of numerous divisions and proliferations, continuous evolution, epigenetics, genomics and microenvironment changes, which collectively give rise to heterogeneity (76). Heterogeneity represents a significant challenge in the treatment of breast cancer, manifesting as diverse histological subtypes, varying treatment sensitivities, and disparate clinical outcomes among patients. ST offers a valuable avenue for elucidating the genetic and phenotypic distinctions of tumor cells, analyzing the spatial heterogeneity of gene expression across diverse tissue regions, and providing a multi-dimensional and comprehensive perspective for investigating tumor heterogeneity (77–79).
The majority of studies have integrated ST with pathological annotation and deconvolution, utilizing diverse methodologies to accurately delineate the location, enrichment characteristics and differentially expressed genes (DEGs) of breast cancer subtypes, thereby elucidating the spatial distribution and transcriptional heterogeneity of breast cancer (80, 81). In 2013, the heterogeneity of 31 gene expression in breast cancer tissues was first observed using ISS technology. On this basis, Svedlund et al. conducted a combined molecular and morphological analysis of the expression of 91 genes in breast cancer tissues, revealing the intratumoral heterogeneity associated with breast cancer subtypes (82). Single-cell transcriptome sequencing of 26 primary breast cancers and spatial transcriptome sequencing of six breast cancers were performed to construct single-cell and spatial maps of breast cancers, revealing cellular heterogeneity in recurrent tumors and identifying novel PD-L1/PD-L2+ macrophage populations associated with clinical outcomes. Through this integrative analysis, the researchers were able to classify breast cancers into nine ‘ecotypes’ with distinct cellular composition and clinical outcomes, suggesting that the use of ST is statistically significant in classifying breast cancers and predicting response to therapy (83). A paucity of zinc finger protein 689 (ZNF689) has been demonstrated to facilitate LINE-1 reverse transcription, thereby exacerbating genomic instability and resulting in elevated intratumoral heterogeneity (ITH) (84).
The combined analysis of ST and single-cell omics allows for the consideration of both resolution and gene-level flux, which is an essential method for elucidating tumor heterogeneity (85). Studies have found that ER+ breast cancer has four heterogeneous groups, namely estrogen-responsive group, proliferative group, hypoxia-inducible group and inflammation-related group. The proliferative compartment has been identified as a key factor in the invasiveness of luminal B breast cancer (86). Liu et al. not only determined the heterogeneity of characteristics, origin and function of different cell populations in breast cancer tissues, but also employed the spatial distribution map of tumor subsets to describe the different stromal cell types in different tissue regions, providing a comprehensive method for describing the heterogeneity and structure of breast cancer (87). In luminal epithelial cells, Kohler et al. observed 11 distinct breast cancer clusters and differentiation trajectories originating from keratin 15+ (K15+) progenitor cells, indicating a robust correlation between normal K15+ primordium and basal-like breast cancer (88).
Copy number variation (CNV) can be used to reveal precancerous lesions (89). Lu et al. employed a combination of FISH and LCM with Smart-3SEQ to ascertain that HER2-amplified ductal carcinoma in situ (DCIS) has multiple subclones with different CNVs (90).
In addition, the spatial heterogeneity of metabolic structure is also an important part of tumor heterogeneity. ST studies have shown that glycolysis, profibrotic states, and metabolic reprogramming based on lipid metabolism play an important role in the extensive heterogeneity of breast cancer (91–93).
The high cost of ST limits its clinical application. To address this issue, some more cost-effective new methods have emerged. The deep learning algorithms ST-Net and BrST-Net can link the cell morphological characteristics and gene expression in ST, predict gene expression using HE tissue case images, and identify the heterogeneity of gene expression in tumors (94, 95). Based on the optimization and verification of a single deep convolutional neural network, the expression-morphology (EMO) analysis of the transcriptomics coverage area was established to predict gene transcription differences and proliferation markers in different regions of H&E staining images. The new method is reliable, economical and scalable (96).
Table 3 provides an overview of studies on the use of ST to study heterogeneity in breast cancer.
4.2 Tumor microenvironment
Intratumoral heterogeneity emphasizes the difference of cells or cell subsets in TME, which can reshape the occurrence, development, invasion, metastasis and drug sensitivity of tumors (97). TME is constituted of a heterogeneous population of cells, including tumor cells and surrounding immune cells, interstitial cells, and active mediators (98). It can be divided into two distinct categories: an immune microenvironment, which is dominated by immune cells, and a non-immune microenvironment, which is dominated by fibroblasts. Currently, cancer research and treatment have changed from a cancer-centric model to a tumor microenvironment-centric model. ST technology is instrumental in elucidating the diversity and intricacies of cell populations within the breast cancer microenvironment, and has become a powerful tool for a comprehensive investigation of the occurrence, progression and escape mechanisms of breast cancer (99–101).
Through spatial analysis of samples from HER2+ breast cancer patients before and after neoadjuvant HER2-targeted therapy, it was found that treatment resulted in decreased activation of HER2-dependent signaling, increased immune infiltration, and that increased expression of CD45 at the time of treatment more highly accurately predicted the presence of a pathological complete remission (pCR) in patients, findings which reveal the role of ST in predicting treatment response and guide clinical management (102).
Tumor-infiltrating macrophages play a crucial role in tumor progression, and their accumulation in TME is associated with poor prognosis and drug resistance in a multitude of cancers. In TME, macrophages exhibit different phenotypic alterations, leading to heterogeneous immune activation (103). At present, the treatment of macrophages focuses on inhibiting their recruitment or reprogramming their phenotype from M2-like to M1-like. The development of ST has improved the understanding of the macrophage body, phenotype and functional plasticity, and provided a new way for macrophage targeted therapy (104). In order to study the determinants of functional heterogeneity of TME macrophages, researchers employed an in vitro organ-type TME (oTME) model, an in vivo mouse model, and human samples of breast cancer. The study found a subset of F4/80highSca-1+ self-renewing macrophages mediated by Notch4 and maintained by type I interferon (IFN), which is closely related to the growth and metastasis of breast cancer (105). In breast cancer cell-infiltrating lymph nodes, macrophages are abundant in pathways that promote tumor growth, NF-κB and NOD-like receptor signaling pathways, and elucidate the differentiation trajectory of macrophages from active chemokine production to active lymphocyte activation (106). Tzeng et al. found that LSM1, a part of the cytoplasmic protein complex Lsm1-7-Pat1, affects the energy metabolism of TME macrophages and breast cancer metastasis, suggesting its potential as a diagnostic and prognostic marker (107).
Dysregulation of lipid metabolism is also a feature of TME, especially excessive free fatty acids (FFA). Studies have found that lipid-associated macrophages (LAMs), a unique macrophage subset, can highly express macrophage markers, lipid metabolism genes and lipid receptors, have M2 macrophage phagocytosis characteristics, and can enhance lipid accumulation (108). In addition, macrophages exhibit strong tryptophan metabolic activity in breast cancer, which can promote the polarization of M1 macrophages and predict the efficacy of immunotherapy (109).
The spatial distribution of tumor infiltrating lymphocytes (TILs) can be a predictor of the prognosis of breast cancer and the response to systemic treatment. Romanens et al. employed ST-FFE to analyze the immune cell subsets in the stromal and intraepithelial regions of TNBC. Their findings indicated that, in comparison to stromal cells, intraepithelial T cells and B cells exhibited reduced diversity and stronger clonality (110).
Cancer-associated fibroblasts (CAFs) are another diverse cell population indispensable for remodeling TME. Studies have identified four distinct functional fibroblast subsets and characterized their spatial distribution patterns. Among these, iCAFs are enriched in breast cancer and play a pivotal role in promoting an immunosuppressive microenvironment, which is altered in accordance with the status of breast cancer (111, 112). In a study by Kang et al., the impact of homologous recombination defect (HRD) on TME development at varying scales was investigated. The findings revealed that HRD alters the dominant phenotype of CAFs from mCAFs to iCAFs and is associated with immune cell function (113).
Table 4 presents a summary of studies on the applications of ST to breast cancer microenvironment.
4.3 Drug resistance
Drug resistance is the main obstacle to the treatment of breast cancer, which is related to tumor heterogeneity. Using ST to detect drug resistance-related cell subsets and explore their unique gene expression patterns is essential for predicting and reversing drug resistance. Basal epithelial subsets are located in the matrix and have the characteristics of drug resistance, and the up-regulated ACTN1 can promote chemotherapy resistance (114). MCU, as a diagnostic biomarker for breast cancer, is associated with advanced clinical status and a low overall survival rate. Inactivation of this protein has been shown to increase sensitivity to specific small molecule drugs (115).
The objective of immunotherapy is to reinstate a normal anti-tumor immune response by reactivating and sustaining the tumor-immune cycle, thereby controlling and clarifying tumors. Immunotherapy includes immune checkpoint inhibitors, therapeutic antibodies, and small molecule inhibitors. In recent years, immunotherapy has also shown robust anti-tumor activity in breast cancer, bringing long-term clinical benefits to some patients (116). Unfortunately, due to the loss of immune response and the emergence of drug resistance, this treatment is only effective for some patients. ST can facilitate a deeper understanding of the immunobiology of TME and identify biomarkers for predicting the efficacy of immunotherapy in patients (117). Tashireva et al. observed that PD-L1-positive tumors lack direct contact between PD1 receptor-expressing cells and PD-L1 ligand-expressing cells, and the lack of this specific immune response may be the reason for the low sensitivity of TNBC patients to immune checkpoint inhibitors (ICIs) (118). Indeed, only a single subgroup of TNBC exhibits responsiveness to ICIs. Investigations have examined the correlation between the spatial organization of immune cells in TNBC and T cell evasion, identifying that deposits of collagen-10, enhanced glycolysis, and activation of TGF-β/VEGF pathways are associated with resistance to anti-PD1 (119). TNBC patients with high Treg infiltration-related scores lack the stimulation of immune activation pathways, which results in resistance to immunotherapy against PD-1 (120).
OTUD4/CD73 proteolytic axis therapy represents a promising avenue for the treatment of immunosuppressive TNBC. The ubiquitination of CD73 by OTUD4 can counteract the deubiquitination of CD73 by TRIM21, resulting in the stability of CD73 and the prevention of a tumor immune response. ST80, a novel inhibitor that destroys the proteolytic interaction between OTUD4 and CD73, has the potential to enhance the efficacy of PD-L1 treatment in TNBC (121).
IMM2902 is an anti-CD47/HER2 bispecific antibody that has been developed for the treatment of patients with trastuzumab-resistant breast cancer. Furthermore, IMM2902 can also stimulate macrophages to produce CXCL9 and CXCL10, thereby increasing the level of T cells and NK cells at the tumor site (122).
Antibody-drug conjugates (ADCs) are a type of targeted therapy for cancer that connect antibodies to drugs through precisely designed chemical ligation methods, showing excellent efficacy in the treatment of breast cancer. ADCs hinder the progression of cell cycle. Cetuximab-CDK inhibitors can use the overexpression of epidermal growth factor receptor and the disordered cell cycle in invasive and drug-resistant tumors to target and effectively deliver drugs targeting cell cycle to basal-like/TNBCs, bringing new hope for the treatment of drug-resistant breast cancer patients (123).
A summary of studies applying ST to breast cancer drug resistance is presented in Table 5.
5 New hope for individualized precise treatment
Breast cancer has a variety of subtypes, which are different in genetics, histology and clinical features, and there are potential variations between different subtypes (124, 125). The treatment of breast cancer is dependent on a number of factors, including the specific subtype of cancer, the extent to which it has spread from the breast to the lymph nodes (stage II or III) or other parts of the body (stage IV), the patient’s response to previous treatment, the molecular composition of tumor cells, and the state of the patients’ immune system (126). As research progresses, the treatment of breast cancer has completed the transformation from single therapy to comprehensive therapy, and its treatment methods are also emerging in an endless stream, including surgery, chemotherapy, targeted therapy, endocrine therapy, immunotherapy and so on (127). Nevertheless, the question of how to ascertain the most appropriate treatment for each patient, taking into account their individual circumstances, has remained a significant challenge in the research and development of anti-tumor drugs and clinical treatment (128).
The ST sequencing analysis of breast cancer patients can better determine the differences in tumor heterogeneity and TME in patients with different disease stages and molecular typing, and can enable a comprehensive investigation of the diverse cell types and their molecular mechanisms in vivo. On this basis, targeted markers and new therapeutic targets can be determined, thereby enabling the personalized and precise treatment of cancer patients (129). For instance, some researchers have proposed a model that can reproduce the clinical treatment response. By comparing patients with pCR and no response progression (pNR), they identified differences in transcription procedures in tumors, stroma and immune infiltration, providing new ideas for the improvement and personalization of TNBC treatment methods (130).
In addition, it is important to identify the most efficacious method of eliciting anti-tumor immune responses, thereby providing patients with a more robust and effective immunotherapy. Shiao SL et al. analyzed TNBC biopsy tissue by integrating scRNA-seq and ST, and found that patients exhibiting resistance to immunotherapy exhibited a lack of immune infiltration both before and after treatment. Conversely, patients who responded to immunotherapy demonstrated two distinct patterns: one with a pre-existing immune response and the other with an immune response only evident following αPD1 and radiotherapy (131). At the same time, research has revealed that quiescent cancer cells (QCCs) in primary TNBC establish a micro-network with damaged dendritic cells, which can resist the attack of T cells, reduce immune infiltration, and demonstrate heightened tumorigenic potential (132).
In conclusion, as shown in Figure 3, the integration of ST into basic and translational research can facilitate the development of new drugs by identifying potential therapeutic targets, reveal promising biomarkers to monitor therapeutic effects and guide treatment decisions, and predict unique sensitive drugs for different patients, which is of great significance for promoting individualized precise treatment of breast cancer.
Figure 3. Application of ST in the field of breast cancer research. (Created with BioRender.com).
The application of ST in BC treatment provides new perspectives for understanding the TME and individualized treatment, but whether the results are statistically significant or not relies on multifactorial analysis and clinical validation. Currently, ST is still in the exploratory stage, but some preliminary studies have demonstrated its potential in breast cancer and provided statistical significance to some extent. At present, the high cost of ST technology itself and the complexity of data processing also limit its widespread application in clinical settings. Therefore, future studies need to further optimize the technology, expand the sample size and combine it with other histological data (e.g. proteomics, metabolomics, etc.) for a comprehensive analysis to ensure the statistical significance and usefulness of its clinical application.
6 Conclusion and discussion
ST is a pioneering technology used to investigate the heterogeneity of tumors and their microenvironment. It can identify and differentiate the expression of functional genes in specific spatial locations, elucidate the mechanisms, potential drug targets, and novel biomarkers associated with the origin and progression of tumor cell heterogeneity, and investigate the interactions between diverse tumor cell types and immune cells, interstitial components (133–135). ST combines the advantages of traditional transcriptomics and spatial histology to capture the distributional features of genes in tissues without disrupting the tissue structure (136). Key advances in this technology include high-throughput spatially localized RNA sequence analysis, precise mapping of transcript distribution through high-resolution microscopy and imaging, and large-scale spatial transcriptome data acquisition through novel microarray platforms (137).
Initial ST technologies relied on combining tissue slices with microarray platforms by immobilizing RNA on tissue slices and analyzing them using high-throughput sequencing techniques. However, this approach had limited spatial resolution. As the technology has evolved, ST with single-cell resolution has gradually emerged, allowing precise localization of gene expression in individual cells by integrating high-resolution microscopy and genomics approaches (138, 139). In recent years, innovations based on new platforms such as nanopore technology, microfluidic chips and fluorescence imaging have led to more efficient and precise data acquisition of spatial transcriptomes and greatly improved spatial resolution (140, 141). In addition, advances in data processing technologies, such as the application of deep learning and machine learning algorithms, have helped to better analyze complex spatial transcriptome data, revealing cellular interactions and spatial heterogeneity. Nevertheless, despite the rapid advancement of this field, certain limitations persist.
Firstly, the current ST technology is unable to reach the single-cell level in spatial resolution, and the efficiency of gene detection is also low (142). The imaging-based ST technology is subject to limitations, including optical crowding, which restricts the number of genes that can be distinguished. Furthermore, the ST technology of in situ capture and sequencing encounters challenges in achieving genuine separation and analysis between single cells, due to factors such as varying capture point resolutions and potential transcript diffusion. The current commercial ST technology has not yet achieved single-cell resolution, and is therefore reliant on bioinformatics algorithms to predict cell types in specific regions (143).
Secondly, the sample data set obtained by ST is not comprehensive enough, which may be due to the differences in sections from different directions and levels. The majority of current ST techniques are based on thin tissue sections, which obtain the expression of genes in two-dimensional space, and still cannot achieve a truly three-dimensional ST. The primary technical challenge lies in developing methods to make the tissue transparent and image the signal within it, or to introduce an imaging probe into the tissue for reaction.
Furthermore, the economic costs and operational complexity also limit the wide application of ST. On the one hand, data processing is a significant obstacle. The complexity, dimensionality, and scale of ST data require the support of powerful computational platforms and algorithms to extract meaningful biological information from large amounts of data. Existing data processing tools and algorithms often struggle to effectively deal with noise, low signal-to-noise ratios, and non-homogeneity of spatial information in spatial data, which limits the effectiveness and feasibility of ST techniques in a wide range of applications (144, 145). On the other hand, ST is costly, especially in the stages of sample preparation, data acquisition and subsequent analysis. High-quality spatial transcriptome experiments usually require expensive equipment (e.g. high-resolution microscopes) as well as the consumption of a large number of reagents, which makes the popularization of the technology face financial barriers (146). Finally, the problem of preservation of tissue samples is also a bottleneck limiting the development of ST technology. Most tissue samples need to be processed and analyzed within a short period of time, and prolonged cryopreservation can lead to RNA degradation, which affects the quality of the data (147).
Finally, ST technology is still in the stage of continuous development, and although it has made significant progress in the field of basic research, its maturity is still low. Currently, the application of spatial transcriptome technology is mainly focused on gene expression analysis at the cellular and tissue level, and has not yet been widely used in clinical practice (148). The shortcomings of this technology include the lack of a unified standardized and regulated process, poor data comparability between different platforms, and complex and demanding experimental operations. In particular, no uniform operational specifications have been developed for sample preparation, data collection, and analysis, leading to problems with the reproducibility and accuracy of the technology. In addition, although ST has been significantly improved in terms of resolution and coverage, there is still a need to improve the spatial resolution, especially the precise expression localization at the single-cell level, in order to better reveal cellular heterogeneity and microenvironmental changes. These limitations directly affect the results of ST. In terms of data processing, the complex amount of data and the high-dimensional nature of spatial information may lead to the inability of researchers to accurately capture certain subtle biological features, affecting the interpretation of results (149). Cost constraints, on the other hand, make it difficult for many small and medium-sized laboratories or clinical institutions to afford large-scale ST experiments, further limiting the popularity and scope of application of the technology. Tissue preservation problems may prevent the timely processing of certain samples with valuable biological information, thus affecting the comprehensiveness and accuracy of the study.
Nevertheless, there are multiple innovative techniques that have the potential to overcome these challenges in the future. Firstly, improvements in imaging technologies, especially advances in super-resolution microimaging, can increase the spatial resolution and thus enhance the sensitivity and accuracy of ST techniques (150). Secondly, the application of new computational methods, such as deep learning and artificial intelligence algorithms, is expected to effectively deal with the complexity in ST data, reduce the noise, extract valuable biological information, and help researchers to mine more potential disease biomarkers from the huge amount of data. For the cost issue, the development of automated platforms can improve the efficiency of experiments and reduce the cost of a single experiment, thus promoting the popularity of ST technology. In addition, innovations in molecular preservation technologies may solve the problem of RNA degradation in tissue preservation and improve the reproducibility and stability of spatial transcriptome data (151).
Future research should focus on the following areas: first, improving the efficiency of spatial resolution and data processing, especially in multi-scale ST, how to process cellular to tissue-level data simultaneously; second, developing standardized and automated experimental processes applicable to clinical samples to reduce experimental costs and improve data reproducibility; third, exploring the combination of ST with other high-throughput genomics technologies (e.g., single-cell genomics, metabolomics) to build a more integrated multi-omics platform to comprehensively reveal biological processes (152–154).
The potential applications of ST in many types of cancer are beginning to emerge, particularly in elucidating the TME and cell-cell interactions. In addition to breast cancer, ST has demonstrated great research value in several cancer types, including lung cancer, colorectal cancer, gastric cancer and brain tumors (155–157). In lung cancer, ST has provided important insights into the TME, macrophage reprogramming, immune checkpoints and immunosuppressive changes in brain metastases, providing potential prognostic biomarkers and therapeutic strategies targeting immune and fibrotic pathways (158, 159). In colorectal cancer, on the other hand, ST has helped to investigate genetic and epigenetic co-evolution in colorectal cancer, revealing the influence of chromatin modifications and transcription factors on tumor progression and metastasis (160, 161). The cellular composition and interactions in the TME have been described in detail, identifying key immune cell and fibroblast subpopulations that play important roles in tumor progression, metastasis and prognosis, providing new targets for precision therapy (162). In gastric cancer, ST can help to study the biological properties of gastrointestinal tract tumors at different tissue levels, reveal the gene expression differences between tumors and adjacent normal tissues, and provide important information for early diagnosis and prediction of tumor recurrence (163, 164). And in brain tumors (e.g. glioblastoma), ST can help to deeply analyze the cellular heterogeneity in different regions of the tumor, and explore how the tumor cells interact with the surrounding normal neural tissue, so as to discover new therapeutic targets (165).
In addition to cancer, ST has shown immense potential in the study of other diseases (166–168). For example, in neurodegenerative diseases such as Alzheimer’s and Parkinson’s, ST can help study gene expression patterns in different regions of the brain, revealing the spatial heterogeneity of neurodegeneration and key biomarkers (169). In cardiovascular disease research, ST can analyze cell-to-cell interactions in heart tissue, uncovering the spatiotemporal dynamics of cell recovery after myocardial injury, thus providing new therapeutic strategies for heart repair and regeneration (170, 171). Furthermore, in the study of immune-related diseases and metabolic disorders such as diabetes, ST allows for a detailed analysis of cellular composition and gene expression changes in diseased tissues, offering crucial insights into disease mechanisms and precision therapies.
To translate ST technology into clinical practice, the next key step is to achieve its high-throughput and cost-effective standardized application (172). This can be accomplished by developing an ST-based diagnostic platform and integrating it with existing clinical pathology evaluation methods, thereby enhancing the accuracy of early disease diagnosis. Additionally, for research on specific diseases such as cancer and neurodegenerative diseases, greater emphasis should be placed on integrating ST with clinical trials to assess its potential in personalized treatment. Through these practical measures, ST will be able to overcome current limitations and advance the development of precision medicine.
In summary, ST offers the potential to elucidate the spatial information of gene expression with different precision, thereby opening avenues for investigating tumor heterogeneity and TME, and overcoming tumor drug resistance. It is of great importance to identify new biomarkers, study new therapeutic targets and improve the prediction of tumor prognosis. The wide application of ST can improve our comprehension of breast cancer and facilitate the implementation of personalized and precise therapeutic strategies for patients, which will have a profound impact on the field of breast cancer research.
Author contributions
JA: Conceptualization, Data curation, Visualization, Writing – original draft. YL: Data curation, Visualization, Writing – review & editing. YXC: Data curation, Writing – review & editing. YLC: Data curation, Writing – review & editing. ZZ: Writing – review & editing. JC: Writing – review & editing. CP: Data curation, Funding acquisition, Writing – review & editing. RH: Data curation, Writing – review & editing. FP: Conceptualization, Funding acquisition, Supervision, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. The study was supported by National Natural Science Foundation of Science and Technology Department of Sichuan Province (no.2023NSFSC1928; 2023NSFSC1992), Project of State Administration of Traditional Chinese Medicine of China (no.ZYYCXTD-D-202209), Project of Sichuan Province Pharmaceutical Administration (no.2022C001), the Open Research Fund of State Key Laboratory of Southwestern Chinese Medicine Resources (no.SKLTCM202401) and the Fundamental Research Funds for the central Universities.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2024.1499301/full#supplementary-material
Abbreviations
AI, Artificial intelligence; ADCs, antibody-drug conjugates; CAFs, cancer-associated fibroblasts; CNN, convolutional neural network; CNV, copy number variation; DCIS, ductal carcinoma in situ; DEGs, differentially expressed genes; DNB, DNA nanospheres; DSP, digital spatial profiling; ECM, extracellular matrix; EMO, expression-morphology; FFA, free fatty acids; EMT, epithelial-mesenchymal transition; FFA, free fatty acid; FISH, fromfluorescently-labeled RNA sequences; GEO-seq, geographical position sequencing; HDST, igh definition spatial transcriptomics; HRD, homologous recombination defect; ICIs, immune checkpoint inhibitors; IFN, interferon; IMPC, invasive micropapillary carcinomal; ISH, in situ hybridization; ISS, in situ sequencing; ITH, intratumoral heterogeneity; K15+, keratin 15+; LAMs, lipid-associated macrophages; LCM, laser capture microdissection; Lock-Seq, lock-probe in situ sequencing; NGS, next-generation sequencing; OCT, optimal cutting temperature compound; OGC, osteoclast-like giant cells; oTME, organ-type TME; pCR, pathological complete response; PDX, patient-derived xenograft; pNR, non-response progression; QCCs, quiescent cancer cells; RNA-seq, transcriptomics sequencing; scRNA-seq, single-cell transcriptomics sequencing; smFISH, single-molecule RNA fluorescence in situ hybridization; ST, spatial transcriptomics; TIVA, transcriptomics in vivo analysis; TME, tumor microenvironment; tomo-seq, RNA tomography; ZNF689, zinc finger protein 689.
References
1. Mosele MF, Westphalen CB, Stenzinger A, Barlesi F, Bayle A, Bièche I, et al. Recommendations for the use of next-generation sequencing (Ngs) for patients with advanced cancer in 2024: A report from the esmo precision medicine working group. Ann Oncol. (2024) 35:588–606. doi: 10.1016/j.annonc.2024.04.005
2. Han H, Nutiu R, Moffat J, Blencowe BJ. Snapshot: high-throughput sequencing applications. Cell. (2011) 146:1044,.e1–2. doi: 10.1016/j.cell.2011.09.002
3. Wang Z, Gerstein M, Snyder M. Rna-seq: A revolutionary tool for transcriptomics. Nat Rev Genet. (2009) 10:57–63. doi: 10.1038/nrg2484
4. Stark R, Grzelak M, Hadfield J. Rna sequencing: the teenage years. Nat Rev Genet. (2019) 20:631–56. doi: 10.1038/s41576-019-0150-2
5. Tang F, Barbacioru C, Wang Y, Nordman E, Lee C, Xu N, et al. Mrna-seq whole-transcriptome analysis of a single cell. Nat Methods. (2009) 6:377–82. doi: 10.1038/nmeth.1315
6. Liu J, Zhang Q, Zhao W, Guo J, Wong YK, Zhang C, et al. Single-cell rna sequencing reveals the effects of hederasaponin C in the treatment of diabetic nephropathy. Acta Materia Med. (2023) 2:449–65. doi: 10.15212/AMM-2023-0031
7. Ståhl PL, Salmén F, Vickovic S, Lundmark A, Navarro JF, Magnusson J, et al. Visualization and analysis of gene expression in tissue sections by spatial transcriptomics. Science. (2016) 353:78–82. doi: 10.1126/science.aaf2403
8. Moor AE, Itzkovitz S. Spatial transcriptomics: paving the way for tissue-level systems biology. Curr Opin Biotechnol. (2017) 46:126–33. doi: 10.1016/j.copbio.2017.02.004
9. Marx V. Method of the year: spatially resolved transcriptomics. Nat Methods. (2021) 18:9–14. doi: 10.1038/s41592-020-01033-y
10. Bahcecioglu G, Basara G, Ellis BW, Ren X, Zorlutuna P. Breast cancer models: engineering the tumor microenvironment. Acta Biomater. (2020) 106:1–21. doi: 10.1016/j.actbio.2020.02.006
11. Carlberg K, Korotkova M, Larsson L, Catrina AI, Ståhl PL, Malmström V. Exploring inflammatory signatures in arthritic joint biopsies with spatial transcriptomics. Sci Rep. (2019) 9:18975. doi: 10.1038/s41598-019-55441-y
12. Nolan E, Lindeman GJ, Visvader JE. Deciphering breast cancer: from biology to the clinic. Cell. (2023) 186:1708–28. doi: 10.1016/j.cell.2023.01.040
13. Gall JG, Pardue ML. Formation and detection of rna-DNA hybrid molecules in cytological preparations. Proc Natl Acad Sci U.S.A. (1969) 63:378–83. doi: 10.1073/pnas.63.2.378
14. Kwon S. Single-molecule fluorescence in situ hybridization: quantitative imaging of single rna molecules. BMB Rep. (2013) 46:65–72. doi: 10.5483/bmbrep.2013.46.2.016
15. Combs PA, Eisen MB. Sequencing mrna from cryo-sliced drosophila embryos to determine genome-wide spatial patterns of gene expression. PloS One. (2013) 8:e71820. doi: 10.1371/journal.pone.0071820
16. Junker JP, Noël ES, Guryev V, Peterson KA, Shah G, Huisken J, et al. Genome-wide rna tomography in the zebrafish embryo. Cell. (2014) 159:662–75. doi: 10.1016/j.cell.2014.09.038
17. Simone NL, Bonner RF, Gillespie JW, Emmert-Buck MR, Liotta LA. Laser-capture microdissection: opening the microscopic frontier to molecular analysis. Trends Genet. (1998) 14:272–6. doi: 10.1016/s0168-9525(98)01489-9
18. Rao BH, Souček P, Hlaváč V. Laser capture microdissection: A gear for pancreatic cancer research. Int J Mol Sci. (2022) 23:14566. doi: 10.3390/ijms232314566
19. Zeng H, Huang J, Ren J, Wang CK, Tang Z, Zhou H, et al. Spatially resolved single-cell translatomics at molecular resolution. Science. (2023) 380:eadd3067. doi: 10.1126/science.add3067
20. Nichterwitz S, Chen G, Aguila Benitez J, Yilmaz M, Storvall H, Cao M, et al. Laser capture microscopy coupled with smart-seq2 for precise spatial transcriptomic profiling. Nat Commun. (2016) 7:12139. doi: 10.1038/ncomms12139
21. Chen J, Suo S, Tam PP, Han JJ, Peng G, Jing N. Spatial transcriptomic analysis of cryosectioned tissue samples with geo-seq. Nat Protoc. (2017) 12:566–80. doi: 10.1038/nprot.2017.003
22. Lovatt D, Ruble BK, Lee J, Dueck H, Kim TK, Fisher S, et al. Transcriptome in vivo analysis (Tiva) of spatially defined single cells in live tissue. Nat Methods. (2014) 11:190–6. doi: 10.1038/nmeth.2804
23. Medaglia C, Giladi A, Stoler-Barak L, De Giovanni M, Salame TM, Biram A, et al. Spatial reconstruction of immune niches by combining photoactivatable reporters and scrna-seq. Science. (2017) 358:1622–6. doi: 10.1126/science.aao4277
24. Boisset JC, Vivié J, Grün D, Muraro MJ, Lyubimova A, van Oudenaarden A. Mapping the physical network of cellular interactions. Nat Methods. (2018) 15:547–53. doi: 10.1038/s41592-018-0009-z
25. Buongiorno-Nardelli M, Amaldi F. Autoradiographic detection of molecular hybrids between rna and DNA in tissue sections. Nature. (1970) 225:946–8. doi: 10.1038/225946a0
26. Femino AM, Fay FS, Fogarty K, Singer RH. Visualization of single rna transcripts in situ. Science. (1998) 280:585–90. doi: 10.1126/science.280.5363.585
27. Wang F, Flanagan J, Su N, Wang LC, Bui S, Nielson A, et al. Rnascope: A novel in situ rna analysis platform for formalin-fixed, paraffin-embedded tissues. J Mol Diagn. (2012) 14:22–9. doi: 10.1016/j.jmoldx.2011.08.002
28. Chen KH, Boettiger AN, Moffitt JR, Wang S, Zhuang X. Rna imaging. Spatially resolved, highly multiplexed rna profiling in single cells. Science. (2015) 348:aaa6090. doi: 10.1126/science.aaa6090
29. Eng CL, Lawson M, Zhu Q, Dries R, Koulena N, Takei Y, et al. Transcriptome-scale super-resolved imaging in tissues by rna seqfish. Nature. (2019) 568:235–9. doi: 10.1038/s41586-019-1049-y
30. Lubeck E, Coskun AF, Zhiyentayev T, Ahmad M, Cai L. Single-cell in situ rna profiling by sequential hybridization. Nat Methods. (2014) 11:360–1. doi: 10.1038/nmeth.2892
31. Weinstein JA, Regev A, Zhang F. DNA microscopy: optics-free spatio-genetic imaging by a stand-alone chemical reaction. Cell. (2019) 178:229–41.e16. doi: 10.1016/j.cell.2019.05.019
32. Weinstein J. Refining the lens of DNA microscopy. Nat Comput Sci. (2024) 4:94–5. doi: 10.1038/s43588-024-00601-y
33. Ke R, Mignardi M, Pacureanu A, Svedlund J, Botling J, Wählby C, et al. In situ sequencing for rna analysis in preserved tissue and cells. Nat Methods. (2013) 10:857–60. doi: 10.1038/nmeth.2563
34. Alon S, Goodwin DR, Sinha A, Wassie AT, Chen F, Daugharthy ER, et al. Expansion sequencing: spatially precise in situ transcriptomics in intact biological systems. Science. (2021) 371:eaax2656. doi: 10.1126/science.aax2656
35. Lee JH, Daugharthy ER, Scheiman J, Kalhor R, Ferrante TC, Terry R, et al. Fluorescent in situ sequencing (Fisseq) of rna for gene expression profiling in intact cells and tissues. Nat Protoc. (2015) 10:442–58. doi: 10.1038/nprot.2014.191
36. Wang X, Allen WE, Wright MA, Sylwestrak EL, Samusik N, Vesuna S, et al. Three-dimensional intact-tissue sequencing of single-cell transcriptional states. Science. (2018) 361:eaat5691. doi: 10.1126/science.aat5691
37. Chen X, Sun YC, Church GM, Lee JH, Zador AM. Efficient in situ barcode sequencing using padlock probe-based baristaseq. Nucleic Acids Res. (2018) 46:e22. doi: 10.1093/nar/gkx1206
38. Fürth D, Hatini V, Lee JH. In situ transcriptome accessibility sequencing (Insta-seq). bioRxiv [Preprint]. (2019). Available online at: https://www.biorxiv.org/content/10.1101/722819v2 (accessed August 06, 2019).
39. Gyllborg D, Langseth CM, Qian X, Choi E, Salas SM, Hilscher MM, et al. Correction to 'Hybridization-based in situ sequencing (Hybiss) for spatially resolved transcriptomics in human and mouse brain tissue'. Nucleic Acids Res. (2021) 49:9002. doi: 10.1093/nar/gkab630
40. Jemt A, Salmén F, Lundmark A, Mollbrink A, Fernández Navarro J, Ståhl PL, et al. An automated approach to prepare tissue-derived spatially barcoded rna-sequencing libraries. Sci Rep. (2016) 6:37137. doi: 10.1038/srep37137
41. Wang Y, Ma S, Ruzzo WL. Spatial modeling of prostate cancer metabolic gene expression reveals extensive heterogeneity and selective vulnerabilities. Sci Rep. (2020) 10:3490. doi: 10.1038/s41598-020-60384-w
42. Liu Y, Yang M, Deng Y, Su G, Enninful A, Guo CC, et al. High-spatial-resolution multi-omics sequencing via deterministic barcoding in tissue. Cell. (2020) 183:1665–81.e18. doi: 10.1016/j.cell.2020.10.026
43. Rodriques SG, Stickels RR, Goeva A, Martin CA, Murray E, Vanderburg CR, et al. Slide-seq: A scalable technology for measuring genome-wide expression at high spatial resolution. Science. (2019) 363:1463–7. doi: 10.1126/science.aaw1219
44. Vickovic S, Eraslan G, Salmén F, Klughammer J, Stenbeck L, Schapiro D, et al. High-definition spatial transcriptomics for in situ tissue profiling. Nat Methods. (2019) 16:987–90. doi: 10.1038/s41592-019-0548-y
45. Cho CS, Xi J, Si Y, Park SR, Hsu JE, Kim M, et al. Microscopic examination of spatial transcriptome using seq-scope. Cell. (2021) 184:3559–72.e22. doi: 10.1016/j.cell.2021.05.010
46. Srivatsan SR, Regier MC, Barkan E, Franks JM, Packer JS, Grosjean P, et al. Embryo-scale, single-cell spatial transcriptomics. Science. (2021) 373:111–7. doi: 10.1126/science.abb9536
47. Chen A, Liao S, Cheng M, Ma K, Wu L, Lai Y, et al. Spatiotemporal transcriptomic atlas of mouse organogenesis using DNA nanoball-patterned arrays. Cell. (2022) 185:1777–92.e21. doi: 10.1016/j.cell.2022.04.003
48. Satija R, Farrell JA, Gennert D, Schier AF, Regev A. Spatial reconstruction of single-cell gene expression data. Nat Biotechnol. (2015) 33:495–502. doi: 10.1038/nbt.3192
49. Butler A, Hoffman P, Smibert P, Papalexi E, Satija R. Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat Biotechnol. (2018) 36:411–20. doi: 10.1038/nbt.4096
50. Hao Y, Stuart T, Kowalski MH, Choudhary S, Hoffman P, Hartman A, et al. Dictionary learning for integrative, multimodal and scalable single-cell analysis. Nat Biotechnol. (2024) 42:293–304. doi: 10.1038/s41587-023-01767-y
51. Gao R, Yuan X, Ma Y, Wei T, Johnston L, Shao Y, et al. Harnessing tme depicted by histological images to improve cancer prognosis through a deep learning system. Cell Rep Med. (2024) 5:101536. doi: 10.1016/j.xcrm.2024.101536
52. Yoosuf N, Navarro JF, Salmén F, Ståhl PL, Daub CO. Identification and transfer of spatial transcriptomics signatures for cancer diagnosis. Breast Cancer Res. (2020) 22:6. doi: 10.1186/s13058-019-1242-9
53. He S, Jin Y, Nazaret A, Shi L, Chen X, Rampersaud S, et al. Starfysh integrates spatial transcriptomic and histologic data to reveal heterogeneous tumor-immune hubs. Nat Biotechnol. (2024). doi: 10.1038/s41587-024-02173-8
54. Xu C, Jin X, Wei S, Wang P, Luo M, Xu Z, et al. Deepst: identifying spatial domains in spatial transcriptomics by deep learning. Nucleic Acids Res. (2022) 50:e131. doi: 10.1093/nar/gkac901
55. Tanevski J, Flores ROR, Gabor A, Schapiro D, Saez-Rodriguez J. Explainable multiview framework for dissecting spatial relationships from highly multiplexed data. Genome Biol. (2022) 23:97. doi: 10.1186/s13059-022-02663-5
56. Liu N, Bhuva DD, Mohamed A, Bokelund M, Kulasinghe A, Tan CW, et al. Standr: spatial transcriptomic analysis for geomx dsp data. Nucleic Acids Res. (2024) 52:e2. doi: 10.1093/nar/gkad1026
57. de Haan K, Zhang Y, Zuckerman JE, Liu T, Sisk AE, Diaz MFP, et al. Deep learning-based transformation of H&E stained tissues into special stains. Nat Commun. (2021) 12:4884. doi: 10.1038/s41467-021-25221-2
58. Saldanha OL, Quirke P, West NP, James JA, Loughrey MB, Grabsch HI, et al. Swarm learning for decentralized artificial intelligence in cancer histopathology. Nat Med. (2022) 28:1232–9. doi: 10.1038/s41591-022-01768-5
59. Shamai G, Livne A, Polónia A, Sabo E, Cretu A, Bar-Sela G, et al. Deep learning-based image analysis predicts pd-L1 status from H&E-stained histopathology images in breast cancer. Nat Commun. (2022) 13:6753. doi: 10.1038/s41467-022-34275-9
60. Zubair A, Chapple RH, Natarajan S, Wright WC, Pan M, Lee HM, et al. Cell type identification in spatial transcriptomics data can be improved by leveraging cell-type-informative paired tissue images using a bayesian probabilistic model. Nucleic Acids Res. (2022) 50:e80. doi: 10.1093/nar/gkac320
61. Baccin C, Al-Sabah J, Velten L, Helbling PM, Grünschläger F, Hernández-Malmierca P, et al. Combined single-cell and spatial transcriptomics reveal the molecular, cellular and spatial bone marrow niche organization. Nat Cell Biol. (2020) 22:38–48. doi: 10.1038/s41556-019-0439-6
62. Hafemeister C, Satija R. Normalization and variance stabilization of single-cell rna-seq data using regularized negative binomial regression. Genome Biol. (2019) 20:296. doi: 10.1186/s13059-019-1874-1
63. Wang Y, Liu B, Zhao G, Lee Y, Buzdin A, Mu X, et al. Spatial transcriptomics: technologies, applications and experimental considerations. Genomics. (2023) 115:110671. doi: 10.1016/j.ygeno.2023.110671
64. Xu H, Fu H, Long Y, Ang KS, Sethi R, Chong K, et al. Unsupervised spatially embedded deep representation of spatial transcriptomics. Genome Med. (2024) 16:12. doi: 10.1186/s13073-024-01283-x
65. Li Y, Luo Y. Stdgcn: spatial transcriptomic cell-type deconvolution using graph convolutional networks. Genome Biol. (2024) 25:206. doi: 10.1186/s13059-024-03353-0
66. Marconato L, Palla G, Yamauchi KA, Virshup I, Heidari E, Treis T, et al. Spatialdata: an open and universal data framework for spatial omics. Nat Methods. (2024). doi: 10.1038/s41592-024-02212-x
67. Shi H, He Y, Zhou Y, Huang J, Maher K, Wang B, et al. Spatial atlas of the mouse central nervous system at molecular resolution. Nature. (2023) 622:552–61. doi: 10.1038/s41586-023-06569-5
68. Moncada R, Barkley D, Wagner F, Chiodin M, Devlin JC, Baron M, et al. Integrating microarray-based spatial transcriptomics and single-cell rna-seq reveals tissue architecture in pancreatic ductal adenocarcinomas. Nat Biotechnol. (2020) 38:333–42. doi: 10.1038/s41587-019-0392-8
69. Qi J, Sun H, Zhang Y, Wang Z, Xun Z, Li Z, et al. Single-cell and spatial analysis reveal interaction of fap(+) fibroblasts and spp1(+) macrophages in colorectal cancer. Nat Commun. (2022) 13:1742. doi: 10.1038/s41467-022-29366-6
70. Grün D, van Oudenaarden A. Design and analysis of single-cell sequencing experiments. Cell. (2015) 163:799–810. doi: 10.1016/j.cell.2015.10.039
71. Longo SK, Guo MG, Ji AL, Khavari PA. Integrating single-cell and spatial transcriptomics to elucidate intercellular tissue dynamics. Nat Rev Genet. (2021) 22:627–44. doi: 10.1038/s41576-021-00370-8
72. Pelka K, Hofree M, Chen JH, Sarkizova S, Pirl JD, Jorgji V, et al. Spatially organized multicellular immune hubs in human colorectal cancer. Cell. (2021) 184:4734–52.e20. doi: 10.1016/j.cell.2021.08.003
73. Suo C, Dann E, Goh I, Jardine L, Kleshchevnikov V, Park JE, et al. Mapping the developing human immune system across organs. Science. (2022) 376:eabo0510. doi: 10.1126/science.abo0510
74. Turajlic S, Sottoriva A, Graham T, Swanton C. Author correction: resolving genetic heterogeneity in cancer. Nat Rev Genet. (2020) 21:65. doi: 10.1038/s41576-019-0188-1
75. Wang J, Wei J, Inuzuka H. Aging and cancer hallmarks as therapeutic targets. Acta Materia Med. (2023) 2:281–4. doi: 10.15212/AMM-2023-0018
76. Polyak K. Heterogeneity in breast cancer. J Clin Invest. (2011) 121:3786–8. doi: 10.1172/jci60534
77. Marusyk A, Almendro V, Polyak K. Intra-tumour heterogeneity: A looking glass for cancer? Nat Rev Cancer. (2012) 12:323–34. doi: 10.1038/nrc3261
78. Zardavas D, Irrthum A, Swanton C, Piccart M. Clinical management of breast cancer heterogeneity. Nat Rev Clin Oncol. (2015) 12:381–94. doi: 10.1038/nrclinonc.2015.73
79. Zeng Q, Mousa M, Nadukkandy AS, Franssens L, Alnaqbi H, Alshamsi FY, et al. Understanding tumour endothelial cell heterogeneity and function from single-cell omics. Nat Rev Cancer. (2023) 23:544–64. doi: 10.1038/s41568-023-00591-5
80. Bassiouni R, Idowu MO, Gibbs LD, Robila V, Grizzard PJ, Webb MG, et al. Spatial transcriptomic analysis of a diverse patient cohort reveals a conserved architecture in triple-negative breast cancer. Cancer Res. (2023) 83:34–48. doi: 10.1158/0008-5472.Can-22-2682
81. Coutant A, Cockenpot V, Muller L, Degletagne C, Pommier R, Tonon L, et al. Spatial transcriptomics reveal pitfalls and opportunities for the detection of rare high-plasticity breast cancer subtypes. Lab Invest. (2023) 103:100258. doi: 10.1016/j.labinv.2023.100258
82. Svedlund J, Strell C, Qian X, Zilkens KJC, Tobin NP, Bergh J, et al. Generation of in situ sequencing based oncomaps to spatially resolve gene expression profiles of diagnostic and prognostic markers in breast cancer. EBioMedicine. (2019) 48:212–23. doi: 10.1016/j.ebiom.2019.09.009
83. Wu SZ, Al-Eryani G, Roden DL, Junankar S, Harvey K, Andersson A, et al. A single-cell and spatially resolved atlas of human breast cancers. Nat Genet. (2021) 53:1334–47. doi: 10.1038/s41588-021-00911-1
84. Ge LP, Jin X, Ma D, Wang ZY, Liu CL, Zhou CZ, et al. Znf689 deficiency promotes intratumor heterogeneity and immunotherapy resistance in triple-negative breast cancer. Cell Res. (2024) 34:58–75. doi: 10.1038/s41422-023-00909-w
85. Janesick A, Shelansky R, Gottscho AD, Wagner F, Williams SR, Rouault M, et al. High resolution mapping of the tumor microenvironment using integrated single-cell, spatial and in situ analysis. Nat Commun. (2023) 14:8353. doi: 10.1038/s41467-023-43458-x
86. Yoshitake R, Mori H, Ha D, Wu X, Wang J, Wang X, et al. Molecular features of luminal breast cancer defined through spatial and single-cell transcriptomics. Clin Transl Med. (2024) 14:e1548. doi: 10.1002/ctm2.1548
87. Liu SQ, Gao ZJ, Wu J, Zheng HM, Li B, Sun S, et al. Single-cell and spatially resolved analysis uncovers cell heterogeneity of breast cancer. J Hematol Oncol. (2022) 15:19. doi: 10.1186/s13045-022-01236-0
88. Kohler KT, Goldhammer N, Demharter S, Pfisterer U, Khodosevich K, Rønnov-Jessen L, et al. Ductal keratin 15(+) luminal progenitors in normal breast exhibit a basal-like breast cancer transcriptomic signature. NPJ Breast Cancer. (2022) 8:81. doi: 10.1038/s41523-022-00444-8
89. Abdalla M, Tran-Thanh D, Moreno J, Iakovlev V, Nair R, Kanwar N, et al. Mapping genomic and transcriptomic alterations spatially in epithelial cells adjacent to human breast carcinoma. Nat Commun. (2017) 8:1245. doi: 10.1038/s41467-017-01357-y
90. Lu P, Foley J, Zhu C, McNamara K, Sirinukunwattana K, Vennam S, et al. Transcriptome and genome evolution during her2-amplified breast neoplasia. Breast Cancer Res. (2021) 23:73. doi: 10.1186/s13058-021-01451-6
91. Henon C, Vibert J, Eychenne T, Gruel N, Colmet-Daage L, Ngo C, et al. Single-cell multiomics profiling reveals heterogeneous transcriptional programs and microenvironment in dsrcts. Cell Rep Med. (2024) 5:101582. doi: 10.1016/j.xcrm.2024.101582
92. Liu YM, Ge JY, Chen YF, Liu T, Chen L, Liu CC, et al. Combined single-cell and spatial transcriptomics reveal the metabolic evolvement of breast cancer during early dissemination. Adv Sci (Weinh). (2023) 10:e2205395. doi: 10.1002/advs.202205395
93. Lv J, Shi Q, Han Y, Li W, Liu H, Zhang J, et al. Spatial transcriptomics reveals gene expression characteristics in invasive micropapillary carcinoma of the breast. Cell Death Dis. (2021) 12:1095. doi: 10.1038/s41419-021-04380-6
94. He B, Bergenstråhle L, Stenbeck L, Abid A, Andersson A, Borg Å, et al. Integrating spatial gene expression and breast tumour morphology via deep learning. Nat BioMed Eng. (2020) 4:827–34. doi: 10.1038/s41551-020-0578-x
95. Rahaman MM, Millar EKA, Meijering E. Breast cancer histopathology image-based gene expression prediction using spatial transcriptomics data and deep learning. Sci Rep. (2023) 13:13604. doi: 10.1038/s41598-023-40219-0
96. Wang Y, Kartasalo K, Weitz P, Ács B, Valkonen M, Larsson C, et al. Predicting molecular phenotypes from histopathology images: A transcriptome-wide expression-morphology analysis in breast cancer. Cancer Res. (2021) 81:5115–26. doi: 10.1158/0008-5472.Can-21-0482
97. Jin MZ, Jin WL. The updated landscape of tumor microenvironment and drug repurposing. Signal Transduct Target Ther. (2020) 5:166. doi: 10.1038/s41392-020-00280-x
98. Han M, Li S, Fan H, An J, Peng C, Peng F. Regulated cell death in glioma: promising targets for natural small-molecule compounds. Front Oncol. (2024) 14:1273841. doi: 10.3389/fonc.2024.1273841
99. Bergholtz H, Carter JM, Cesano A, Cheang MCU, Church SE, Divakar P, et al. Best practices for spatial profiling for breast cancer research with the geomx(®) digital spatial profiler. Cancers (Basel). (2021) 13:4456. doi: 10.3390/cancers13174456
100. Schlam I, Church SE, Hether TD, Chaldekas K, Hudson BM, White AM, et al. The tumor immune microenvironment of primary and metastatic her2- positive breast cancers utilizing gene expression and spatial proteomic profiling. J Transl Med. (2021) 19:480. doi: 10.1186/s12967-021-03113-9
101. Wu Q, Liu Z, Gao Z, Luo Y, Li F, Yang C, et al. Klf5 inhibition potentiates anti-pd1 efficacy by enhancing cd8(+) T-cell-dependent antitumor immunity. Theranostics. (2023) 13:1381–400. doi: 10.7150/thno.82182
102. Chen J, Larsson L, Swarbrick A, Lundeberg J. Spatial landscapes of cancers: insights and opportunities. Nat Rev Clin Oncol. (2024) 21:660–74. doi: 10.1038/s41571-024-00926-7
103. Mao X, Zhou D, Lin K, Zhang B, Gao J, Ling F, et al. Single-cell and spatial transcriptome analyses revealed cell heterogeneity and immune environment alternations in metastatic axillary lymph nodes in breast cancer. Cancer Immunol Immunother. (2023) 72:679–95. doi: 10.1007/s00262-022-03278-2
104. Nasir I, McGuinness C, Poh AR, Ernst M, Darcy PK, Britt KL. Tumor macrophage functional heterogeneity can inform the development of novel cancer therapies. Trends Immunol. (2023) 44:971–85. doi: 10.1016/j.it.2023.10.007
105. Ben-Chetrit N, Niu X, Sotelo J, Swett AD, Rajasekhar VK, Jiao MS, et al. Breast cancer macrophage heterogeneity and self-renewal are determined by spatial localization. bioRxiv. (2023). doi: 10.1101/2023.10.24.563749
106. Xu K, Wang R, Chen Q, Liu Y, Li X, Mao L, et al. Microenvironment components and spatially resolved single-cell transcriptome atlas of breast cancer metastatic axillary lymph nodes. Acta Biochim Biophys Sin (Shanghai). (2022) 54:1336–48. doi: 10.3724/abbs.2022131
107. Tzeng YT, Hsiao JH, Chu PY, Tseng LM, Hou MF, Tsang YL, et al. The role of lsm1 in breast cancer: shaping metabolism and tumor-associated macrophage infiltration. Pharmacol Res. (2023) 198:107008. doi: 10.1016/j.phrs.2023.107008
108. Liu Z, Gao Z, Li B, Li J, Ou Y, Yu X, et al. Lipid-associated macrophages in the tumor-adipose microenvironment facilitate breast cancer progression. Oncoimmunology. (2022) 11:2085432. doi: 10.1080/2162402x.2022.2085432
109. Xue L, Wang C, Qian Y, Zhu W, Liu L, Yang X, et al. Tryptophan metabolism regulates inflammatory macrophage polarization as a predictive factor for breast cancer immunotherapy. Int Immunopharmacol. (2023) 125:111196. doi: 10.1016/j.intimp.2023.111196
110. Romanens L, Chaskar P, Marcone R, Ryser S, Tille JC, Genolet R, et al. Clonal expansion of intra-epithelial T cells in breast cancer revealed by spatial transcriptomics. Int J Cancer. (2023) 153:1568–78. doi: 10.1002/ijc.34620
111. Croizer H, Mhaidly R, Kieffer Y, Gentric G, Djerroudi L, Leclere R, et al. Deciphering the spatial landscape and plasticity of immunosuppressive fibroblasts in breast cancer. Nat Commun. (2024) 15:2806. doi: 10.1038/s41467-024-47068-z
112. Ma C, Yang C, Peng A, Sun T, Ji X, Mi J, et al. Pan-cancer spatially resolved single-cell analysis reveals the crosstalk between cancer-associated fibroblasts and tumor microenvironment. Mol Cancer. (2023) 22:170. doi: 10.1186/s12943-023-01876-x
113. Kang K, Wu Y, Han C, Wang L, Wang Z, Zhao A. Homologous recombination deficiency in triple-negative breast cancer: multi-scale transcriptomics reveals distinct tumor microenvironments and limitations in predicting immunotherapy response. Comput Biol Med. (2023) 158:106836. doi: 10.1016/j.compbiomed.2023.106836
114. Inayatullah M, Mahesh A, Turnbull AK, Dixon JM, Natrajan R, Tiwari VK. Basal-epithelial subpopulations underlie and predict chemotherapy resistance in triple-negative breast cancer. EMBO Mol Med. (2024) 16:823–53. doi: 10.1038/s44321-024-00050-0
115. Li CJ, Tzeng YT, Hsiao JH, Tseng LM, Hsu TS, Chu PY. Spatial and single-cell explorations uncover prognostic significance and immunological functions of mitochondrial calcium uniporter in breast cancer. Cancer Cell Int. (2024) 24:140. doi: 10.1186/s12935-024-03327-z
116. Emens LA, Romero PJ, Anderson AC, Bruno TC, Capitini CM, Collyar D, et al. Challenges and opportunities in cancer immunotherapy: A society for immunotherapy of cancer (Sitc) strategic vision. J Immunother Cancer. (2024) 12:e009063. doi: 10.1136/jitc-2024-009063
117. Gao ZJ, Fang Z, Yuan JP, Sun SR, Li B. Integrative multi-omics analyses unravel the immunological implication and prognostic significance of cxcl12 in breast cancer. Front Immunol. (2023) 14:1188351. doi: 10.3389/fimmu.2023.1188351
118. Tashireva LA, Kalinchuk AY, Gerashchenko TS, Menyailo M, Khozyainova A, Denisov EV, et al. Spatial profile of tumor microenvironment in pd-L1-negative and pd-L1-positive triple-negative breast cancer. Int J Mol Sci. (2023) 24:1433. doi: 10.3390/ijms24021433
119. Hammerl D, Martens JWM, Timmermans M, Smid M, Trapman-Jansen AM, Foekens R, et al. Spatial immunophenotypes predict response to anti-pd1 treatment and capture distinct paths of T cell evasion in triple negative breast cancer. Nat Commun. (2021) 12:5668. doi: 10.1038/s41467-021-25962-0
120. Huang P, Zhou X, Zheng M, Yu Y, Jin G, Zhang S. Regulatory T cells are associated with the tumor immune microenvironment and immunotherapy response in triple-negative breast cancer. Front Immunol. (2023) 14:1263537. doi: 10.3389/fimmu.2023.1263537
121. Zhu Y, Banerjee A, Xie P, Ivanov AA, Uddin A, Jiao Q, et al. Pharmacological suppression of the otud4/cd73 proteolytic axis revives antitumor immunity against immune-suppressive breast cancers. J Clin Invest. (2024) 134:e176390. doi: 10.1172/jci176390
122. Zhang B, Shi J, Shi X, Xu X, Gao L, Li S, et al. Development and evaluation of a human cd47/her2 bispecific antibody for trastuzumab-resistant breast cancer immunotherapy. Drug Resist Update. (2024) 74:101068. doi: 10.1016/j.drup.2024.101068
123. Cheung A, Chenoweth AM, Johansson A, Laddach R, Guppy N, Trendell J, et al. Anti-egfr antibody-drug conjugate carrying an inhibitor targeting cdk restricts triple-negative breast cancer growth. Clin Cancer Res. (2024) 30:3298–315. doi: 10.1158/1078-0432.Ccr-23-3110
124. Anderson WF, Rosenberg PS, Prat A, Perou CM, Sherman ME. How many etiological subtypes of breast cancer: two, three, four, or more? J Natl Cancer Inst. (2014) 106:dju165. doi: 10.1093/jnci/dju165
125. Peng F, Xiong L, Tang H, Peng C, Chen J. Regulation of epithelial-mesenchymal transition through micrornas: clinical and biological significance of micrornas in breast cancer. Tumour Biol. (2016) 37:14463–77. doi: 10.1007/s13277-016-5334-1
126. Guiu S, Michiels S, André F, Cortes J, Denkert C, Di Leo A, et al. Molecular subclasses of breast cancer: how do we define them? The impakt 2012 working group statement. Ann Oncol. (2012) 23:2997–3006. doi: 10.1093/annonc/mds586
127. An J, Peng C, Xie X, Peng F. New advances in targeted therapy of her2-negative breast cancer. Front Oncol. (2022) 12:828438. doi: 10.3389/fonc.2022.828438
128. Chen P, Zhang X, Ding R, Yang L, Lyu X, Zeng J, et al. Patient-derived organoids can guide personalized-therapies for patients with advanced breast cancer. Adv Sci (Weinh). (2021) 8:e2101176. doi: 10.1002/advs.202101176
129. Kumar T, Nee K, Wei R, He S, Nguyen QH, Bai S, et al. A spatially resolved single-cell genomic atlas of the adult human breast. Nature. (2023) 620:181–91. doi: 10.1038/s41586-023-06252-9
130. Donati B, Reggiani F, Torricelli F, Santandrea G, Rossi T, Bisagni A, et al. Spatial distribution of immune cells drives resistance to neoadjuvant chemotherapy in triple-negative breast cancer. Cancer Immunol Res. (2024) 12:120–34. doi: 10.1158/2326-6066.Cir-23-0076
131. Shiao SL, Gouin KH 3rd, Ing N, Ho A, Basho R, Shah A, et al. Single-cell and spatial profiling identify three response trajectories to pembrolizumab and radiation therapy in triple negative breast cancer. Cancer Cell. (2024) 42:70–84.e8. doi: 10.1016/j.ccell.2023.12.012
132. Baldominos P, Barbera-Mourelle A, Barreiro O, Huang Y, Wight A, Cho JW, et al. Quiescent cancer cells resist T cell attack by forming an immunosuppressive niche. Cell. (2022) 185:1694–708.e19. doi: 10.1016/j.cell.2022.03.033
133. Greenwald AC, Darnell NG, Hoefflin R, Simkin D, Mount CW, Gonzalez Castro LN, et al. Integrative spatial analysis reveals a multi-layered organization of glioblastoma. Cell. (2024) 187:2485–501.e26. doi: 10.1016/j.cell.2024.03.029
134. Lein E, Borm LE, Linnarsson S. The promise of spatial transcriptomics for neuroscience in the era of molecular cell typing. Science. (2017) 358:64–9. doi: 10.1126/science.aan6827
135. Liu H, Gao J, Feng M, Cheng J, Tang Y, Cao Q, et al. Integrative molecular and spatial analysis reveals evolutionary dynamics and tumor-immune interplay of in situ and invasive acral melanoma. Cancer Cell. (2024) 42:1067–85.e11. doi: 10.1016/j.ccell.2024.04.012
136. Rao A, Barkley D, França GS, Yanai I. Exploring tissue architecture using spatial transcriptomics. Nature. (2021) 596:211–20. doi: 10.1038/s41586-021-03634-9
137. Zormpas E, Queen R, Comber A, Cockell SJ. Mapping the transcriptome: realizing the full potential of spatial data analysis. Cell. (2023) 186:5677–89. doi: 10.1016/j.cell.2023.11.003
138. Schott M, León-Periñán D, Splendiani E, Strenger L, Licha JR, Pentimalli TM, et al. Open-st: high-resolution spatial transcriptomics in 3d. Cell. (2024) 187:3953–72.e26. doi: 10.1016/j.cell.2024.05.055
139. Zhou Z, Zhong Y, Zhang Z, Ren X. Spatial transcriptomics deconvolution at single-cell resolution using redeconve. Nat Commun. (2023) 14:7930. doi: 10.1038/s41467-023-43600-9
140. Fazal FM, Han S, Parker KR, Kaewsapsak P, Xu J, Boettiger AN, et al. Atlas of subcellular rna localization revealed by apex-seq. Cell. (2019) 178:473–90.e26. doi: 10.1016/j.cell.2019.05.027
141. Zhu J, Pang K, Hu B, He R, Wang N, Jiang Z, et al. Custom microfluidic chip design enables cost-effective three-dimensional spatiotemporal transcriptomics with a wide field of view. Nat Genet. (2024) 56:2259–70. doi: 10.1038/s41588-024-01906-4
142. Liao J, Lu X, Shao X, Zhu L, Fan X. Uncovering an organ's molecular architecture at single-cell resolution by spatially resolved transcriptomics. Trends Biotechnol. (2021) 39:43–58. doi: 10.1016/j.tibtech.2020.05.006
143. Dries R, Chen J, Del Rossi N, Khan MM, Sistig A, Yuan GC. Advances in spatial transcriptomic data analysis. Genome Res. (2021) 31:1706–18. doi: 10.1101/gr.275224.121
144. Qian X, Harris KD, Hauling T, Nicoloutsopoulos D, Muñoz-ManChado AB, Skene N, et al. Probabilistic cell typing enables fine mapping of closely related cell types in situ. Nat Methods. (2020) 17:101–6. doi: 10.1038/s41592-019-0631-4
145. Zappia L, Phipson B, Oshlack A. Exploring the single-cell rna-seq analysis landscape with the scrna-tools database. PloS Comput Biol. (2018) 14:e1006245. doi: 10.1371/journal.pcbi.1006245
146. Fang S, Chen B, Zhang Y, Sun H, Liu L, Liu S, et al. Computational approaches and challenges in spatial transcriptomics. Genomics Proteomics Bioinf. (2023) 21:24–47. doi: 10.1016/j.gpb.2022.10.001
147. Lebrigand K, Bergenstråhle J, Thrane K, Mollbrink A, Meletis K, Barbry P, et al. The spatial landscape of gene expression isoforms in tissue sections. Nucleic Acids Res. (2023) 51:e47. doi: 10.1093/nar/gkad169
148. Cao J, Li C, Cui Z, Deng S, Lei T, Liu W, et al. Spatial transcriptomics: A powerful tool in disease understanding and drug discovery. Theranostics. (2024) 14:2946–68. doi: 10.7150/thno.95908
149. Qiao C, Li D, Liu Y, Zhang S, Liu K, Liu C, et al. Rationalized deep learning super-resolution microscopy for sustained live imaging of rapid subcellular processes. Nat Biotechnol. (2023) 41:367–77. doi: 10.1038/s41587-022-01471-3
150. Danishuddin, Khan S, Kim JJ. Spatial transcriptomics data and analytical methods: an updated perspective. Drug Discovery Today. (2024) 29:103889. doi: 10.1016/j.drudis.2024.103889
151. Dannhorn A, Kazanc E, Flint L, Guo F, Carter A, Hall AR, et al. Morphological and molecular preservation through universal preparation of fresh-frozen tissue samples for multimodal imaging workflows. Nat Protoc. (2024) 19:2685–711. doi: 10.1038/s41596-024-00987-z
152. Ben-Chetrit N, Niu X, Swett AD, Sotelo J, Jiao MS, Stewart CM, et al. Integration of whole transcriptome spatial profiling with protein markers. Nat Biotechnol. (2023) 41:788–93. doi: 10.1038/s41587-022-01536-3
153. Fan Y, Andrusivová Ž, Wu Y, Chai C, Larsson L, He M, et al. Expansion spatial transcriptomics. Nat Methods. (2023) 20:1179–82. doi: 10.1038/s41592-023-01911-1
154. Moffitt JR, Lundberg E, Heyn H. The emerging landscape of spatial profiling technologies. Nat Rev Genet. (2022) 23:741–59. doi: 10.1038/s41576-022-00515-3
155. Barkley D, Moncada R, Pour M, Liberman DA, Dryg I, Werba G, et al. Cancer cell states recur across tumor types and form specific interactions with the tumor microenvironment. Nat Genet. (2022) 54:1192–201. doi: 10.1038/s41588-022-01141-9
156. Lewis SM, Asselin-Labat ML, Nguyen Q, Berthelet J, Tan X, Wimmer VC, et al. Spatial omics and multiplexed imaging to explore cancer biology. Nat Methods. (2021) 18:997–1012. doi: 10.1038/s41592-021-01203-6
157. Zhong Z, Hou J, Yao Z, Dong L, Liu F, Yue J, et al. Domain generalization enables general cancer cell annotation in single-cell and spatial transcriptomics. Nat Commun. (2024) 15:1929. doi: 10.1038/s41467-024-46413-6
158. De Zuani M, Xue H, Park JS, Dentro SC, Seferbekova Z, Tessier J, et al. Single-cell and spatial transcriptomics analysis of non-small cell lung cancer. Nat Commun. (2024) 15:4388. doi: 10.1038/s41467-024-48700-8
159. Zhang Q, Abdo R, Iosef C, Kaneko T, Cecchini M, Han VK, et al. The spatial transcriptomic landscape of non-small cell lung cancer brain metastasis. Nat Commun. (2022) 13:5983. doi: 10.1038/s41467-022-33365-y
160. Heide T, Househam J, Cresswell GD, Spiteri I, Lynn C, Mossner M, et al. The co-evolution of the genome and epigenome in colorectal cancer. Nature. (2022) 611:733–43. doi: 10.1038/s41586-022-05202-1
161. Yang S, Zhang D, Sun Q, Nie H, Zhang Y, Wang X, et al. Single-cell and spatial transcriptome profiling identifies the transcription factor bhlhe40 as a driver of emt in metastatic colorectal cancer. Cancer Res. (2024) 84:2202–17. doi: 10.1158/0008-5472.Can-23-3264
162. Wang F, Long J, Li L, Wu ZX, Da TT, Wang XQ, et al. Single-cell and spatial transcriptome analysis reveals the cellular heterogeneity of liver metastatic colorectal cancer. Sci Adv. (2023) 9:eadf5464. doi: 10.1126/sciadv.adf5464
163. Sun C, Wang A, Zhou Y, Chen P, Wang X, Huang J, et al. Spatially resolved multi-omics highlights cell-specific metabolic remodeling and interactions in gastric cancer. Nat Commun. (2023) 14:2692. doi: 10.1038/s41467-023-38360-5
164. Sundar R, Liu DH, Hutchins GG, Slaney HL, Silva AN, Oosting J, et al. Spatial profiling of gastric cancer patient-matched primary and locoregional metastases reveals principles of tumour dissemination. Gut. (2021) 70:1823–32. doi: 10.1136/gutjnl-2020-320805
165. Ravi VM, Will P, Kueckelhaus J, Sun N, Joseph K, Salié H, et al. Spatially resolved multi-omics deciphers bidirectional tumor-host interdependence in glioblastoma. Cancer Cell. (2022) 40:639–55.e13. doi: 10.1016/j.ccell.2022.05.009
166. Jain S, Eadon MT. Spatial transcriptomics in health and disease. Nat Rev Nephrol. (2024) 20:659–71. doi: 10.1038/s41581-024-00841-1
167. Park HE, Jo SH, Lee RH, Macks CP, Ku T, Park J, et al. Spatial transcriptomics: technical aspects of recent developments and their applications in neuroscience and cancer research. Adv Sci (Weinh). (2023) 10:e2206939. doi: 10.1002/advs.202206939
168. Wang WJ, Chu LX, He LY, Zhang MJ, Dang KT, Gao C, et al. Spatial transcriptomics: recent developments and insights in respiratory research. Mil Med Res. (2023) 10:38. doi: 10.1186/s40779-023-00471-x
169. Chen WT, Lu A, Craessaerts K, Pavie B, Sala Frigerio C, Corthout N, et al. Spatial transcriptomics and in situ sequencing to study alzheimer's disease. Cell. (2020) 182:976–91.e19. doi: 10.1016/j.cell.2020.06.038
170. Farah EN, Hu RK, Kern C, Zhang Q, Lu TY, Ma Q, et al. Spatially organized cellular communities form the developing human heart. Nature. (2024) 627:854–64. doi: 10.1038/s41586-024-07171-z
171. Sun J, Singh P, Shami A, Kluza E, Pan M, Djordjevic D, et al. Spatial transcriptional mapping reveals site-specific pathways underlying human atherosclerotic plaque rupture. J Am Coll Cardiol. (2023) 81:2213–27. doi: 10.1016/j.jacc.2023.04.008
Keywords: spatial transcriptomics, tumor microenvironment, individualized precise treatment, breast cancer, heterogeneity
Citation: An J, Lu Y, Chen Y, Chen Y, Zhou Z, Chen J, Peng C, Huang R and Peng F (2024) Spatial transcriptomics in breast cancer: providing insight into tumor heterogeneity and promoting individualized therapy. Front. Immunol. 15:1499301. doi: 10.3389/fimmu.2024.1499301
Received: 20 September 2024; Accepted: 05 December 2024;
Published: 19 December 2024.
Edited by:
Bramanandam Manavathi, University of Hyderabad, IndiaReviewed by:
Yuquan Chen, Monash University, AustraliaXiaofei Liu, Affiliated Hospital of Shandong University of Traditional Chinese Medicine, China
Copyright © 2024 An, Lu, Chen, Chen, Zhou, Chen, Peng, Huang and Peng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Cheng Peng, cGVuZ2NoZW5nY2hlbmdkdUAxMjYuY29t; Ruizhen Huang, SHVhbmdSdWktemhlbjQxQGZveG1haWwuY29t; Fu Peng, cGVuZ2ZAc2N1LmVkdS5jbg==