Skip to main content

Selective treatment pressure in colon cancer drives the molecular profile of resistant circulating tumor cell clones


The characterization of circulating tumor cells (CTCs) holds promises for precision medicine because these cells are an important clinical indicator of treatment efficacy. We established the first and still only nine permanent colon CTC lines from peripheral blood samples of a patient with metastatic colon cancer collected at different time points during treatment and cancer progression. The study objectives were (i) to compare the gene expression profiles of these CTC lines, and (ii) to determine the main features acquired during treatment. The number of upregulated genes was higher in the CTC lines obtained after treatment, indicating that they acquired properties to escape treatment pressure. Among these upregulated genes, some are involved in the mTOR and PI3K/AKT signaling pathways. Moreover, cytidine deaminase expression was significantly increased in the CTC lines obtained after failure of the first- and second-line 5-fluorouracile-based treatments, suggesting that these CTCs can eliminate this specific drug and resist to therapy. Several enzymes involved in xenobiotic metabolism also were upregulated after treatment, suggesting the activation of detoxification mechanisms in response to chemotherapy. Finally, the significant higher expression of aldolase B in four of the six CTC lines obtained after treatment withdrawal and cancer progression indicated that these clones originated from liver metastases. In conclusion, these CTC lines generated at different time points during treatment of metastatic colon cancer in a single patient are characterized by the deregulation of different genes that promote (i) drug resistance, (ii) xenobiotic and energy metabolism, and (iii) stem cell properties and plasticity.

Main text

Liquid biopsy-based strategies may transform oncology in the near future [1]. Many studies analyzed circulating tumor cells (CTCs) as a real-time liquid biopsy in epithelial tumors (i.e. breast, prostate, and colon cancer) [2] and showed that their study contributes to the prognostic evaluation, patient stratification, and real-time monitoring of treatment efficacy, as well as to the identification of therapeutic targets and tumor resistance mechanisms. However, the low CTC concentration in blood is a crucial limiting factor for the identification of CTCs that may initiate a new tumor at a distant site (i.e. “metastasis-initiator CTCs” or “metastasis-competent CTCs”). Moreover, the establishment and long-term maintenance of in vitro CTC cultures are a major challenge, achieved only by few groups [3, 4].

We established the first nine permanent CTC lines from a 57-year-old patient with metastatic colon cancer (MCC). At diagnosis, he had unresectable widespread MCC with abdominal and mediastinal lymph node invasion and liver metastases. The patient was first treated with the 5-fluorouracile (5-FU)-irinotecan (FOLFIRI) combination and bevacizumab (5 cycles). The first biological progression was observed after the fifth FOLFIRI cycle, and second-line treatment was initiated with 5-FU-oxaliplatin (FOLXFOX) and bevacizumab, but was stopped after the third cycle because of disease progression (clinical and morphological). The patient died about 6 months after diagnosis due to cancer progression. We derived these CTC lines from blood samples collected before treatment initiation (CTC-MCC-41) [5], after the first- and second-line treatments (CTC-MCC-41.4), and 1 week before the patient’s death (CTC-MCC-41.5 [A-G]) [6]. Their phenotypic and molecular characterization indicated that they present common traits and all display epithelial-to-mesenchymal plasticity and stem cell-like characteristics [6, 7].

Here, we obtained and compared the transcriptomic profiles of these CTC lines to identify expression changes linked to tumor progression and treatment pressure, with the aim of identifying biomarkers and mechanisms involved in colon cancer progression.

Results and discussion

CTC selection during cancer progression and treatment pressure

