Selective treatment pressure in colon cancer drives the molecular profile of resistant circulating tumor cell clones
Molecular Cancer volume 20, Article number: 30 (2021)
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.
Liquid biopsy-based strategies may transform oncology in the near future . Many studies analyzed circulating tumor cells (CTCs) as a real-time liquid biopsy in epithelial tumors (i.e. breast, prostate, and colon cancer)  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) , 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]) . 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 . All our CTC lines harbored wild type PIK3CA and AKT , but these signaling pathways were deregulated. Some studies demonstrated that PI3K/mTOR pathway inhibitors could be used in primary and metastatic colorectal cancer . 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 .
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 . 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 .
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 . Analysis of CDA expression in colorectal cancer using publicly available data  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.
As CDA is secreted in the extracellular compartment (https://www.uniprot.org/uniprot/P32320) 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 220.127.116.11) 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 . Colorectal cancer is known to preferentially metastasize (∼70% patients) to the liver , 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 . 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 (https://www-ncbi-nlm-nih-gov.proxy.insermbiblio.inist.fr/geo) through the accession numbers: GSM1198032, GSM1198034, GSM1198038, GSM1198047, GSM1198050, GSM1198030, GSM1198035, GSM1198037, GSM1198046, GSM1198048, GSM1198023, GSM1198031, GSM1198033, GSM1198042, GSM1198052, GSM557089, GSM557091, GSM557093, GSM557099, GSM557107.
Circulating tumor cell
Epithelial to mesenchymal transition
Alix-Panabières C. The future of liquid biopsy. Nature. 2020;579(7800):S9.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Patanaphan V, Salazar OM. Colorectal cancer: metastatic patterns and prognosis. South Med J. 1993;86(1):38–41.
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.
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).
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.
CAP declare to have had an honorarium from Menarini. The other authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary materials and methods.
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.
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.
About this article
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). https://doi.org/10.1186/s12943-021-01326-6