- Open Access
Single-cell and spatial analyses of cancer cells: toward elucidating the molecular mechanisms of clonal evolution and drug resistance acquisition
Inflammation and Regeneration volume 41, Article number: 22 (2021)
Even within a single type of cancer, cells of various types exist and play interrelated roles. Each of the individual cells resides in a distinct microenvironment and behaves differently. Such heterogeneity is the most cumbersome nature of cancers, which is occasionally uncountable when effective prevention or total elimination of cancers is attempted. To understand the heterogeneous nature of each cell, the use of conventional methods for the analysis of “bulk” cells is insufficient. Although some methods are high-throughput and compressive regarding the genes being detected, the obtained data would be from the cell mass, and the average of a large number of the component cells would no longer be measured. Single-cell analysis, which has developed rapidly in recent years, is causing a drastic change. Genome, transcriptome, and epigenome analyses at single-cell resolution currently target cancer cells, cancer-associated fibroblasts, endothelial cells of vessels, and circulating and infiltrating immune cells. In fact, surprisingly diverse features of clonal evolution of cancer cells, during the development of cancer or acquisition of drug resistance, accompanied by corresponding gene expression changes in the circumstantial stromal cells, appeared in recent single-cell analyses. Based on the obtained novel insights, better optimal drug selection and new drug administration sequences were started. Even a remaining concern of the single cell analyses is being addressed. Until very recently, it was impossible to obtain positional information of cells in cancer via single-cell analysis because such information is lost during preparation of single-cell suspensions. A new method, collectively called spatial transcriptome (ST) analysis, has been developed and rapidly applied to various clinical specimens. In this review, we first outline the recent achievements of single-cell cancer analysis in analyzing the molecular basis underlying the acquisition of drug resistance, particularly focusing on the latest anti-epidermal growth factor receptor tyrosine kinase inhibitor, osimertinib. Further, we review the currently available ST analysis methods and introduce our recent attempts regarding the respective topics.
A detailed understanding of cellular diversity and constituting ecosystems of cancers is expected to serve as a foundation for cancer eradication and control. Particularly, heterogeneity of cancer cells derived from clonal evolution during anticancer drug treatment is directly related to drug resistance and thus is crucial for treatment success. The selection of the optimal drug and drug administration procedure is not always based on molecular evidence, considering the possible emergence of drug-resistant cells. To this end, conventional omics analysis approaches such as RNA-seq, chromatin immunoprecipitation sequencing (ChIP-seq), and whole-genome/exome sequencing (WGS/WES) have limited power, although they may be comprehensive regarding the covered genes and genome. These methods usually use the “bulk” for the material. Inherently, they measure the average value of a large number of component cells; therefore, they cannot detect the behavior of minor cells or the presence of heterogeneous cellular populations therein. For this purpose, an individual cell analysis, so-called single-cell analysis, is needed. Indeed, several papers have been reported. For example, Kim et al. combined single-nucleus RNA-seq, single-nucleus DNA-seq (for copy number profiling), and bulk exome sequencing .
Despite the substantial advantages of single-cell analysis, the current single-cell analytical methods have been argued against due to several drawbacks. The most significant drawback is that the positional information of various cells is lost when single cells are dissociated from the tissue. With this disadvantage, it is sometimes difficult to grasp the whole picture of the interconnected structure of the precancerous, cancer, drug-resistant, and transient cells of those stages. Spatial transcriptome (ST) analysis was recently developed to address this flaw. Here, we review recent single-cell analysis approaches for a better understanding of cancer heterogeneity. A particular focus will be placed on the molecular mechanisms underlying drug resistance. Additionally, we introduce ongoing ST analysis-based studies aimed at further exploring the spatial heterogeneity of cancer cells. Finally, perspectives on the integration of these two approaches are discussed.
Single-cell sequencing for intratumor heterogeneity of cancer cells
Single-cell analysis for understanding intratumor heterogeneity and clonal evolution of cancer cells
During cancer progression, tumor cells proliferate with an accumulation of genomic mutations. Single-cell genome sequencing is potentially the most direct method for analyzing intratumor heterogeneity and cancer evolution. However, it is practically impossible to obtain precise single-cell genome data simply because an individual human cell has only two DNA copies. In addition, substantial sequencing depth is needed to cover the entire human genome. Single-cell copy number analysis has been employed in several studies to elucidate cancer genomic diversity instead of in-depth whole genome sequencing. From relatively shallow genomic sequence data, an evolutionary model for each critical phase of cancer, including clonal proliferation, invasion, relapse, and metastasis, was successfully constructed [2,3,4,5]. It was also shown to be possible to estimate copy number aberrations and clonal structure from single-cell RNA sequencing (scRNA-seq) data, based on the assumption that largescale chromosomal aberrations would affect the transcriptome features of the genes residing therein. Several computational methods, such as inferCNV  and CopyKAT , have been developed for this purpose.
Furthermore, epigenomic and transcriptomic features change during cancer progression, providing a gene expression base for malignant phenotypes, such as dedifferentiation and epithelial-mesenchymal transition (EMT). scRNA-seq analysis is the most powerful tool for inferring the transition statuses of cancers. For this purpose, a computational method called pseudo-time analysis was employed. In glioblastoma, for example, transcriptomic heterogeneity and their transition during the differentiation of cancers have been analyzed. In 2014, Patel et al. reported scRNA-seq analysis using 430 primary glioblastomas cells. They classified cells at various stages as cells having stem-like to fully differentiated signatures . In recent studies, a much larger number of cells have been analyzed. These studies have collectively delineated a complicated lineage hierarchy and differentiation paths of glioblastoma cells and their transcriptional regulators [8,9,10].
Single-cell analysis for characterizing microenvironmental non-tumor cells
There have been a large number of single-cell studies focusing on non-tumor cells and their microenvironment. The first study of such an approach, designated as a “cancer ecosystem” study, was conducted for metastatic melanoma , in which 4645 cells, including cancer cells and non-malignant cells, were analyzed. Another study conducted by Lambrechts et al., for example, reported a single-cell transcriptome catalog of stromal cells in lung cancer patients . They analyzed approximately 100,000 cells (52,698 cells for the initial catalog and 40,250 cells for validation) and identified 52 stromal cell clusters, including immune cells, endothelial cells, and fibroblasts. They characterized each type of stromal cell and found that their distribution was associated with tumor characteristics and patient prognosis.
In recent years, tumor-infiltrating lymphocytes (TILs) have received much attention. Their omics features might be associated with patient prognosis and efficacy of immunotherapies, such as immune checkpoint inhibitors (ICIs) . Numerous studies have appeared to be related to TILs, such as CD8+ T cells, regulatory T cells, dendritic cells, and tumor-associated macrophages (TAMs). For example, a research group at Stanford University reported immune cell profiling using scRNA-seq before and after anti-PD-1 treatment for basal cell carcinoma patients . They found that pre-existing TILs were neither activated nor expanded after PD-1 blockade; instead, novel T cell clones appeared to undergo sequential reactions. GAPFREE2 consortium, including the authors of this manuscript, also conducted scRNA-seq using clinical samples from gastrointestinal cancer patients to characterize the transcriptome features of a series of cells constructing cancers, primarily focusing on immune cells  (Fig. 1a). In this study, we employed 10x Genomics Chromium scRNA-seq system, which is based on the micro-droplet system . Briefly, cells resuspended in reagent for reverse transcription, Gel beads for barcoding and oil for making droplet are loaded in the chip. After the reverse transcription, each cell has unique cell barcode, 10x BC, and each mRNA has unique molecular identifier, UMI, thus we can process pooled samples and make library for sequencing. With those 10xBC and UMI, we can identify expression level of each cell (Fig. 1b). We characterized several groups of regulatory T cells that infiltrated cancer tissues using scRNA-seq (Fig. 1c). In this study, we also performed protein level analysis using cytometry by time of flight (CyTOF), as well as immunohistochemistry (IHC) concurrently with scRNA-seq. Basically, mass cytometry, including CyTOF, is a combination of flow cytometry and elemental mass spectrometry. It enabled us to analyze cell with much more parameters in the quantitative way. In CyTOF analysis, cells are stained with markers simultaneously. Though FACS analysis use fluorescence, CyTOF use heavy metal tags . After staining, cells are processed in the machine (Fig. 1d). We compared immune cell compositions detected via scRNA-seq and CyTOF using PBMCs and biopsy/surgically dissected cancer tissues. Some unexpected cells, such as plasma cells, could only be characterized in scRNA-seq datasets because CyTOF could only measure the levels of the preset markers. In contrast, CyTOF is known to have advantage when we have specific target for analysis . It could more clearly classify some cell types, such as NK cells, due to differences in the expression of cell-type markers at the mRNA and protein levels. Their results demonstrated that scRNA-seq could be an even more powerful tool when used in combination with CyTOF by mutually complementing weak points [17, 18]. We summarize the character of scRNA-seq and CyTOF in Table 1.
Single-cell sequencing unveils molecular mechanisms of anticancer drug resistance
Anticancer drugs are used for surgically inoperable cancer tissues that have already developed metastasis, severe invasion, or recurrence. Molecular targeting drugs that target signaling molecules specifically expressed in cancer cells have been administered in clinical practice in addition to traditional cytotoxic anticancer drugs, leading to improved clinical outcomes. However, there are still not a few cases in which the cancers become drug-resistant and recur after a certain period of post-drug treatment. One of the most significant barriers to the complete elimination of cancer is the emergence of drug-resistant cells during cancer evolution with a complicated microenvironment in cancer tissues. Several important studies have been performed to elucidate the drug resistance mechanism of cancer cells using single-cell sequencing technologies and have revealed the mechanisms of resistance to anticancer drugs, such as molecular targeting drugs, conventional cytotoxic drugs, and ICIs.
Single-cell sequencing presented the hypothesis for mechanisms of drug resistance
In 2018, Ho et al. analyzed the mechanism of BRAF inhibitors in melanoma using single-cell RNA-seq analysis and K/clustering evaluation (SAKE), which can process the data analysis from QC to enrichment analysis and provide detailed gene expression profiles. They compared SAKE with known scRNA-seq data analysis tools, SC3, Seurat, SINCERA, CIDR, and RaceID, and demonstrated its accuracy and robustness. They applied SAKE to reveal the mechanism of resistance to BRAF inhibitors using the Fluidigm C1 system, 10x Genomics Chromium, and bulk RNA-seq and indicated that a rare population of resistant cells to BRAF-i exists ab initio .
Sharma et al. focused on the resistance mechanism of both phenotypically homogeneous and heterogeneous cell lineages using patient-derived primary oral squamous cell carcinoma (OSCC) cell lines. They showed that cells with high heterogeneity have an advantage in the selection of intratumor heterogeneity. In contrast, cells with high homogeneity utilize tumor evolution via epigenetic changes, and the stem cell factors SOX2 and SOX9 play a role in this drug-induced adaptation .
Kim et al. conducted single-cell DNA/RNA analysis of breast cancer specimens before and after chemotherapy to elucidate that cancer cell heterogeneity is a cause of drug tolerance. Samples from 20 triple-negative breast cancer (TNBC) patients undergoing neoadjuvant chemotherapy (NAC) showed pre-existing and drug-induced selection of NAC-resistant tumor cells . Schnepp et al. also focused on chemotherapy resistance; using two prostate cancer cell line models, a docetaxel-sensitive variant and a docetaxel-resistant variant, they showed that NP1 was a key regulator .
Jerby-Arnon et al. reported that there are multiple mechanisms of resistance to an ICI. They conducted scRNA-seq and bulk RNA-seq of 33 melanoma tumors. Based on the results, they demonstrated that their predictive biomarkers and combination therapy are effective .
Collectively, it has been revealed that both of the following mechanisms are possible and have been observed: (i) the primary mechanism to acquire drug resistance is a series of genomic mutations and their clonal expansion at the same time as the clonal evolution for further fitting the niche; (ii) the tumor-forming cancer cells have a certain diversity, originally including already drug-resistant cells; and (iii) gene expression was reprogrammed after responding to chemotherapy and changed to a newly acquired expression profile.
Revealing the mechanism of resistance to epidermal growth factor receptor tyrosine kinase inhibitors using the single-cell technology
Our study focused on the mechanism of resistance acquisition to epidermal growth factor receptor tyrosine kinase inhibitors (EGFR-TKIs) in non-small cell lung cancers (NSCLCs). Resistance to EGFR-TKIs in tumors is known to be complex. The most well-elucidated somatic mutations in the EGFR genes are detected in 10–30% of NSCLC cases. These mutations promote downstream pro-survival and anti-apoptotic signals [23, 24]. Although most patients with EGFR mutations show dramatic responses to first- and second-generation EGFR-TKIs, such as gefitinib and erlotinib, and afatinib and dacomitinib, respectively, the majority of patients subsequently acquire resistance within 2 years . Osimertinib, a third-generation EGFR-TKI, was recently developed to target cells harboring a T790M resistance mutation. However, even for osimertinib, the emergence of resistance mutations has already been reported [26, 27]. We considered that resistance acquisition, which may occur in a small population of cells, can be identified by single-cell multiomics analysis. Our group, in collaboration with Kobayashi group at National Cancer Center of Japan, attempted to identify key gene expression changes during the process of resistance acquisition, first using cell lines and then clinical specimens through scRNA-seq and single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq). ATAC-seq, assay for transposase-accessible chromatin using sequencing, was first reported in 2013. Using ATAC-seq, we can know regulatory region by inserting sequence adaptors into accessible region of the genomes. One of the advances of ATAC-seq is it requires only low sample input. Single-cell ATAC-seq was reported in 2015 [28, 29]. First, microfluidic based system was used for scATAC-seq, thus it can process only limited number of cells. However, currently, we can use 10x Genomics for scATAC-seq and analyze much more number of cells . In scATAC-seq using 10x Genomics, cells are processed to isolate nuclei. The condition for nuclei isolation depends on the sample conditions. We show examples of processing of PBMC and cell line (Fig. 2a, left bottom). After nuclei isolations, samples are processed with transposase Tn5. Next, samples are encapsulated with beads with barcode using Chromium. We make libraries for sequencing. Detailed information for ATAC-seq is described in Buenrosto et al. . We also summarized the character of scATAC-seq in Table 1 . Schematic representation of the study design is shown in Fig. 2b–e. Collectively, we believe molecular features at multiple omics layers at single cell resolution would provide new insights into the mechanisms of EGFR-TKIs resistance (Kashima et al. submitted).
Overcoming the limitation of recent single-cell technologies: single-cell multiome analysis
Until very recently, for the scRNA-seq and scATAC-seq analyses, the sample was prepared independently due to technical limitations. Therefore, it was impossible to associate epigenomic and transcriptomic features for each single cell directly. Recently, a novel method called GEX-ATAC sequencing has been developed. In this method, it is possible to conduct both scRNA-seq and scATAC-seq simultaneously for a given single cell. The representative commercial kit was released from 10x Genomics. Although most of the workflow for library preparation is similar to that of previous scRNA-seq or scATAC seq, GEX-ATAC employs several unique processes. First, GEX-ATAC utilizes nuclei as a material (Fig. 3a). Therefore, the transcriptome data represented by GEX-ATAC is from nuclear mRNA; thus, it is not always consistent with cytosolic mRNA (Fig. 3b). Second, as shown in Fig. 3c, each gel bead has two types of primers: one is a primer for scATAC-seq, and the other is for scRNA-seq (Fig. 3c, left). The primer for scRNA-seq consisted of Illumina Truseq Read1, 16 nt of 10x Barcodes, working as a cell identifier, unique molecular identifier (UMI), and poly(dT)VN (Fig. 3c, middle). The other primer for scATAC-seq consisted of the Illumina P5 sequence, 16 nt of 10x Barcodes, and spacer (Fig. 3c, right). After library construction and sequencing, the sequence data were separated into scRNA-seq or scATAC-seq and each individual cell based on these primer sequences (10x Genomics User guide CG000338 RevA). A technical advantage of GEX-ATAC is that it can be started with frozen samples, which may broaden the sample selection subjected to the analysis (10x Genomics website: https://www.10xgenomics.com/products/single-cell-multiome-atac-plus-gene-expression).
Overcoming the limitation of the recent single-cell technologies: ST analysis
As described above, in a substantial number of cases for various cancer species, analyses have been performed to characterize the diverse behavior of individual cells in a given cancer. However, it has been argued that a novel approach is needed to quantify single cells that retain pathological spatial information. The field of ST analysis has emerged to address this issue. Recently, this field has expanded rapidly, and new methods have been developed one after another. These large numbers of methodologies have fundamentally different approaches, with advantages and disadvantages unique to each method and many issues to be resolved, such as dependence on sensitivity, effort, and organizational type. None of the single methods can use all of these important parameters. In this chapter, we outline the available ST analysis techniques in four groups according to the differences in their approaches: (1) microdissection-based, (2) in situ hybridization-based (ISH), (3) in situ sequencing-based (ISS), and (4) in situ capture-based (ISC). Table 1 presents a comparison of the methods. Further, we present an example of ST analysis using Visium from 10x Genomics, a commercial version of ISC-based ST analysis.
The most robust method for obtaining spatial gene expression information is to perform gene expression profiling of the targeted areas visually specified from the tissues. In laser-capture microdissection (LCM), after physically excising the target area with a laser beam under a microscope, RNA is extracted in vitro, and the gene expression information is acquired [33, 34]. In 2017, Chen et al. introduced an extended version of the LCM protocol named geographic position sequencing (Geo-Seq) that combines LCM and scRNA-seq technology . They obtained the transcriptome of a small tissue region equivalent to 10 cells while maintaining positional information. While it is possible to completely extract the different types of RNA by resectioning individual regions, such as analyzing feasible isoforms, the throughput is limited due to the time and effort of the experimental procedure. It is difficult to grasp the entire image of the tissue because the analysis area is quite limited.
A method that physically resects and acquires localized spatial gene expression data and LCM thinly slices tissue samples and prepares individual sequencing libraries from each slice. In 2014, RNA tomography (tomo-seq) was performed, in which zebrafish embryos were sliced along the body axis at a thickness of 18 μm, and RNA was extracted from each frozen section (~ 100 sections) to obtain gene expression information. They then computationally reconstructed the 3D transcriptional profile of the embryo by overlapping RNA-seq information from all frozen sections .
The transcriptome in vivo analysis (TIVA) was published in 2014 to perform ST analysis of any region using photonics . This method introduces a light-activated polyA oligo called a TIVA-tag into living cells in the tissue. This oligo has a cell transmembrane peptide that allows it to enter the cytoplasm. By irradiating the laser beam to any targeted cell, the TIVA-tag is activated inside the cell and hybridizes with the cellular mRNA. Thereafter, the TIVA-tag-mRNA hybrid was purified from the selected cells with streptavidin, and the captured mRNA was analyzed by RNA-seq. Although it is an excellent method using optical engineering to achieve single-cell resolution, it is currently limited to living cells.
In 2017, Medaglia et al. developed NICHE-seq as a method that used photoactivation to determine the cellular and molecular composition of the microenvironment and demonstrated the experiment using mice expressing photoactivatable GFP . In this method, it is possible to analyze thousands of cells in a microenvironment, but the exact location of the cells in the photoactivated region remains unknown. In 2018, Boisset et al. reported ProximID focusing on the actual physical interaction between cells to construct cellular networks, although the location of cells in the microenvironment is still unknown . This method creates a cell structural unit containing two or three interacting cells by dissociating the tissue under mild conditions so that the small interaction structure is maintained. Then, by performing scRNA-seq on cell structural units collected manually, it shows that we can find a new cell-cell interaction even if we do not know the constituent cell types in advance. However, the location information at the tissue level is unknown because of the low throughput due to manual work and the small amount of cell masse.
In medicine and molecular biology, the term “in situ” refers to “the original place in the organism where the cells belong.” Instead of the abovementioned approach of resecting individual cells (or parts) in the tissue and extracting RNA molecules, ISH-based methods are available by directly visualizing RNA where it should be. ISH uses a probe complementary to the nucleic acid and hybridizes it in a cell. This technology has existed since the 1960s and has been used to visualize gene expression since the late 1970s . The single-molecule RNA fluorescence in situ hybridization (smFISH) method has achieved high sensitivity by hybridizing multiple short fluorescent-labeled probes to different regions of the transcript . It is still of general use because of its high intracellular spatial resolution and the feasibility of quantitative measurement of transcripts. The intracellular location of RNA molecules can provide important information on biological implications as a snapshot of transcription. However, only a few genes can be targeted at once because of the limit of spectral overlapping in standard microscopy. Cai et al. reported the seqFISH method in 2014 . This method involves hybridization of the fluorescent-labeled probes (single-stranded DNA) to intracellular RNA and imaging, after which the probes are stripped with DNase I and hybridized again with different fluorescent-labeled probes. By repeating this cycle, the mRNA-type information is converted into color information, and the gene is identified by the combination pattern. Theoretically, by repeating N times with F types of fluorescent dyes, it is possible to identify the RNA types of F to the Nth power. However, in practice, errors and noise are introduced from the measurements, such as gaps in positional information, when repeating measurements. In 2015, Zhuang et al. identified 1001 genes using the multiplexed error-robust FISH method. The imaging algorithm and probe design were optimized to correct the errors and restore the true signal from the data, including errors . In 2019, Cai et al. reported that seqFISH successfully identified approximately 10,000 genes by adopting a similar concept (seqFISH+) . However, the throughput and adaptation to large tissues are concerns because both methods require optical resolution at the cellular level for imaging. When simultaneously detecting a larger number of different transcripts, ISH-based methods requiring a known target sequence had significant overlapping spectrum restrictions. With recent advancements in ingenious sequencing strategies, it is possible to obtain almost all transcriptome readouts.
In the ISH-based methods described above, the spatial gene expression in tissues or cells was detected by barcoding the fluorescent materials and probes to “see” the expressed genes in the cells. In contrast, ISS-based methods attempt to detect nucleotide sequences by directly “reading” them. In 2013, the first ISS approach using a padlock probe (PLP) targeting a known gene was reported [45, 46]. This method first reverse-transcribed mRNA in the tissue to synthesize cDNA. Hybridizing PLPs, micrometer-sized RCA products (RCPs) are obtained by a targeted amplification method called rolling circle amplification (RCA). The genes were then identified by decoding the RCP nucleotide sequence using the sequencing-by-ligation (SBL) method. In 2017, CARTANA was established to commercialize padlock-based ISS technology, with commercialization planned by 10x Genomics. In 2020, the HybISS method was presented as an improved version of the padlock-based ISS technology, with SBL changing to sequencing-by-hybridization (SBH) . They modified the probe design, which allows for a new barcoding system via SBH chemistry for an increased signal-to-noise ratio. In a method called barcode in situ targeted sequencing (BaristaSeq) reported in 2018, the read length was extended by cross-linking RCPs with extracellular matrices and using the sequencing-by-synthesis (SBS) method used in Illumina sequencers for base decoding to achieve longer read lengths, up to 15 bases . STARmap (spatially resolved transcript amplicon readout mapping) reported by Stanford University in 2018 is another approach using PLP [45, 46]. In this method, it is possible to bypass the reverse transcription step by using a set of primers and PLPs that hybridize to a specific RNA, called NSAIL. In the following step, the RCPs are amplified using the RCA method. Incorporated by the amine-modified bases during RCA, the RCPs were fixed in the location. After removing proteins and lipids to increase tissue permeability, sequencing was performed using the modified SBL method. The PLPs contain a 5-mer barcode that labels each gene, and the amplified DNA can be identified by the sequence of the 5-mer barcode depending on the gene. The possible analysis in 3D is the most distinctive point of this method. According to the paper, the resolution is still 100–150 μm thick, but it is expected to be further improved by increasing the 5-mer of the barcode. The methods described above are based on known sequences. In 2014, a method named fluorescent in situ RNA sequencing (FISSEQ) was reported, which enables the capture of unspecified RNAs . This method involves RNA reverse transcription using random primers. Incorporation of aminoallyl dUTPs makes it possible to cross-link aminoallyl dUTPs with BS (PEG)9 and allow RCPs to be anchored to the cell protein substrate. Finally, the SBL method determines the sequences of reading lengths of 30 bases.
Although ISH and ISS enabled the identification of molecules at the cellular level resolution, they require large systems. In contrast, ISC-based methods, an approach that conducts ex-situ sequencing after capturing RNA in situ, may facilitate gene analysis while maintaining the spatial organization of the tissue relatively easily.
The spatial transcriptomics technology was developed in 2016 at the Karolinska Institute, Sweden, as the first technology to adopt this approach . In this method, the oligo dT primer to trap poly A-RNA was pre-coated onto the glass slides on which the tissue sections were attached. Each primer was preloaded with a barcode to identify the location of the coordinates on the glass slide. Specifically, the slide area was subdivided into spots with a diameter of 100 μm, with different barcode-labeled capture primer sticks on each spot. Upon attaching the tissue to a slide and disrupting the cell membrane, mRNA in the tissue was captured by nearby primers. Adding reverse transcriptase to synthesize cDNA makes it possible to synthesize a gene library that reflects its localization. Sequences were determined by next-generation sequencing. Because the determined sequence contains the coordinate information, the detailed expression information of the coordinates can be obtained. The 10x Genomics company acquired, developed, and improved this ST technology to 55 μm/diameter resolution and commercialized it under the name “10x Visium” at the end of 2018. Instead of using oligos, barcoded beads named pucks were laid on glass slides by a method called slide-seq to achieve a concept similar to the ST technique by a group at the Massachusetts Institute of Technology . Because the beads have a diameter of 10 μm and a resolution higher than that of the STs described above, it is possible to acquire information at the single-cell level. However, unlike ST, it should be noted that the current version has limited sensitivity and needs the support of scRNA-seq data for proper cell-type mapping. In addition, tissue images can be obtained from adjacent sections but not from the same sections that provide RNA data. Shortly after presenting slide-seq, a method using smaller barcode beads named high-definition spatial transcriptomics (HDST) was reported by the group that developed ST . This bead with a diameter of 2 μm has a higher resolution, but the proper mapping of cell types still requires supporting scRNA-seq data. One method that can be adapted to formalin-fixed paraffin-embedded (FFPE) samples is digital spatial profiling (DSP) . In this method, FFPE samples are stained with indexing oligonucleotides covalently attached to primary antibodies or mRNA hybridization probes with a UV-photocleavable (PC) linker. The user manually selects regions of interest (ROIs) followed by excitation with UV light. UV cleaves oligo-off antibodies in the ROI. Photocleaved oligos are collected and counted digitally using nCounter or a next-generation sequencer readout. A comparison of each ST seq technology is shown in Table 2.
An example of STseq applications
Next, we profiled the spatial gene expression of breast cancer specimens via STseq using the Visium platform of 10x Genomics (Pleasanton, CA, USA), which is a commercial version of ISC-based ST analysis . Figure 4 shows the Visium results. Clustering of expression information by the K-means method showed that cancer cells were classified into clusters 1–3, confirming the spatial heterogeneity of breast cancer cells. The cancer microenvironment is classified into four clusters. Interestingly, depending on the spatial differences in cancers, the microenvironment also differed. Even among morphologically identical cancer cells, there is heterogeneity in their expression profiles, suggesting that each cancer cell forms its own microenvironment.
Cutting-edge Tools of STseq informatic analysis
As the experimental instrument technology developed, bioinformatic tools for spatial transcriptome are up coming. We introduce the paper using two different tools for STseq. Zhang et al. revealed the spatial trajectory of NSCLC samples from the data of 4 LUAD and 8 LUSC, lung squamous carcinoma using Visium dataset . In addition to reveal the heterogeneity using Seurat , they conducted analysis using spatial-pseudo-time and revealed the spatial evolutional trajectory. First, to analyze transient gene expression along spatial trajectory, they used “SPATA,” a framework for spatial transcriptome analysis, . It can analyze DEGs in spatial context. With SPATA analysis, they showed the invasion direction and continuous features (Fig. 5). They also used another tool called “stLEARN” to analyzed pseudo-time trajectory. With stLEARN, we can reveal the connections within subclusters (Fig. 5) . From the result, they showed possible spatial evolution in NSCLC. Though the tool for STseq is still limited, we believe we can get much detailed biological mechanism using those technologies.
We first reviewed a recent approach to better understand the molecular mechanisms underlying drug resistance. We introduced our recent single-cell analysis of osimertinib-resistant lung cancer cells. In the latter section, we reviewed the currently available ST analysis methods.
Single-cell analysis technologies allow for a bottom-up view of the organs and tissues from the perspective of cells. As a result, we can define various aspects, such as development, cancer heterogeneity, and infectious diseases, with a new and entirely different view. scRNA-seq revealed intratumor and intertumor heterogeneity. Recently, an increasing number of clinical samples, including anticancer therapy-resistant tumors, have been analyzed in combination with other single-cell level omics technologies. These single-cell multiomics analyses are expected to provide new findings.
Conversely, once the cellular architecture is destroyed during single-cell suspension preparations, the intuitive information viewed top-down is lost. STseq is a powerful method that can provide gene expression information without compromising the positional relationship between cells. Using 10x Genomics’ Visium, which is already commercially available, it is now possible to measure gene expression comprehensively and quantitatively with relative ease. However, the resolution of this method is insufficient.
The spatial concept within a cell is important to delve into the microenvironment within individual cells. It is possible to estimate the critical information related to the actual biological implications from the intracellular localization of the RNA molecule. Although ISS using PLPs, such as FishSeq and HybISS, may identify the genes at the subcellular level, the number of genes to be identified is limited due to the large size of RCPs and the short sequences to be decoded. In the future, it will be necessary to use both technologies in combination for analysis. We expect rapid progression of single-cell spatial omics technology that can capture time and space.
Availability of data and materials
Chromatin immunoprecipitation sequencing
Single-cell RNA sequencing
Single-cell RNA-seq analysis and K/clustering evaluation
Single-cell assay for transposase-accessible chromatin using sequencing
Non-small cell lung cancers
Epidermal growth factor receptor tyrosine kinase inhibitor
Unique molecular identifier
Cytometry by time of flight
Kim C, Gao R, Sei E, Brandt R, Hartman J, Hatschek T, et al. Chemoresistance evolution in triple-negative breast cancer delineated by single-cell sequencing. Cell. 2018;173:879–893.e13 https://doi.org/10.1016/j.cell.2018.03.041.
Wang Y, Waters J, Leung ML, Unruh A, Roh W, Shi X, et al. Clonal evolution in breast cancer revealed by single nucleus genome sequencing. Nature. 2014;512(7513):155–60 https://doi.org/10.1038/nature13600.
Gao R, Davis A, McDonald TO, Sei E, Shi X, Wang Y, et al. Punctuated copy number evolution and clonal stasis in triple-negative breast cancer. Nat Genet. 2016;48(10):1119–30 https://doi.org/10.1038/ng.3641.
Leung ML, Davis A, Gao R, Casasent A, Wang Y, Sei E, et al. Single-cell DNA sequencing reveals a late-dissemination model in metastatic colorectal cancer. Genome Res. 2017;27(8):1287–99 https://doi.org/10.1101/gr.209973.116.
Casasent AK, Schalck A, Gao R, Sei E, Long A, Pangburn W, et al. Multiclonal invasion in breast tumors identified by topographic single cell sequencing. Cell. 2018;172:205–217.e12 https://doi.org/10.1016/j.cell.2017.12.007.
Tickle T, Tirosh I, Georgescu C, Brown M, Haas B. inferCNV of the Trinity CTAT Project. Cambridge: Klarman Cell Observatory, Broad Institute of MIT and Harvard; 2019. https://github.com/broadinstitute/inferCNV
Gao R, Bai S, Henderson YC, Lin Y, Schalck A, Yan Y, et al. Delineating copy number and clonal substructure in human tumors from single-cell transcriptomes. Nat Biotechnol. 2021;39(5):599–608 https://doi.org/10.1038/s41587-020-00795-2.
Couturier CP, Ayyadhury S, Le PU, Nadaf J, Monlong J, Riva G, et al. Single-cell RNA-seq reveals that glioblastoma recapitulates a normal neurodevelopmental hierarchy. Nat Commun. 2020;11(1):3406 https://doi.org/10.1038/s41467-020-17186-5.
Neftel C, Laffy J, Filbin MG, Hara T, Shore ME, Rahme GJ, et al. An integrative model of cellular states, plasticity, and genetics for glioblastoma. Cell. 2019;178:835–849.e21 https://doi.org/10.1016/j.cell.2019.06.024.
Castellan M, Guarnieri A, Fujimura A, Zanconato F, Battilana G, Panciera T, et al. Single-cell analyses reveal YAP/TAZ as regulators of stemness and cell plasticity in glioblastoma. Nat Cancer. 2020;2(2):174–88 https://doi.org/10.1038/s43018-020-00150-z.
Tirosh I, Izar B, Prakadan SM, Wadsworth MH, Treacy D, Trombetta JJ, et al. Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq. Science. 2016;352(6282):189–96 https://doi.org/10.1126/science.aad0501.
Lambrechts D, Wauters E, Boeckx B, Aibar S, Nittner D, Burton O, et al. Phenotype molding of stromal cells in the lung tumor microenvironment. Nat Med. 2018;24(8):1277–89 https://doi.org/10.1038/s41591-018-0096-5.
Yofe I, Dahan R, Amit I. Single-cell genomic approaches for developing the next generation of immunotherapies. Nat Med. 2020;26(2):171–7 https://doi.org/10.1038/s41591-019-0736-4.
Yost KE, Satpathy AT, Wells DK, Qi Y, Wang C, Kageyama R, et al. Clonal replacement of tumor-specific T cells following PD-1 blockade. Nat Med. 2019;25(8):1251–9 https://doi.org/10.1038/s41591-019-0522-3.
Kashima Y, Togashi Y, Fukuoka S, Kamada T, Irie T, Suzuki A, et al. Potentiality of multiple modalities for single-cell analyses to evaluate the tumor microenvironment in clinical specimens. Sci Rep. 2021;11(1):341 https://doi.org/10.1038/s41598-020-79385-w.
Zheng GXY, Terry JM, Belgrader P, Ryvkin P, Bent ZW, Wilson R, et al. Massively parallel digital transcriptional profiling of single cells. Nat Commun. 2017;8(1):14049 https://doi.org/10.1038/ncomms14049.
Spitzer MH, Nolan GP. Mass cytometry: single cells, many features. Cell. 2016;165(4):780–91 https://doi.org/10.1016/j.cell.2016.04.019.
Bjornson ZB, Nolan GP, Fantl WJ. Single-cell mass cytometry for analysis of immune system functional states. Curr Opin Immunol. 2013;25(4):484–94 https://doi.org/10.1016/j.coi.2013.07.004.
Ho Y-J, Anaparthy N, Molik D, Mathew G, Aicher T, Patel A, et al. Single-cell RNA-seq analysis identifies markers of resistance to targeted BRAF inhibitors in melanoma cell populations. Genome Res. 2018;28(9):1353–63 https://doi.org/10.1101/gr.234062.117.
Sharma A, Cao EY, Kumar V, Zhang X, Leong HS, Wong AML, et al. Longitudinal single-cell RNA sequencing of patient-derived primary cells reveals drug-induced infidelity in stem cell hierarchy. Nat Commun. 2018;9(1):4931 https://doi.org/10.1038/s41467-018-07261-3.
Schnepp PM, Shelley G, Dai J, Wakim N, Jiang H, Mizokami A, et al. Single-cell transcriptomics analysis identifies nuclear protein 1 as a regulator of docetaxel resistance in prostate cancer cells. Mol Cancer Res. 2020;18(9):1290–301 https://doi.org/10.1158/1541-7786.MCR-20-0051.
Jerby-Arnon L, Shah P, Cuoco MS, Rodman C, Su M-J, Melms JC, et al. A cancer cell program promotes T cell exclusion and resistance to checkpoint blockade. Cell. 2018;175:984–997.e24 https://doi.org/10.1016/j.cell.2018.09.006.
Citri A, Yarden Y. EGF-ERBB signalling: towards the systems level. Nat Rev Mol Cell Biol. 2006;7(7):505–16 https://doi.org/10.1038/nrm1962.
Nguyen K-SH, Kobayashi S, Costa DB. Acquired resistance to epidermal growth factor receptor tyrosine kinase inhibitors in non-small-cell lung cancers dependent on the epidermal growth factor receptor pathway. Clin Lung Cancer. 2009;10(4):281–9 https://doi.org/10.3816/CLC.2009.n.039.
Gazdar AF. Activating and resistance mutations of EGFR in non-small-cell lung cancer: role in clinical response to EGFR tyrosine kinase inhibitors. Oncogene. 2009;28(S1):S24–31 https://doi.org/10.1038/onc.2009.198.
Thress KS, Paweletz CP, Felip E, Cho BC, Stetson D, Dougherty B, et al. Acquired EGFR C797S mutation mediates resistance to AZD9291 in non-small cell lung cancer harboring EGFR T790M. Nat Med. 2015;21(6):560–2 https://doi.org/10.1038/nm.3854.
Yang Z, Yang N, Ou Q, Xiang Y, Jiang T, Wu X, et al. Investigating novel resistance mechanisms to third-generation EGFR tyrosine kinase inhibitor osimertinib in non-small cell lung cancer patients. Clin Cancer Res. 2018;24(13):3097–107 https://doi.org/10.1158/1078-0432.CCR-17-2310.
Cusanovich DA, Daza R, Adey A, Pliner HA, Christiansen L, Gunderson KL, et al. Multiplex single-cell profiling of chromatin accessibility by combinatorial cellular indexing. Science. 2015;348:910–4 https://doi.org/10.1126/science.aab1601.
Buenrostro JD, Wu B, Litzenburger UM, Ruff D, Gonzales ML, Snyder MP, et al. Single-cell chromatin accessibility reveals principles of regulatory variation. Nature. 2015;523(7561):486–90 https://doi.org/10.1038/nature14590.
Satpathy AT, Granja JM, Yost KE, Qi Y, Meschi F, McDermott GP, et al. Massively parallel single-cell chromatin landscapes of human immune cell development and intratumoral T cell exhaustion. Nat Biotechnol. 2019;37(8):925–36 https://doi.org/10.1038/s41587-019-0206-z.
Buenrostro JD, Giresi PG, Zaba LC, Chang HY, Greenleaf WJ. Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position. Nat Methods. 2013;10(12):1213–8 https://doi.org/10.1038/nmeth.2688.
Pott S, Lieb JD. Single-cell ATAC-seq : strength in numbers. Genome Biol. 2015:1–4 https://doi.org/10.1186/s13059-015-0737-7.
Bonner RF, Emmert-Buck M, Cole K, Pohida T, Chuaqui R, Goldstein S, et al. Laser capture microdissection: molecular analysis of tissue. Science. 1997;278:1481,1483 https://doi.org/10.1126/science.278.5342.1481.
Emmert-Buck MR, Bonner RF, Smith PD, Chuaqui RF, Zhuang Z, Goldstein SR, et al. Laser capture microdissection. Science. 1996;274(5289):998–1001 https://doi.org/10.1126/science.274.5289.998.
Chen J, Suo S, Tam PP, Han J-DJ, Peng G, Jing N. Spatial transcriptomic analysis of cryosectioned tissue samples with Geo-seq. Nat Protoc. 2017;12(3):566–80 https://doi.org/10.1038/nprot.2017.003.
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(3):662–75 https://doi.org/10.1016/j.cell.2014.09.038.
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(2):190–6 https://doi.org/10.1038/nmeth.2804.
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 https://doi.org/10.1126/science.aao4277.
Boisset J-C, Vivié J, Grün D, Muraro MJ, Lyubimova A, van Oudenaarden A. Mapping the physical network of cellular interactions. Nat Methods. 2018;15(7):547–53 https://doi.org/10.1038/s41592-018-0009-z.
Capco DG, Jeffery WR. Differential distribution of poly(A)-containing RNA in the embryonic cells of oncopeltus fasciatus. Dev Biol. 1978;67(1):137–51 https://doi.org/10.1016/0012-1606(78)90305-6.
Raj A, van den Bogaard P, Rifkin SA, van Oudenaarden A, Tyagi S. Imaging individual mRNA molecules using multiple singly labeled probes. Nat Methods. 2008;5(10):877–9 https://doi.org/10.1038/nmeth.1253.
Lubeck E, Coskun AF, Zhiyentayev T, Ahmad M, Cai L. Single-cell in situ RNA profiling by sequential hybridization. Nat Methods. 2014;11(4):360–1 https://doi.org/10.1038/nmeth.2892.
Chen KH, Boettiger AN, Moffitt JR, Wang S, Zhuang X. Spatially resolved, highly multiplexed RNA profiling in single cells. Science. 2015;348:aaa6090 https://doi.org/10.1126/science.aaa6090.
Eng C-HL, 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(7751):235–9 https://doi.org/10.1038/s41586-019-1049-y.
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(9):857–60 https://doi.org/10.1038/nmeth.2563.
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 https://doi.org/10.1126/science.aat5691.
Gyllborg D, Langseth CM, Qian X, Choi E, Salas SM, Hilscher MM, et al. Hybridization-based in situ sequencing (HybISS) for spatially resolved transcriptomics in human and mouse brain tissue. Nucleic Acids Res. 2020;48(19):e112 https://doi.org/10.1093/nar/gkaa792.
Chen X, Sun Y-C, Church GM, Lee JH, Zador AM. Efficient in situ barcode sequencing using padlock probe-based BaristaSeq. Nucleic Acids Res. 2018;46(4):e22 https://doi.org/10.1093/nar/gkx1206.
Lee JH, Daugharthy ER, Scheiman J, Kalhor R, Yang JL, Ferrante TC, et al. Highly multiplexed subcellular RNA sequencing in situ. Science. 2014;343(6177):1360–3 https://doi.org/10.1126/science.1250212.
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 https://doi.org/10.1126/science.aaf2403.
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 https://doi.org/10.1126/science.aaw1219.
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(10):987–90 https://doi.org/10.1038/s41592-019-0548-y.
Merritt CR, Ong GT, Church SE, Barker K, Danaher P, Geiss G, et al. Multiplex digital spatial profiling of proteins and RNA in fixed tissue. Nat Biotechnol. 2020;38(5):586–99 https://doi.org/10.1038/s41587-020-0472-9.
Nagasawa S, Kuze Y, Maeda I, Kojima Y, Motoyoshi A, Onishi T, et al. Genomic profiling reveals heterogeneous populations of ductal carcinoma in situ of the breast. Commun Biol. 2021;4(1):438 https://doi.org/10.1038/s42003-021-01959-9.
Zhang L, Mao S, Yao M, Chao N, Yang Y, Ni Y, Song T, Liu Z, Yang Y, Li W. Spatial transcriptome sequencing revealed spatial trajectory in the non-small cell lung carcinoma. https://doi.org/10.1101/2021.04.26.441394.
Stuart T, Butler A, Hoffman P, Hafemeister C, Papalexi E, Mauck WM, et al. Comprehensive integration of single-cell data. Cell. 2019;177:1888–1902.e21 https://doi.org/10.1016/j.cell.2019.05.031.
Kueckelhaus J, von Ehr J, Ravi VM, Will P, Joseph K, Beck J, Hofmann UG, Delev D, Schnell O, Heiland DH. Inferring spatially transient gene expression pattern from spatial transcriptomic studies. https://doi.org/10.1101/2020.10.20.346544.
Pham D, Tan X, Xu J, LF Grice, Lam PY, Raghubar A, Vukovic J, Ruitenberg MJ, Nguyendoi Q. stLearn: integrating spatial location, tissue morphology and gene expression to find cell types, cell-cell interactions and spatial trajectories within undissociated tissues. https://doi.org/10.1101/2020.05.31.125658.
We thank and S.S. Kobayashi and A. Ohashi for their suggestions and advice for manuscript preparation.
This work was supported by JSPS KAKENHI Grant Number 19K16879 (to Y.K.).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Nagasawa, S., Kashima, Y., Suzuki, A. et al. Single-cell and spatial analyses of cancer cells: toward elucidating the molecular mechanisms of clonal evolution and drug resistance acquisition. Inflamm Regener 41, 22 (2021). https://doi.org/10.1186/s41232-021-00170-x
- Single-cell RNA-seq
- Anticancer drug resistance
- Single-cell multiome analysis
- Spatial transcriptome analysis