The analysis of the transcriptomic profiles of the nine CTC lines using Affymetrix HG-U133P microarray chips clearly separated them in four distinct groups: CTC-MCC-41, CTC-MCC-41.4, CTC-MCC-41.5 [ABFG], and CTC-MCC-41.5 [CDE] (Fig. 1a). We then compared the transcriptomic profiles of the pre-treatment CTC-MCC-41 line and of the other eight CTC lines (CTC-MCC-41.4 and CTC-MCC-41.5[A-G]) to understand the treatment impact on CTC clonal selection (Fig. 1b and Additional file 2: Tables S1–2). The significantly higher number of upregulated genes in the eight cell lines obtained after treatment initiation (Fig. 1c) suggests that they acquired properties to adapt and resist to treatment. For instance, genes involved in the mTOR and PI3K/AKT signaling cascades, which are implicated in cancer development by coordinating cell growth, survival and proliferation, and in resistance to chemotherapy, were upregulated in the post-treatment CTC lines (Fig. 1d). Colon cancer-specific mortality is higher in patients with tumors harboring mutated PIK3CA than wild type PIK3CA [8]. All our CTC lines harbored wild type PIK3CA and AKT [6], but these signaling pathways were deregulated. Some studies demonstrated that PI3K/mTOR pathway inhibitors could be used in primary and metastatic colorectal cancer [9]. We recently reported that CTC-MCC-41 cells also respond to mTOR and AKT inhibitors, suggesting these therapies are effective even in the absence of mutations [10].

Fig. 1

Transcriptomic analysis of the nine CTC lines obtained from a patient with metastatic colon cancer at different times during disease progression. a Unsupervised 3D PCA representing the gene expression patterns of the CTC line obtained before treatment (n = 3 samples) and the eight CTC lines obtained after treatment initiation (n = 10 samples); b Hierarchical clustering of the CTC lines based on the differentially expressed transcripts (DETs) in CTC-MCC-41 (before treatment) and the other eight lines. The color intensity indicates the transcript expression level (red for upregulated transcripts and blue for downregulated transcripts); c Volcano plots showing the distribution of gene expression fold changes and p-values using the TAC software. Transcripts upregulated in CTC-MCC-41 are indicated in green, and transcripts upregulated in the other eight CTC lines are in red. d Top enriched pathways for the DETs (ANOVA analysis), ranked in increasing order (Fisher’s exact test p-value) identified by Ingenuity® Pathway Analysis, by comparing CTC-MCC-41 (before treatment) and the other eight lines. The “EIF2 signaling”, “mTOR and “PI3K/AKT signaling” pathways were significantly deregulated. Red: genes upregulated in CTC-MCC-41 compared with the other eight CTC lines; green: genes upregulated in the eight post-treatment CTC lines; white: no overlap with the data set; Orange - log (p-value); e Volcano plots showing the distribution of gene expression fold changes and p-values using the TAC software, with transcripts upregulated in the CTC-MCC-41.5 [ABFG] subgroup (green) and in the CTC-MCC-41.5 [CDE] subgroup (red). f Top enriched pathways for the DETs based on the ANOVA analysis, ranked in increasing order (Fisher’s exact test p-values) identified by Ingenuity® Pathway Analysis, in the comparison of the seven CTC-MCC-41.5 lines obtained just before the patient’s death. Many “metabolism” pathways were significantly deregulated, including “xenobiotic metabolism signaling”. Abbreviations: TAC, Transcriptome Analysis Console; DET, differentially expressed transcripts

Conversely, we did not observe many differences between the CTC-MCC-41.4 line, obtained after the last treatment, and the seven CTC-MCC-41.5 [A-G] lines obtained before the patient’s death (Additional file 2: Table S3). This suggests that, without drug pressure, the rapid cancer worsening was not link anymore to a clonal evolution but is due to natural disease progression with the replication of the CTC clones already selected under treatment. Indeed, the seven last CTC lines seem to be already present in the pooled CTC-MCC-41.4 line and they have been selected in the in vitro culture from the last blood sample. However, we could clearly segregate the last seven CTC cell lines (CTC-MCC-41.5 [A-G]), in two groups, [ABFG] and [CDE] (Fig. 1a), with different gene expression profiles (Fig. 1e). Most of the significantly deregulated pathways were involved in metabolism signaling (Fig. 1f), including xenobiotic metabolism. This suggests that detoxification mechanisms were induced upon exposure to anti-cancer drugs, as indicated by the deregulation of the irinotecan/SN38 pathway specifically in the [CDE] group. Furthermore, lipid metabolism upregulation appeared to be more represented in the “CDE signature”. Lipid metabolism is a key function on the basis of the enrichment of different signaling cascades leading to energy metabolism. Since Warburg’s work, metabolic reprogramming is one of the main hallmarks of cancer cells and plays a critical role in the continued tumor growth and progression and is driven by a complex interplay between the tumor mutational landscape, epigenetic modifications, and microenvironmental influences [11]. This topic is actively studied and a high-throughput metabolic-based assay was developed for rapid detection of rare metabolically active disseminated tumor cells in pleural effusion of lung cancer [12].

