Research | Open | Published:
Spatio-temporal mutation profiles of case-matched colorectal carcinomas and their metastases reveal unique de novo mutations in metachronous lung metastases by targeted next generation sequencing
Molecular Cancervolume 15, Article number: 63 (2016)
Targeted next generation sequencing (tNGS) has become part of molecular pathology diagnostics for determining RAS mutation status in colorectal cancer (CRC) patients as predictive tool for decision on EGFR-targeted therapy. Here, we investigated mutation profiles of case-matched tissue specimens throughout the disease course of CRC, to further specify RAS-status dynamics and to identify de novo mutations associated with distant metastases.
Case-matched formalin-fixed and paraffin-embedded (FFPE) resection specimens (n = 70; primary tumours, synchronous and/or metachronous liver and/or lung metastases) of 14 CRC cases were subjected to microdissection of normal colonic epithelial, primary and metastatic tumour cells, their DNA extraction and an adapted library protocol for limited DNA using the 48 gene TruSeq Amplicon Cancer PanelTM, MiSeq sequencing and data analyses (Illumina).
By tNGS primary tumours were RAS wildtype in 5/14 and mutated in 9/14 (8/9 KRAS exon 2; 1/9 NRAS Exon 3) of cases. RAS mutation status was maintained in case-matched metastases throughout the disease course, albeit with altered allele frequencies. Case-matched analyses further identified a maximum of three sequence variants (mainly in APC, KRAS, NRAS, TP53) shared by all tumour specimens throughout the disease course per individual case. In addition, further case-matched de novo mutations were detected in synchronous and/or metachronous liver and/or lung metastases (e.g. in APC, ATM, FBXW7, FGFR3, GNAQ, KIT, PIK3CA, PTEN, SMAD4, SMO, STK11, TP53, VHL). Moreover, several de novo mutations were more frequent in synchronous (e.g. ATM, KIT, PIK3CA, SMAD4) or metachronous (e.g. FBXW7, SMO, STK11) lung metastases. Finally, some de novo mutations occurred only in metachronous lung metastases (CDKN2A, FGFR2, GNAS, JAK3, SRC).
Together, this study employs an adapted FFPE-based tNGS approach to confirm conservation of RAS mutation status in primary and metastatic tissue specimens of CRC patients. Moreover, it identifies genes preferentially mutated de novo in late disease stages of metachronous CRC lung metastases, several of which might be actionable by targeted therapies.
Predictive molecular pathology mutation testing in selected genes is a routine diagnostic application in several epithelial tumour entities. This includes extended RAS testing (KRAS and NRAS exons 2,3,4) in advanced colorectal cancer (CRC) patients, which functions as well-established predictive biomarker for resistance to EGFR-targeted therapy (e.g. cetuximab or panitumumab) [1–5]. In addition, BRAF mutation testing is recommended for molecular grading of undifferentiated CRCs or as supportive tool for further molecular classification of microsatellite-instable CRCs .
Whilst dideoxy-Sanger sequencing, pyrosequencing or quantitative polymerase chain reaction-based assays were routinely used in the past, targeted next generation sequencing (tNGS) has by now evolved as a robust, time- and cost-efficient technique to accommodate the growing demand for mutation profiling in molecular pathology laboratories. Indeed, there are increasing numbers of reports on the general reliability and applicability of tNGS for mutation analyses of formalin-fixed and paraffin-embedded (FFPE) tissue specimens and/or the generation of tumour-entity specific tNGS gene panels, using different platforms and library preparation approaches [7–12]. In addition, several investigators compared primary colon and/or rectal tumours and metastases by NGS approaches in cohorts of 13 cases , 15 cases , 18 cases , 20 cases , 24 cases , 34 cases  or >400 cases . These studies were based on fresh-frozen tissue specimens  or assessment of FFPE tissue specimens by amplicon-based semi-conductor NGS technology [13, 16, 17, 19], or evaluation of fresh-frozen tissue specimens with >70 % tumour cell content by whole exome sequencing [14, 18]. Reference to the otherwise mostly small tissue specimens and limited tumour cells of CRC liver or lung metastases in NGS performance and data interpretation is still sparse. Moreover, few of the previous studies focused on the individual clinico-pathological and molecular characteristics of the investigated cases. One study  performed tNGS data analysis of 13 matched pairs of the primary tumour plus each one liver metastasis, reporting a 78 % match of mutations. In another study using semi-conductor NGS technology, mutations in for example APC, KRAS, FBXW7, PIK3CA, BRAF, SMAD4 were concordant between primary CRCs and their metastases, whereas 4 of 24 cases also showed de novo mutations in SYNE1, CTNNB1, TP53 and PTEN . Similar findings were found in another study, which examined 17 paired primary and mainly synchronous metastatic CRC tissue specimens in a >400 CRC cohort . In two whole exome sequencing based studies [14, 18], mutation profiles of key CRC associated genes were similar in primary and metastatic tumours by 53 %  and 57 % , with additional de novo mutations in metastases of 47 %  and 43 % . No correlations of case-matched tumour specimens to clinico-pathological and molecular characteristics were made.
Taken together, there is hence still a lack of studies which address the question of mutation profiles of CRC tumour progression based on case-matched tissue specimen analyses by another NGS platform using limited numbers of microdissected tumour cells, and case-matched data analyses reflecting the individual clinico-pathological and molecular characteristics of CRCs.
In this study, we therefore established a sequencing-by-synthesis targeted NGS approach for analysis of 48 genes in limited numbers of tumour cells microdissected from different types of FFPE tissue specimens. We applied this approach to study mutation profiles of 14 CRC cases with 70 matched tissue specimens from primary tumours and corresponding synchronous and metachronous liver and/or lung metastases. We took into account clinical (treatment regimens, disease course, e.g. synchronous and metachronous metastases), pathological (e.g. primary tumour location, histotype) and molecular (e.g. microsatellite-status, CpG island methylator phenotype) parameters. By this approach, we confirm conservation of predictive RAS mutation status throughout the disease course and identify de novo mutations associated with either synchronous or metachronous lung metastases.
Patients and tissue specimens
For establishing the novel tNGS approach, matched Fresh-Frozen (FF) and Formalin-fixed and Paraffin-embedded (FFPE) tissue specimens of 3 non-small cell lung, breast and colorectal carcinomas were selected as technical testing cohort. For further validation, FFPE tissue specimens of two breast cancers (n = 4 histologic lesions) and two CRC cases (n = 7, including normal epithelium, primary tumour, liver and lung metastases) were selected as technical validation cohort (Additional file 1: Text S1, Table S1).
For analysis of case-matched disease course tissue specimens of CRC patients, FFPE tissue specimens were obtained from 14 CRC patients with synchronous and/or metachronous liver and/or lung metastases (total 70 tissue specimens) undergoing surgery at the Medical Center - University of Freiburg, Germany. In case of multiple lesions at one surgically addressed disease time point (e.g. two lesions in distinct liver segments or two lesions of lung metastases in the right and left lung), these were analyzed individually. Processed tissue specimens hence included normal colonic epithelium dissected from the tumour-free resection margins, primary tumours derived from central tumour areas of resection specimens as well as metastatic liver or lung lesions of resection specimens. All tissue specimens were re-classified and marked by a qualified pathologist (LL) for microdissection and tumour cell content on newly prepared H&E sections serial to those for microdissection.
Definition of synchronous versus metachronous metastases was according to resection with > 6 months post primary tumour resection. Clinico-pathological parameters - including treatment regimens - of the cases are given in Fig. 1 and Table 1. Note that rectal carcinomas receiving neoadjuvant chemotherapy prior to resection of the primary tumors (4/5 cases) are hence classified with “ypTNM”. Note also that different treatment regimens were given to each case throughout the individual complex disease courses. Therefore, unless otherwise specified, data was analyzed in truly case-matched approach addressing each individual patient.
The use of tissue specimens had been approved by the local ethics Institution (#251/04, #66/07 and #191/13, Ethik-Kommission, University of Freiburg, Germany).
Microdissection, DNA isolation and DNA quality control
For fresh-frozen tissue specimens, 10 to 20 μm serial sections were cut from each tissue specimen for manual microdissection of normal or tumour cells under morphological control, respective microscope using fine needles. DNA isolation was with the QIAamp Micro Kit, according to the manufacturer’s protocol (Qiagen, Hilden, Germany).
For FFPE tissue specimens, 5 to 10 μm serial sections were cut from each tissue specimen for manual microdissection of normal or tumour cells under morphological control, respective microscope using fine needles. From the testing and validation tissue specimens, DNA isolation was with the QIAamp FFPE DNA Kit, according to the manufacturer`s protocol (Qiagen, Hilden, Germany). For the disease course CRC tissue specimens, DNA was isolated using the Allprep DNA/RNA FFPE Kit according to the manufacturer’s instructions (Qiagen, Hilden, Germany).
For each sample, DNA concentration and purity was determined using the Nanodrop-8000 (ND1000, Peqlab, Erlangen, Germany) and the DNA quality was measured in triplicate using the “FFPE QC Kit” according to the manufacturer`s protocols (Illumina, San Diego, USA). For the latter, the average difference of Cycle-treshold (Ct) values between the DNA samples and a positive DNA control were calculated as relative “QC value”, with “<2” suggested as a valid DNA sample for library preparation by the manufacturer.
MSI and CIMP analyses of CRC tissue specimens
For MSI testing, a previously published multiplex PCR protocol was used with slight modifications . The CpG island methylator phenotype (CIMP) status was determined by pyrosequencing addressing methylation of five CIMP-related genes (RUNX3, CACNA1G, EPM2AIP1/MLH1, NEUROG1,CRABP1) with modifications to previously published protocols [21–23].
Library preparation and sequencing
Initially, libraries of the technical testing cohort were prepared with the TruSeq Amplicon Cancer Panel (TSACP), covering 225 amplicons of 48 cancer associated genes, following exactly the manufacturer’s protocol (Illumina, San Diego, USA). Since this did not yield reliable libraries for the different types of FFPE tissue specimens (Additional file 1: Table S1A), modifications to the library protocol were made to the first hybridization step with oligo pools and to subsequent PCR cycling parameters (see Additional file 1: Table S1B). The comparison of the performance of these different library protocols is presented for the technical testing and validation cohorts (Additional file 1: Figure S1).
The selection of specific library protocols according to DNA amount and quality for the 70 DNA samples of the CRC study cohort is provided together with the tumour cell content in Additional file 1: Table S2).
Libraries were examined for quantity and quality using the Bioanalyzer 2100 system (Agilent Technologies, Santa Clara, USA) and subjected to clean up, normalization and pooling according to the manufacturer (Illumina, San Diego, USA). MiSeq sequencing runs were performed with each 16–20 pooled libraries using the MiSeq Reagent Kit v2 (300 cycle) (Illumina, San Diego, USA) and paired-end sequencing with 2 × 151 bp.
Targeted NGS data analysis and validation of detected variants
Established software tools associated with MiSeq technology were used, including Real Time Analysis, MiSeq Reporter analysis and VariantStudio (all Illumina, San Diego, USA). For variant filtering, a 10 % cut-off of allele frequency, with exclusion of synonymous variants and those at a population frequency >5 % (i.e. common sequence variants in KDR p.Q472H, KIT p.M541L and TP53 p.P72R are not reported in text, tables and figures) were used. The Integrative Genomics Viewer (IGV; Broad Institute, UK) [24, 25] was used to visualize read alignments and presence of variants against the reference genome (hg39).
Data sets of normal colonic epithelium were subtracted from those of case-matched tumour (primary, metastases), yielding somatic sequence variants exclusively. De novo mutations in metastases were obtained by subtracting the primary tumour data sets from those of case-matched metastases. All detected sequence variants are provided in Additional file 2: Table S4.
Validation of selected somatic variants was done for 27 DNA samples of 8 CRC cases by dideoxy sequencing according to SOP-driven molecular pathology protocols or for TP53 as described previously  (Additional file 1: Table S3). All primer sequences are available upon request.
Establishment of library preparation from DNA samples derived from different origin tissue specimens
We first tested the standard tNGS library protocol for different types of matched fresh-frozen (FF) and FFPE tissue specimens from a technical testing cohort of 3 pairs of non-small-cell lung (NSCLC), colorectal (CRC) and breast carcinomas (Additional file 1: Figure S1A, Table S1A). To allow tNGS analysis also from low input and poorer DNA quality, we adapted the library protocol (oligo hybridization, PCR cycling conditions) using DNA samples derived from a technical validation cohort of FFPE tissue specimens of two CRC cases (n = 7, including normal epithelium, primary tumour, liver or lung metastases) and two breast carcinomas (n = 4, different histologic lesions) (Additional file 1: Figure S1B, Table S1B).
Thereby, modification of the library protocol allowed analysis of different types of FFPE tissue specimens, including those with limited DNA amounts from few microdissected tumour cells. Hence, this also enables analysis of limited viable tumour cells from CRC tissue specimens post neoadjuvant therapy or from liver metastases with marked necrosis.
Targeted NGS of CRC cases with longitudinal disease course tissue specimens
We next applied the new tNGS approach to profile the tissue specimens obtained throughout the individual disease course of 14 CRC patients, each case presenting with specific clinico-pathological parameters, which were included in tNGS data evaluation and interpretation. From each case, normal epithelium, primary tumour and at least two distant metastases were selected for mutation profiling, resulting in n = 70 DNA samples in total (Fig. 1, Table 1, Additional file 1: Table S2).
As expected from the size and histology of the lesions of primary tumour and metastatic sites (Fig. 2a), the DNA samples derived from the microdissected tumour cells showed a wide range of DNA concentrations (5 to 452 ng/μl, reflecting the maximum DNA of enriched tumour cells from the partially small tissue specimens of metastatic lesions; mean tumour cell content of enriched samples = 69 %, with 4 DNA samples from cases #1, #6, #13 showing only 10–30 % tumour cell enrichment; Additional file 1: Table S2) and quality scores (QC range 0.1 to 3, with a QC value >2 in 24/70 DNA samples), therefore requiring adjusted tNGS library preparation according to the newly established protocols (Additional file 1: Table S2). Thereby, all 70 DNA samples yielded sufficient libraries, resulting in mean amplicon coverage of 5767x (range 1003x to 9517x) upon sequencing (Fig. 2b).
With a cut-off of 10 % allele frequency for variant calling, the most frequently mutated genes, respective sequence variants occurring in at least 2 tissue specimens of primary tumour and/or liver and/or lung metastases were: TP53 (10/14 cases), APC (8/14 cases), KRAS (8/14 cases), SMAD4 (2/14 cases) as well as FGFR3 (1/14 cases) and PIK3CA (1/14 cases) and NRAS (1/14 cases) (Table 2). See also Additional file 2: Table S4 for all detected sequence vairants.
RAS mutation status is maintained in case-matched primary CRCs and their metastatic sites
We next examined the RAS status throughout the disease course of the individual CRC cases.
Of the 14 cases, 8/14 (57 %) exhibited KRAS (exon 2, codon 12 or 13) mutations in the primary tumour and 1/14 (7 %) a NRAS (exon 3, codon 61) mutation (Table 3). These RAS mutations were detectable throughout the disease course in each resected tumour lesion in 7/9 (78 %) mutated cases. Thereby, the allele frequency of the RAS mutations varied between primary tumour, liver and/or lung metastases. In two cases, the allele frequencies were below 10 % in the synchronous liver metastasis (9.7 %, case #5) and in one of the two synchronous lung metastastic lesions (5.6 %, case #13). In cases #3 and #6, the KRAS p.G12V mutation was not detected in the metachronous lung and liver metastasis even at a cut-off level of >1 % allele frequency. The drop of allele frequency in mainly small metastatic lesions most likely relates to the limited or only single distributed (e.g. case #6 liver metastasis, Fig. 2a) tumour cells available for microdissection.
Finally, none of the 5 cases with wildtype RAS in the primary tumour subsequently developed RAS mutations in the EGFR-targeted therapy relevant “hotspot” codons (Table 3).
Mutation profiles of CRC patients throughout the disease course
To determine mutation profiles of tumour sites, respective locations of disease progression, we first examined the NGS data sets of all 70 DNA samples in a case-mixed manner.
This revealed tumour site associated shared sequence variants in 1) APC, PDGFRA, PTEN for primary colon and/or rectal tumours, 2) ATM, TP53 for synchronous liver metastases, 3) FGFR3 and TP53 for metachronous liver metastases, 4) ABL1, APC, ATM, ERBB4 and SMAD4 for synchronous lung metastases, and 5) APC and TP53 for metachronous lung metastases (Table 4).
Hence, this case-mixed analysis suggests some gene sequence variants specific for liver or lung metastases, but gives no information on whether or not these were already present in e.g. the primary tumour of individual CRC patients.
Case-matched mutation profiles throughout disease course reveal distinct de novo mutations in synchronous versus metachronous lung metastases
Since the clinico-pathological parameters and disease course are individual for each CRC patient (e.g. neoadjuvant treatment in 4/5 rectal tumours as well as adjuvant treatment regimens), we next analyzed the case-matched mutation profiles, specifically focusing on the late(st) stages of lung metastasis.
As shown by the Venn diagrams for each CRC case in Fig. 3, generally few distinct mutations (maximum of 3, including APC, KRAS, NRAS, PIK3CA, SMAD4 and/or TP53) were shared between all case-matched DNA samples (primary tumour, synchronous and/or metachronous liver and/or lung metastases) of an individual CRC case (Table 2). Taking into account the individual clinico-pathological parameters, there was no association of mutation profiles to primary tumour site (left colon, right colon, rectum) or histology, tumour stage, grading or molecular subtype (Fig. 3).
However, interestingly numerous case-matched de novo mutations occurred in liver (cases #1, 2, 6, 7, 8, 10, 13) and lung (cases #1, 2, 3, 6, 7, 8, 9, 10, 11, 12) metastases, which were not present in the case-matched primary tumours (Fig. 3): Case-matched analysis and comparison of synchronous (n = 12 DNA samples of 11 CRC cases) versus metachronous (n = 6 DNA samples of 5 CRC cases) liver metastases yielded 4/12 (33 %) synchronous versus 4/6 (67 %) metachronous liver metastases with unique de novo mutations. Thereby, the de novo mutation profiles ranged from 1–90 and 2–20 variants in synchronous as compared to metachronous liver metastases, including mutations in FBXW7, FGFR3, GNAQ and PTEN occurring in more than 2 cases with synchronous liver metastases (Figs. 3 and 4a).
In contrast, case-matched analysis and comparison of synchronous (n = 8 of 5 CRC cases) and metachronous (n = 16 of 10 CRC cases) lung metastases revealed a broad mutation profile in both synchronous (range of variants: 21–90) and metachronous (range of variants: 1–98) lung metastases (Figs. 3 and 4b). Thereby, synchronous lung metastasis more frequently showed mutations in ATM, KIT, PIK3CA and SMAD4 (each with 3/5 cases as compared to 1/5 cases with metachronous liver metastases). Moreover, metachronous lung metastases showed exclusive mutations in CDKN2A, FGFR2, GNAS, JAK3 and SRC (each 1/10 cases as compared to 0/5 cases with synchronous lung metastases) as well as frequent mutations in FBXW7 (4/10 cases; synchronous: 2/5 cases), SMO (3/10 cases; synchronous: 1/5 cases) and STK11 (3/10 cases; synchronous: 1/5 cases).
As above, also the identified de novo mutations in lung metastases were not associated with clinico-pathological parameters. Thereby, neither the number of de novo mutations and/or the genes affected by de novo mutations were directly linked to specific prior treatment regimens, as summarized in Table 5.
Hence, these data allow case-matched and tumour site location specific identification of gene mutations involved in disease progression of individual CRC cases.
Multiple lesions of lung, but not liver metastases resected at one time point may differ in the mutation profile
Finally, de novo mutations in each two separate topographic lesions of liver or lung metastases of one resection in 3 CRC cases (#3, #13, #14) were examined (Fig. 5):
Exactly the same mutation profile (mutations in APC, NRAS, TP53) were detected in two lesions (liver segments 5 and 7) of the liver metastasis of case #14.
In contrast, in two lesions (upper, lower lobes) of the synchronous lung metastasis of case #13 shared mutations occurred in ABL1, APC, ATM and ERBB4, but this was accompanied by 39 (“lower lobe left”) and 90 (“upper lobe left”) additional variants for each separate lesion. Similarly, in two lesions (lung segments 7 and 8) of the metachronous lung metastasis of case #3 shared mutations were detected in APC (n = 2), FGFR2 (n = 1) and TP53 (n = 1), which was associated with 70 (“segment 7”) and 98 (“segment 8”) additional variants. Thereby, the different mutation profiles of each the two lung metastatic lesions per resection were seen for the mucinous (case #13) and the tubular/NOS (case #3) CRC case (Fig. 5).
This study comprehensively evaluated the feasibility and outcome of targeted next generation sequencing of case-matched colorectal cancer (CRC) Formalin-fixed and Paraffin-embedded tissue specimens in a close morphological context.
The study describes a novel protocol which allows targeted next generation sequencing (tNGS) even from low input and/or poorer quality DNA derived from microdissected tumour cells of archival tissue specimens using a sequencing-by-synthesis tNGS approach.
Moreover, the present study shows that mutation profiling of CRC, respective investigation of “driver mutations”, intratumoural heterogeneity and tumour evolution need to be performed in a fully case-matched (i.e. taking into account individual treatment regimens, including e.g. neoadjuvant therapy of rectal cancers prior to resection of the primary tumour) and morphology-controlled manner: The spectrum of apparent “driver mutations” of e.g. lung metastasis were different upon evaluation of all lung metastases in a case-mixed manner as compared to analysis of true “de novo” mutations of lung metastases in a case-matched manner. In addition, the present study reveals an obvious difference of the mutation profile in (case-matched) synchronous and metachronous liver or lung metastases, suggesting tumour heterogeneity and distinct tumour evolution. Hence, the investigative approach of the CRC cohort here focused on those mutations possibly accounting for tumour progression of individual patients more than 2 years post resection of the primary tumour.
The CRC cohort here only included 14 CRC patients, but a total of 70 tissue specimens obtained during case-matched follow-up, respective disease progression, especially resected lung metastases, were investigated in a case-matched manner also for tNGS data analysis and interpretation. Although this appears a very small number of cases, this is in fact similar in other CRC screening studies. For example, a large study of initially 2000 CRC patients, then selected 468 cases for NGS and of these only 19 cases had matched metastases, with 12/19 cases presenting distant metastases . In other studies, 13 cases with mainly liver metastases ; 24 cases with either one (n = 7), two (n = 11) or three (n = 6) local or distant metastatic lesions (lymph node, liver, lung, brain) ; 16 rectal carcinomas with lymph node (n = 13/16), liver (n = 1/16) or lung (n = 2/16)  were included. All of these studies were performed by semi-conductor NGS technology as compared to the present study using a sequencing-by-synthesis tNGS approach. Few studies commented on the morphological aspects of the (microdissected?) tumour specimens in case when allele frequencies drop, which may merely be due to few (enriched?) tumour cells from CRC tissue specimens (e.g. especially in neoadjuvant treated primary rectal tumours or highly necrotic liver metastases). Clearly, statistical sub-group analyses in such a small cohort are not feasible to address all of these clinico-pathological factors. Hence, the present study reports on and presents case-matched tNGS results with associated clinico-pathological data and morphological observations.
In fact, our present study addresses the morphological aspect from different angles: First, by precise microdissection under a microscope using fine needles we clearly eliminate e.g. lymphocyte-rich or necrotic areas from DNA sample preparation. Second, we provide the tumour cell content in terms of the percent of enriched (tumour) cell nuclei, not in terms of tissue area. Third, we report a feasible approach that is applicable to routine molecular pathology diagnostics – the initial diagnostic biopsy or resection specimen is used (not a separate tissue specimen obtained just for molecular analysis by an invasive approach, which is frequently impossible in advanced disease stages of CRC and other solid tumours) and all data interpretation is made with respect to the morphology, as specifically shown for the selected multiple distant metastatic tumour lesions of single resection specimens.
In this setting, we first report on the findings of the RAS mutation status by tNGS in the CRC cohort, with 9/14 cases showing known KRAS or NRAS mutations. Thereby, mutation status was maintained throughout the disease course albeit with altered allele frequencies. In fact, 2/9 cases showed a RAS mutation in the primary tumour, which was apparently “lost” at first sight in their metastases. Careful evaluation of the morphology clearly showed only few tumour cells available for DNA extraction in serial sections of the tissue specimens, which without morphological controlled enrichment might have been diluted in “non-tumour cell (and necrotic) background”. However, a morphology-adapted tNGS cut-off to lower allele frequencies then indeed showed RAS mutations at <5 % in the 2/5 cases. Similarly, observations were reported in a previous study for KRAS and other gene mutations , however no reference was made as to the morphology of these lesions, which in liver and lung tend to be highly necrotic with a paucity of tumour cells for enrichment by microdissection.
Hence, the present data validates mutations in genes (e.g. APC, KRAS, NRAS, TP53) known to be involved in (sporadic) CRC development and throughout progression [13, 16, 17, 27] with close reference to the clinico-pathological features. In fact, these key mutations, here referring to the exact sequence variant, were not directly associated with the primary tumour location (left or right colon, rectum) and histology, tumour stage, grading or molecular sub-type (Microsatellite-(in)stability; CpG-Island methylator phenotype) in this small cohort. Moreover, in the present study only a maximum of 3 sequence variants were shared between all CRC tumour sites throughout disease progression upon case-matched analyses, whereas broader and new mutation profiles were seen especially upon tumour progression to (synchronous or metachronous) lung metastases. The limited overlap of shared sequence variants between primary CRC and multiple metastatic lesions was not reported before. However, also others showed the same few genes to be similarly affected in primary tumours and associated distant metastases [13, 15, 17, 19].
In addition to the above key genes affected, de novo mutations occurring in liver or lung metastases of case-matched tumours in at least 3 cases were in APC, ATM, FBXW7, FGFR3, GNAQ, KIT, PIK3CA, PTEN, SMAD4, SMO, STK11, TP53 and VHL. Thereby, synchronous vs metachronous as well as liver vs lung metastases displayed distinct gene patterns suggesting distinct tumour evolution. Of these, several mutations may be actionable through targeted therapies interfering with the affected signaling pathways, e.g. FBXW7 (Additional file 1: Figure S2) (this study and [14, 16, 17, 28–30]) and its functional consequences on mTOR and PTEN signaling [31–33] and mTOR inhibitors . These data suggest that when looking at advanced disease stages of CRCs, each patient needs to be individually addressed due to the case-specific prior surgical and clinical treatment regimens and associated case-specific molecular progression. Hence, sequential biopsies – if clinically possible – or accompanying liquid biopsies may eventually evolve as supplementary tool to SOP-driven molecular tissue based analyses for defining individual treatment options.
Together, our present study demonstrates a novel approach for targeted next generation sequencing of routine diagnostic FFPE tissue specimens and thereby again highlights the importance of morphological-controlled tissue specimens. By this, the study identifies only few shared mutations, respective exact sequence variants throughout the disease course of individual colorectal cancer patients, including KRAS and NRAS. Furthermore, the study shows case-matched de novo altered mutation profiles in synchronous vs metachrounous liver and/or lung metastasis. Some of the affected genes - e.g. FBXW7 and PTEN may thereby be further explored for targeted therapy, including validation of findings in larger prospective case series.
Douillard JY, Siena S, Cassidy J, et al. Final results from PRIME: randomized phase III study of panitumumab with FOLFOX4 for first-line treatment of metastatic colorectal cancer. Ann Oncol. 2014;25:1346–55.
Wong NA, Gonzalez D, Salto-Tellez M, et al. RAS testing of colorectal carcinoma—a guidance document from the Association of Clinical Pathologists Molecular Pathology and Diagnostics Group. J Clin Pathol. 2014;67:751–7.
Schwartzberg LS, Rivera F, Karthaus M, et al. PEAK: a randomized, multicenter phase II study of panitumumab plus modified fluorouracil, leucovorin, and oxaliplatin (mFOLFOX6) or bevacizumab plus mFOLFOX6 in patients with previously untreated, unresectable, wild-type KRAS exon 2 metastatic colorectal cancer. J Clin Oncol. 2014;32:2240–7.
Ashraf N, Kothari N, Kim R. Predictive biomarkers for anti-epidermal growth factor receptor therapy: beyond KRAS testing. J Natl Compr Canc Netw. 2014;12:1433–42.
Modest DP, Stintzing S, von Weikersthal LF, et al. Impact of Subsequent Therapies on Outcome of the FIRE-3/AIO KRK0306 Trial: First-Line Therapy With FOLFIRI Plus Cetuximab or Bevacizumab in Patients With KRAS Wild-Type Tumors in Metastatic Colorectal Cancer. J Clin Oncol. 2015;33:3718–26.
Gonsalves WI, Mahoney MR, Sargent DJ, et al. Patient and tumor characteristics and BRAF and KRAS mutations in colon cancer, NCCTG/Alliance N0147. J Natl Cancer Inst. 2014;106(7).
Singh RR, Patel KP, Routbort MJ, et al. Clinical validation of a next-generation sequencing screen for mutational hotspots in 46 cancer-related genes. J Mol Diagn. 2013;15:607–22.
Spencer DH, Sehn JK, Abel HJ, Watson MA, Pfeifer JD, Duncavage EJ. Comparison of clinical targeted next-generation sequence data from formalin-fixed and fresh-frozen tissue specimens. J Mol Diagn. 2013;15:623–33.
Endris V, Penzel R, Warth A, et al. Molecular diagnostic profiling of lung cancer specimens with a semiconductor-based massive parallel sequencing approach: feasibility, costs, and performance compared with conventional sequencing. J Mol Diagn. 2013;15:765–75.
Betge J, Kerr G, Miersch T, et al. Amplicon sequencing of colorectal cancer: variant calling in frozen and formalin-fixed samples. PLoS One. 2015;10:e0127146.
Hagemann IS, Devarakonda S, Lockwood CM, et al. Clinical next-generation sequencing in patients with non-small cell lung cancer. Cancer. 2015;121:631–9.
Cheng DT, Mitchell TN, Zehir A, et al. Memorial Sloan Kettering-Integrated Mutation Profiling of Actionable Cancer Targets (MSK-IMPACT): A Hybridization Capture-Based Next-Generation Sequencing Clinical Assay for Solid Tumour Molecular Oncology. J Mol Diagn. 2015;17:251–64.
Vignot S, Lefebvre C, Frampton GM, et al. Comparative analysis of primary tumour and matched metastases in colorectal cancer patients: evaluation of concordance between genomic and transcriptional profiles. Eur J Cancer. 2015;51:791–9.
Lee SY, Haq F, Kim D, et al. Comparative genomic analysis of primary and synchronous metastatic colorectal cancers. PLoS One. 2014;9:e90459.
Tan IB, Malik S, Ramnarayanan K, 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;6:32.
Crumley SM, Pepper KL, Phan AT, Olsen RJ, Schwartz MR, Portier BP. Next-Generation Sequencing of Matched Primary and Metastatic Rectal Adenocarcinomas Demonstrates Minimal Mutation Gain and Concordance to Colonic Adenocarcinomas. Arch Pathol Lab Med. 2016;140(6):529–35.
Jesinghaus M, Wolf T, Pfarr N, et al. Distinctive Spatiotemporal Stability of Somatic Mutations in Metastasized Microsatellite-stable Colorectal Cancer. Am J Surg Pathol. 2015;39:1140–7.
Lim B, Mun J, Kim JH, et al. Genome-wide mutation profiles of colorectal tumors and associated liver metastases at the exome and transcriptome levels. Oncotarget. 2015;6:22179–90.
Kim R, Schell MJ, Teer JK, Greenawalt DM, Yang M, Yeatman TJ. Co-evolution of somatic variation in primary and metastatic colorectal cancer may expand biopsy indications in the molecular era. PLoS One. 2015;10:e0126670.
Berg KD, Glaser CL, Thompson RE, Hamilton SR, Griffin CA, Eshleman JR. Detection of microsatellite instability by fluorescence multiplex polymerase chain reaction. J Mol Diagn. 2000;2:20–8.
Weisenberger DJ, Siegmund KD, Campan M, et al. CpG island methylator phenotype underlies sporadic microsatellite instability and is tightly associated with BRAF mutation in colorectal cancer. Nat Genet. 2006;38:787–93.
Ogino S, Kawasaki T, Kirkner GJ, Kraft P, Loda M, Fuchs CS. Evaluation of markers for CpG island methylator phenotype (CIMP) in colorectal cancer by a large population-based sample. J Mol Diagn. 2007;9:305–14.
Zlobec I, Bihl M, Foerster A, Rufle A, Lugli A. Comprehensive analysis of CpG island methylator phenotype (CIMP)-high, −low, and -negative colorectal cancers based on protein marker expression and molecular features. J Pathol. 2011;225:336–43.
Robinson JT, Thorvaldsdóttir H, Winckler W, et al. Integrative genomics viewer. Nat Biotechnol. 2011;29:24–6.
Thorvaldsdóttir H, Robinson JT, Mesirov JP. Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2013;14:178–92.
Fichter CD, Herz C, Münch C, Opitz OG, Werner M, Lassmann S. Occurrence of multipolar mitoses and association with Aurora-A/-B kinases and p53 mutations in aneuploid esophageal carcinoma cells. BMC Cell Biol. 2011;12:13.
Carethers JM, Jung BH. Genetics and Genetic Biomarkers in Sporadic Colorectal Cancer. Gastroenterology. 2015;149:1177–90.
Xie T, Cho YB, Wang K, et al. Patterns of somatic alterations between matched primary and metastatic colorectal tumors characterized by whole-genome sequencing. Genomics. 2014;104:234–41.
Mlecnik B, Bindea G, Kirilovsky A, et al. The tumor microenvironment and Immunoscore are critical determinants of dissemination to distant metastasis. Sci Transl Med. 2016;8:327ra26.
Malapelle U, Pisapia P, Sgariglia R, et al. Less frequently mutated genes in colorectal cancer: evidences from next-generation sequencing of 653 routine cases. J Clin Pathol. 2016.
Wang L, Ye X, Liu Y, Wei W, Wang Z. Aberrant regulation of FBW7 in cancer. Oncotarget. 2014;5:2000–15.
Li L, Sarver AL, Khatri R, et al. Sequential expression of miR-182 and miR-503 cooperatively targets FBXW7, contributing to the malignant transformation of colon adenoma to adenocarcinoma. J Pathol. 2014;234:488–501.
Mao JH, Kim IJ, Wu D, et al. FBXW7 targets mTOR for degradation and cooperates with PTEN in tumor suppression. Science. 2008;321:1499–502.
Jardim DL, Wheler JJ, Hess K, et al. FBXW7 mutations in patients with advanced cancers: clinical and molecular characteristics and outcomes with mTOR inhibitors. PLoS One. 2014;9:e89388.
The authors thank Ms. Marion Kunz, Ms. Bianca Riedel and Ms. Xenia Ungefug (Institute of Surgical Pathology, Medical Center, University of Freiburg, Germany) for expert technical assistance. The study was supported by the German Consortium for Translational Cancer Research (DKTK, grant to SL, MW) and in part by the Mushett Family Foundation (Chester, NJ, USA; grant to SL, MW).
This study was funded by the German Cancer Consortium (DKTK, grants to MW and SL) and in part by the Mushett Family Foundation (Chester, New Jersey, USA).
Availability of data and material
The datasets during and/or analyzed during the current study available from the corresponding author on reasonable request.
VK analyzed NGS and ddSeq data, ALG reevaluated and retrieved the cohort and tissue specimens, performed DNA extraction, established CIMP analyses, VKO and ALG drafted the manuscript, LL performed histological analyses, RF, FM, SW, UH and BP provided clinical expertise, MW approved histological data, SL designed and supervised the study, established the NGS workflow and analyzed NGS data, compiled graphs, figures and wrote the manuscript. All authors have read and agreed to the final version of the manuscript.
The authors declare that they have no ompeting interests.
Consent for publication
Ethics approval and consent to participate
The use of archival tissue specimens had been approved by the local ethics Institution (#251/04, #66/07 and #191/13, Ethik-Kommission, University of Freiburg, Germany).
Text S1. Establishment of library preparation from DNA samples derived from different origin tissue specimens. Adaptations of the tNGS library protocol. Table S1. Modifications to tNGS library protocol of low input or poorer quality DNA improves sequencing performance. Table S2. Details of tissue specimens and associated DNA and library processing data. Table S3. Validation of variants detected by tNGS by dideoxy sequencing. Figure S1. Quantity and quality of tNGS libraries of the standard and adapted library protocols. Figure S2. De novo FBXW7 sequence variants detected in colorectal carcinoma liver and/or lung metastases. (DOCX 556 kb)
Summary file of all detected sequence variants. (XLSX 69 kb)