Skip to main content
  • Research article
  • Open access
  • Published:

Comparative sequence analysis of patient-matched primary colorectal cancer, metastatic, and recurrent metastatic tumors after adjuvant FOLFOX chemotherapy



In the era of genome-guided personalized cancer treatment, we must understand chemotherapy-induced genomic changes in tumors. This study evaluated whether adjuvant FOLFOX chemotherapy modifies the mutational profile of recurrent colorectal cancer (CRC).


Whole exome sequencing was performed on samples from primary CRC tumors, untreated metastatic tumors, and recurrent tumors following adjuvant FOLFOX chemotherapy. The samples were resected from four patients.


The number of mutations or the mutation spectrum in individual patients was nearly identical. Copy number variants persisted regardless of FOLFOX therapy administration. The genomic signature of oxaliplatin exposure (G > T/C > A, T > A/A > T) was not enriched after FOLFOX chemotherapy. Overlapping single nucleotide variants (SNVs) and indels remained in 26–65% of the patient-matched tumor samples. One patient harbored an AKT1 E17K mutation in the recurrent tumor, whereas PIK3CA E542K and E88Q mutations were detected in the primary and untreated metastatic tumor samples. Genes related to intracellular Ca2+ homeostasis were enriched among the genes uniquely mutated after FOLFOX chemotherapy.


We found that the mutation rates, mutation spectrum, and copy number variants were nearly identical regardless of the administration of FOLFOX therapy in the four CRC cases. The mutational discordance between the patient-matched tumor samples is likely caused by tumor heterogeneity and chemotherapy-induced clonal selection. These findings might be useful as pilot data for larger studies to clarify the changes in the mutational landscape induced by adjuvant FOLFOX chemotherapy.

Peer Review reports


Colorectal cancer (CRC) is the third most common type of cancer and the fourth leading cause of cancer death worldwide [1]. Combination chemotherapy with cytotoxic and molecular targeting agents has prolonged the survival time of patients with metastatic CRC. In the era of genome-based personalized cancer treatment, an assessment of the molecular profile of individual tumors is necessary to guide the selection of appropriate therapy methods. For example, activating KRAS and NRAS mutations are negative predictive markers for the effectiveness of the anti-epidermal growth factor receptor (EGFR) antibodies panitumumab and cetuximab [2, 3]. Intensive chemotherapy is recommended for patients with the BRAF V600E mutation because this mutation is a strong prognostic factor for poor survival [4, 5]. Moreover, several genetic alterations that are potential prognostic and predictive biomarkers or therapeutic targets have been explored. Extensive data sets of the mutational profiles of CRC have been generated [6], and large collaborations have created gene expression-based classifications that predict patient outcomes [7].

However, systemic chemotherapies could alter the mutational landscape of several cancers [8, 9]. A previous exome sequencing study revealed that mutagenic chemotherapy regimens, such as adjuvant chemotherapy with the DNA-alkylating-like agent temozolomide to treat glioma, can induce new mutations and cause the malignant progression of recurrent tumors [9].

FOLFOX is a combination chemotherapy regimen that consists of leucovorin-modulated 5-fluorouracil (5-FU) and oxaliplatin (L-OHP), which are commonly used worldwide as standard adjuvant chemotherapies for curatively resected stage III and IV CRCs [10]. L-OHP is a third-generation platinum (Pt)-containing antitumor compound that induces DNA damage associated with intra- and inter-strand cross-links (Pt-GG adducts) [11,12,13]. Previously, in vitro studies have demonstrated the mutagenic activity of L-OHP [14]. Therefore, adjuvant FOLFOX chemotherapy has the potential to alter the mutational profiles of recurrent cancers so that they differ from those of primary CRC tumors. Our previous report showed that the mutational status of predictive biomarker genes for the effectiveness of anti-EGFR-antibodies was not altered by FOLFOX therapy [15]. However, the influence of FOLFOX therapy on exome-wide mutational profiles has not been reported previously.

This study used whole exome sequencing to compare gene alteration profiles of recurrent cancers after adjuvant FOLFOX chemotherapy in patient-matched primary CRC and metastatic tumor samples prior to chemotherapy.


Patient selection

We reviewed the clinical records of patients with CRC who had been treated with adjuvant FOLFOX chemotherapy after curative resection at the National Cancer Center Hospital East (Kashiwa, Japan). From January 2006 to December 2009, 156 patients were treated with adjuvant FOLFOX at our institution, and 66 patients developed recurrent tumors during or after adjuvant FOLFOX chemotherapy. Of these patients, 26 underwent curative resection of recurrent tumors. We selected four CRC patients for whom tumor specimens from primary tumors, metastases resected prior to FOLFOX chemotherapy (pre-FOLFOX metastasis), and recurrent metastatic tumors in the same organ after adjuvant FOLFOX (post-FOLFOX metastasis) were available. Informed consent to use tissue specimens for this study was obtained from all patients, and the tissue samples were provided by the National Cancer Center Biobank, Japan. The institutional review board at the National Cancer Center approved the study protocol. This study was performed according to the Epidemiological Study Guidelines of the Ministry of Health, Labor, and Welfare of Japan. We disclosed the study design on the National Cancer Center website and gave the relatives of deceased patients the opportunity to decline participation.

DNA samples

We obtained matched primary CRC, pre-FOLFOX metastasis, post-FOLFOX metastasis, and normal colorectal tissue samples from four patients. The normal colorectal tissues were collected from surgical specimens of the primary tumors. All tissue samples were formalin-fixed, paraffin-embedded (FFPE) specimens. DNA samples were obtained from macroscopically dissected FFPE specimens cut into 10-μm-thick sections. Genomic DNA was extracted using EZ1 Advanced XL and EZ1 DNA Tissue Kits (Qiagen, Hilden, Germany) according to the manufacturer’s instructions [16]. Nucleic acid yields were determined using a NanoDrop 2000 (Thermo Fisher Scientific, Waltham, MA, USA), and the quality of genomic DNA was examined using a Quant-iT picoGreen dsDNA (Life Technologies, Carlsbad, CA, USA) assay kit.