Cytidine deaminase as a drug resistance biomarker

Besides deregulation of the irinotecan/SN38 pathway in the [CDE] group, the gene encoding cytidine deaminase (CDA) was deregulated in the post-treatment CTC lines (Additional file 2: Table S2). This finding in a patient treated with 5-FU, a pyrimidine analogue, is interesting because CDA maintains the cellular pyrimidine pool by catalyzing the hydrolytic deamination of cytidine and deoxycytidine to uridine and deoxyuridine. Some studies showed that in patients treated with direct cytidine analogues, such as gemcitabine and cytosine arabinoside, CDA overexpression might be a marker of resistance [13]. Analysis of CDA expression in colorectal cancer using publicly available data [14] indicated that it was significantly downregulated in colon adenocarcinoma compared with normal colon (Fig. 2a). However, CDA was strongly upregulated in the CTC-MCC41.4 cell line obtained directly after failure of the first- and second-line 5-FU-based treatments (RT-qPCR analysis in Fig. 2b) and also, to a lower extent, in the CTC lines (CTC-MCC-41.5 A-G) obtained just before death, compared with the pre-treatment CTC-MCC-41 line.

Fig. 2

Major features of the aggressive colon cancer CTC clones selected after treatment. a Boxplot showing cytidine deaminase (CDA) relative expression in normal colon (n = 41) and colon adenocarcinoma (n = 286) samples using UALCAN (p-value = 1.45E-06). b RT-qPCR analysis of CDA gene expression in the nine CTC lines (mean values of triplicate experiments with standard deviation). c ELISA measurement of CDA protein level (ng/ml) in the conditioned medium from the nine CTC lines (mean values from triplicate experiments with standard deviation). Aldolase B (ALDOB) expression analysis (d) by RT-qPCR in the nine CTC lines [*significant upregulation in the CTC-MCC-41.5 [ABFG] cell lines compared with all the other CTC lines (p-value = 0.0134), and only with the CTC-MCC-41.5 [CDE] lines (p-value = 0.0466); **significantly upregulation compared with all the other CTC lines (p-value = 0.0017)], and e in tissue samples (normal colon, normal liver, primary colorectal tumor, liver metastases of colon cancer) (n = 5/tissue). f Representation of the cancer hallmarks of the CTC lines resistant to treatment: (i) Stemness and Plasticity: metastasis-competent potential, (ii) Drug resistance: different molecules and pathways involved in treatment resistance and drug detoxification, such as the CDA enzyme implicated in pyrimidine analogue metabolism (e.g. for 5-FU), (iii) Metabolism: xenobiotic metabolism linked to therapy resistance, energy metabolism linked to mitochondrial activity, and ALDOB activity that gave information also on the liver metastasis origin of the CTC41.5[ABFG] lines. Abbreviation: EMT, epithelial-mesenchymal transition

