Comparative sequencing analysis reveals high genomic concordance between matched primary and metastatic colorectal cancer lesions
- A Rose Brannon1,
- Efsevia Vakiani1,
- Brooke E Sylvester2,
- Sasinya N Scott1,
- Gregory McDermott1, 7,
- Ronak H Shah1,
- Krishan Kania2,
- Agnes Viale5,
- Dayna M Oschwald6,
- Vladimir Vacic6,
- Anne-Katrin Emde6,
- Andrea Cercek3,
- Rona Yaeger3,
- Nancy E Kemeny3,
- Leonard B Saltz3,
- Jinru Shia1,
- Michael I D’Angelica4,
- Martin R Weiser4,
- David B Solit2, 3, 7Email author and
- Michael F Berger1, 2, 7Email author
© Brannon et al.; licensee BioMed Central Ltd. 2014
Received: 30 April 2014
Accepted: 22 August 2014
Published: 28 August 2014
Colorectal cancer is the second leading cause of cancer death in the United States, with over 50,000 deaths estimated in 2014. Molecular profiling for somatic mutations that predict absence of response to anti-EGFR therapy has become standard practice in the treatment of metastatic colorectal cancer; however, the quantity and type of tissue available for testing is frequently limited. Further, the degree to which the primary tumor is a faithful representation of metastatic disease has been questioned. As next-generation sequencing technology becomes more widely available for clinical use and additional molecularly targeted agents are considered as treatment options in colorectal cancer, it is important to characterize the extent of tumor heterogeneity between primary and metastatic tumors.
We performed deep coverage, targeted next-generation sequencing of 230 key cancer-associated genes for 69 matched primary and metastatic tumors and normal tissue. Mutation profiles were 100% concordant for KRAS, NRAS, and BRAF, and were highly concordant for recurrent alterations in colorectal cancer. Additionally, whole genome sequencing of four patient trios did not reveal any additional site-specific targetable alterations.
Colorectal cancer primary tumors and metastases exhibit high genomic concordance. As current clinical practices in colorectal cancer revolve around KRAS, NRAS, and BRAF mutation status, diagnostic sequencing of either primary or metastatic tissue as available is acceptable for most patients. Additionally, consistency between targeted sequencing and whole genome sequencing results suggests that targeted sequencing may be a suitable strategy for clinical diagnostic applications.
Precision oncology relies on the accurate characterization of targetable oncogenic mutations present at the time of metastatic disease. However, it is often challenging to obtain biopsies of metastatic tumors, and it is still preferable to use the least invasive screening methods possible. Emerging evidence of both inter- and intra-tumor genetic heterogeneity in several solid tumor types raises concerns that molecular profiling of primary tumors may not be representative of metastatic disease -. In colorectal cancer (CRC), comparative lesion sequencing of a small number of cases found a high degree of concordance between primary tumors and metastases . In contrast, a recent study of 21 patients using next generation sequencing reported a high degree of mutational discordance between primary and metastatic samples .
We previously showed that when analysis was performed on the invasive compartment of primary tumors, KRAS, NRAS, BRAF, and TP53 mutations were highly concordant between primary and metastatic tumors . This study provided a preliminary indication that the use of archived primary tumor for molecular profiling may be suitable for clinical decision making in metastatic CRC. However, this conclusion was based on the analysis of only a small number of genes by mass-spectrometry based genotyping and Sanger sequencing.
To determine the extent of additional, clinically relevant genetic heterogeneity, we extended this analysis by performing high coverage, next generation sequencing analysis of 230 cancer-associated genes. Specifically, we performed targeted sequencing on primary, metastatic, and normal tissue from 69 colorectal cancer patients. We found that there was a high degree of concordance with regard to early occurring and recurrent mutations. KRAS, NRAS, and BRAF mutations were always identical in both the primary and metastatic tumors. Whole genome sequencing of two concordant and two discordant patient sets upheld the targeted sequencing results and revealed few additional recurrent mutations. In sum, these data suggest that for current clinical practices, either primary or metastatic tissue can be selected for testing, and targeted sequencing of key cancer genes is a suitable strategy for identifying clinically actionable alterations.
Clinical characteristics of CRC cases subjected to targeted sequencing
N = 69
Median age, years
Primary tumor location
Stage at presentation
Mutation profiles are highly concordant between primary and metastatic tumors
Despite the overall high concordance in mutation profiles we observed, we chose to further investigate the remaining discordance. To determine whether any of this discordance could be explained by intra-tumor heterogeneity, we sequenced 97 additional samples from spatially separate regions of the primary tumors (n = 62) and metastases (n = 37) from 22 patients, encompassing 46 discordant mutations. All 97 samples were derived from formalin fixed paraffin embedded (FFPE) tissue and reviewed for morphology and tumor cellularity. Sequencing of multiple regions and samples resolved 17/46 discordant mutations, including 12/22 (55%) mutations that were originally detected only in the metastasis but were subsequently found to be subclonal in the primary tumor (Additional file 1: Table S1). Altogether these results suggest that the small proportion of discordant mutations we observed may itself be an overestimate of the true discordance when accounting for intra-tumor heterogeneity.
Private oncogenic alterations are occasionally detected
Convergent phenotypic evolution was observed in two patients harboring distinct mutations in individual genes, one involving independent mutations in TP53 (R248Q and Y163*) and the other involving separate hotspot mutations in PIK3CA (E542K and E545K, Additional file 2: Figure S2). Three additional private events were observed in PIK3CA including two E545K mutations specific to the primary tumor and one in-frame deletion at N107 specific to the metastasis. Private mutations of unknown significance were also found in PIK3CD, PIK3CG, PIK3C2G, PIK3R1, and PTEN. Four of these PI3K pathway events were found in other regions, indicative of subclonality. These results suggest that despite the overall high level of genomic concordance between primary and metastatic CRC, heterogeneity in potentially actionable genes, such as those within the PI3K pathway, is present in at least a subset of patients. This finding could have therapeutic implications given the ongoing evaluation of PI3K inhibitors in clinical trials as well as the potential benefit of aspirin therapy in CRC patients harboring PIK3CA mutations .
Clinical correlates of mutational concordance
We next examined whether mutational concordance was correlated with major clinical characteristics (Additional file 3: Table S2). We observed no significant difference in the number of mutations or degree of mutational concordance for cases with different primary tumor locations (right colon, left colon, or rectum) or for cases with different time intervals between primary and metastasis resections (concurrent versus subsequent; Additional file 4: Table S3). We also observed no disease-specific survival differences for patients with concordant and discordant mutation profiles (Additional file 2: Figure S5). However, among patients whose tumors were concurrently resected, those that did not receive prior treatment were more likely to harbor discordant mutations (22/28, 79%) compared to those that received prior therapy (11/24, 46%; chi-square P = 0.01. Patients without prior therapy also had a higher total number of mutations, although the difference was not statistically significant (6.6 vs. 5.8, P = 0.1). These differences may be due to either general tumor debulking and/or decreased tumor heterogeneity from effective drug treatment . In accordance with this hypothesis, the primary tumor size in patients receiving prior treatment was slightly smaller than in chemonaive patients (4.2 vs. 5.2 cm, T test P = 0.02). Additionally, pre-treated patients with concurrently resected tumors were less likely to harbor primary-only mutations than patients with subsequent resections where only the metastasis received treatment (6/24, 25% vs. 7/11, 64%; chi-square P = 0.03).
Whole genome analysis reveals consistent mutation patterns
While our capture-based assay encompasses all well-established targetable or actionable genes in CRC, this approach by definition would fail to identify discordant mutations not included in the assay design. To determine whether our results were representative of the level of mutational concordance genome-wide, we performed whole genome sequencing (WGS) of two cases where the primary tumor and metastasis were concordant by IMPACT analysis (patients 14 and 54) and two cases that harbored discordant mutations (patients 3 and 19, with one and three discordant mutations by IMPACT, respectively). Given the slight disparity in mutations for pre-treated versus chemonaive tumors, we chose the patients such that within each group, one patient’s tumors were chemonaive (patients 3 and 54) and the other patient’s tumors had received prior treatment (patients 14 and 19). Primary and metastatic tumors were sequenced to >80× coverage, and matched normal controls were sequenced to >40× coverage. All SNVs and indels detected at higher than 5% allele frequency by IMPACT were detected in our WGS analysis, and no additional mutations in the 230 genes were identified, providing independent validation of the results produced by both methods.
In contrast, WGS of the paired tumors from the two patients who exhibited 100% concordance by IMPACT analysis (patients 14 and 54) showed a high degree of similarity, with 80% and 83% shared mutations overall. As above, the level of concordance increased even further upon restricting the analysis to non-synonymous protein-coding mutations (87% and 86% shared) and loss-of-function mutations (91% and 91% shared). For patient 54, a single discordant mutation in RUNXT1 was the only discordant alteration present in COSMIC, and no discordant mutations were nonsense or splice-site mutations. Similarly none of the discordant mutations in patient 14 are currently represented in the COSMIC database. Case 14 was notable in that analysis by IMPACT revealed no mutations in both the beta-catenin and RAS pathways, although the tumor harbored an amplification of CDK8, which has be shown to dysregulate the WNT/beta-catenin pathway . No additional beta-catenin pathway alterations were identified by WGS; however, two nonsense mutations in RASA1 (C372* and R679*) were detected in the both primary and metastatic samples. RASA1, a negative regulator of RAS, was altered in 2% of CRC TCGA samples, with mutations in this gene mutually exclusive of KRAS and NRAS mutations. This result prompted the addition of capture probes for the RASA1 gene in our current version of IMPACT and highlights the utility of whole genome sequencing to identify rare alterations in driver pathways such as RAS in tumors pan-negative for the more commonly mutated genes.
We performed deep sequencing of 230 cancer-associated genes in 69 primary CRC tumors and matched metastases to define the mutational concordance of these genes in primary and metastatic tumors. We identified a pattern of genomic alterations consistent with prior studies of CRC, including frequent mutations in the APC, TP53, and KRAS genes ,. Mutations in TP53 were more prevalent in our series than TCGA. TCGA analyzed only primary tumors, the majority of which were derived from patients with stage 1 to 3 disease, whereas all patients in our cohort by definition had metastatic disease. This difference is likely the basis for both the higher prevalence of TP53 mutation and the lower prevalence of NRAS mutations in our cohort versus TCGA. The higher prevalence of certain mutations observed from IMPACT may also be the result of greater sensitivity for mutation detection afforded by deeper sequence coverage (mean coverage of 692X). Although we did not observe significant differences in the degree of heterogeneity between those tumors that were resected concurrently or subsequently, future studies are required to determine the extent to which systemic treatments result in preferential selection of specific mutational events.
The only major clinical feature that seemed to correlate with mutational concordance was whether the tumors had received prior treatment. In patients where tumors were concurrently resected, chemonaive tumors were more likely to have private mutations than pre-treated tumors. Similarly, patients with pre-treated metastatic tumors subsequently resected were more likely to harbor primary-only mutations than patients with pre-treated concurrent resections. These findings raise the possibility that prior treatment results in a decrease of apparent tumor heterogeneity, though larger studies are needed to further explore this hypothesis.
Overall, our exon capture data indicate a high degree of concordance in mutational profiles, especially when considering events that occur early in colorectal carcinogenesis. Further, multiple examples of discordance involved parallel evolution demonstrated by independent distinct mutations in the same gene . A subsequent whole genome analysis of four paired primary and metastatic samples further suggested that concordance by IMPACT is a good surrogate of genome-wide mutational concordance. Notably, mutations in KRAS, NRAS, and BRAF were 100% concordant between primary tumors and metastases in our cohort, which is consistent with prior studies ,-. These results are also consistent with the recommendation that molecular testing of the primary tumor is appropriate in most clinical scenarios. Exceptions would include patients with a history of multiple primary tumors or polyps and patients in which interval drug treatment may result in clonal selection of clinically actionable mutations . An example of the latter would be treatment with anti-EGFR therapy, which has been shown to promote the selection of RAS mutant subclones in tumors that are otherwise KRAS, NRAS, and BRAF wild-type ,. In primary tumors lacking KRAS, NRAS, or BRAF mutations, we did identify occult alterations in the EGFR/RAS pathway that in some patients were private to the metastatic tumor. Discordant mutations were also occasionally observed in components of the PI3K pathway, which could be clinically important as selective inhibitors of the pathway are being actively studied in patients with CRC.
Finally, we found that targeted sequencing was sufficient to identify the most clinically actionable alterations found by WGS. This result suggests that an updated IMPACT test may be a viable time- and cost-saving alternative to WGS for molecular profiling of CRC patients.
Inter- and intra-tumor heterogeneity, a growing concern in the molecular diagnostics field, is apparent in CRC. However, in patients who have not been treated with anti-EGFR therapies, the current clinically actionable genes, KRAS, NRAS, and BRAF, are 100% concordant between primary and metastatic tissues. As the mutational status of these genes guides current clinical practice, diagnostic testing from either tissue site as available is appropriate in most clinical scenarios. Additionally, targeted sequencing is becoming a more common practice in research and clinical settings, and our results demonstrate the clinical utility of this approach both through comparison of primary and metastatic tissue and similarities of results to whole genome sequencing.
With Institutional Review Board approval (WA0129-12) and compliance with the Helsinki Declaration, we analyzed 69 matched trios (normal, primary, and metastatic tissues) from patients undergoing resection at our institution where frozen tissue was available. The set was enriched for patients with stage IV disease at diagnosis to decrease potential discordant alterations caused by time. All specimens were reviewed for histological verification of a colorectal adenocarcinoma diagnosis and to ensure greater than 50% tumor content. Macrodissection was performed on specimens with less than 50% viable tumor to minimize stromal contamination. Normal DNA was obtained from normal colon tissue located at least 15 cm away from the tumor. Normal and tumor DNA were extracted from shaved sections cut from frozen tissue blocks using the Qiagen DNeasy Blood & Tissue Kit. In two cases where the frozen section contained an adenoma, invasive regions of the primary were macrodissected from formalin-fixed paraffin-embedded sections, and DNA was extracted using the Qiagen DNeasy Blood & Tissue Kit modified for deparaffinization and these sequencing results were used in place of the original frozen data. Additional regions for determination of private mutations were selected from available FFPE tissue blocks and similarly macrodissected and prepared.
The microsatellite instability (MSI) status was determined for each case using a 5-microsatellite marker (BAT25, BAT26, D17S250, D2S123, and D5S346) genotyping platform according to a standard protocol . Fluorophore-labeled primers were designed (Applied Biosystems) targeting the five loci. All microsatellite loci were amplified for matched normal and tumor DNA in a multiplex polymerase chain reaction (PCR) and submitted for genotyping to the MSKCC Genomics Core. Microsatellite marker stability was analyzed using Peak Scanner™ software. MSI status was categorized as microsatellite stable if all markers were stable, MSI-low if <30% of markers were unstable, and MSI-high if ≥30% of markers were unstable.
IMPACT targeted sequencing
Library preparation and sequence analysis is as previously described . Briefly, 100 to 500 ng DNA from frozen or formalin fixed paraffin embedded (FFPE) tissue was prepared using NEBNext DNA Library Prep Kit for Illumina with Kapa HiFi DNA Polymerase for PCR steps. A total of 100 ng resulting library was pooled in sets of 12 to 24 samples for capture with custom Nimblegen probes. Each pool was sequenced in a single lane of an Illumina HiSeq 2000. Resulting fastq files were aligned according to best practice with BWA , GATK, Picard , and Samtools . Mutations and indels were called with Mutect  and GATK SomaticIndelDetector, respectively, and annotated with Oncotator  and Cosmic v65 . Only non-synonymous mutations above 5% allele frequency in exonic regions were retained. For primary or metastasis-specific mutations, the allele counts at the corresponding site in the alternate sample were examined. If the mutant allele was present in at least 3 total reads and at least 2% of all reads, the event was labeled as `both’ primary and metastasis by comparative analysis. Copy number analysis was performed using average read depth from GATK, loess-normalized for GC content, and compared to diploid normal.
Cell lines and culture
Human embryonic kidney 293H cells were maintained in DME-HG medium supplemented with 10% FBS, 2 mM glutamine, and 50 units/mL each of penicillin and streptomycin. MAP2K1 mutations were generated from the MEK1-GFP plasmid (Addgene, 14746) using the QuickChange Site-Directed Mutagenesis Kit (Stratagene) as recommended. All mutant plasmids were verified by DNA sequencing. 293H cells were seeded for 70% to 90% confluency at the time of transfection in the culture medium without penicillin and streptomycin overnight. Cells were transiently transfected with wild-type or mutant MAP2K1 DNA using the Lipofectamine 2000 Transfection Reagent as recommended.
Western blot analysis
At 24 h post transfection, cells were collected and lysed in 1% NP-40 lysis buffer and processed for immunoblotting as previously described . Rabbit polyclonal antibodies recognizing phosphorylated Erk1/2 (Thr202/Tyr204), Erk1/2, phosphorylated p90RSK (Thr359/Ser363), and phosphorylated p90RSK (Ser380) were obtained from Cell Signaling. Rabbit monoclonal antibodies recognizing RSK1/2/3, GFP, and GAPDH were obtained from Cell Signaling. After incubation with horseradish peroxidase-conjugated secondary antibodies, proteins were detected by chemiluminescence (SuperSignal West Dura Chemiluminescent Substrate, Thermo Scientific) and visualized using the Fuji LAS-4000 (GE Life Sciences).
Fluorescence in situ hybridization
Fluorescence in situ hybridization was performed by the MSKCC Molecular Cytogenetics Core using BAC clone RP11-339 F13 and PAC clone RP5-1091E12 spanning the EGFR locus in 7p11, both labeled by nick translation with Red dUTP, together with a chromosome 7 centromeric repeat probe (p7t1) labeled with Green dUTP (Enzo Life Sciences, supplied by Abbott Molecular). Initial hybridization showed weak FISH signals against a high level of background autofluorescence. Fresh probe was prepared and new slides were hybridized. FISH signals for the second hybridization were much stronger but the background autofluorescence was still high. Images stacks (9× 0.5 um slices for Red & Green) were captured through the depth of the tissue for recording and analysis. R/G signal ratios were scored in representative tumor fields with a minimum of 30 nuclei, where possible.
Whole genome sequencing
Five micrograms of DNA from patient-matched frozen primary tumor, metastatic tumor, and normal tissue for EV-003, EV-014, EV-019, and EV-054 were sent for whole genome sequencing at the New York Genome Center (NYGC) and Illumina. Median coverage was 87X for tumors and 50X for normals. Mutations were called on BAM files produced by the Illumina CASAVA alignment pipeline, using Strelka . These mutations were then genotyped on BAM files produced by the BWA aligner, using GATK Unified Genotyper ,. Calls with fewer than 3 reads or 5% frequency in either the primary and metastatic tumor were removed, as well as those that were less than five times the frequency in the normal tissue. Remaining calls were annotated for the canonical isoform using Oncotator and COSMIC v65 .
Data are publically available through dbGaP (accession phs000790.v1.p1) and the Memorial Sloan Kettering Cancer Center cBioPortal for Cancer Genomics (; study `Colorectal Adenocarcinoma Triplets’) .
ARB, EV, DBS, and MFB conceived the project, designed and analyzed the experiments, and wrote the manuscript. BES, SNS, and AV performed and supervised laboratory experiments. ARB, SNS, GCM, RHS, KK, DMO, VV, and AKE analyzed next generation sequencing data. EV, AC, RY, NEK, LBS, JS, MID, and MRW collected and analyzed patient samples. All authors read and approved the final manuscript.
Integrated mutation profiling of actionable cancer targets
Single nucleotide variant
The Cancer Genome Atlas
Whole genome sequencing
The authors wish to thank the MSKCC Geoffrey Beene Core and Genomics Core for their assistance in sample preparation and sequencing; the MSKCC Molecular Cytogenetics Core for FISH studies; and the New York Genome Center for whole genome sequencing. Research reported in this publication was supported by the Geoffrey Beene Cancer Research Center, the Farmer Family Foundation, and the National Cancer Institute of the National Institutes of Health under award numbers T32CA160001 and 1K08CA175150-01A1.
- Bedard PL, Hansen AR, Ratain MJ, Siu LL: Tumour heterogeneity in the clinic. Nature. 2013, 501: 355-364. 10.1038/nature12627.PubMedView ArticleGoogle Scholar
- Gerlinger M, Rowan AJ, Horswell S, Larkin J, Endesfelder D, Gronroos E, Martinez P, Matthews N, Stewart A, Tarpey P, Varela I, Phillimore B, Begum S, McDonald NQ, Butler A, Jones D, Raine K, Latimer C, Santos CR, Nohadani M, Eklund AC, Spencer-Dene B, Clark G, Pickering L, Stamp G, Gore M, Szallasi Z, Downward J, Futreal A, Swanton C: Intratumor heterogeneity and branched evolution revealed by multiregion sequencing. N Engl J Med. 2012, 366: 883-892. 10.1056/NEJMoa1113205.PubMedView ArticleGoogle Scholar
- Yachida S, Jones S, Bozic I, Antal T, Leary R, Fu B, Kamiyama M, Hruban RH, Eshleman JR, Nowak MA, Velculescu VE, Kinzler KW, Vogelstein B, Iacobuzio-Donahue CA: Distant metastasis occurs late during the genetic evolution of pancreatic cancer. Nature. 2010, 467: 1114-1117. 10.1038/nature09515.PubMedPubMed CentralView ArticleGoogle Scholar
- Jones S, Chen WD, Parmigiani G, Diehl F, Beerenwinkel N, Antal T, Traulsen A, Nowak MA, Siegel C, Velculescu VE, Kinzler KW, Vogelstein B, Willis J, Markowitz SD: Comparative lesion sequencing provides insights into tumor evolution. Proc Natl Acad Sci U S A. 2008, 105: 4283-4288. 10.1073/pnas.0712345105.PubMedPubMed CentralView ArticleGoogle Scholar
- Vermaat JS, Nijman IJ, Koudijs MJ, Gerritse FL, Scherer SJ, Mokry M, Roessingh WM, Lansu N, de Bruijn E, van Hillegersberg R, van Diest PJ, Cuppen E, Voest EE: 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-699. 10.1158/1078-0432.CCR-11-1965.PubMedView ArticleGoogle Scholar
- Vakiani E, Janakiraman M, Shen R, Sinha R, Zeng Z, Shia J, Cercek A, Kemeny N, D’Angelica M, Viale A, Heguy A, Paty P, Chan TA, Saltz LB, Weiser M, Solit DB: Comparative genomic analysis of primary versus metastatic colorectal carcinomas. J Clin Oncol. 2012, 30: 2956-2962. 10.1200/JCO.2011.38.2994.PubMedPubMed CentralView ArticleGoogle Scholar
- Comprehensive molecular characterization of human colon and rectal cancer. Nature. 2012, 487: 330-337. 10.1038/nature11252.
- Fearon ER, Vogelstein B: A genetic model for colorectal tumorigenesis. Cell. 1990, 61: 759-767. 10.1016/0092-8674(90)90186-I.PubMedView ArticleGoogle Scholar
- Domingo E, Church DN, Sieber O, Ramamoorthy R, Yanagisawa Y, Johnstone E, Davidson B, Kerr DJ, Tomlinson IP, Midgley R: Evaluation of PIK3CA mutation as a predictor of benefit from nonsteroidal anti-inflammatory drug therapy in colorectal cancer. J Clin Oncol. 2013, 31: 4297-4305. 10.1200/JCO.2013.50.0322.PubMedView ArticleGoogle Scholar
- Lavi O, Greene JM, Levy D, Gottesman MM: The role of cell density and intratumoral heterogeneity in multidrug resistance. Cancer Res. 2013, 73: 7168-7175. 10.1158/0008-5472.CAN-13-1768.PubMedView ArticleGoogle Scholar
- Peradziryi H, Kaplan NA, Podleschny M, Liu X, Wehner P, Borchers A, Tolwinski NS: PTK7/Otk interacts with Wnts and inhibits canonical Wnt signalling. EMBO J. 2011, 30: 3729-3740. 10.1038/emboj.2011.236.PubMedPubMed CentralView ArticleGoogle Scholar
- Firestein R, Bass AJ, Kim SY, Dunn IF, Silver SJ, Guney I, Freed E, Ligon AH, Vena N, Ogino S, Chheda MG, Tamayo P, Finn S, Shrestha Y, Boehm JS, Jain S, Bojarski E, Mermel C, Barretina J, Chan JA: CDK8 is a colorectal cancer oncogene that regulates beta-catenin activity. Nature. 2008, 455: 547-551. 10.1038/nature07179.PubMedPubMed CentralView ArticleGoogle Scholar
- Voss MH, Hakimi AA, Pham CG, Brannon AR, Chen YB, Cunha LF, Akin O, Liu H, Takeda S, Scott SN, Socci ND, Viale A, Schultz N, Sander C, Reuter VE, Russo P, Cheng EH, Motzer RJ, Berger MF, Hsieh JJ: Tumor Genetic Analyses of Patients with Metastatic Renal Cell Carcinoma and Extended Benefit from mTOR Inhibitor Therapy. Clin Cancer Res. 2014, 20: 1955-1964. 10.1158/1078-0432.CCR-13-2345.PubMedPubMed CentralView ArticleGoogle Scholar
- Italiano A, Hostein I, Soubeyran I, Fabas T, Benchimol D, Evrard S, Gugenheim J, Becouarn Y, Brunet R, Fonck M, Francois E, Saint-Pual MC, Pedeutour F: KRAS and BRAF mutational status in primary colorectal tumors and related metastatic sites: biological and clinical implications. Ann Surg Oncol. 2010, 17: 1429-1434. 10.1245/s10434-009-0864-z.PubMedView ArticleGoogle Scholar
- Knijn N, Mekenkamp LJ, Klomp M, Vink-Borger ME, Tol J, Teerenstra S, Meijer JW, Tebar M, Riemersma S, van Krieken JH, Punt CJ, Nagtegaal ID: KRAS mutation analysis: a comparison between primary tumours and matched liver metastases in 305 colorectal cancer patients. Br J Cancer. 2011, 104: 1020-1026. 10.1038/bjc.2011.26.PubMedPubMed CentralView ArticleGoogle Scholar
- Santini D, Loupakis F, Vincenzi B, Floriani I, Stasi I, Canestrari E, Rulli E, Maltese PE, Andreoni F, Masi G, Graziano F, Baldi GG, Salvatore L, Russo A, Perrone G, Tommasino MR, Magnani M, Falcone A, Tonini G, Ruzzo A: High concordance of KRAS status between primary colorectal tumors and related metastatic sites: implications for clinical practice. Oncologist. 2008, 13: 1270-1275. 10.1634/theoncologist.2008-0181.PubMedView ArticleGoogle Scholar
- Diaz LA, Williams RT, Wu J, Kinde I, Hecht JR, Berlin J, Allen B, Bozic I, Reiter JG, Nowak MA, Kinzler KW, Oliner KS, Vogelstein B: The molecular evolution of acquired resistance to targeted EGFR blockade in colorectal cancers. Nature. 2012, 486: 537-540.PubMedPubMed CentralGoogle Scholar
- Misale S, Yaeger R, Hobor S, Scala E, Janakiraman M, Liska D, Valtorta E, Schiavo R, Buscarino M, Siravegna G, Bencardino K, Cercek A, Chen CT, Veronese S, Zanon C, Sartore-Bianchi A, Gambacorta M, Gallicchio M, Vakiani E, Boscaro V: Emergence of KRAS mutations and acquired resistance to anti-EGFR therapy in colorectal cancer. Nature. 2012, 486: 532-536.PubMedPubMed CentralGoogle Scholar
- Boland CR, Thibodeau SN, Hamilton SR, Sidransky D, Eshleman JR, Burt RW, Meltzer SJ, Rodriguez-Bigas MA, Fodde R, Ranzani GN, Srivastava S: A National Cancer Institute Workshop on Microsatellite Instability for cancer detection and familial predisposition: development of international criteria for the determination of microsatellite instability in colorectal cancer. Cancer Res. 1998, 58: 5248-5257.PubMedGoogle Scholar
- Won HH, Scott SN, Brannon AR, Shah RH, Berger MF: Detecting somatic genetic alterations in tumor specimens by exon capture and massively parallel sequencing. J Vis Exp. 2013, 80: e50710-PubMedGoogle Scholar
- Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25: 1754-1760. 10.1093/bioinformatics/btp324.PubMedPubMed CentralView ArticleGoogle Scholar
- Picard. , [http://broadinstitute.github.io/picard/]
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009, 25: 2078-2079. 10.1093/bioinformatics/btp352.PubMedPubMed CentralView ArticleGoogle Scholar
- Cibulskis K, Lawrence MS, Carter SL, Sivachenko A, Jaffe D, Sougnez C, Gabriel S, Meyerson M, Lander ES, Getz G: Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples. Nat Biotechnol. 2013, 31: 213-219. 10.1038/nbt.2514.PubMedView ArticleGoogle Scholar
- Oncotator. , [http://www.broadinstitute.org/oncotator]
- Forbes SA, Bindal N, Bamford S, Cole C, Kok CY, Beare D, Jia M, Shepherd R, Leung K, Menzies A, Butler AP, Teague JW, Campbell PJ, Stratton MR, Futreal PA: COSMIC: mining complete cancer genomes in the Catalogue of Somatic Mutations in Cancer. Nucleic Acids Res. 2011, 39: D945-D950. 10.1093/nar/gkq929.PubMedPubMed CentralView ArticleGoogle Scholar
- Pratilas CA, Hanrahan AJ, Halilovic E, Persaud Y, Soh J, Chitale D, Shigematsu H, Yamamoto H, Sawai A, Janakiraman M, Taylor BS, Pao W, Toyooka S, Ladanyi M, Gazdar A, Rosen N, Solit DB: Genetic predictors of MEK dependence in non-small cell lung cancer. Cancer Res. 2008, 68: 9375-9383. 10.1158/0008-5472.CAN-08-2223.PubMedPubMed CentralView ArticleGoogle Scholar
- Saunders CT, Wong WS, Swamy S, Becq J, Murray LJ, Cheetham RK: Strelka: accurate somatic small-variant calling from sequenced tumor-normal sample pairs. Bioinformatics. 2012, 28: 1811-1817. 10.1093/bioinformatics/bts271.PubMedView ArticleGoogle Scholar
- McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, DePristo MA: The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010, 20: 1297-1303. 10.1101/gr.107524.110.PubMedPubMed CentralView ArticleGoogle Scholar
- cBioPortal for Cancer Genomics. , [http://cbioportal.org]
- Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, Jacobsen A, Byrne CJ, Heuer ML, Larsson E, Antipin Y, Reva B, Goldberg AP, Sander C, Schultz N: The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012, 2: 401-404. 10.1158/2159-8290.CD-12-0095.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.