Whole exome sequencing and variant calling

Using genomic DNA from tumors and matched normal samples, we performed exome capture sequencing. Using an Agilent SureSelect Human All Exome V5 + UTRs kit (Agilent Technologies, Santa Clara, CA, USA), whole exome sequencing was performed using an Illumina HiSeq 2000 system (Illumina, San Diego, CA, USA) to generate 100-bp paired-end sequencing reads according to the manufacturer’s instructions. Burrows-Wheeler Aligner (BWA, [17] was used to align the sequencing reads to the human reference genome (hg19). The Genome Analysis ToolKit version 1.6 (GATK, was used for the local realignment and score recalibration of the sequencing reads [18]. We employed Picard ( for the basic processing and management of the sequencing data. To reduce the false-positive rate, the following filtering criteria were applied: (i) GATK confidence score [18] ≥ 50; (ii) number of forward and reverse reads ≥1; and (iii) variants present in at least 10% of the reads. All mutations detected in the paired non-tumor colon tissues were excluded from our analysis. We also excluded alterations present in dbSNP151, the 1000 Genomes Project, and in-house Japanese exomes derived from 299 normal tissues in our previous studies, with the aim of identifying tumor-specific variants. We also performed a visual inspection to filter out false-positive variants.

All mutations with clinical inference were annotated using ANNOVAR [19].

Gene ontology analysis

A Gene Ontology (GO) analysis was performed using the Database for Annotation, Visualization, and Integrated Discovery (DAVID, [20]. Adjusted P-values less than 0.05 were considered significant. The GO analysis was performed for nonsense mutations, small insertions/deletions (indels), and missense mutations that were predicted as “probably damaging” and “possibly damaging” by PolyPhen2 ( [21].

Copy number variant analysis

The log ratio of the depth of coverage between the tumor and normal colorectal tissues was calculated using the GATK-depth of coverage tool. Then, copy number variant (CNV) segments were identified from the log ratio of the depth of coverage using the Exome CNV R package [22]. Log ratios of the depth of coverage that were greater than two were considered indicative of significant copy number amplification.

Statistical analysis

The Wilcoxon signed-rank test was used to evaluate the difference in the number of genetic alterations between primary, pre-FOLFOX metastatic, and post-FOLFOX metastatic tumors. Increases in specific mutation types among the post-FOLFOX unique mutations were also assessed using this method. Microsoft Office Excel 2013 (Microsoft Corporation, Redmond, WA, USA) was used to perform all the statistical analyses.


Patient characteristics and clinical courses

Table 1 shows the patient characteristics. There were two male and two female participants with a median age of 68 years. The primary tumor sites were the colon in one patient and the rectum in three patients. After curative resection of their primary and metastatic tumors, all patients were treated with a modified FOLFOX6 (mFOLFOX6) regimen including an L-OHP dose of 85 mg/m2 administered every 14 days; 12 treatment cycles were planned [23]. Post-FOLFOX metastasis developed during (cases 1 and 2) or after (cases 3 and 4) adjuvant chemotherapy. Histopathological analyses diagnosed the primary CRC tumors as well-differentiated adenocarcinoma (cases 2 and 4) and moderately differentiated adenocarcinoma (cases 1 and 3). All metastatic tumors exhibited histology concordant with the corresponding primary colorectal adenocarcinoma.

Table 1 Patient characteristics

Additional file 1 Figure S1 summarizes the clinical courses of these four patients. In case 1, adjuvant mFOLFOX6 was initiated after the colectomy for primary sigmoid colon cancer and hepatic resection was performed for the synchronous colorectal metastases. After three cycles of mFOLFOX6, recurrence in the remnant liver was examined by computed tomography (CT) imaging. In case 2, the patient underwent high anterior resection and liver metastasectomy. Early recurrence in the liver was identified by CT imaging after three cycles of adjuvant mFOLFOX6, which was continued according to the clinician’s discretion. In total, the patient received nine cycles of mFOLFOX6 before resection of the liver recurrence. In case 3, the patient was diagnosed with stage IIA rectal adenocarcinoma, received lower anterior resection without chemoradiotherapy, and was followed without adjuvant chemotherapy. Liver metastasis was diagnosed 14 months after the first operation. Adjuvant mFOLFOX6 was administered following liver metastasectomy, and adjuvant therapy was discontinued at 11 cycles as a result of intolerable peripheral sensory neuropathy. Liver recurrence was identified four months after the end of adjuvant chemotherapy. In case 4, the patient underwent lower anterior resection and lung metastasectomy for metastatic diseases. Nine months after completing the planned adjuvant chemotherapy, lung recurrence was identified by CT imaging.

A median of 10 mFOLFOX6 cycles was reported in this study (range, 4–12 cycles), and the median disease-free survival was 218.5 days (range, 97–556 days).

Comparison of the mutation rates and analysis of the mutation spectrum

Whole exome sequencing was performed to investigate the profile of somatic alterations in all tumor samples. The tumors were sequenced to an average 124-fold coverage (range, 90–155), enabling the sensitive detection of single nucleotide variants (SNVs) and indels to a 10% variant frequency. To identify germline mutations, we sequenced the paired non-tumor colon tissues in addition to the tumor tissues from each patient. Additional file 2: Table S1 shows a summary of the coverage details.

The mutation rate of the somatic SNVs and indels in each tumor sample ranged from 0.90 to 3.26/Mb, with a median of 1.66/Mb. The rates were consistent with those of non-hypermutated CRC cases reported by The Cancer Genome Atlas (TCGA) [6]. Although the mutation rates increased slightly after FOLFOX administration, there were no significant differences between the matched primary CRC, pre-FOLFOX metastatic, and post-FOLFOX metastatic samples (P > 0.05, Fig. 1a).

Fig. 1
figure 1

a Comparison of mutation rates: Box-and-whisker plot showing the range, median, and 25th and 75th percentile levels of the number of mutations in all primary, pre-FOLFOX metastatic, and post-FOLFOX metastatic tumor samples. The difference in the number of mutations was not significant among the primary, pre-FOLFOX metastatic, and post-FOLFOX metastatic tumor samples across all cases (P > 0.05, Wilcoxon signed-rank test) b Comparison of mutation spectra: Stacked bar plot showing the proportion of mutations accounted for by each of the six mutation types in all primary, pre-FOLFOX metastatic, and post-FOLFOX metastatic tumor samples for each case

Next, the mutation spectrum of each tumor was investigated. All tumor samples harbored a predominance of C > T/G > A transitions (range, 45–70%), which is reported as a mutational signature of CRC [6, 24] (Fig. 1b). C > A/G > T and T > A/A > T transversions were reported previously as the most common mutations produced by L-OHP in cultured cells [14]. In addition, C > A mutations were particularly prominent in a CpC context (equivalently, G > T in a GpG context) as a result of intrastrand cross-links induced by cisplatin, which is a Pt-containing antitumor compound similar to L-OHP [25]. Therefore, we evaluated whether these platinum exposure signature mutations were enriched among the post-FOLFOX unique mutations. Figure 2 shows the mutation spectra of the primary CRC, pre-FOLFOX, and post-FOLFOX unique mutations. There was no significant increase in C > A/G > T (primary vs post-FOLFOX metastasis, P = 0.27; pre- vs post-FOLFOX metastasis, P = 1.00) or T > A/A > T (primary vs post-FOLFOX metastasis, P = 1.00; pre- vs post-FOLFOX metastasis, P = 0.14) transversions among the post-FOLFOX unique mutations. The incidence of C > A mutations in a CpC context was not increased in post-FOLFOX metastasis (primary vs post-FOLFOX metastasis, P = 0.18; pre- vs post-FOLFOX metastasis, P = 0.18). In addition, there were no significant differences in other mutation fractions (Additional file 3: Figure S2).

Fig. 2
figure 2

Comparison of the mutation spectrum among unique mutations: The relative frequencies of 96 trinucleotides are shown for unique mutations in the primary, pre-FOLFOX metastatic, and post-FOLFOX metastatic tumor samples. The presented number of mutations is the sum of the data recorded for the four patients. For specific mutation types, a Wilcoxon signed-rank test was used, and significant enrichment in the platinum exposure signature, C > A in a CpC context (equivalently, G > T in a GpG context), was not observed

Analysis of CNVs

To evaluate whether the CNV was affected by FOLFOX treatment, the CNVs in post-FOLFOX metastatic tumors were compared to those in primary tumors and pre-FOLFOX metastatic tumors. Our comparative analysis revealed that CNVs persisted regardless of the administration of FOLFOX therapy in any tumor sample in cases 1, 2, and 4 (Fig. 3). In case 3, we observed focal amplification of the 7q21, 10q22, and 10q23 chromosomal regions, which persisted regardless of FOLFOX therapy administration. However, the amplification of some genes, such as SEMA3E, SEMA3A, PCLO, AK055932, and BX647900, was observed only in pre- and post-FOLFOX metastasis and not in the primary tumor in case 3 (Additional file 4: Table S2).

Fig. 3
figure 3

Comparison of copy number alterations: The CNV patterns were similar among primary (blue), pre-FOLFOX metastatic (purple), and post-FOLFOX metastatic (red) tumor samples in all cases. The x-axis coordinates represent positions along the genome; the vertical bars indicate the borders between chromosomes. The y-axis represents the chromosomal copy number in the tumor compared to normal tissues. A log ratio greater than two was considered significant copy number amplification

Overlap of SNVs and indels between patient-matched tumor samples

The overlap of the detected mutations between the primary, pre-FOLFOX metastatic, and post-FOLFOX metastatic tumor samples in the individual cases was investigated. Additional file 5: Table S3 lists the details of all detected gene mutations. Of the gene mutations detected in each post-FOLFOX metastatic sample, 112 (82%) in case 1, 114 (46%) in case 2, 71 (73%) in case 3, and 125 (86%) in case 4 were shared in the matched tumor samples, indicating that 14–54% of the mutations were post-FOLFOX unique mutations (Fig. 4a).

Fig. 4
figure 4

a Overlap of SNVs and indels: All somatic mutations and indels were classified into three categories: mutations present in all tumor regions (common, blue), mutations shared in more than two but not all regions (shared, green), and mutations in one tumor region (unique, yellow) in each case. The distributions of all somatic mutations in each tumor sample are shown in a stacked bar plot. The number of analyzed mutations is displayed on the right side of each bar. b Comparison of Mut-driver gene mutations: Mutations in 14 Mut-driver genes were detected in all tumor samples, and eight of these were classified as common mutations. A discordant mutation in the PI3K pathway before and after FOLFOX administration was observed in case 2

Next, we focused on Mut-driver genes, as advocated by Vogelstein et al. [26]; these genes include tumor suppressor genes for which at least 20% of the mutations caused truncation of the gene product and oncogenes for which at least 20% of the missense mutations occurred at a single position along the polypeptide chain. High mutational concordance of these genes, including frequent mutation in CRC genes such as APC, KRAS, FBXW7, TP53, and PIK3CA, was observed between the matched tumor samples in the individual cases (Fig. 4b). However, in case 2, a lack of PIK3CA mutations was found in the post-FOLFOX metastatic samples, although the PIK3CA E542K and E88Q mutations were detected in both the primary tumor and pre-FOLFOX metastatic samples. By contrast, AKT1 E17K (C > T) mutations were found only in the post-FOLFOX metastasis samples in case 2.

After FOLFOX administration, we also identified the gain or loss of some mutations that were predicted as functionally important variants by Polyphen2 [21]. Additional file 6: Table S4 lists the details of these mutations.

Gene ontology analysis of post-FOLFOX unique mutations

Post-FOLFOX unique mutations in recurrent tumors may reflect the mechanism of chemo-resistance because recurrent tumors are thought to develop through chemotherapy-induced selective pressure. A GO analysis was performed to identify the molecular functions enriched in post-FOLFOX metastatic samples. For this analysis, we selected gene mutations that predict a possible impact of an amino acid substitution on the structure and function of a human protein (the analyzed gene lists are shown in Additional file 7: Table S5). Although no identical gene mutations were detected in more than two cases, four significant functional clusters were identified (Table 2). The genes annotated to the GO term “calcium ion binding” were detected in all cases, and genes annotated to “calcium ion transport” were detected in three cases. These GO terms were not enriched among the pre-FOLFOX unique gene mutations or among common mutations (Additional file 8: Table S6).

Table 2 Significant functional clusters among post-FOLFOX unique mutations


In this study, we performed whole exome sequencing on primary colorectal, metastatic, and recurrent tumor samples after adjuvant FOLFOX therapy to evaluate the influence of this cytotoxic chemotherapy regimen on the mutational profile of recurrent CRC. Recently, a wide range of genomic alterations have displayed associations with cancer behavior, and most of these alterations are typically found in coding regions [27]. Therefore, whole exome sequencing is a reasonable strategy for identifying clinically actionable alterations induced by adjuvant FOLFOX chemotherapy in recurrent CRC. Furthermore, one of the most important aims of this study was to evaluate whether adjuvant FOLFOX chemotherapy produces genetic alterations in recurrent CRC. Thus, we excluded variants that could possibly be regarded as germline alterations from our analysis by using the sequencing results from patient-matched normal tissues and previous SNP databases.

Our data showed that the mutation rates and mutation spectrum were nearly identical between the primary CRC, pre-FOLFOX metastatic, and post-FOLFOX metastatic samples. Although the mutagenic activity of L-OHP was demonstrated in cultured cells [14], an enrichment in the L-OHP signature was not observed in post-FOLFOX metastatic tumors. Differences between the results of that study and our results may be due to differences in the drug concentration used in that in vitro study and the concentration used in clinical practice. A significant dose-dependent increase in mutation frequency was observed when CHO-K1 cells were exposed to 10–40 μM L-OHP [14]. However, the Cmax of L-OHP has been reported as approximately 3.6 μM in the clinical dose setting with mFOLFOX6 therapy (85 mg/m2) [28], which might explain the differences between the in vitro findings and those of the present study. A recent study revealed that some recurrent gliomas were hypermutated and harbored driver mutations in the RB and Akt-mTOR pathways that bore the signature of temozolomide-induced mutagenesis after adjuvant temozolomide chemotherapy. Recurrent gliomas that showed evidence of temozolomide-induced hypermutation underwent malignant progression to high-grade tumors with a poorer prognosis [9]. By contrast, our study suggests that FOLFOX is a safe regimen that lacks the potential risk of inducing new driver mutations and malignant progression, unlike adjuvant chemotherapy involving temozolomide to treat glioma.

The patterns of CNVs were nearly identical before and after FOLFOX chemotherapy in all the cases. The focal amplification of the 7q21, 10q22, and 10q23 chromosomal regions in the post-FOLFOX metastasis samples was also detected in the primary or pre-FOLFOX metastasis samples in case 3. However, we revealed that some gene amplifications were observed only in the pre- and post-FOLFOX metastasis samples and not in the primary tumor samples in case 3. Among these genes, we focused on SEMA3E, which has an expression level reported to be positively correlated with increased metastasis in ovarian, melanoma, and colon cancers [29]. Our results suggest that the acquisition of SEMA3E gene amplification might relate to the metastatic potential in CRC and that CNVs can change during the process of tumor progression.

One of the limitations of our study is its small sample size due to the strict patient selection criteria. Several studies have reported that differences in mutational profiles depend on the metastatic site [30,31,32]. Therefore, it was necessary to analyze patient- and organ-matched pre-FOLFOX metastasis and post-FOLFOX metastasis samples. However, such samples are rarely found in daily clinical practice. Our analysis used valuable samples from four CRC cases and accurately assessed the influence of adjuvant FOLFOX chemotherapy on the mutational profile of CRC. The other limitation is the possibility that our cases may have had primary resistance to FOLFOX therapy. This resistance may have been why we failed to observe a chemotherapy-induced signature in the recurrent tumor samples. Our findings should be further validated in a large cohort that includes multiple patients with various backgrounds, such as variations in the duration of FOLFOX chemotherapy or the timing of recurrence.

It is also possible that there was simply not enough time to observe a chemotherapy-induced signature in the recurrent tumor samples that we used. Twelve-cycle FOLFOX chemotherapy was planned as the full adjuvant chemotherapy course, but it was often discontinued due to tumor relapse or intolerable adverse events. Our data reflect a correlation between adjuvant FOLFOX chemotherapy and the mutation profile of CRC in clinical practice, but further investigation is necessary regarding whether long-term treatment with FOLFOX chemotherapy, such as chemotherapy for unresectable CRC, can induce additional mutations. A previous study suggested that a circulating tumor DNA analysis could be used to reliably monitor tumor dynamics in subjects with cancer who were undergoing surgery or chemotherapy for CRC [33]. This non-invasive approach may be a more suitable way to test residual tumor cells and analyze temporal changes in mutation profiles.

Previous comparative genomic studies have reported varying degrees of divergence in the genomic profiles of primary CRC and matched metastases. Substantial mutational divergence between paired primary tumor and metastasis samples was reported in a whole exome sequencing analysis with a low sequencing depth [34] and in shallow-targeted sequencing analysis variants with a variant allele frequency greater than 20% [35]. By contrast, a recent deep-targeted sequencing analysis [36, 37] and whole exome sequencing study with an average sequencing depth of > 100-fold [38] showed a high degree of mutational concordance, generally 50–80%, between primary CRC and metastases. These studies also suggest an exceedingly high level of concordance between primary CRC and metastases in genetic alterations that occur early in colorectal carcinogenesis, such as alterations in APC, KRAS, NRAS, and BRAF. In our study, we achieved an average coverage depth of 124-fold, enabling the sensitive detection of SNVs and indels to a 10% variant frequency. Our results are comparable to the results of previous high-resolution sequencing analyses. Therefore, it is more accurate to state that the majority of mutations are shared between primary CRC and metastatic tumors in at least some cases. Furthermore, our finding confirms our previous results [15] showing that both primary tumors and their subsequent metastases could be valid sources of DNA for patient selection before commencing anti-EGFR therapy.

Brannon et al. demonstrated that some of the mutational discordance could be explained by spatial heterogeneity in the patient tumors. They also suggest chemotherapy-derived clonal selection based on the finding that pre-treated tumors were less likely to have unique mutations than chemo-naive tumors in patients with concurrently resected tumors [36]. The increase in unique mutations after FOLFOX chemotherapy in our cases is also likely explained by intra-tumor heterogeneity and clonal selection rather than by the mutagenicity of L-OHP. According to the Big Bang model of human colorectal tumor growth [39], numerous heterogeneous subclones predominantly expand through tumor evolution. After treatment-derived clonal selection via methods such as surgery and cytotoxic chemotherapy, clones acquiring a fitness advantage in the environment become dominant in residual tumors. Our findings raise the possibility that surgery and adjuvant FOLFOX chemotherapy change the clonal composition, resulting in differences in genetic alterations before and after FOLFOX administration. This hypothesis may explain the mutational discordance in case 2. The AKT1 E17K mutation can activate the PI3K/AKT/mTOR pathway in a similar manner as PIK3CA alterations and promotes carcinogenesis through increasing cell proliferation or survival [40] via the PI3K pathway. AKT1 mutant clones, which were found in the minority of primary tumor samples, might sufficiently increase via clonal selection to detectable levels following surgery and adjuvant chemotherapy. The high rate of unique mutations in post-FOLFOX metastatic samples (54% in case 2) also indicates that there are differences in the clonal composition of these tumors compared to those of primary tumor and pre-FOLFOX metastatic samples. However, because of the very limited number of cases and the heterogeneity of the tumors, further validation is necessary. Furthermore, in a recent report, Angelova M et al. demonstrated that the immune system also influences tumor clonal composition and tumor evolution during the metastatic process [41]. In future studies, we must focus on clonal selection induced by not only chemotherapeutic agents but also by the immune system.

In addition, we reported that some genes that might be related to drug resistance were gained or lost during chemotherapy. Previous studies have indicated that CIAPIN1 is involved in the development of multiple drug resistance (MDR) [42]. PTPRJ is expressed in CRC cells, and it is reported that the sustained inhibition of PTPRJ increased cell resistance to 5-fluorouracil (5-FU)-induced apoptosis [43]. Wang et al. indicated that the TopBP1 expression level was related to the prognosis of non-small cell lung cancer patients treated with platinum-based chemotherapy [44]. The CIAPIN1 R132W and PTPRJ L738 V mutations were identified among the post-FOLFOX unique mutations in case 2. Conversely, in case 2, the TOPBP1 S630 L mutation was identified in primary and pre-FOLFOX metastasis but not in post-FOLFOX metastasis. These mutations were non-synonymous and were predicted as functionally important variants by ANNOVAR [19] in addition to Polyphen2 [21]. Thus, these mutations may affect the function of proteins related to drug resistance. However, the true function of these mutations is unclear, and further investigation is necessary.

The results of the GO analysis showed an enrichment of genes that are annotated as “calcium ion transport” among genes uniquely mutated after FOLFOX. For example, TRPM7 encodes a calcium-permeant ion channel that is notable for its inherent serine/threonine kinase activity [45]. The product of CACNA1E is the alpha 1E subunit of R-type voltage-dependent calcium channels [46]. The calcium-selective channel encoded by TRPV5 is activated by a low internal calcium level [47]. SLC24A6 encodes a family of potassium-dependent sodium/calcium exchangers that maintain cellular calcium homeostasis [48]. Previous in vitro studies have indicated that there is a relationship between intracellular Ca2+ homeostasis and the P-glycoprotein-dependent MDR phenotype [49], which is considered one of the mechanisms underlying resistance to L-OHP [50]. Although further investigation of the true function of these mutations is necessary in a larger cohort, alterations in genes involved in intracellular Ca2+ homeostasis may be related to resistance to FOLFOX therapy with the development of MDR.


In conclusion, our data showed that the mutation rates, mutation spectra, or CNVs were nearly identical between the primary tumor, pre-FOLFOX metastatic, and post-FOLFOX metastatic samples in four CRC cases. We found that some gene mutations that might be related to FOLFOX resistance were gained or lost during chemotherapy, and the inter-tumor discordance of the mutational profiles suggests the existence of intra-tumor heterogeneity and the induction of clonal selection as a result of response to the FOLFOX chemotherapy. Our findings might be useful as pilot data for a larger study to clarify the changes in the mutational landscape induced by adjuvant FOLFOX chemotherapy.





Burrows-Wheeler Aligner


Copy number variant


Colorectal cancer


Computed tomography


Database for Annotation, Visualization, and Integrated Discovery


Epidermal growth factor receptor


Formalin-fixed, paraffin-embedded


Genome Analysis ToolKit


Gene Ontology


Small insertions/deletions




Multi-drug resistance




Single nucleotide variant


  1. DeSantis CE, Lin CC, Mariotto AB, Siegel RL, Stein KD, Kramer JL, et al. Cancer treatment and survivorship statistics, 2014. CA Cancer J Clin. 2014;64:252–71.

    Article  Google Scholar 

  2. Douillard JY, Oliner KS, Siena S, Tabernero J, Burkes R, Barugel M, et al. Panitumumab-FOLFOX4 treatment and RAS mutations in colorectal cancer. N Engl J Med. 2013;369:1023–34.

    Article  CAS  Google Scholar 

  3. Van Cutsem E, Lenz HJ, Kohne CH, Heinemann V, Tejpar S, Melezinek I, et al. Fluorouracil, leucovorin, and irinotecan plus cetuximab treatment and RAS mutations in colorectal cancer. J Clin Oncol. 2015;33:692–700.

    Article  Google Scholar 

  4. Loupakis F, Cremolini C, Masi G, Lonardi S, Zagonel V, Salvatore L, et al. Initial therapy with FOLFOXIRI and bevacizumab for metastatic colorectal cancer. N Engl J Med. 2014;371:1609–18.

    Article  Google Scholar 

  5. Yuan ZX, Wang XY, Qin QY, Chen DF, Zhong QH, Wang L, et al. The prognostic role of BRAF mutation in metastatic colorectal cancer receiving anti-EGFR monoclonal antibodies: a meta-analysis. PLoS One. 2013;8:e65995.

    Article  CAS  Google Scholar 

  6. The Cancer Genome Atlas Network. Comprehensive molecular characterization of human colon and rectal cancer. Nature. 2012;487:330–7.

    Article  Google Scholar 

  7. Guinney J, Dienstmann R, Wang X, de Reynies A, Schlicker A, Soneson C, et al. The consensus molecular subtypes of colorectal cancer. Nat Med. 2015;21:1350–6.

    Article  CAS  Google Scholar 

  8. Ding L, Ley TJ, Larson DE, Miller CA, Koboldt DC, Welch JS, et al. Clonal evolution in relapsed acute myeloid leukaemia revealed by whole-genome sequencing. Nature. 2012;481:506–10.

    Article  CAS  Google Scholar 

  9. Johnson BE, Mazor T, Hong C, Barnes M, Aihara K, McLean CY, et al. Mutational analysis reveals the origin and therapy-driven evolution of recurrent glioma. Science. 2014;343:189–93.

    Article  CAS  Google Scholar 

  10. Benson AB 3rd, Venook AP, Bekaii-Saab T, Chan E, Chen YJ, Cooper HS, et al. Colon cancer, version 3.2014. J Natl Compr Cancer Netw. 2014;12:1028–59.

    Article  CAS  Google Scholar 

  11. Woynarowski JM, Faivre S, Herzig MC, Arnett B, Chapman WG, Trevino AV, et al. Oxaliplatin-induced damage of cellular DNA. Mol Pharmacol. 2000;58:920–7.

    Article  CAS  Google Scholar 

  12. Hah SS, Sumbad RA, de Vere White RW, Turteltaub KW, Henderson PT. Characterization of oxaliplatin-DNA adduct formation in DNA and differentiation of cancer cell drug sensitivity at microdose concentrations. Chem Res Toxicol. 2007;20:1745–51.

    Article  CAS  Google Scholar 

  13. Sharma S, Gong P, Temple B, Bhattacharyya D, Dokholyan NV, Chaney SG. Molecular dynamic simulations of cisplatin- and oxaliplatin-d (GG) intrastand cross-links reveal differences in their conformational dynamics. J Mol Biol. 2007;373:1123–40.

    Article  CAS  Google Scholar 

  14. Silva MJ, Costa P, Dias A, Valente M, Louro H, Boavida MG. Comparative analysis of the mutagenic activity of oxaliplatin and cisplatin in the Hprt gene of CHO cells. Environ Mo Mutagen. 2005;46:104–15.

    Article  CAS  Google Scholar 

  15. Kawamoto Y, Tsuchihara K, Yoshino T, Ogasawara N, Kojima M, Takahashi M, et al. KRAS mutations in primary tumours and post-FOLFOX metastatic lesions in cases of colorectal cancer. Br J Cancer. 2012;107:340–4.

    Article  CAS  Google Scholar 

  16. Bando H, Yoshino T, Tsuchihara K, Ogasawara N, Fuse N, Kojima T, et al. KRAS mutations detected by the amplification refractory mutation system-scorpion assays strongly correlate with therapeutic effect of cetuximab. Br J Cancer. 2011;105:403–6.

    Article  CAS  Google Scholar 

  17. Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009;25:1754–60.

    Article  CAS  Google Scholar 

  18. DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011;43:491–8.

    Article  CAS  Google Scholar 

  19. Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38:e164.

    Article  Google Scholar 

  20. Huang d W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.

    Article  CAS  Google Scholar 

  21. Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7:248–9.

    Article  CAS  Google Scholar 

  22. Sathirapongsasuti JF, Lee H, Horst BA, Brunner G, Cochran AJ, Binder S, et al. Exome sequencing-based copy-number variation and loss of heterozygosity detection: ExomeCNV. Bioinformatics. 2011;27:2648–54.

    Article  CAS  Google Scholar 

  23. Allegra CJ, Yothers G, O'Connell MJ, Sharif S, Colangelo LH, Lopa SH, et al. Initial safety report of NSABP C-08: a randomized phase III study of modified FOLFOX6 with or without bevacizumab for the adjuvant treatment of patients with stage II or III colon cancer. J Clin Oncol. 2009;27:3385–90.

    Article  CAS  Google Scholar 

  24. Wood LD, Parsons DW, Jones S, Lin J, Sjoblom T, Leary RJ, et al. The genomic landscapes of human breast and colorectal cancers. Science. 2007;318:1108–13.

    Article  CAS  Google Scholar 

  25. Meier B, Cooke SL, Weiss J, Bailly AP, Alexandrov LB, Marshall J, et al. C. Elegans whole-genome sequencing reveals mutational signatures related to carcinogens and DNA repair deficiency. Genome Res. 2014;24:1624–36.

    Article  CAS  Google Scholar 

  26. Vogelstein B, Papadopoulos N, Velculescu VE, Zhou S, Diaz LA Jr, Kinzler KW. Cancer genome landscapes. Science. 2013;339:1546–58.

    Article  CAS  Google Scholar 

  27. Choi M, Scholl UI, Ji W, Liu T, Tikhonova IR, Zumbo P, et al. Genetic diagnosis by whole exome capture and massively parallel DNA sequencing. Proc Natl Acad Sci U S A. 2009;106:19096–101.

    Article  CAS  Google Scholar 

  28. Ehrsson H, Wallin I, Yachnin J. Pharmacokinetics of oxaliplatin in humans. Med Oncol. 2002;19:261–5.

    Article  CAS  Google Scholar 

  29. Casazza A, Finisguerra V, Capparuccia L, Camperi A, Swiercz JM, Rizzolio S, et al. Sema3E-Plexin D1 signaling drives human cancer cell invasiveness and metastatic spreading in mice. J Clin Invest. 2010;120:2684–98.

    Article  CAS  Google Scholar 

  30. Kim MJ, Lee HS, Kim JH, Kim YJ, Kwon JH, Lee JO, et al. Different metastatic pattern according to the KRAS mutational status and site-specific discordance of KRAS status in patients with colorectal cancer. BMC Cancer. 2012;12:347.

    Article  CAS  Google Scholar 

  31. Baas JM, Krens LL, Guchelaar HJ, Morreau H, Gelderblom H. Concordance of predictive markers for EGFR inhibitors in primary tumors and metastases in colorectal cancer: a review. Oncologist. 2011;16:1239–49.

    Article  CAS  Google Scholar 

  32. Kopetz S, Overman MJ, Chen K, Lucio-Eterovic AK, Fogelman BKKDR, Dasari A, et al. Mutation and copy number discordance in primary versus metastatic colorectal cancer (mCRC). J Clin Oncol. 2014;32(Suppl 5s):3509.

    Article  Google Scholar 

  33. Diehl F, Schmidt K, Choti MA, Romans K, Goodman S, Li M, et al. Circulating mutant DNA to assess tumor dynamics. Nat Med. 2008;14:985–90.

    Article  CAS  Google Scholar 

  34. Lee SY, Haq F, Kim D, Jun C, Jo HJ, Ahn SM, et al. Comparative genomic analysis of primary and synchronous metastatic colorectal cancers. PLoS One. 2014;9:e90459.

    Article  Google Scholar 

  35. Vermaat JS, Nijman IJ, Koudijs MJ, Gerritse FL, Scherer SJ, Mokry M, et al. Primary colorectal cancers and their subsequent hepatic metastases are genetically different: implications for selection of patients for targeted treatment. Clin Cancer Res. 2012;18:688–99.

    Article  CAS  Google Scholar 

  36. Brannon AR, Vakiani E, Sylvester BE, Scott SN, McDermott G, Shah RH, et al. Comparative sequencing analysis reveals high genomic concordance between matched primary and metastatic colorectal cancer lesions. Genome Biol. 2014;15:454–7.

    Article  Google Scholar 

  37. Tan IB, Malik S, Ramnarayanan K, McPherson JR, Ho DL, Suzuki Y, et al. High-depth sequencing of over 750 genes supports linear progression of primary tumors and metastases in most patients with liver-limited metastatic colorectal cancer. Genome Biol. 2015;16:32.

    Article  Google Scholar 

  38. Lim B, Mun J, Kim JH, Kim CW, Roh SA, Cho DH, et al. Genome-wide mutation profiles of colorectal tumors and associated liver metastases at the exome and transcriptome levels. Oncotarget. 2015;6:22179–90.

    PubMed  PubMed Central  Google Scholar 

  39. Sottoriva A, Kang H, Ma Z, Graham TA. A big bang model of human colorectal tumor growth. Nat Genet. 2015;47:209–16.

    Article  CAS  Google Scholar 

  40. Cully M, You H, Levine AJ, Mak TW. Beyond PTEN mutations: the PI3K pathway as an integrator of multiple inputs during tumorigenesis. Nat Rev Cancer. 2006;6:184–92.

    Article  CAS  Google Scholar 

  41. Angelova M, Mlecnik B, Vasaturo A, Bindea G, Fredriksen T, Lafontaine L, Buttard B, Morgand E, Bruni D, Jouret-Mourin A, Hubert C, Kartheuser A, Humblet Y, Ceccarelli M, Syed N, Marincola FM, Bedognetti D, Van den Eynde M, Galon J. Evolution of metastases in space and time under immune selection. Cell. 2018;18:751–65.

    Article  Google Scholar 

  42. Zhang YF, Li XH, Shi YQ, Wu YY, Li N, He Q, et al. CIAPIN1 confers multidrug resistance through up-regulation of MDR-1 and Bcl-L in LoVo/Adr cells and is independent of p53. Oncol Rep. 2011;25:1091–8.

    CAS  PubMed  Google Scholar 

  43. Yan CM, Zhao YL, Cai HY, Miao GY, Ma W. Blockage of PTPRJ promotes cell growth and resistance to 5-FU through activation of JAK1/STAT3 in the cervical carcinoma cell line C33A. Oncol Rep. 2015;33:1737–44.

    Article  CAS  Google Scholar 

  44. Wang LR, He LJ, Wang Y, Li YY, Lou Y, Zhang GB, Li Y. Chen. Correlation between BRCA1 and TopBP1 protein expression and clinical outcome of non-small cell lung cancer treated with platinum-based chemotherapy. J Cancer Chemother Pharmacol. 2015;76:163–70.

    Article  CAS  Google Scholar 

  45. Runnels LW, Yue L, Clapham DE. TRP-PLIK, a bifunctional protein with kinase and ion channel activities. Science. 2001;291:1043–7.

    Article  CAS  Google Scholar 

  46. Vajna R, Schramm M, Pereverzev A, Arnhold S, Grabsch H, Klockner U, et al. New isoform of the neuronal Ca2+ channel alpha1E subunit in islets of Langerhans and kidney--distribution of voltage-gated Ca2+ channel alpha1 subunits in cell lines and tissues. Eur J Biochem. 1998;257:274–85.

    Article  CAS  Google Scholar 

  47. Hoenderop JG, van der Kemp AW, Hartog A, van de Graaf SF, van Os CH, Willems PH, et al. Molecular identification of the apical Ca2+ channel in 1, 25-dihydroxyvitamin D3-responsive epithelia. J Biol Chem. 1999;274:8375–8.

    Article  CAS  Google Scholar 

  48. Cai X, Lytton J. The cation/ca(2+) exchanger superfamily: phylogenetic analysis and structural implications. Mol Biol Evol. 2004;21:1692–703.

    Article  CAS  Google Scholar 

  49. Sulova Z, Seres M, Barancik M, Gibalova L, Uhrik B, Polekova L, et al. Does any relationship exist between P-glycoprotein-mediated multidrug resistance and intracellular calcium homeostasis. Gen Physiol Biophys 2009;28 Spec No Focus:F89–F95.

  50. Martinez-Balibrea E, Martinez-Cardus A, Gines A, Ruiz de Porras V, Moutinho C, Layos L, et al. Tumor related molecular mechanisms of Oxaliplatin resistance. Mol Cancer Ther. 2015;14:1767–76.

    Article  CAS  Google Scholar 

Download references


The National Cancer Center Biobank is supported by the National Cancer Center Research and Development Fund, Japan.


This study was supported by research funding from the National Cancer Center. The funding body helped in material collection and data analyses.

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article and its additional files.

Author information

Authors and Affiliations



KT and WO conceived the study with advice from SM and TY. KH, WO, SM, and RY performed the data analysis. YK (Kawamoto), HB, TY, and AO provided clinical data and helped collect the tumor tissues. KH and WO wrote the manuscript with advice from SM, YK (Kawamoto), HB, RY, SY, TY, YK (Komatsu), AO, NS, and KT. All authors have read and approved the manuscript.

Corresponding author

Correspondence to Wataru Okamoto.

Ethics declarations

Ethics approval and consent to participate

The study protocol was approved by the Institutional Review Board of the National Cancer Center. Samples of the four patients admitted before June 2011 and whose written informed consents could not be obtained were used with the approval of the Institutional Review Board, which privileged the opt-out for the patients, according to the Ethical Guidelines for Epidemiological Research issued by Ministry of Education, Culture, Sports, Science and Technology, and Ministry of Health, Labour and Welfare of Japan.

Consent for publication

Informed consent to publish identifying patient details was obtained from all participants included in this study.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional files

Additional file 1:

Figure S1. Clinical courses of the four cases. The blue arrows indicate the day of surgery, and the red arrows indicate the day on which the recurrent tumors were diagnosed during or after adjuvant FOLFOX therapy. The arrowheads (light blue) indicate the number of FOLFOX treatments. Disease-free survival (DFS) is calculated from the time of the final operation until post-FOLFOX recurrence (TIF 259 kb)

Additional file 2:

Table S1. Coverage details. Summary of the sequencing coverage depth in each sample. (XLSX 10 kb)

Additional file 3:

Figure S2. Comparison of the mutation spectrum among unique mutations. Relative mutation frequencies are shown for unique mutations in primary, pre-FOLFOX metastatic, and post-FOLFOX metastatic tumor samples. The presented number of mutations is the sum of data recorded in the four patients. There were no significant differences in all mutation fractions (Wilcoxon signed-rank test). (TIF 119 kb)

Additional file 4:

Table S2. Focal amplification in case 3. Summary of the gene amplification in case 3. (XLSX 10 kb)

Additional file 5:

Table S3. Detailed information for all detected gene mutations. List of all gene mutations detected in our analysis. (XLSX 86 kb)

Additional file 6:

Table S4. The gain or loss of the mutations that were predicted as functionally important. (XLSX 13 kb)

Additional file 7:

Table S5. Post-FOLFOX unique gene mutations used for the GO analysis. List of selected gene mutations for the GO analysis. (XLSX 16 kb)

Additional file 8:

Table S6. Significant functional clusters in each mutation category. List of the Gene Ontology terms enriched among unique gene mutations before FOLFOX administration or among common mutations. (XLSX 10 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Harada, K., Okamoto, W., Mimaki, S. et al. Comparative sequence analysis of patient-matched primary colorectal cancer, metastatic, and recurrent metastatic tumors after adjuvant FOLFOX chemotherapy. BMC Cancer 19, 255 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: