Research and application of single-cell sequencing in tumor heterogeneity and drug resistance of circulating tumor cells
Biomarker Research volume 8, Article number: 60 (2020)
Malignant tumor is a largely harmful disease worldwide. The cure rate of malignant tumors increases with the continuous discovery of anti-tumor drugs and the optimisation of chemotherapy options. However, drug resistance of tumor cells remains a massive obstacle in the treatment of anti-tumor drugs. The heterogeneity of malignant tumors makes studying it further difficult for us. In recent years, using single-cell sequencing technology to study and analyse circulating tumor cells can avoid the interference of tumor heterogeneity and provide a new perspective for us to understand tumor drug resistance.
Malignant tumor is a common disease whose incidence increases yearly. It has become a considerable threat to human health. In recent years, the five-year survival rate of patients has substantially improved with the widespread use of anti-tumor drugs. However, the vast majority of tumors have drug resistance, which is a large barrier to treatment . Heterogeneity of tumor cells is considered an important cause of drug resistance. In the past, we studied drug resistance in tumor cells through high-throughput sequencing based on numerous mixed cell samples, which ignored the heterogeneity of tumor cells and resulted in the dilution of the genetic characteristics of low-abundance but functionally essential cells such as circulating tumor cells (CTC). In recent years, the research methods for CTC have become more diversified. Therefore, single-cell sequencing of CTC can analyse the information of a single cell genome, transcriptome and epigenetic group, which reduces the interference of tumor heterogeneity  and provides a new perspective for understanding the drug resistance of tumors.
Relationship between tumor cell heterogeneity and drug resistance
Scientists currently believe that two mechanisms lead to drug resistance in tumors: inherent drug resistance and acquired drug resistance (Fig. 1). Inherent drug resistance is the existing drug resistance before the use of anti-tumor drugs. Acquired drug resistance occurs during or after treatment. Inherent drug resistance may arise from rare pre-existing subclones, whereas acquired drug resistance is an acquired new mutation . After multiple divisions and proliferation of tumor cells, their progeny cells show inconsistencies in genomic and biological characteristics, and this inconsistency makes several biological characteristics of tumor cells different, which is called tumor heterogeneity. Tumor heterogeneity can be divided into inter-tumor heterogeneity and intra-tumor heterogeneity. The current research focus is on heterogeneity within the tumor. Heterogeneity in the tumor includes spatial heterogeneity and temporal heterogeneity (Fig. 2). In tumors, different cellular clones at different spatial sites lead to spatial heterogeneity. Tumor cells change with time, which is the temporal heterogeneity of tumor cells . tumor cells also affect the stroma, immune cells and other cells, which constitutes the heterogeneity of the tumor microenvironment (TME) . Numerous studies have shown that tumor heterogeneity is an important cause of drug resistance in tumor cells . For example, cells with strong drug resistance will gradually replace cells sensitive to drugs with the progress of chemotherapy . Thus, we need to have a deeper understanding of tumor heterogeneity.
Value of single-cell sequencing in the study of tumor cell heterogeneity
Important information such as mutation status, epigenetic status and related protein expression levels of tumor cells may be expressed only in few cells or even in a single cell . Heterogeneity is ignored if mixed tumor cells are used for analysis. Studying the drug resistance of tumor cells at the single-cell level is important. Single-cell sequencing technology refers to a technique for sequencing the genome and transcriptome at the single-cell level. Compared with previous sequencing methods, it can perform a more thorough analysis of healthy cells and tumor cells . It can also identify previously unknown cell types [10, 11]. Thus, it can better reveal the heterogeneity of tumor cells at the cellular and molecular levels. Using single-cell sequencing technology to study the heterogeneity of tumor cells has been widely practiced in malignant tumors such as breast cancer, melanoma and lung cancer [12,13,14,15,16].
Single-cell sequencing of CTC and drug resistance
As mentioned earlier, the heterogeneity of tumor cells, especially the transcriptome information, including time and space limitation, is likely to change constantly. Hence, studying the heterogeneity of tumor cells dynamically can better explore the problem of tumor drug resistance . The main approaches of obtaining tissue specimens are tissue biopsy and cell puncture , which is an invasive procedure with the risk of tumor spread, especially in patients with advanced cancer and multiple metastases . Moreover, researchers may not be able to acquire sufficient experimental standard quality tissue specimens for various reasons. CTC, a type of tumor cell that is separated from the primary focus or metastasis of solid tumors and enters the peripheral blood circulation, has gradually come into people’s field of vision. Studies have confirmed that CTC has characteristics similar to those of tissue cells at the single-cell sequencing level . Extracting tumor cells from peripheral blood has the advantage of being less invasive. In the state of peripheral blood circulation, various types of tumors are relatively uniform, which can more accurately reflect the temporal and spatial heterogeneity of tumor cells . Therefore, increasing studies are conducted on tumor resistance through CTC. Scholars study choriocarcinoma , colorectal cancer [23,24,25,26], liver cancer , lung cancer [28, 29], breast cancer [30, 31], lymphoma , gastric cancer  and other malignant tumors through CTC to explore the drug resistance of tumor cells. However, most current CTC analyses are based on CTC epithelial biomarkers, and epithelial biomarkers may not be expressed in several tumor types . Analyses have not been studied at the genomic level, thus lacking individual CTC characteristics or genetic characteristics of patients . CTC accounts for a small proportion of blood, and only a few dozen CTCs can be obtained at a single time. The genetic material isolated from a single cell is minimal. Applying traditional sequencing methods to this kind of research is challenging. The minimal raw materials of cells hinders the analysis of tumor cells. Although the whole genome amplification (WGA) technology has been developed for decades, the previous WGA technology has a greater difficulty in achieving high-quality genome amplification . The application of single-cell sequencing to analyse CTC can study the information at the gene level in more detail to explore the relevant mechanisms of drug resistance in tumor cells further, accurately predict the effect of anti-tumor drugs and formulate more effective drug regimens .
Single-cell sequencing of CTC to explore the mechanism of drug resistance in tumor cells
In recent years, scientists studied CTC through single-cell sequencing technology and discovered new mechanisms of drug resistance of tumor cells. Miyamoto et al.  isolated 133 CTC from blood samples of patients with prostate cancer and performed single-cell RNA sequencing on CTC, original tumor samples and samples of prostate cancer cell lines. Previous studies showed that the androgen receptor (AR) pathway is a first-line therapeutic target for prostate cancer. However, Miyamoto found that AR point mutations associated with altered signalling are not common in patients with castration-resistant prostate cancer (CRPC). tumor heterogeneity exists between patients such that different CTCs have different or more mRNA splicing changes. Compared with metastatic prostate cancer without enzalutamide (an AR inhibitor), non-classical Wnt signalling is considerably enriched in prostate CTC treated with enzalutamide during the progression of radiological imaging or prostate-specific antigen. CTCs with low glucocorticoid receptor (GR) expression levels are also enriched in non-classical Wnt signals in patients with enzalutamide progression but not in CTCs with high GR levels. Studies have shown that non-classical Wnt signalling may be a mechanism of CRPC resistance.
Anaplastic lymphoma kinase (ALK) gene, first discovered in a subtype of anaplastic large cell lymphoma, is responsible for encoding a receptor tyrosine kinase (RTK) called ALK. The common pathogenic mutation in the ALK gene is gene rearrangement, which can lead to tumors. ALK gene fusion often occurs in a subset of patients with non-small-cell lung cancer (NSCLC); thus, ALK inhibitors are useful in the treatment of such patients. However, such patients will inevitably become resistant to ALK inhibitors, leading to treatment failure. Pailler et al.  studied drug resistance mutations in CTCs isolated from 17 patients with ALK-rearranged NSCLC (14 were resistant to crizotinib, and 3 were resistant to loratinib) and studied more than 48 cancer-related genes and 14 ALK mutation regions at the single-cell level. Genomic heterogeneity was found in CTC from patients with crizotinib-resistant ALK rearrangement. In the ALK-independent pathway, the kirsten rat sarcoma viral oncogene and TP53 pathways play a significant role. In a patient resistant to loratinib, Pailler identified two ALK multiple mutations, which are a ‘targeted’ resistance mechanisms. These mutations are likely the mechanism of resistance of patients with ALK-rearrangement and NSCLC to ALK inhibitors.
According to Miyamoto’s sequencing results, Schissler  studied the mechanism of prostate cancer resistance from another perspective and argued that Miyamoto’s single-cell sequencing results lack a functional explanation, and the approach to sequencing does not fully apply to small populations of cells such as CTC. Schissler et al. designed an experimental model to simulate transcriptome dynamics by analysing the cell–cell statistical distance of aggregation within a biomolecular pathway to identify differentially expressed pathways associated with drug resistance. They designed a new aggregation method, cell-centric statistics method, and proved its effectiveness in predicting a single CTC’s drug resistance. They found that five types of pathways are substantially overexpressed in the CTC of patients with enzalutamide resistance. Syndecan-4-mediated signal transduction pathways play an important role in the resistance mechanism.
Most breast cancers express estrogen receptor (ER); thus, targeted therapy has become the first choice for the treatment of ER-positive breast cancer . Selective ER modulators or decomposers can be used to target the ER pathway, and aromatase inhibitors (AI) and other drugs that cause estrogen deficiency can be used for treatment. Studies have shown that estrogen receptor 1 (ESR1) mutations are known resistance mechanisms of tumor cells to AI or gonadotropin-releasing hormone analogues. Franken et al.  performed CTC single-cell sequencing on 46 patients with metastatic breast cancer and found that mutations in ESR1 (including known hot spot mutations and new mutations) are only present in patients receiving estrogen deprivation therapy (EDT). By contrast, no mutations were detected in patients not receiving EDT or other types of therapy. The researchers believe that the newly discovered mutations could lead to resistance to targeted drugs because the newly discovered mutations are found only in patients receiving targeted therapy and affect highly conserved amino acids. Hong et al.  believed in an interaction between genetic and non-genetic factors in the development of resistance to endocrine therapy in breast cancer. They used single-cell RNA sequencing to analyse tumor cells and defined a subpopulation of cells that may be pre-resistant to endocrine therapy. This subpopulation is substantially expressed in CTC. Transcription reprogramming and copy number changes to achieve complete drug resistance also reveal a multi-step drug resistance mechanism that interacts with genetic and non-genetic factors.
Researchers  found that the emergence of KRAS mutant clones is the second mechanism of colon cancer resistance through the single-cell analysis of CTC in patients with colon cancer. Researchers  found another acquired epidermal growth factor receptor (EGFR) extracellular region mutation (S492R) in the CTC of patients with colorectal cancer, which can prevent cetuximab (an EGFR blocking antibody), enabling its resistance to cetuximab.
Many studies fully demonstrated that using single-cell sequencing technology to analyse CTC can better study the mechanism of tumor drug resistance.
Single-cell sequencing of CTC to predict drug resistance in tumor cells
Scientists currently believe that drug resistance of tumor cells is a complex biological behaviour with continuous evolution and dynamic change in tumor development. In this process, multiple factors affect the drug resistance of tumor cells with dynamic changes over time. The use of single-cell sequencing to monitor the abnormal signal pathway of CTC can predict drug resistance and realise individualisation and precise treatment.
Malignant tumor cells involve a series of genomic changes, including copy number variations (CNVs), single nucleotide variations (SNVs) and insertions/deletions (indels). In lung cancer, the timely detection of mutations that lead to drug resistance is essential for selecting the appropriate treatment for patients. Lung cancer with EGFR mutation is initially well treated with tyrosine kinase inhibitor , but drug resistance is inevitable with drugs such as Gefitinib and Erlotinib. This resistance can be overcome in combination with other drugs . Therefore, the EGFR expression profile of patients with NSCLC must be monitored in time to capture the dynamic time changes of the disease for treatment. Ni et al.  effectively solved the problem of insufficient DNA content in CTC by using multiple annealing circular loop amplification (MALBAC) and achieved genome-wide amplification. Thus, CNV and SNV can be accurately measured in a single cell. Mutations in the EGFR gene were found in the CTC of one patient, and PIK3CA mutations were detected in the CTC of this patient through single-cell genome sequencing of 11 patients with NSCLC. Studies have shown that PIK3CA mutation is associated with erlotinib resistance . The patient developed rapidly after one month of treatment with erlotinib. Ni believed that the single-cell sequencing of CTC could effectively predict the drug resistance of anti-tumor drugs and facilitate the further guidance of drug use. Tan et al.  also performed EGFR mutation detection on a single CTC in seven patients with NSCLC, such as L858R and T790M. Compared with the matched tumor biopsy results, the results showed the right consistency between the two. Thus, single-cell analysis can provide a more accurate disease map, which is crucial for the timely detection of mutations leading to drug resistance and the selection of appropriate treatment options for patients. Similar conclusion can also be made in single-cell sequencing of CTC in prostate cancer. A single-cell research of CTCs from blood samples of patients with localised high-risk prostate cancer showed 202,241 SNVs and 137,407 indels. The researcher explored the effect of the SNVs identified in CTCs on drug response targets and found that nine genetic variations are associated with the response to docetaxel, and 48 SNVs influence drug response for 24 known cancer drugs. A gene set enrichment analysis was carried out to study the correlation between CTC-shared copy number alterations (CNAs) and pathways. The results showed some gene amplification affecting DNA damage repair pathways, which could result in chemotherapy resistance .
For small-cell lung cancer, no convincing treatment plan other than chemotherapy currently exists. However, its drug resistance also leads to the failure of many treatments, and several patients even develop drug resistance after the first use of chemotherapy drugs. Clinically, patients with small-cell lung cancer can be divided into a chemotherapy-resistant group and a chemotherapy-sensitive group according to the response within 3 months after initial chemotherapy. Tumor cell adhesion is poor in small-cell lung cancer. The expression characteristics of CTC may be higher than those of other cancer types. Su et al.  performed whole genome sequencing detection of CTC in 10 patients with small-cell lung cancer who received standard chemotherapy regimens of etoposide + platinum (cisplatin, neda platinum or carboplatin). Amongst the 91 CTCs selected, the CNA spectra of patients’ CTC were found highly consistent, which also confirmed that CTC could be used as an ideal sample for clinical correlation analysis of patients with cancer. Based on the sequencing of a single CTC, the researchers established a CNA score to analyse the correlation between the CNA data and the efficacy and survival of first-line chemotherapy. Su et al. used the CNA score as a standard to predict the clinical subtypes of patients, correctly predicting 20 out of 25 patients with chemotherapy resistance and 15 out of 16 patients with chemotherapy sensitivity. Similarly, Carter et al.  detected CNA on 88 CTC cells isolated from 13 patients with small-cell lung cancer and generated a CNA-based classifier, which was then tested on another 18 patients. After verifying 112 CTC cells from 18 patients, the accuracy was found to be 83.3%. Studies have sufficiently demonstrated that the CTC single-cell sequencing technology can be used to determine whether patients are sensitive or resistant to chemotherapy drugs, providing theoretical support for the realisation of accurate treatment.
Single-cell genome sequencing can provide genetic heterogeneity and cell lineage information, whereas single-cell transcriptome sequencing can more dynamically represent the sum of all RNAs produced by a species or a specific cell under a functional state, which can better define the current table of cells type. Cheng et al.  used hydro-seq technology to conduct transcriptome sequencing of CTC and selected 666 CTCs from 21 patients with breast cancer. The drug targets of breast cancer, including ER, AR and human epidermal growth factor receptor 2 were successfully detected from CTC, and the full transcriptome analysis of CTCs was effectively realised to study the heterogeneity of tumor cells further. This finding also proves that CTC single-cell transcriptome sequencing is an effective method to analyse the molecular characteristics of tumor cells and can be used for treatment selection and patient monitoring.
Single-cell sequencing of CTC to investigate the association with drug resistance and TME
As mentioned previously, tumor heterogeneity is an important theoretical basis for drug resistance. Tumor heterogeneity is not only affected by tumor cells themselves but also closely related to infiltrating immune cells , endothelial cells forming blood vessels , cancer-associated fibroblast (CAF)  and extracellular matrix . These internal and external environments that are closely related to the development of tumor and tumor cells together constitute the TME. Therefore, the study of TME has become an important means to explore the drug resistance of tumors. Considerable research supported the viewpoint that TME could influence the therapy response by different mechanisms. A recent study reported that TME fibroblasts can be differentiated into CAFs by tumor-derived exosomes, which promote drug-resistant phenotypes . Extracellular vesicle-mediated intercellular communication has also been proven to play a crucial role in the response to anti-tumor treatments . In addition, tumor-associated macrophages and other stromal cells from TME could affect tumor resistance . In recent years, studies showed that the detection of CTC could reflect its relationship with TME and the heterogeneity of TME [62, 63]. This technique applied in CTC showed the potential to explore the association with drug resistance and TME with the development of single-cell sequencing. Heather et al.  applied single cell RNA sequencing to identify two populations of CTCs amongst 1707 cells from 14 patients with breast cancer. According to transcriptome analysis, CTCs could be divided into two subgroups: one is rich in transcripts that are characterised by estrogen responsiveness and increased proliferation, and the other is rich in transcripts characterised by reduced proliferation and epithelial–mesenchymal transition (EMT). A cell–cell communication tool was used to investigate the association between CTCs and TME. The results showed that CTCs with increased EMT markers have a more diverse pattern of interaction. The gene expression pathway activation analysis showed that a CTC population with EMT characteristics has a lower likelihood of activating the apoptotic pathway, which could lead to drug resistance. Heather’s study indicates that the single-cell sequencing of CTC can be used to further explore the tumor microenvironment and tumor drug resistance.
The single-cell sequencing of CTC is a burgeoning technology, but it is not yet fully mature and has several flaws. For example, fewer nucleotide feedstocks may be more likely to cause biological noise than conventional sequencing. The latest linear amplification via transposon insertion method  can address this problem. Another limitation is that very few materials are available. For example, in single-cell transcriptome sequencing of CTC, each cell contains an average of only about 10 pique total RNA, in which the mRNA is only about 0.1 pique; people think that CTC is not found in the blood of all patients with cancer, especially patients with early cancer ; hence, the single-cell sequencing of CTC is still not used on a large scale. Moreover, whether CTC can fully express all tumor information remains controversial, and studies found no apparent relationship between tumor tissue and CTC . To what extent CTC represents the corresponding tumor remains controversial, such that people  proposed combining single-cell sequencing with traditional sequencing technology and further analysing the results to reduce the error. The single-cell sequencing of CTC will inevitably lead to a partial loss of spatial information of tumor cells. The spatial location of interactions between various tumor cells and cells in the TME is also critical to the development of tumors. Single-cell sequencing is necessary for our future research on tumor resistance to avoid the spatial information loss caused by CTC.
In addition, the bioinformatics and computational methods of single-cell sequencing remain unable to match the corresponding data fully, and the new single-cell sequencing data make the existing analysis tools appear impractical. At present, single-cell sequencing technology has also produced unprecedented data types. For example, two new sequencing algorithms [69, 70], which consider the spatial analysis of variance of malignant tumors, have been proposed recently. Therefore, the expanding data are also a considerable challenge for single-cell analysis.
Although the single-cell sequencing of CTC remains immature, this cross-era technology has its natural advantages , and we believe that several problems will be gradually solved with the continuous optimisation of genome-wide, transcriptome amplification methods and the rapid development of bioinformatics.
The single-cell sequencing of CTC contributes to the study of genetic heterogeneity and drug resistance of tumor cells. The single-cell sequencing of CTC certainly provides more powerful tools to understand many of the remaining mysteries of cancer cell resistance.
Availability of data and materials
Anaplastic lymphoma kinase
Anaplastic large cell lymphoma
Analysis of variance
Castration-resistant prostate cancer
Circulating tumor cell
Copy number alteration
Copy number variations
Estrogen Deprivation Therapy
Epidermal growth factor receptor
Estrogen receptor 1
Gonadotropin-releasing hormone analogs
Human epidermal growth factor receptor 2
Kirsten rat sarcoma viral oncogene
Linear amplification via transposon insertion
Multiple annealing circular loop amplification
Non-small cell lung cancer
Prostate specific antigen
Receptor tyrosine kinase
Single nucleotide variations
- TKI :
Tyrosine kinase inhibitor
Whole genome sequencing
Whole genome amplification
Ren X, Kang B, Zhang Z. Understanding tumor ecosystems by single-cell sequencing: promises and limitations. Genome Biol. 2018;19(1):211.
Huang L, Ma F, Chapman A, et al. Single-cell whole genome amplification and sequencing: methodology and applications. Annu Rev Genomics Hum Genet. 2015;16:79–102.
Yardley DA. Drug resistance and the role of combination chemotherapy in improving patient outcomes. Int J Breast Cancer. 2013;2013:137414.
Valkenburg KC, de Groot AE, Pienta KJ. Targeting the tumor stroma to improve cancer therapy. Nat Rev Clin Oncol. 2018;15:366–81.
Maley CC, Aktipis A, Graham TA, et al. Classifying the evolutionary and ecological features of neoplasms. Nat Rev Cancer. 2017;17(10):605–19.
Lee JY, Yoon JK, Kim B, et al. Tumor evolution and intratumor heterogeneity of an epithelial ovarian cancer investigated using next-generation sequencing. BMC Cancer. 2015;15:85.
Lynch M, et al. Genetic drift, selection and the evolution of the mutation rate. Nat Rev Genet. 2016;17:704–14.
Gerlinger M, Horswell S, Larkin J, et al. Genomic architecture and evolution of clear cell renal cell carcinomas defined by multiregion sequencing. Nat Genet. 2014;46(3):225–33.
Van den Bos H, Bakker B, Spierings DCJ, Lansdorp PM, Foijer F. Single-cell sequencing to quantify genomic integrity in cancer. Int J Biochem Cell Biol. 2018;94:146–50.
Tang Q, Iyer S, Lobbardi R, et al. Dissecting hematopoietic and renal cell heterogeneity in adult zebrafish at single-cell resolution using RNA sequencing. J Exp Med. 2017;214(10):2875–87.
Edil BH, Luo W, Li M. Genomic and single cell sequencing facilitate the dissection of heterogeneity of pancreatic tumors. BMC Med. 2020;18(1):177.
Kim C, Gao R, Sei E, et al. Chemoresistance evolution in triple-negative breast cancer delineated by single-cell sequencing. Cell. 2018;173(4):879–893.e13.
Tirosh I, Izar B, Prakadan SM, et al. Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq. Science. 2016;352(6282):189–96.
Huss WJ, Hu Q, Glenn ST, et al. Comparison of SureSelect and Nextera exome capture performance in single-cell sequencing. Hum Hered. 2018;83(3):153–62.
Lee MC, Lopez-Diaz FJ, Khan SY, et al. Single-cell analyses of transcriptional heterogeneity during drug tolerance transition in cancer cells by RNA sequencing. Proc Natl Acad Sci U S A. 2014;111(44):E4726–35.
Kim KT, Lee HW, Lee HO, et al. Single-cell mRNA sequencing identifies subclonal heterogeneity in anti-cancer drug responses of lung adenocarcinoma cells. Genome Biol. 2015;16(1):127.
Roukos DH. Spatiotemporal diversification of intrapatient genomic clones and early drug development concepts realize the roadmap of precision cancer medicine. Drug Discov Today. 2017;22(8):1148–64.
Schmitt MW, Loeb LA, Salk JJ. The influence of subclonal resistance mutation resistance mutations on targeted cancer therapy. Nat Rev Clin Oncol. 2016;13(6):335–47.
Janku F. Tumor heterogeneity in the clinic: is it a real problem? Ther Adv Med Oncol. 2014;6(2):43–51.
Jiang W, Wang H, Cui Y, et al. Polymer nanofiber-based microchips for EGFR mutation analysis of circulating tumor cells in lung adenocarcinoma. Int J Nanomedicine. 2018;13:1633–42.
Calabuig-Fariñas S, Jantus-Lewintre E, Herreros-Pomares A, Camps C. Circulating tumor cells versus circulating tumor DNA in lung cancer-which one will win? Transl Lung Cancer Res. 2016;5(5):466–82.
He W, Hou M, Zhang H, et al. Clinical significance of circulating tumor cells in predicting disease progression and chemotherapy resistance in patients with gestational choriocarcinoma. Int J Cancer. 2019;144(6):1421–31.
Matsusaka S, Hanna DL, Ning Y, et al. Epidermal growth factor receptor mRNA expression: a potential molecular escape mechanism from regorafenib. Cancer Sci. 2020;111(2):441–50.
Troncarelli Flores BC, Souza E, Silva V, Ali Abdallah E, et al. Molecular and kinetic analyses of circulating tumor cells as predictive markers of treatment response in locally advanced rectal cancer patients. Cells. 2019;8(7):641.
Wills B, Gorse E, Lee V. Role of liquid biopsies in colorectal cancer. Curr Probl Cancer. 2018;42(6):593–600.
Sun Q, Liu Y, Liu B, Liu Y. Use of liquid biopsy in monitoring colorectal cancer progression shows strong clinical correlation. Am J Med Sci. 2018;355(3):220–7.
Zhou J, Huang A, Yang XR. Liquid biopsy and its potential for management of hepatocellular carcinoma. J Gastrointest Cancer. 2016;47(2):157–67.
Lei L, Wang WX, Zhu YC, et al. Potential mechanism of primary resistance to icotinib in patients with advanced non-small cell lung cancer harboring uncommon mutant epidermal growth factor receptor: a multi-center study. Cancer Sci. 2020;111(2):679–86.
Lei L, Wang WX, Zhu YC, et al. Real-world efficacy and potential mechanism of resistance of icotinib in Asian advanced non-small cell lung cancer with EGFR uncommon mutations: a multi-center study. Cancer Med. 2020;9(1):12–8.
Wang Q, Guldner IH, Golomb SM, et al. Single-cell profiling guided combinatorial immunotherapy for fast-evolving CDK4/6 inhibitor-resistant HER2-positive breast cancer. Nat Commun. 2019;10(1):3817.
Jordan NV, Bardia A, Wittner BS, et al. HER2 expression identifies dynamic functional states within circulating breast cancer cells. Nature. 2016;537(7618):102–6.
Alcaide M, Rushton C, Morin RD. Ultrasensitive detection of circulating tumor DNA in lymphoma via targeted hybridization capture and deep sequencing of barcoded libraries. Methods Mol Biol. 1956;2019:383–435.
Li Y, Zhang X, Liu D, et al. Evolutionary expression of HER2 conferred by chromosome aneuploidy on circulating gastric cancer cells contributes to developing targeted and chemotherapeutic resistance. Clin Cancer Res. 2018;24(21):5261–71.
Joosse SA, Pantel K. Biologic challenges in the detection of circulating tumor cells. Cancer Res. 2013;73(1):8–11.
Wu CP, Wu P, Zhao HF, Liu WL, Li WP. Clinical applications of and challenges in single-cell analysis of circulating tumor cells. DNA Cell Biol. 2018;37(2):78–89.
Zhang X, Marjani SL, Hu Z, Weissman SM, Pan X, Wu S. Single-cell sequencing for precise cancer research: progress and prospects. Cancer Res. 2016;76(6):1305–12.
Puram SV, Tirosh I, Parikh AS, et al. Single-cell transcriptomic analysis of primary and metastatic tumor ecosystems in head and neck cancer. Cell. 2017;171(7):1611–1624.e24.
Miyamoto DT, Zheng Y, Wittner BS, et al. RNA-Seq of single prostate CTCs implicates noncanonical Wnt signaling in antiandrogen resistance. Science. 2015;349(6254):1351–6.
Pailler E, Faugeroux V, Oulhen M, et al. Acquired resistance mutations to ALK inhibitors identified by single circulating tumor cell sequencing in ALK-rearranged non-small-cell lung cancer. Clin Cancer Res. 2019;25(22):6671–82.
Schissler AG, Li Q, Chen JL, et al. Analysis of aggregated cell-cell statistical distances within pathways unveils therapeutic-resistance mechanisms in circulating tumor cells. Bioinformatics. 2016;32(12):i80–9.
Huang B, Warner M, Gustafsson JÅ. Estrogen receptors in breast carcinogenesis and endocrine therapy. Mol Cell Endocrinol. 2015;418(Pt 3):240–4.
Franken A, Honisch E, Reinhardt F, et al. Detection of ESR1 mutations in single circulating tumor cells on estrogen deprivation therapy but not in primary tumors from metastatic luminal breast cancer patients. J Mol Diagn. 2020;22(1):111–21.
Hong SP, Chan TE, Lombardo Y, et al. Single-cell transcriptomics reveals multi-step adaptations to endocrine therapy. Nat Commun. 2019;10(1):3840.
Misale S, Yaeger R, Hobor S, et al. Emergence of KRAS mutations and acquired resistance to anti-EGFR therapy in colorectal cancer. Nature. 2012;486(7404):532–6.
Montagut C, Dalmases A, Bellosillo B, et al. Identification of a mutation in the extracellular domain of the epidermal growth factor receptor conferring cetuximab resistance in colorectal cancer [published correction appears in Nat Med. 2012 Sep;18(9):1445. Somasekar, Seshagiri [corrected to Seshagiri, Somasekar]]. Nat Med. 2012;18(2):221–3.
Ettinger DS, Wood DE, Akerley W, et al. Non-small cell lung cancer, version 6.2015. J Natl Compr Cancer Netw. 2015;13(5):515–24.
Yan D, Parker RE, Wang X, et al. MERTK promotes resistance to irreversible EGFR tyrosine kinase inhibitors in non-small cell lung cancers expressing wild-type EGFR family members. Clin Cancer Res. 2018;24(24):6523–35.
Ni X, Zhuo M, Su Z, et al. Reproducible copy number variation patterns among single circulating tumor cells of lung cancer patients. Proc Natl Acad Sci U S A. 2013;110(52):21083–8.
Sequist LV, Waltman BA, Dias-Santagata D, et al. Genotypic and histological evolution of lung cancers acquiring resistance to EGFR inhibitors. Sci Transl Med. 2011;3(75):75ra26.
Tan SJ, Yeo T, Sukhatme SA, Kong SL, Lim WT, Lim CT. Personalized treatment through detection and monitoring of genetic aberrations in single circulating tumor cells. Adv Exp Med Biol. 2017;994:255–73.
Rangel-Pozzo A, Liu S, Wajnberg G, et al. Genomic analysis of localized high-risk prostate cancer circulating tumor cells at the single-cell level. Cells. 2020;9(8):1863. https://doi.org/10.3390/cells9081863 PMID: 32784507; PMCID: PMC7466090.
Su Z, Wang Z, Ni X, et al. Inferring the evolution and progression of small-cell lung cancer by single-cell sequencing of circulating tumor cells. Clin Cancer Res. 2019;25(16):5049–60.
Carter L, Rothwell DG, Mesquita B, et al. Molecular analysis of circulating tumor cells identifies distinct copy-number profiles in patients with chemosensitive and chemorefractory small-cell lung cancer. Nat Med. 2017;23(1):114–9.
Cheng YH, Chen YC, Lin E, et al. Hydro-Seq enables contamination-free high-throughput single-cell RNA-sequencing for circulating tumor cells. Nat Commun. 2019;10(1):2163.
Zhang Y, Song J, Zhao Z, et al. Single-cell transcriptome analysis reveals tumor immune microenvironment heterogenicity and granulocytes enrichment in colorectal cancer liver metastases. Cancer Lett. 2020;470:84–94. https://doi.org/10.1016/j.canlet.2019.10.016 Epub 2019 Oct 11. PMID: 31610266.
Norton KA, Jin K, Popel AS. Modeling triple-negative breast cancer heterogeneity: effects of stromal macrophages, fibroblasts and tumor vasculature. J Theor Biol. 2018;452:56–68. https://doi.org/10.1016/j.jtbi.2018.05.003 Epub 2018 May 8. PMID: 29750999; PMCID: PMC6127870.
Neuzillet C, Tijeras-Raballand A, Ragulan C, et al. Inter- and intra-tumoral heterogeneity in cancer-associated fibroblasts of human pancreatic ductal adenocarcinoma. J Pathol. 2019;248(1):51–65. https://doi.org/10.1002/path.5224 Epub 2019 Feb 22. PMID: 30575030; PMCID: PMC6492001.
Le Large TY, Mantini G, Meijer LL, et al. Microdissected pancreatic cancer proteomes reveal tumor heterogeneity and therapeutic targets. JCI Insight. 2020;5(15):e138290. https://doi.org/10.1172/jci.insight.138290 PMID: 32634123; PMCID: PMC7455080.
Ringuette Goulet C, Bernard G, Tremblay S, et al. Exosomes induce fibroblast differentiation into cancer-associated fibroblasts through TGFβ signaling. Mol Cancer Res. 2018;16(7):1196–204. https://doi.org/10.1158/1541-7786.MCR-17-0784 Epub 2018 Apr 10. PMID: 29636362.
Zheng P, Chen L, Yuan X, et al. Exosomal transfer of tumor-associated macrophage-derived miR-21 confers cisplatin resistance in gastric cancer cells. J Exp Clin Cancer Res. 2017;36(1):53.
Wang T, Nasser MI, Shen J, et al. Functions of exosomes in the triangular relationship between the tumor, inflammation, and immunity in the tumor microenvironment. J Immunol Res. 2019;2019:1–10.
Hamilton G, Rath B, Klameth L, et al. Small cell lung cancer: recruitment of macrophages by circulating tumor cells. Oncoimmunology. 2015;5:e1093277.
Wei C, Yang C, Wang S, et al. Crosstalk between cancer cells and tumor associated macrophages is required for mesenchymal circulating tumor cell-mediated colorectal cancer metastasis. Mol Cancer. 2019;18(1):64. https://doi.org/10.1186/s12943-019-0976-4 PMID: 30927925; PMCID: PMC6441214.
Brechbuhl HM, Vinod-Paul K, Gillen AE, et al. Analysis of circulating breast cancer cell heterogeneity and interactions with peripheral blood mononuclear cells. Mol Carcinog. 2020;59(10):1129–39. https://doi.org/10.1002/mc.23242 Epub 2020 Aug 21. PMID: 32822091.
Chen C, Xing D, Tan L, et al. Single-cell whole-genome analyses by linear amplification via transposon insertion (LIANTI). Science. 2017;356(6334):189–94.
Tewes M, Kasimir-Bauer S, Welt A, Schuler M, Kimmig R, Aktas B. Detection of disseminated tumor cells in bone marrow and circulating tumor cells in blood of patients with early-stage male breast cancer. J Cancer Res Clin Oncol. 2015;141(1):87–92.
Raimondi C, Nicolazzo C, Gradilone A, et al. Circulating tumor cells: exploring intratumor heterogeneity of colorectal cancer. Cancer Biol Ther. 2014;15(5):496–503.
Salehi S, Steif A, Roth A, Aparicio S, Bouchard-Côté A, Shah SP. ddClone: joint statistical inference of clonal populations from single cell and bulk tumor sequencing data. Genome Biol. 2017;18(1):44.
Edsgärd D, Johnsson P, Sandberg R. Identification of spatial expression trends in single-cell gene expression data. Nat Methods. 2018;15(5):339–42.
Svensson V, Teichmann SA, Stegle O. SpatialDE: identification of spatially variable genes. Nat Methods. 2018;15(5):343–6.
Zhang J, Wang W, Huang J, Wang X, Zeng Y. How far is single-cell sequencing from clinical application? [published online ahead of print, 2020 Jul 5]. Clin Transl Med. 2020;10(3):e117.
This study was supported by Jiangsu Provincial Key Research and Development Special Fund (BE2015666), Jiangsu Innovative team leading talent fund (CXTDC2016006), Jiangsu Natural Science Foundation (BK20171304), Jiangsu six high peak talent fund (WSW-205), and Jiangsu 333 talent fund (BRA2016140).
Jiangsu Provincial Key Research and Development Special Fund (BE2015666), Jiangsu Innovative team leading talent fund (CXTDC2016006), Jiangsu Natural Science Foundation (BK20171304), Jiangsu six high peak talent fund (WSW-205), and Jiangsu 333 talent fund (BRA2016140).
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
Dai, Z., Gu, Xy., Xiang, Sy. et al. Research and application of single-cell sequencing in tumor heterogeneity and drug resistance of circulating tumor cells. Biomark Res 8, 60 (2020). https://doi.org/10.1186/s40364-020-00240-1