As CDA is secreted in the extracellular compartment ( and can be detected in blood, we quantified CDA concentration by ELISA in conditioned medium from the nine CTC-MCC lines. As observed for the CDA gene, CDA protein level was highest in CTC-MCC-41.4 cells (Fig. 2c). These findings suggest that in this patient, CDA was directly produced and secreted by resistant and aggressive CTCs in response to 5-FU-based chemotherapy. Therefore, CDA might represent a candidate plasmatic biomarker to monitor 5-FU efficacy and resistance development.

Aldolase B as a marker to identify CTCs released only by liver metastases

Comparison of the signatures of the two CTC-MCC-41.5 sub-groups ([ABFG] and [CDE]) (Additional file 2: Tables S4–5) highlighted ALDOB upregulation in the [ABFG] group (Fig. 2d). This gene encodes the aldolase B enzyme (fructose-bisphosphate aldolase B or liver-type aldolase), one of three isoenzymes (A, B, and C) of the class I fructose 1,6-bisphosphate aldolase enzyme (EC that plays a key role in glycolysis and gluconeogenesis. This enzyme is preferentially expressed in the liver and at lower extend in the kidney and the small intestine. ALDOB downregulation correlates with poor overall survival in liver and gastric cancer [15, 16], whereas ALDOB overexpression in colorectal cancer has been associated with poor overall survival and epithelial-to-mesenchymal transition promotion [17]. Colorectal cancer is known to preferentially metastasize (70% patients) to the liver [18], which is the main organ for glycogenesis and gluconeogenesis. In addition, this specific patient had liver metastases. Comparison of ALDOB expression in different tissues using publicly available data showed that: (i) ALDOB was strongly expressed in normal liver compared with normal colon, (ii) ALDOB expression was comparable in primary colon adenocarcinoma and normal colon samples, and (iii) ALDOB was upregulated in colon cancer liver metastases compared with primary tumors (Fig. 2e). Bu et al. showed that during liver colonization, colon cancer cells undergo metabolic reprogramming by upregulating ALDOB [19]. This enhances fructose metabolism and promotes the growth of colon cancer liver metastases. These data suggest that the [ABFG] lines were derived from CTCs released by liver metastases. It is thought that CTCs are released by the primary tumor and/or metastases; however, to our knowledge, this is the first time that this could be demonstrated in a patient.


From Nowell’s work in 1976, it is acknowledged that cancer is an evolutionary process and that treatments exert selective pressures that drive the tumor cell evolution, favoring the appearance of resistant clones. The present results confirm that CTC line profiling is a relevant approach to study clonal selection during disease progression and to discover new CTC biomarkers for monitoring treatment response. The transcriptomic analysis of the CTC lines obtained after treatments (CTC-MCC-41.4 and CTC-MCC-41.5[A-G]) showed the progressive deregulation of genes involved in cancer aggressiveness hallmarks: (i) drug resistance, with the upregulation of molecules and pathways implicated in drug detoxification (e.g. CDA and the irinotecan/SN38 pathway), (ii) metabolism changes, including upregulation of genes implicated in the xenobiotic metabolism, linked to therapy resistance, and energy metabolism (mitochondrial activity to provide energy to hyper-metabolic cells), and the ALDOB gene linked to liver metabolism, and (iii) stemness and plasticity, highlighting their metastasis-competent potential (Fig. 2f). This work is, however, specific of this cancer patient and the biomarkers found and highlighted in these CTC lines need now to be validated on CTCs and plasma from independent patients with colorectal cancer.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

The 20 publicly available Affymetrix datasets on normal colon (n = 5 samples), primary tumor (n = 5 samples), liver metastases (n = 5 samples) and normal liver (n = 5 samples) are accessible at the gene expression Omnibus (GEO) repository ( through the accession numbers: GSM1198032, GSM1198034, GSM1198038, GSM1198047, GSM1198050, GSM1198030, GSM1198035, GSM1198037, GSM1198046, GSM1198048, GSM1198023, GSM1198031, GSM1198033, GSM1198042, GSM1198052, GSM557089, GSM557091, GSM557093, GSM557099, GSM557107.





Aldolase B


Cytidine deaminase


Circulating tumor cell


Epithelial to mesenchymal transition


  1. 1.

    Alix-Panabières C. The future of liquid biopsy. Nature. 2020;579(7800):S9.

    Article  Google Scholar 

  2. 2.

    Pantel K, Alix-Panabieres C. Liquid biopsy and minimal residual disease - latest advances and implications for cure. Nat Rev Clin Oncol. 2019;16(7):409–24.

    CAS  Article  Google Scholar 

  3. 3.

    Koch C, Kuske A, Joosse SA, Yigit G, Sflomos G, Thaler S, et al. Characterization of circulating breast cancer cells with tumorigenic and metastatic capacity. EMBO Mol Med. 2020;12(9):e11908.

    CAS  Article  Google Scholar 

  4. 4.

    Cortés-Hernández LE, Eslami SZ, Alix-Panabières C. Circulating tumor cell as the functional aspect of liquid biopsy to understand the metastatic cascade in solid cancer. Mol Asp Med. 2020;72:100816.

    Article  Google Scholar 

  5. 5.

    Cayrefourcq L, Mazard T, Joosse S, Solassol J, Ramos J, Assenat E, et al. Establishment and characterization of a cell line from human circulating colon cancer cells. Cancer Res. 2015;75(5):892–901.

    CAS  Article  Google Scholar 

  6. 6.

    Soler A, Cayrefourcq L, Mazard T, Babayan A, Lamy PJ, Assou S, et al. Autologous cell lines from circulating colon cancer cells captured from sequential liquid biopsies as model to study therapy-driven tumor changes. Sci Rep. 2018;8(1):15931.

    Article  Google Scholar 

  7. 7.

    Alix-Panabieres C, Cayrefourcq L, Mazard T, Maudelonde T, Assenat E, Assou S. Molecular portrait of metastasis-competent circulating tumor cells in Colon Cancer reveals the crucial role of genes regulating energy metabolism and DNA repair. Clin Chem. 2017;63(3):700–13.

    CAS  Article  Google Scholar 

  8. 8.

    Ogino S, Nosho K, Kirkner GJ, Shima K, Irahara N, Kure S, et al. PIK3CA mutation is associated with poor prognosis among patients with curatively resected colon cancer. J Clin Oncol. 2009;27(9):1477–84.

    CAS  Article  Google Scholar 

  9. 9.

    Narayanankutty A. PI3K/ Akt/ mTOR pathway as a therapeutic target for colorectal Cancer: a review of preclinical and clinical evidence. Curr Drug Targets. 2019;20(12):1217–26.

    CAS  Article  Google Scholar 

  10. 10.

    Smit DJ, Cayrefourcq L, Haider MT, Hinz N, Pantel K, Alix-Panabières C, et al. High Sensitivity of Circulating Tumor Cells Derived from a Colorectal Cancer Patient for Dual Inhibition with AKT and mTOR Inhibitors. Cells. 2020;9(9):2129.

    CAS  Article  Google Scholar 

  11. 11.

    Peitzsch C, Gorodetska I, Klusa D, Shi Q, Alves TC, Pantel K, et al. Metabolic regulation of prostate cancer heterogeneity and plasticity. Semin Cancer Biol. 2020;S1044-579X(20):30261–3.

    Google Scholar 

  12. 12.

    Tang Y, Wang Z, Li Z, Kim J, Deng Y, Li Y, et al. High-throughput screening of rare metabolically active tumor cells in pleural effusion and peripheral blood of lung cancer patients. Proc Natl Acad Sci. 2017;114(10):2544–9.

    CAS  Article  Google Scholar 

  13. 13.

    Weizman N, Krelin Y, Shabtay-Orbach A, Amit M, Binenbaum Y, Wong RJ, et al. Macrophages mediate gemcitabine resistance of pancreatic adenocarcinoma by upregulating cytidine deaminase. Oncogene. 2014;33(29):3812–9.

    CAS  Article  Google Scholar 

  14. 14.

    Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi B, et al. UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses. Neoplasia (New York, NY). 2017;19(8):649–58.

    CAS  Article  Google Scholar 

  15. 15.

    He J, Jin Y, Chen Y, Yao HB, Xia YJ, Ma YY, et al. Downregulation of ALDOB is associated with poor prognosis of patients with gastric cancer. Onco Targets Ther. 2016;9:6099–109.

    Article  Google Scholar 

  16. 16.

    Tao QF, Yuan SX, Yang F, Yang S, Yang Y, Yuan JH, et al. Aldolase B inhibits metastasis through ten-eleven translocation 1 and serves as a prognostic biomarker in hepatocellular carcinoma. Mol Cancer. 2015;14:170.

    Article  Google Scholar 

  17. 17.

    Li Q, Li Y, Xu J, Wang S, Xu Y, Li X, et al. Aldolase B overexpression is associated with poor prognosis and promotes tumor progression by epithelial-Mesenchymal transition in colorectal adenocarcinoma. Cell Physiol Biochem. 2017;42(1):397–406.

    CAS  Article  Google Scholar 

  18. 18.

    Patanaphan V, Salazar OM. Colorectal cancer: metastatic patterns and prognosis. South Med J. 1993;86(1):38–41.

    CAS  Article  Google Scholar 

  19. 19.

    Bu P, Chen KY, Xiang K, Johnson C, Crown SB, Rakhilin N, et al. Aldolase B-Mediated Fructose Metabolism Drives Metabolic Reprogramming of Colon Cancer Liver Metastasis. Cell Metab. 2018;27(6):1249–62.e4.

    CAS  Article  Google Scholar 

Download references


The authors thank Dr. Elisabetta Andermarcher for assistance with her comments and proofreading that greatly improved the manuscript.


This work was supported by the National Cancer Institute (INCA) and the General Direction for Caregiving (DGOS) for patient recruitment and analysis in the COLOSPOT study (NCT01596790), by SIRIC Montpellier Cancer Grant INCa_Inserm_DGOS_12553 for staff salary as well as by the MAVA Foundation, the CNRS and the ANR project TRANSCAN (ANR-18-CE35–0009).

Author information




LC and SA analyzed and interpreted the data. LC, SA and CAP were major contributors in the conception of the study and in the writing of the manuscript. All authors red and approved the final manuscript.

Corresponding authors

Correspondence to Said Assou or Catherine Alix-Panabières.

Ethics declarations

Ethics approval and consent to participate

The patient gave his consent to participate in the COLOSPOT study to investigate CTCs in metastatic colorectal cancer. The bioethics committee gave its approval (ID-RCB 2011-A01130–41).

Consent for publication

The content of this manuscript has not been previously published and is not under consideration for publication elsewhere.

Competing interests

CAP declare to have had an honorarium from Menarini. The other authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Supplementary materials and methods.

Additional file 2: Table S1.

List of DETs upregulated in CTC-MCC-41 vs all the other CTC lines. Table S2. List of DETs downregulated in CTC-MCC-41 vs all the other CTC lines. Table S3. List of DETs upregulated in all CTC-MCC-41.5 cell lines vs the CTC-MCC-41.4 cell line. Table S4. List of DETs upregulated in the [ABFG] subgroup vs the [CDE] subgroup of the CTC-MCC-41.5 cell lines. Table S5. List of DETs downregulated in the [ABFG] subgroup vs the [CDE] subgroup of the CTC-MCC-41.5 cell lines. Table S6. Primer sequences.

Additional file 3: Figure S1.

Box-plots showing the expression level of a set of differentially expressed transcripts (DETs) in the CTC-MCC-41, CTC-MCC-41.4, and CTC-MCC-41.5 [ABFG] and [CDE] cell lines. Abbreviations: CTC-BT, CTC line derived before treatment initiation; CTC-AT, CTC lines derived after treatment initiation. Figure S2. RT-qPCR validation of the microarray data. Abbreviations: CTC-BT, CTC line derived before treatment initiation; CTC-AT, CTC lines derived after treatment initiation.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Cayrefourcq, L., Thomas, F., Mazard, T. et al. Selective treatment pressure in colon cancer drives the molecular profile of resistant circulating tumor cell clones. Mol Cancer 20, 30 (2021).

Download citation


  • Circulating tumor cells
  • Gene expression
  • Clonal evolution
  • Colon cancer
  • CDA