Next Article in Journal
Exosomal microRNA in Pancreatic Cancer Diagnosis, Prognosis, and Treatment: From Bench to Bedside
Previous Article in Journal
Prognostic and Therapeutic Role of Angiogenic Microenvironment in Thyroid Cancer
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 13
Issue 11
10.3390/cancers13112776
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Somatic Mutations in Solid Metastatic Pan-Cancer Patients

by
Moom R. Roosan
1,†,
Isa Mambetsariev
2,†,
Rebecca Pharaon
2,
Jeremy Fricke
2,
Angel R. Baroz
2,
Joseph Chao
2,
Chen Chen
3,
Mohd W. Nasser
4,
Ramakanth Chirravuri-Venkata
4,
Maneesh Jain
4,
Lynette Smith
5,
Susan E. Yost
2,
Karen L. Reckamp
2,6,
Raju Pillai
7,
Leonidas Arvanitis
7,
Michelle Afkhami
7,
Edward W. Wang
2,
Vincent Chung
2,
Mihaela Cristea
2,
Marwan Fakih
2,
Marianna Koczywas
2,
Erminia Massarelli
2,
Joanne Mortimer
2,
Yuan Yuan
2,
Surinder K. Batra
4,
Sumanta Pal
2,*,‡ and
Ravi Salgia
2,*,‡add Show full author list remove Hide full author list
1
School of Pharmacy, Chapman University, Irvine, CA 92618, USA
2
Department of Medical Oncology & Therapeutics Research, City of Hope, Duarte, CA 91010, USA
3
Applied AI and Data Science, City of Hope, Duarte, CA 91010, USA
4
Department of Biochemistry and Molecular Biology, University of Nebraska Medical Center, Omaha, NE 68198, USA
5
Department of Biostatistics, University of Nebraska Medical Center, Omaha, NE 68198, USA
6
Cedars-Sinai Medical Center, Department of Medicine, Division of Medical Oncology, Los Angeles, CA 90048, USA
7
Department of Pathology, City of Hope, Duarte, CA 91010, USA
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work and should be considered co-first authors.
‡
These authors should be considered co-senior authors.
Cancers 2021, 13(11), 2776; https://doi.org/10.3390/cancers13112776
Submission received: 15 April 2021 / Revised: 27 May 2021 / Accepted: 28 May 2021 / Published: 3 June 2021
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Cancer metastasis significantly contributes to cancer-related mortality. Our retrospective cohort study aimed to evaluate the mutational landscape of seven solid metastatic tumors and mutational effects on survival using a single molecular testing panel. Additionally, we assessed the treatments used in advanced cancer. We identified somatic mutations that were mutually exclusive in seven gene pairs. Among them, somatic mutations in APC and CDKN2A showed an opposite effect on overall survival (OS). Longer OS was associated with metastatic cases diagnosed post-2015. Progression-free survival was associated with the use of targeted treatments. Our results highlight complex interactions of mutational landscape with a single molecular test, time of metastatic diagnosis, and the impact of targeted therapy usage on survival using a pan-cancer cohort.

Abstract

Metastasis continues to be the primary cause of all cancer-related deaths despite the recent advancements in cancer treatments. To evaluate the role of mutations in overall survival (OS) and treatment outcomes, we analyzed 957 metastatic patients with seven major cancer types who had available molecular testing results with a FoundationOne CDx® panel. The most prevalent genes with somatic mutations were TP53, KRAS, APC, and LRP1B. In this analysis, these genes had mutation frequencies higher than in publicly available datasets. We identified that the somatic mutations were seven mutually exclusive gene pairs and an additional fifty-two co-occurring gene pairs. Mutations in the mutually exclusive gene pair APC and CDKN2A showed an opposite effect on the overall survival. However, patients with CDKN2A mutations showed significantly shorter OS (HR: 1.72, 95% CI: 1.34–2.21, p < 0.001) after adjusting for cancer type, age at diagnosis, and sex. Five-year post metastatic diagnosis survival analysis showed a significant improvement in OS (median survival 28 and 43 months in pre-2015 and post-2015 metastatic diagnosis, respectively, p = 0.00021) based on the year of metastatic diagnosis. Although the use of targeted therapies after metastatic diagnosis prolonged OS, the benefit was not statistically significant. However, longer five-year progression-free survival (PFS) was significantly associated with targeted therapy use (median 10.9 months (CI: 9.7–11.9 months) compared to 9.1 months (CI: 8.1–10.1 months) for non-targeted therapy, respectively, p = 0.0029). Our results provide a clinically relevant overview of the complex molecular landscape and survival mechanisms in metastatic solid cancers.

1. Introduction

Despite recent advances in targeted therapies, metastatic cancer continues to be one of the main causes of morbidity and mortality worldwide. Cancer results from the accumulation of mutations. The phenotypic features and survival outcomes of cancer can be linked to the mutational clone that enables it to proliferate, invade, and overcome its immune environment [1]. These mutational clones lead to heterogeneity, contributing to the challenges in selecting targeted treatment. A better understanding of the mutational landscape and personalized treatment is necessary. The advent of precision medicine has led to the identification of actionable driver mutations regardless of tumor site or origin [2]. The presence of this Darwinian dynamic drives cancer evolution in the process of natural selection of “driver” somatic mutations that give rise to persisting epigenetic changes responsible for proliferation and tumorigenesis [3,4]. The identification and the classification of these mutational drivers has unlocked the potential to therapeutically inhibit this process of perpetual autonomous expansion through targeted therapy [5,6,7,8,9].
The major NCI-MATCH study showed that molecular profiling was successful in 93% of cases, and an actionable alteration was identified in almost 40% of patients [2]. These results affirm the promise of precision medicine to transform clinical approaches from an organ-focused model to a comprehensive mutational signature model of cancer care. However, the application of targeted therapy across solid tumors has not been previously described in a real-world scenario. In this study, we have assessed and characterized the genomic and therapeutic landscape of a cohort of solid metastatic tumors at a single academic site. We performed a retrospective cohort analysis of 957 patients from seven major cancer groups who received next-generation sequencing (NGS) to investigate the relationship between alterations and overall survival (OS). The utilization of a single panel, FoundationOne CDx®, allowed us to independently evaluate the heterogeneity of our cohort, as well as to understand the alteration frequencies and their involvement in cancer diagnosis, progression, and OS.

2. Materials and Methods

2.1. Patient Characteristics

Patients (N = 1001) with solid tumors who underwent a commercial NGS molecular testing panel (Foundation Medicine, Cambridge, MA, USA) at City of Hope (COH) between 2010 and 2018 were evaluated. FoundationOne CDx® assay is a clinically and analytically valid somatic mutation detection assay approved by FDA in November 2017 for all solid advanced tumors. Using their proprietary DNA extraction method from formalin-fixed paraffin-embedded (FFPE) tumor tissue samples, FoundationOne CDx® can detect single nucleotide variants (SNVs), insertion-deletions (indels), and copy number alterations (CNAs) in 324 genes, as well as select rearrangements in 36 genes, microsatellite instability (MSI), and tumor mutational burden (TMB). Unlike hotspot tests, FoundationOne CDx® offers comprehensive genomic profiling of the included genes. Therefore, the risk of selection bias of actionable mutation in this study is minimal with this assay for the solid cancer types analyzed. The COH patients with solid cancer with available FoundationOne CDx® testing reports were eligible for evaluation in this study. The inclusion criteria for the current study include metastatic disease of solid cancers and availability of FoundationOne CDx® test results after metastatic diagnosis. Cancer types that had less than 15 samples or samples that had an unknown date of metastatic disease diagnosis were excluded from the analysis. The majority of the specimens were collected from the primary sites as shown in Supplemental Table S1. Based on the inclusion-exclusion criteria, 957 patients were eligible for the study (Supplemental Figure S1). The clinical data collected from the electronic health record (EHR) included sex, age at diagnosis, date of metastatic diagnosis, race, stage, vital status at last contact, date of death/last contact, histology, cancer type, specimen type, date of molecular testing sample collection and results, first-line treatments, progression-free survival (PFS) and OS. The study was approved by the COH Institutional Review Board (IRB #18038) and the Department of Health and Human Services. The study was conducted according to the standards of Good Clinical Practice, the Declaration of Helsinki, and the US Department of Health and Human Services. All patients were evaluated for their written consent or exemption of consent according to requirements of the IRB and ethics committees, and all recorded data for this study were deidentified.

2.2. Statistical Analysis

Patient demographics and somatic mutations were summarized using descriptive statistics. Wilcoxon rank-sum test was used to compare the frequencies of mutations in two independent metastatic cohorts, specifically the COH pan-cancer and the University of Michigan 2017 metastatic cancer cohorts [10]. We evaluated mutation frequencies, co-occurrence, mutual exclusivity, and prognostic values. Cox regression models were used for univariate and multivariate analysis of OS endpoints with complete observations for clinical variables considered. For the mutation-based survival analysis, somatic mutations in 27 genes with 95 or more occurrences in the cohort were considered, and p ˂ 0.0017 (Bonferroni corrected for multiple testing) was considered to be significant for survival analyses using the Kaplan–Meier log-rank test. Mutations in >35 patient cases were considered for both co-occurrence and mutual exclusivity analysis. For multiple testing corrections, a false discovery rate (FDR) < 0.01 was considered to be statistically significant for analyses of co-occurrence and mutual exclusivity. For all other survival analyses, p < 0.05 was deemed to be significant. For pairwise OS comparison of APC and CDKN2A mutation statuses, Log-Rank test with Bonferroni adjusted p value was used for the significance test. Time to test was calculated in days between the date of metastatic diagnosis and FoundationOne CDx® report date. Student t-test was used to evaluate times to test between pre-2015 and 2015 or later metastatic diagnosis. OS days were calculated from the date of metastatic diagnosis to the date of last follow-up, while PFS days were calculated from the date of metastatic diagnosis to the date of progression after first-line interventions. Data regarding PFS was obtained from the electronic medical record based on the primary oncologist’s notes and assessment of progression. We also compared our cohort to the mutation case frequencies from the Genomic Data Commons (GDC) and 500 adult patients with solid metastatic tumors [10,11].
The R packages used were ComplexHeatmap, GenVisR for heatmap analyses, discover for mutation co-occurrence and mutual exclusivity, and survival and survminer for OS analysis [12,13,14,15,16,17]. All analyses were performed using R version 3.6.2 [18].

3. Results

3.1. Clinical and Genomic Features

Of the 957 patients included in the analysis, 51.9% were female, and 69.7% were Caucasian. The median age at diagnosis was 60 years (IQR: 50–69 years) (Table 1). Seven major cancer diagnosis were included-breast (9.1%), colorectal (33.6%), gastrointestinal (GI; 12.9%), genitourinary (GU; 12.3%), gynecological (GYN; 6.2%), head and neck (H&N; 2%), and thoracic (23.9%) cancers. Colon adenocarcinoma (274, 28.6%) and lung adenocarcinoma (171, 17.9%) were the most common cancer subtypes in the cohort. A total of 390 (40.6%) patients harbored at least one targetable somatic mutation in actionable genes—AKT1, ATM, ALK, BRAF, BRCA1, BRCA2, EGFR, ERBB2, ESR1, FGF3, FGFR1, FGFR2, FGFR3, MET, NTRK1, NTRK2, NTRK3, PIK3CA, ROS1, and RET (NCCN colon, esophagus, head and neck, kidney, lung, ovary, prostate, urinary tract, and uterus guidelines; accessed on 13 December 2020) [19]. The most frequently altered gene in the COH pan-cancer cohort was TP53 (57.3%), followed by APC (31%) and KRAS (27%). In the COH cohort, TP53 alterations were the most prevalent in all cancer types except GU cancers (Figure 1, Supplemental Table S2). Colorectal cancer had the highest frequencies of TP53, APC, and KRAS mutations (76.4%, 73.9%, and 47.8%, respectively). H&N cancer had the highest frequencies of LRP1B (31.6%), and breast cancer had the highest frequency of PIK3CA and MYC (32.2% and 28.7%, respectively). Additionally, the COH pan-cancer cohort and the University of Michigan 2017 metastatic cancer cohort were significantly different within the top eleven common mutation frequencies (Figure 1B, Wilcoxon rank-sum test p 0.016). TP53, KRAS, and APC mutation frequencies were higher in the COH compared to the University of Michigan cohort [10]. In comparison with the GDC cancer types, TP53 was the most frequently altered mutation (35%). However, TP53 mutations occurred at a much higher rate in the COH cohort than the University of Michigan metastatic cohort (49%) and other primary tumor cohorts of TCGA (41%), TCGA 27.0 (36.6%), GENIE 8.1 (39%), and PCAWG (36.9%) [20,21,22,23]. The second most prevalent altered gene in our cohort was APC (31%), which was over 20% higher than the next highest dataset when compared to the University of Michigan metastatic cohort (0.07%) and mixed cohorts from TCGA (8%), TCGA 27.0 (8.1%), GENIE 8.1 (10.3%), and PCAWG (3.4%) [20,21,22,23]. The next most prevalent altered gene was KRAS (27%), which was 12% higher than TCGA (7%), TCGA 27.0 (9.8%), GENIE 8.1 (15.1%), and PCAWG (11.1%) [20,21,22,23].

3.2. Co-Occurrence and Mutual-Exclusivity in Pan-Cancer Cohort

There were 147 somatic mutations detected with more than 35 mutation cases within the cohort. Seven gene-pairs were found to be mutually exclusive and 51 gene-pairs to be co-occurring (Supplemental Table S3). Of the mutually exclusive pairs, KRAS was mutually exclusive with EGFR, ERBB2, and VHL. APC was mutually exclusive with CCDN1 and CDKN2A. Of these seven pairs, a groupwise mutual exclusivity analysis with eight genes (KRAS, EGFR, ERBB2, VHL, RB1, APC, CDKN2B, CCND1) were also mutually exclusive (p < 0.0001, Figure 2). Alternatively, ARFRP1, ASXL1, AURKA, BCL2L1, GNAS, SRC, TOP1, ZNF217 were the most commonly co-occurring, with each co-occurring with seven other genes (FDR < 0.01). Of the actionable genes, BRCA2, ERBB2, and MYC co-occurred with CDK8, CDK12, and RUNX1T1, respectively.

3.3. Survival Analysis

The COH metastatic solid cancer patients had a median PFS and OS of 10 months and 27 months, respectively. Out of 957 patients, 831 patients progressed despite the first-line interventions with pharmacotherapy, surgery, or radiation and 44 patients (21 colorectal, 4 GI, 7 GU, 2 gynecological, 10 thoracic cancer patients) did not progress. The remaining 82 patients had an unknown status of progression. The presence of APC (HR: 0.68, 95% CI: 0.54–0.85, p < 0.0001) and CDKN2A (HR: 2, 95% CI: 1.6–2.6, p < 0.0001) was significant predictors of the longer and shorter OS, respectively. However, only CDKN2A status remained significant (HR: 1.72, 95% CI: 1.34–2.21, p < 0.001, Supplemental Figure S2) after adjusting for age at metastatic diagnosis, gender, and cancer type in the multivariate Cox regression models. Due to the small number of patients with stable disease or remission, PFS was not significantly associated with any somatic mutation status. Patients with APC+/CDKN2A- metastatic cancer survived significantly longer than the patients who are APC-/CDKN2A+ (median survival 52 months compared to 20 months, p < 0.0001, Figure 3). Although APC and CDKN2A co-mutations are rare and have contrasting effects on OS, patients with APC mutation showed statistically nonsignificant longer OS in the presence of CDKN2A (median survival 34 months with APC+/CDKN2A+ mutations compared to 20 months with APC-/CDKN2A+ mutations).

3.4. Effect of Molecular Testing and First-Line Treatments

The median time from initial metastatic diagnosis to molecular testing was 265 days (IQR: 75–639 days) with a significantly accelerated time to testing since 2015 (mean 200.9 days compared to 890.6 days prior to 2015, p < 0.0001). The median time to reporting from sample collection was 115 days (IQR 41–401 days) and was significantly associated with time to testing (Pearson’s correlation: 0.46; p < 0.0001). Five-year post metastatic diagnosis survival analysis showed a significant improvement in OS (median survival 28 and 43 months in pre-2015 and post-2015 metastatic diagnosis, respectively, p 0.00021) by the year of metastatic diagnosis potentially due to faster time to testing and availability of targeted therapy based on the test results (Supplemental Figure S3).
First-line treatment information was available for 941 (98.3%) patients, with 387 (40.4%) patients receiving targeted therapies as a first-line treatment after metastatic diagnosis. Of the known first-line treatments received by these patients, 5-FU, leucovorin, oxaliplatin, bevacizumab, irinotecan, capecitabine, and carboplatin were among the most frequently used chemotherapies across all seven cancer types (Figure 4). These drugs were most commonly used in colorectal and GI cancers. Carboplatin and pemetrexed regimen were most frequently administered in thoracic cancer. Although patients who received a targeted therapy as a first-line treatment after metastatic diagnosis tended to have a longer 5-year OS (45 months compared to 37 months for patients who received traditional therapies), the survival difference was not statistically significant (p = 0.072, Figure 5A). However, PFS was significantly longer for patients who received targeted therapies compared to patients who received traditional therapies such as chemotherapies, surgery, or radiation (median PFS 10.9 months [CI: 9.7–11.9 months] compared to 9.1 months [CI: 8.1–10.1 months], respectively, p 0.0029, Figure 5B). Use of targeted therapies as first-line treatment remained significantly associated with PFS (HR 0.81; 95% CI: 0.69–0.94; p 0.007) after adjusting for patient’s cancer type, age at diagnosis, and sex.

4. Discussion

Treatment of metastatic cancer continues to be challenging, resulting in high mortality among all cancer patients [24]. Here, we analyzed 957 metastatic solid cancer patients treated at the COH. In comparison with other published pan-cancer datasets, including The Cancer Genome Atlas (TCGA), the Pan-Cancer Analysis of Whole Genomes (PCAWG), and the recent data releases for the Genomics Evidence Neoplasia Information Exchange (GENIE) 8.1 and TCGA 27.0, there was significant variation in top somatic mutation frequencies of TP53, APC, KRAS, and LRP1B genes [20,21,22,23]. One possible explanation for the differences in frequency is the metastatic status of our cohort compared to the mixed cancer patient population from TCGA 27.0 and GENIE. Further, there were only 957 patients in our COH cohort, while the most recent data releases for TCGA 27.0 and GENIE 8.1 were for 84,392 and 95,918 primary tumor samples, respectively [21,22,25]. Also, we analyzed breast, colorectal, gastrointestinal, genitourinary, gynecological, head and neck, and thoracic major cancer types. Colorectal (322, 33.6%) and thoracic (229, 23.9%) accounted for 57.5% of all cancer diagnosis in our cohort, and TP53, APC, KRAS, and LRP1B mutation rates are higher in colorectal cancers and non-small cell lung cancers (NSCLC) [26,27,28,29]. APC gene mutations are known to be associated with chromosomal instability, however, similar to our findings several studies have identified that mutated APC was associated with improved survival as compared to wild-type APC in colorectal cancer [30,31,32,33,34].
Mutual exclusivity and co-occurrence among oncogenic driver mutations warrant careful discussion when considering targeted therapies and treatment plans. While biology drives mutual exclusivity, most co-occurrences are most likely by chance [14]. Previous studies have reported the prevalence of concomitant driver mutations among the different patient populations. In our study, the seven mutually exclusive gene pairs with FDR < 0.01 are KRAS/EGFR, CDKN2A/APC, KRAS/ERBB2, VHL/TP53, RB1/CDKN2B, VHL/KRAS, and CCND1/APC. A well-known pair of genes that validates our mutual exclusion data is KRAS with mutated EGFR loci, most commonly seen in lung cancer [35]. The EGFR family of receptors EGFR and HER2 (ERBB2) act as upstream regulators of KRAS in the Ras/Raf/MEK/ERK pathway. Concurrent mutations in EGFR or ERBB2 and KRAS are not selected for and may be considered redundant in tumor cells, thus selective pressures may account for the mutual exclusivity of EGFR/ERBB2 and KRAS [36,37,38]. However, recent publications have shown small cohorts of patients that harbor the rare EGFR/KRAS co-mutations. Zhuang et al. reported a comprehensive study analyzing concomitant driver gene mutations in 3774 NSCLC Chinese patients harboring EGFR, ALK, ROS1, KRAS, and BRAF alterations. Only 63 patients (1.7%) harbored mutations in two or three of these genes and among these patients, EGFR/KRAS was the most frequent co-alteration (31.7%) [39]. Pesek et al. reported four clinical cases of NSCLC patients with confirmed co-mutations in four of the five clinical cases presented [40]. Additionally, our analysis confirmed that the presence of TP53 mutations was mutually exclusive with VHL mutations. Previously, it was found that mutual exclusivity among the mutations depends on different histological subtypes in renal cell carcinomas (RCCs) [41]. Szymańska et al. reported 77% of clear-cell carcinomas from a multi-center case-control study harbor VHL mutations; in contrast, only 22% of RCCs of other histological types harbor VHL mutations [42]. Upon stratification, clear cell tumors showed mutual exclusivity among TP53 and VHL mutations, whereas non-clear cell tumors had TP53 mutations in the absence of VHL mutations.
We reported the following actionable gene pairs as co-occurring—ERBB2/CDK12, BRCA2/CDK8, and MYC/RUNX1T1 with additional gene pairs, including PREX2/RUNX1T1 and STK11/KEAP1. Somatic mutations that co-occur are common across various cancer types. For example, CDK12 is an essential transcription-associated CDK responsible for DNA damage response and is often found in the ERBB2 amplicon [43,44,45]. Capra et al. reported that CDK12 gene amplification can often co-occur with ERBB2 amplification and may act as an oncogenic driver in HER2-positive breast cancer [46,47,48]. Wilson et al. reported similar co-occurring gene pair mutations among breast cancer subtypes (HER2+, hormone receptor positive [HR+], and triple-negative breast cancer [TNBC]) [49]. Interestingly, the commonly occurring gene pairs MYC/RUNX1T1, PREX2/RUNX1T1, PRKDC/PREX2, and MYC/PREX2 in our study are coamplified genes on the same loci, within the 8q loci, and reported by Wilson et al. in 13% of patients with HR+, HER2- disease and 19% with TNBC [49]. Lastly, our analysis showed STK11/KEAP1 co-occurring mutations (FDR < 0.01). Concurrent mutations in STK11/KEAP1 act as a major driver in primary resistance to PD1 blockade in KRAS-mutated lung adenocarcinoma [50]. Papillon-Cavanagh et al. investigated the utility of STK11/KEAP1 mutations as predictive biomarkers in patients diagnosed with stage IIIB, IIIC, IVA, or IVB non-squamous NSCLC. It was determined that the co-occurring mutations were associated with poor prognosis in both anti-PD-1/PD-L1 treated and chemotherapy-treated populations [51].
Rapid advancements in NGS technology have allowed physicians to understand in real-time the genetic makeup of a patient’s tumor and track these alterations over time [52,53,54]. Groisberg et al. reported the identification of an oncogenic mutation in at least 92% of patients [55]. In a more comprehensive analysis of 10,000 patients, Zehir et al. reported 91% successful sequencing and identified 36.7% of patients with at least one actionable alteration [56]. These rates of actionable alterations were previously estimated to be lower, but recent studies show that actionability in solid tumors can vary between 35–80% [57,58,59,60,61]. This has been compounded by the recent results from the NCI-match study that sequenced 93% of 5954 patients and identified actionable alterations in 37.8% of cases [2]. Our results suggest a similar frequency to other NGS-guided studies that match targeted therapy based on the identification of at least one actionable alteration [62,63]. The advent of novel therapeutics, including new NTRK, RAS, RET, and MET inhibitors has transformed the clinical trial landscape and presented more therapeutic options to patients than ever before [64,65,66,67]. These advancements in genomic sequencing and trial accrual validate precision medicine to deliver improved outcomes in all solid tumors.
It is worth noting that our study has several limitations. First, this is a retrospective cohort study using seven solid metastatic tumor patients from a single center. Therefore, the study results may not be generalizable to other metastatic cohorts. Second, survival may depend on various factors beyond the age of diagnosis, sex, cancer type, and treatments used. Also, treatment options are different depending on cancer types, and the efficacy of treatments may vary significantly based on how soon in the disease progression they were initiated. Therefore, our survival assessment was oversimplified. Third, we did not capture clonal heterogeneity and evaluate the effect of intratumor heterogeneity on OS. Future studies need to address these limitations in the evaluation of somatic mutations in solid metastatic pan-cancer patients.

5. Conclusions

The promise of precision medicine relies on identifying mutations in patients and developing novel therapeutics for these drug targets. Beyond actionable mutations, gene expression and epigenetics data may further elucidate critical pathway activation and identify potential targeted therapies [68,69]. Although gene expression and epigenetics data are not routinely implemented in clinical decision-making, using these in collaboration with mutation data will shed light on the selection of optimal targeted therapies and potentially avoid resistance development. Targeted therapies are becoming more common and the efficacy of these agents has been proven in several solid tumors [67,70,71,72,73,74,75,76]. These agents provide a measurable outcome advantage compared to other therapeutic options, including chemotherapy and immunotherapy—with the longest median OS of 7 years in patients with ALK-mutated metastatic NSCLC, an achievement unmatched by any other therapeutic option [77]. However, the primary barrier to targeted therapeutics remains acquired resistance due to evolutionary competition between resistant and sensitive subclones, which leads to survival of persistent subclones that alter tumor behavior and outcomes [78,79,80,81]. Our results showed the application of targeted therapies improved PFS, potentially owing to early testing and availability of increasing number of targeted therapies. Combinatory targeted therapeutic options are a potential non-toxic avenue for combatting this resistance, and several early and late-stage clinical trials, including NCI #10327, NCI #9466, and NCT03944772, are currently underway to understand their efficacy. The promise of precision medicine is quickly transforming into the reality of precision therapeutics with the discovery of effective KRAS-inhibitors and other drug targets that were previously thought to be undruggable [65,71]. Therefore, the outcomes of patients with solid tumors rely on the successful inhibition of oncogenic mutational signatures that drive disease progression by utilizing informed targeted therapy.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cancers13112776/s1: Figure S1: Study Flow Diagram Showing Metastatic Pan-cancer Patient Selection for Analysis. A total of 1001 patients with available FoundationOne CDx® test results were evaluated. Three patients with unknown cancer types, one patient with skin cancer and forty patients without test results after metastatic diagnosis were excluded, resulting in 957 eligible patients for the study analysis. Figure S2: Forest plot of variables analyzed for overall survival. Hazard ratios of each variable along with the p-value are shown. Mutations in CDKN2A gene was associated with significantly worse survival after adjusting for age at diagnosis, sex, APC mutations and cancer type. Figure S3: Kaplan-Meier overall survival analysis in solid pan-cancer metastatic patients five year after metastatic diagnosis. Patients who had the metastatic diagnosis before 2015 had significantly shorter overall survival compared to those diagnosed in 2015 or later. Table S1: Summary of the specimen site classifications for the various cancer subgroups. Table S2: Prevalence in cases (%) of the top thirty mutated genes across seven cancer types. Table S3: List of all mutually exclusive and co-occurring gene pairs with false discovery rate (FDR) < 0.01.

Author Contributions

M.R.R., Conceptualization, Data Curation, Formal Analysis, Methodology, Software, Validation, Visualization, Writing—Original Draft; I.M., Conceptualization, Data Curation, Formal Analysis, Methodology, Validation, Visualization, Writing—Original Draft; R.P. (Rebecca Pharaon), Data Curation, Formal Analysis, Validation, Writing—Review and Editing; J.F., Data Curation, Formal Analysis, Validation, Writing—Review and Editing; A.R.B., Data Curation, Formal Analysis, Validation, Writing—Review and Editing; J.C., Data Curation, Formal Analysis, Writing—Review and Editing; C.C., Data Curation; Formal Analysis, Writing—Review and Editing; M.W.N., Formal Analysis, Methodology, Writing—Review and Editing; R.C.-V., Formal Analysis, Methodology, Writing—Review and Editing; M.J., Formal Analysis, Methodology, Writing—Review and Editing; L.S., Data Curation, Formal Analysis, Methodology, Writing—Review and Editing; S.E.Y., Data Curation, Formal Analysis, Writing—Review and Editing; K.L.R., Data Curation, Formal Analysis, Validation, Writing—Review and Editing; R.P. (Raju Pillai), Data Curation, Methodology, Validation, Writing—Review and Editing; L.A., Data Curation, Formal Analysis, Validation, Writing—Review and Editing; M.A., Data Curation, Formal Analysis, Validation, Writing—Review and Editing; E.W.W., Data Curation, Formal Analysis, Writing—Review and Editing; V.C., Data Curation, Validation, Writing—Review and Editing; M.C., Data Curation, Validation, Writing—Review and Editing; M.F., Data Curation, Validation, Writing—Review and Editing; M.K., Data Curation, Validation, Writing—Review and Editing; E.M., Data Curation, Validation, Writing—Review and Editing; J.M., Data Curation, Validation, Writing—Review and Editing; Y.Y., Data Curation, Formal Analysis, Validation, Writing—Review and Editing; S.K.B., Formal Analysis, Validation, Writing—Review and Editing; S.P., Conceptualization, Data Curation, Formal Analysis, Validation, Writing—Original Draft; R.S., Conceptualization, Funding, Data Curation, Formal Analysis, Validation, Writing—Original Draft. All authors have read and agreed to the published version of the manuscript.

Funding

The work was supported by the National Cancer Institute of the National Institutes of Health under award numbers P30CA033572, U54CA209978, R01CA247471, and R01CA218545.

Institutional Review Board Statement

The study was approved by the COH Institutional Review Board (IRB #18038) and the Department of Health and Human Services. The study was conducted according to the standards of Good Clinical Practice, the Declaration of Helsinki, and the US Department of Health and Human Services.

Informed Consent Statement

All patients were evaluated for their written consent or exemption of consent according to requirements of the IRB and ethics committees, and all recorded data for this study were deidentified.

Data Availability Statement

All of the pertinent data and codes for the study is available at https://github.com/mumtahena/COH_metastatic_solid (accessed on 4 January 2020).

Acknowledgments

The results shown here are in whole or partly based upon data generated by the TCGA Research Network: https://www.cancer.gov/tcga (accessed on 4 January 2020).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hanahan, D.; Weinberg, R.A. Hallmarks of Cancer: The Next Generation. Cell 2011, 144, 646–674. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Flaherty, K.T.; Gray, R.J.; Chen, A.P.; Li, S.; McShane, L.M.; Patton, D.; Hamilton, S.R.; Williams, P.M.; Iafrate, A.J.; Sklar, J.; et al. Molecular Landscape and Actionable Alterations in a Genomically Guided Cancer Clinical Trial: National Cancer Institute Molecular Analysis for Therapy Choice (NCI-MATCH). J. Clin. Oncol. 2020, 38, 3883–3894. [Google Scholar] [CrossRef] [PubMed]
  3. Flavahan, W.A.; Gaskell, E.; Bernstein, B.E. Epigenetic plasticity and the hallmarks of cancer. Science 2017, 357, eaal2380. [Google Scholar] [CrossRef] [Green Version]
  4. Vaz, M.; Hwang, S.Y.; Kagiampakis, I.; Phallen, J.; Patil, A.; O’Hagan, H.; Murphy, L.; Zahnow, C.A.; Gabrielson, E.; Velculescu, V.E.; et al. Chronic Cigarette Smoke-Induced Epigenomic Changes Precede Sensitization of Bronchial Epithelial Cells to Single-Step Transformation by KRAS Mutations. Cancer Cell 2017, 32, 360–376.e6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Nik-Zainal, S.; Davies, H.; Staaf, J.; Ramakrishna, M.; Glodzik, D.; Zou, X.; Martincorena, I.; Alexandrov, L.B.; Martin, S.; Wedge, D.; et al. Landscape of somatic mutations in 560 breast cancer whole-genome sequences. Nat. Cell Biol. 2016, 534, 47–54. [Google Scholar] [CrossRef] [PubMed]
  6. Petljak, M.; Alexandrov, L.B. Understanding mutagenesis through delineation of mutational signatures in human cancer. Carcinogenesis 2016, 37, 531–540. [Google Scholar] [CrossRef] [Green Version]
  7. Hayward, N.; Wilmott, J.; Waddell, N.; Johansson, P.A.; Field, M.A.; Nones, K.; Patch, A.-M.; Kakavand, H.; Alexandrov, L.B.; Burke, H.; et al. Whole-genome landscapes of major melanoma subtypes. Nat. Cell Biol. 2017, 545, 175–180. [Google Scholar] [CrossRef] [PubMed]
  8. Polak, P.; Kim, J.; Braunstein, L.Z.; Karlic, R.; Haradhavala, N.J.; Tiao, G.; Rosebrock, D.; Livitz, D.; Kübler, K.; Mouw, K.W.; et al. A mutational signature reveals alterations underlying deficient homologous recombination repair in breast cancer. Nat. Genet. 2017, 49, 1476–1486. [Google Scholar] [CrossRef]
  9. Stratton, M.R.; Campbell, P.J.; Futreal, P.A. The cancer genome. Nat. Cell Biol. 2009, 458, 719–724. [Google Scholar] [CrossRef] [Green Version]
  10. Robinson, D.R.; Wu, Y.-M.; Lonigro, R.J.; Vats, P.; Cobain, E.; Everett, J.; Cao, X.; Rabban, E.; Kumar-Sinha, C.; Raymond, V.; et al. Integrative clinical genomics of metastatic cancer. Nat. Cell Biol. 2017, 548, 297–303. [Google Scholar] [CrossRef]
  11. Grossman, R.L.; Heath, A.P.; Ferretti, V.; Varmus, H.E.; Lowy, D.R.; Kibbe, W.A.; Staudt, L.M. Toward a Shared Vision for Cancer Genomic Data. N. Engl. J. Med. 2016, 375, 1109–1112. [Google Scholar] [CrossRef]
  12. Gu, Z.; Eils, R.; Schlesner, M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics 2016, 32, 2847–2849. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Skidmore, Z.; Wagner, A.H.; Lesurf, R.; Campbell, K.M.; Kunisaki, J.; Griffith, O.; Griffith, M. GenVisR: Genomic Visualizations in R. Bioinformatics 2016, 32, 3012–3014. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Canisius, S.; Martens, J.W.M.; Wessels, L.F.A. A novel independence test for somatic alterations in cancer shows that biology drives mutual exclusivity but chance explains most co-occurrence. Genome Biol. 2016, 17, 1–17. [Google Scholar] [CrossRef] [Green Version]
  15. Zhao, H.; Chen, K.-Z.; Hui, B.-G.; Zhang, K.; Yang, F.; Wang, J. Role of circulating tumor DNA in the management of early-stage lung cancer. Thorac. Cancer 2018, 9, 509–515. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Kassambara, A.K.; Biecek, M. P Survminer: Drawing Survival Curves Using ’ggplot2′, R package version 0.4.6; 2019. Available online: https://rpkgs.datanovia.com/survminer/index.html (accessed on 8 August 2020).
  17. Therneau, T. A Package for Survival Analysis in R, R package version 3.1-12; 2020. Available online: https://cran.r-project.org/web/packages/survival/index.html (accessed on 8 July 2020).
  18. Team, R.C. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019. [Google Scholar]
  19. NCCN® Clinical Practical Guidelines in Oncology. Available online: https://www.nccn.org/professionals/physician_gls/default.aspx (accessed on 13 December 2020).
  20. Hoadley, K.A.; Yau, C.; Wolf, D.M.; Cherniack, A.D.; Tamborero, D.; Yang, T.-H.O.; Leiserson, M.D.; Niu, B.; McLellan, M.D.; Uzunangelov, V.; et al. Multiplatform Analysis of 12 Cancer Types Reveals Molecular Classification within and across Tissues of Origin. Cell 2014, 158, 929–944. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Cerami, E.; Gao, J.; Dogrusoz, U.; Gross, B.E.; Sumer, S.O.; Aksoy, B.A.; Jacobsen, A.; Byrne, C.J.; Heuer, M.L.; Larsson, E.; et al. The cBio Cancer Genomics Portal: An Open Platform for Exploring Multidimensional Cancer Genomics Data. Cancer Discov. 2012, 2, 401–404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Gao, J.; Aksoy, B.A.; Dogrusoz, U.; Dresdner, G.; Gross, B.; Sumer, S.O.; Sun, Y.; Jacobsen, A.; Sinha, R.; Larsson, E.; et al. Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the cBioPortal. Sci. Signal. 2013, 6, pl1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Campbell, P.J.; Getz, G.; Korbel, J.O.; Stuart, J.M.; Jennings, J.L.; Stein, L.D.; Perry, M.D.; Nahal-Bose, H.K.; Ouellette, B.F.F.; Li, C.H.; et al. Pan-cancer analysis of whole genomes. Nature 2020, 578, 82–93. [Google Scholar] [CrossRef] [Green Version]
  24. Siegel, R.L.; Miller, K.D.; Jemal, A. Cancer statistics. CA Cancer J. Clin. 2020, 70, 7–30. [Google Scholar] [CrossRef]
  25. Chen, F.; Zhang, Y.; Varambally, S.; Creighton, C.J. Molecular Correlates of Metastasis by Systematic Pan-Cancer Analysis Across the Cancer Genome Atlas. Mol. Cancer Res. 2019, 17, 476–487. [Google Scholar] [CrossRef] [Green Version]
  26. Donehower, L.A.; Soussi, T.; Korkut, A.; Liu, Y.; Schultz, A.; Cardenas, M.; Li, X.; Babur, O.; Hsu, T.-K.; Lichtarge, O.; et al. Integrated Analysis of TP53 Gene and Pathway Alterations in The Cancer Genome Atlas. Cell Rep. 2019, 28, 1370–1384.e5. [Google Scholar] [CrossRef] [Green Version]
  27. Zhang, L.; Shay, J.W. Multiple Roles of APC and its Therapeutic Implications in Colorectal Cancer. J. Natl. Cancer Inst. 2017, 109, djw332. [Google Scholar] [CrossRef] [Green Version]
  28. Prior, I.A.; Lewis, P.D.; Mattos, C. A Comprehensive Survey of Ras Mutations in Cancer. Cancer Res. 2012, 72, 2457–2467. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Chen, H.; Chong, W.; Wu, Q.; Yao, Y.; Mao, M.; Wang, X. Association of LRP1B Mutation with Tumor Mutation Burden and Outcomes in Melanoma and Non-small Cell Lung Cancer Patients Treated with Immune Check-Point Blockades. Front. Immunol. 2019, 10, 1113. [Google Scholar] [CrossRef] [Green Version]
  30. Fodde, R.; Kuipers, J.; Rosenberg, C.; Smits, R.; Kielman, M.; Gaspar, C.; Van Es, J.H.; Breukel, C.; Wiegant, J.; Giles, R.H.; et al. Mutations in the APC tumour suppressor gene cause chromosomal instability. Nat. Cell Biol. 2001, 3, 433–438. [Google Scholar] [CrossRef] [PubMed]
  31. Wang, C.; Ouyang, C.; Cho, M.; Ji, J.; Sandhu, J.; Goel, A.; Kahn, M.; Fakih, M. Wild-type APC Is Associated with Poor Survival in Metastatic Microsatellite Stable Colorectal Cancer. Oncologist 2021, 26, 208–214. [Google Scholar] [CrossRef] [PubMed]
  32. Schell, M.J.; Yang, M.; Teer, J.K.; Lo, F.Y.; Madan, A.; Coppola, D.; Monteiro, A.N.A.; Nebozhyn, M.V.; Yue, B.; Loboda, A.; et al. A multigene mutation classification of 468 colorectal cancers reveals a prognostic role for APC. Nat. Commun. 2016, 7, 11743. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Jorissen, R.N.; Christie, M.; Mouradov, D.; Sakthianandeswaren, A.; Li, S.; Love, C.; Xu, Z.-Z.; Molloy, P.L.; Jones, I.T.; McLaughlin, S.; et al. Wild-type APC predicts poor prognosis in microsatellite-stable proximal colon cancer. Br. J. Cancer 2015, 113, 979–988. [Google Scholar] [CrossRef] [PubMed]
  34. Yaeger, R.; Chatila, W.K.; Lipsyc, M.D.; Hechtman, J.F.; Cercek, A.; Sanchez-Vega, F.; Jayakumaran, G.; Middha, S.; Zehir, A.; Donoghue, M.T.A.; et al. Clinical Sequencing Defines the Genomic Landscape of Metastatic Colorectal Cancer. Cancer Cell 2018, 33, 125–136. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Ochoa, S.; Martínez-Pérez, E.; Zea, D.J.; Molina-Vila, M.A.; Marino-Buslje, C. Co-mutation and exclusion analysis in human tumors, a means for cancer biology studies and treatment design. bioRxiv 2017, 182501. [Google Scholar] [CrossRef] [Green Version]
  36. Bae, N.C.; Chae, M.H.; Lee, M.H.; Kim, K.M.; Lee, E.B.; Kim, C.H.; Park, T.-I.; Han, S.B.; Jheon, S.; Jung, T.H.; et al. EGFR, ERBB2, and KRAS mutations in Korean non-small cell lung cancer patients. Cancer Genet. Cytogenet. 2007, 173, 107–113. [Google Scholar] [CrossRef]
  37. Bradshaw, R.; Chalkley, R.; Biarc, J.; Burlingame, A.L. Receptor tyrosine kinase signaling mechanisms: Devolving TrkA responses with phosphoproteomics. Adv. Biol. Regul. 2013, 53, 87–96. [Google Scholar] [CrossRef] [Green Version]
  38. Fitzgerald, T.L.; Lertpiriyapong, K.; Cocco, L.; Martelli, A.M.; Libra, M.; Candido, S.; Montalto, G.; Cervello, M.; Steelman, L.; Abrams, S.L.; et al. Roles of EGFR and KRAS and their downstream signaling pathways in pancreatic cancer and pancreatic cancer stem cells. Adv. Biol. Regul. 2015, 59, 65–81. [Google Scholar] [CrossRef] [PubMed]
  39. Ju, L.; Han, M.; Zhao, C.; Li, X. EGFR, KRAS and ROS1 variants coexist in a lung adenocarcinoma patient. Lung Cancer 2016, 95, 94–97. [Google Scholar] [CrossRef] [PubMed]
  40. Benesova, L.; Minarik, M.; Jancarikova, D.; Belsanova, B.; Pesek, M. Multiplicity of EGFR and KRAS mutations in non-small cell lung cancer (NSCLC) patients treated with tyrosine kinase inhibitors. Anticancer Res. 2010, 30, 1667–1671. [Google Scholar] [PubMed]
  41. Kandoth, C.; McLellan, M.D.; Vandin, F.; Ye, K.; Niu, B.; Lu, C.; Xie, M.; Zhang, Q.; McMichael, J.F.; Wyczalkowski, M.A.; et al. Mutational landscape and significance across 12 major cancer types. Nature 2013, 502, 333–339. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Szymańska, K.; Moore, L.; Rothman, N.; Chow, W.; Waldman, F.; Jaeger, E.; Waterboer, T.; Foretova, L.; Navrátilová, M.; Janout, V.; et al. TP53, EGFR, and KRAS mutations in relation to VHL inactivation and lifestyle risk factors in renal-cell carcinoma from central and eastern Europe. Cancer Lett. 2010, 293, 92–98. [Google Scholar] [CrossRef] [PubMed]
  43. Ko, T.K.; Kelly, E.; Pines, J. CrkRS: A novel conserved Cdc2-related protein kinase that colocalises with SC35 speckles. J. Cell Sci. 2001, 114, 2591–2603. [Google Scholar] [CrossRef]
  44. Paculova, H.; Kohoutek, J. The emerging roles of CDK12 in tumorigenesis. Cell Div. 2017, 12, 7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Mertins, P.; Cptac, N.; Mani, D.R.; Ruggles, K.V.; Gillette, M.A.; Clauser, K.R.; Wang, P.; Wang, X.; Qiao, J.W.; Cao, S.; et al. Proteogenomics connects somatic mutations to signalling in breast cancer. Nat. Cell Biol. 2016, 534, 55–62. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Capra, M.; Nuciforo, P.G.; Confalonieri, S.; Quarto, M.; Bianchi, M.; Nebuloni, M.; Boldorini, R.; Pallotti, F.; Viale, G.; Gishizky, M.L.; et al. Frequent Alterations in the Expression of Serine/Threonine Kinases in Human Cancers. Cancer Res. 2006, 66, 8147–8154. [Google Scholar] [CrossRef] [Green Version]
  47. Pilarova, K.; Herudek, J.; Blazek, D. CDK12: Cellular functions and therapeutic potential of versatile player in cancer. NAR Cancer 2020, 2, zcaa003. [Google Scholar] [CrossRef] [Green Version]
  48. Sircoulomb, F.; Bekhouche, I.; Finetti, P.; Adélaïde, J.; Ben Hamida, A.; Bonansea, J.; Raynaud, S.; Innocenti, C.; Charafe-Jauffret, E.; Tarpin, C.; et al. Genome profiling of ERBB2-amplified breast cancers. BMC Cancer 2010, 10, 539. [Google Scholar] [CrossRef] [PubMed]
  49. Wilson, T.R.; Udyavar, A.R.; Chang, C.-W.; Spoerke, J.M.; Aimi, J.; Savage, H.M.; Daemen, A.; O’Shaughnessy, J.A.; Bourgon, R.; Lackner, M.R.; et al. Genomic Alterations Associated with Recurrence and TNBC Subtype in High-Risk Early Breast Cancers. Mol. Cancer Res. 2018, 17, 97–108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Skoulidis, F.; Goldberg, M.E.; Greenawalt, D.M.; Hellmann, M.D.; Awad, M.M.; Gainor, J.F.; Schrock, A.B.; Hartmaier, R.J.; Trabucco, S.E.; Gay, L.; et al. STK11/LKB1 Mutations and PD-1 Inhibitor Resistance in KRAS-Mutant Lung Adenocarcinoma. Cancer Discov. 2018, 8, 822–835. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Papillon-Cavanagh, S.; Doshi, P.; Dobrin, R.; Szustakowski, J.; Walsh, A.M. STK11 and KEAP1 mutations as prognostic biomarkers in an observational real-world lung adenocarcinoma cohort. ESMO Open 2020, 5, e000706. [Google Scholar] [CrossRef] [Green Version]
  52. Morganti, S.; Tarantino, P.; Ferraro, E.; D’Amico, P.; Duso, B.A.; Curigliano, G. Next Generation Sequencing (NGS): A Revolutionary Technology in Pharmacogenomics and Personalized Medicine in Cancer. Adv. Exp. Med. Biol. 2019, 1168, 9–30. [Google Scholar] [CrossRef]
  53. Postel, M.; Roosen, A.; Laurent-Puig, P.; Taly, V.; Wang-Renault, S.-F. Droplet-based digital PCR and next generation sequencing for monitoring circulating tumor DNA: A cancer diagnostic perspective. Expert Rev. Mol. Diagn. 2018, 18, 7–17. [Google Scholar] [CrossRef] [PubMed]
  54. Tan, O.; Shrestha, R.; Cunich, M.; Schofield, D. Application of next-generation sequencing to improve cancer management: A review of the clinical effectiveness and cost-effectiveness. Clin. Genet. 2018, 93, 533–544. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Groisberg, R.; Hong, D.S.; Roszik, J.; Janku, F.; Tsimberidou, A.M.; Javle, M.; Meric-Bernstam, F.; Subbiah, V. Clinical Next-Generation Sequencing for Precision Oncology in Rare Cancers. Mol. Cancer Ther. 2018, 17, 1595–1601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Zehir, A.; Benayed, R.; Shah, R.; Syed, A.; Middha, S.; Kim, H.R.; Srinivasan, P.; Gao, J.; Chakravarty, D.; Devlin, S.M.; et al. Mutational landscape of metastatic cancer revealed from prospective clinical sequencing of 10,000 patients. Nat. Med. 2017, 23, 703–713. [Google Scholar] [CrossRef] [PubMed]
  57. Toomey, S.; Carr, A.; Mezynski, M.J.; Elamin, Y.; Rafee, S.; Cremona, M.; Morgan, C.; Madden, S.; Abdul-Jalil, K.I.; Gately, K.; et al. Identification and clinical impact of potentially actionable somatic oncogenic mutations in solid tumor samples. J. Transl. Med. 2020, 18, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Lee, J.W.; Kim, N.K.D.; Lee, S.H.; Cho, H.W.; Ma, Y.; Ju, H.Y.; Yoo, K.H.; Sung, K.W.; Koo, H.H.; Park, W.-Y. Discovery of actionable genetic alterations with targeted panel sequencing in children with relapsed or refractory solid tumors. PLoS ONE 2019, 14, e0224227. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Drilon, A.; Wang, L.; Arcila, M.E.; Balasubramanian, S.; Greenbowe, J.R.; Ross, J.S.; Stephens, P.J.; Lipson, D.; Miller, V.A.; Kris, M.G.; et al. Broad, Hybrid Capture–Based Next-Generation Sequencing Identifies Actionable Genomic Alterations in Lung Adenocarcinomas Otherwise Negative for Such Alterations by Other Genomic Testing Approaches. Clin. Cancer Res. 2015, 21, 3631–3639. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Palmer, G.A.; Miller, V.A.; Curran, J.; Ross, J.S.; Lipson, D.; Yelensky, R.; Stephens, P.; Lancelotta, M.P.; Cronin, M.T. Next-generation sequencing (NGS) to identify actionable genomic changes in common and rare solid tumors: The FMI experience with the initial 50 consecutive patients. J. Clin. Oncol. 2012, 30, 10590. [Google Scholar] [CrossRef]
  61. Miller, V.A.; Ross, J.S.; Wang, K.; Ali, S.M.; Otto, G.; Curran, J.; Palma, N.A.; Yelensky, R.; Downing, S.; Stephens, P.; et al. Use of next-generation sequencing (NGS) to identify actionable genomic alterations (GA) in diverse solid tumor types: The Foundation Medicine (FMI) experience with 2200+ clinical samples. J. Clin. Oncol. 2013, 31, 11020. [Google Scholar] [CrossRef]
  62. Hainsworth, J.D.; Meric-Bernstam, F.; Swanton, C.; Hurwitz, H.; Spigel, D.R.; Sweeney, C.; Burris, H.A.; Bose, R.; Yoo, B.; Stein, A.; et al. Targeted Therapy for Advanced Solid Tumors on the Basis of Molecular Profiles: Results From MyPathway, an Open-Label, Phase IIa Multiple Basket Study. J. Clin. Oncol. 2018, 36, 536–542. [Google Scholar] [CrossRef] [PubMed]
  63. Tsimberidou, A.-M.; Hong, D.S.; Ye, Y.; Cartwright, C.; Wheler, J.J.; Falchook, G.S.; Naing, A.; Fu, S.; Piha-Paul, S.; Janku, F.; et al. Initiative for Molecular Profiling and Advanced Cancer Therapy (IMPACT): An MD Anderson Precision Medicine Study. JCO Precis. Oncol. 2017, 2017, 1–18. [Google Scholar] [CrossRef]
  64. Doebele, R.C.; Drilon, A.; Paz-Ares, L.; Siena, S.; Shaw, A.T.; Farago, A.F.; Blakely, C.M.; Seto, T.; Cho, B.C.; Tosi, D.; et al. Entrectinib in patients with advanced or metastatic NTRK fusion-positive solid tumours: Integrated analysis of three phase 1–2 trials. Lancet Oncol. 2020, 21, 271–282. [Google Scholar] [CrossRef]
  65. Hong, D.S.; Kuo, J.; Sacher, A.G.; Barlesi, F.; Besse, B.; Kuboki, Y.; Dy, G.K.; Dembla, V.; Krauss, J.C.; Burns, T.F.; et al. CodeBreak 100: Phase I study of AMG 510, a novel KRASG12C inhibitor, in patients (pts) with advanced solid tumors other than non-small cell lung cancer (NSCLC) and colorectal cancer (CRC). J. Clin. Oncol. 2020, 38, 3511. [Google Scholar] [CrossRef]
  66. Wolf, J.; Seto, T.; Han, J.-Y.; Reguart, N.; Garon, E.B.; Groen, H.J.; Tan, D.S.; Hida, T.; de Jonge, M.; Orlov, S.V.; et al. Capmatinib inMETExon 14–Mutated orMET-Amplified Non–Small-Cell Lung Cancer. N. Engl. J. Med. 2020, 383, 944–957. [Google Scholar] [CrossRef] [PubMed]
  67. Drilon, A.; Oxnard, G.R.; Tan, D.S.; Loong, H.H.; Johnson, M.; Gainor, J.; McCoach, C.E.; Gautschi, O.; Besse, B.; Cho, B.C.; et al. Efficacy of Selpercatinib in RET Fusion–Positive Non–Small-Cell Lung Cancer. N. Engl. J. Med. 2020, 383, 813–824. [Google Scholar] [CrossRef]
  68. Rahman, M.; MacNeil, S.M.; Jenkins, D.F.; Shrestha, G.; Wyatt, S.R.; McQuerry, J.A.; Piccolo, S.R.; Heiser, L.M.; Gray, J.W.; Johnson, W.E.; et al. Activity of distinct growth factor receptor network components in breast tumors uncovers two biologically relevant subtypes. Genome Med. 2017, 9, 40. [Google Scholar] [CrossRef] [Green Version]
  69. Bennett, R.L.; Licht, J.D. Targeting Epigenetics in Cancer. Annu. Rev. Pharmacol. Toxicol. 2018, 58, 187–207. [Google Scholar] [CrossRef]
  70. Ramalingam, S.S.; Vansteenkiste, J.; Planchard, D.; Cho, B.C.; Gray, J.E.; Ohe, Y.; Zhou, C.; Reungwetwattana, T.; Cheng, Y.; Chewaskulyong, B.; et al. Overall Survival with Osimertinib in Untreated, EGFR-Mutated Advanced NSCLC. N. Engl. J. Med. 2020, 382, 41–50. [Google Scholar] [CrossRef] [PubMed]
  71. Hong, D.S.; Fakih, M.G.; Strickler, J.H.; Desai, J.; Durm, G.A.; Shapiro, G.I.; Falchook, G.S.; Price, T.J.; Sacher, A.; Denlinger, C.S.; et al. KRASG12C Inhibition with Sotorasib in Advanced Solid Tumors. N. Engl. J. Med. 2020, 383, 1207–1217. [Google Scholar] [CrossRef] [PubMed]
  72. Murthy, R.K.; Loi, S.; Okines, A.; Paplomata, E.; Hamilton, E.; Hurvitz, S.A.; Lin, N.U.; Borges, V.; Abramson, V.; Anders, C.; et al. Tucatinib, Trastuzumab, and Capecitabine for HER2-Positive Metastatic Breast Cancer. N. Engl. J. Med. 2020, 382, 597–609. [Google Scholar] [CrossRef] [PubMed]
  73. Shaw, A.T.; Bauer, T.M.; De Marinis, F.; Felip, E.; Goto, Y.; Liu, G.; Mazieres, J.; Kim, D.-W.; Mok, T.; Polli, A.; et al. First-Line Lorlatinib or Crizotinib in Advanced ALK-Positive Lung Cancer. N. Engl. J. Med. 2020, 383, 2018–2029. [Google Scholar] [CrossRef]
  74. De Bono, J.S.; Mateo, J.; Fizazi, K.; Saad, F.; Shore, N.; Sandhu, S.; Chi, K.N.; Sartor, O.; Agarwal, N.; Olmos, D.; et al. Olaparib for Metastatic Castration-Resistant Prostate Cancer. N. Engl. J. Med. 2020, 382, 2091–2102. [Google Scholar] [CrossRef]
  75. Golan, T.; Hammel, P.; Reni, M.; Van Cutsem, E.; Macarulla, T.; Hall, M.J.; Park, J.-O.; Hochhauser, D.; Arnold, D.; Oh, D.-Y.; et al. Maintenance Olaparib for Germline BRCA-Mutated Metastatic Pancreatic Cancer. N. Engl. J. Med. 2019, 381, 317–327. [Google Scholar] [CrossRef] [PubMed]
  76. Kopetz, S.; Grothey, A.; Yaeger, R.; Van Cutsem, E.; Desai, J.; Yoshino, T.; Wasan, H.; Ciardiello, F.; Loupakis, F.; Hong, Y.S.; et al. Encorafenib, Binimetinib, and Cetuximab in BRAF V600E–Mutated Colorectal Cancer. N. Engl. J. Med. 2019, 381, 1632–1643. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Pacheco, J.M.; Gao, D.; Smith, D.; Purcell, T.; Hancock, M.; Bunn, P.; Robin, T.; Liu, A.; Karam, S.; Gaspar, L.; et al. Natural History and Factors Associated with Overall Survival in Stage IV ALK-Rearranged Non–Small Cell Lung Cancer. J. Thorac. Oncol. 2019, 14, 691–700. [Google Scholar] [CrossRef] [PubMed]
  78. Kochanowski, K.; Morinishi, L.; Altschuler, S.J.; Wu, L.F. Drug persistence—From antibiotics to cancer therapies. Curr. Opin. Syst. Biol. 2018, 10, 1–8. [Google Scholar] [CrossRef]
  79. Levin-Reisman, I.; Brauner, A.; Ronin, I.; Balaban, N.Q. Epistasis between antibiotic tolerance, persistence, and resistance mutations. Proc. Natl. Acad. Sci. USA 2019, 116, 14734–14739. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Mizrahi, S.P.; Gefen, O.; Simon, I.; Balaban, N.Q. Persistence to anti-cancer treatments in the stationary to proliferating transition. Cell Cycle 2016, 15, 3442–3453. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Balaban, N.Q.; Helaine, S.; Lewis, K.; Ackermann, M.; Aldridge, B.; Andersson, D.I.; Brynildsen, M.P.; Bumann, D.; Camilli, A.; Collins, J.J.; et al. Definitions and guidelines for research on antibiotic persistence. Nat. Rev. Genet. 2019, 17, 441–448. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Cancers 13 02776 g001 550
Figure 1. Mutation analysis of 957 patients treated at City of Hope. (A) Co-mutation plot of all somatic mutations with a frequency of 6% or more in 957 metastatic cancer patients. Row percentages represent the frequency of gene mutations in the cohort. Each of the column bars represents the number of mutations for each patient. SNV: single nucleotide variants; VUS: variant of unknown significance; GI: gastrointestinal; GU: genitourinary. (B) Comparison of most frequent mutations from the City of Hope 2020 and University of Michigan 2017 cohorts of solid metastatic cancer patients.
Figure 1. Mutation analysis of 957 patients treated at City of Hope. (A) Co-mutation plot of all somatic mutations with a frequency of 6% or more in 957 metastatic cancer patients. Row percentages represent the frequency of gene mutations in the cohort. Each of the column bars represents the number of mutations for each patient. SNV: single nucleotide variants; VUS: variant of unknown significance; GI: gastrointestinal; GU: genitourinary. (B) Comparison of most frequent mutations from the City of Hope 2020 and University of Michigan 2017 cohorts of solid metastatic cancer patients.
Cancers 13 02776 g001
Cancers 13 02776 g002 550
Figure 2. Mutually exclusive and co-occurring genes. Genes with more than 35 positive mutation cases were considered. Chord diagrams of (A) mutually exclusive, and (B) co-occurring gene pairs with a false discovery rate (FDR) less than 0.01. The arcs shown in the scaled chord-diagram between the pairs of genes represent the mutual exclusivity or co-occurrence. There were seven mutually exclusive gene pairs, of which KRAS was mutually exclusive with EGFR, ERBB2, and VHL. APC was mutually exclusive with CCND1 and CDKN2A. Thirty-one co-occurring gene pairs with more than one co-occurring gene pair were identified. ARFRP1, ASXL1, AURKA, BCL2L1, GNAS, SRC, TOP1, ZNF217 were most commonly co-occurring, and each co-occurred with seven other genes (FDR < 0.01).
Figure 2. Mutually exclusive and co-occurring genes. Genes with more than 35 positive mutation cases were considered. Chord diagrams of (A) mutually exclusive, and (B) co-occurring gene pairs with a false discovery rate (FDR) less than 0.01. The arcs shown in the scaled chord-diagram between the pairs of genes represent the mutual exclusivity or co-occurrence. There were seven mutually exclusive gene pairs, of which KRAS was mutually exclusive with EGFR, ERBB2, and VHL. APC was mutually exclusive with CCND1 and CDKN2A. Thirty-one co-occurring gene pairs with more than one co-occurring gene pair were identified. ARFRP1, ASXL1, AURKA, BCL2L1, GNAS, SRC, TOP1, ZNF217 were most commonly co-occurring, and each co-occurred with seven other genes (FDR < 0.01).
Cancers 13 02776 g002
Cancers 13 02776 g003 550
Figure 3. Patient survival analysis in the City of Hope Metastatic Dataset. Five-year Kaplan–Meier overall survival (OS) plots of patients with known (A) APC, (B) CDKN2A, and (C) both APC and CDKN2A mutations. Shorter OS was significantly associated with the presence of CDKN2A mutations, while APC mutations were significantly associated with longer OS. The OS benefit is observed in patients with APC and CDKN2A co-mutation compared to CDKN2A mutations only.
Figure 3. Patient survival analysis in the City of Hope Metastatic Dataset. Five-year Kaplan–Meier overall survival (OS) plots of patients with known (A) APC, (B) CDKN2A, and (C) both APC and CDKN2A mutations. Shorter OS was significantly associated with the presence of CDKN2A mutations, while APC mutations were significantly associated with longer OS. The OS benefit is observed in patients with APC and CDKN2A co-mutation compared to CDKN2A mutations only.
Cancers 13 02776 g003
Cancers 13 02776 g004 550
Figure 4. Frequency of the first-line drug usage for the cohort (N = 957 patients) by cancer type (bottom), the proportion of patients in each cancer type treated with each therapy (indicated by color on heat map), and the number of patients who received each treatment (right). The 30 most frequently administered first-line anticancer drugs are displayed. Top row box plots represent the median, upper, and lower quartiles of the number of drugs used in each cancer type, and whiskers represent the limits of the distribution.
Figure 4. Frequency of the first-line drug usage for the cohort (N = 957 patients) by cancer type (bottom), the proportion of patients in each cancer type treated with each therapy (indicated by color on heat map), and the number of patients who received each treatment (right). The 30 most frequently administered first-line anticancer drugs are displayed. Top row box plots represent the median, upper, and lower quartiles of the number of drugs used in each cancer type, and whiskers represent the limits of the distribution.
Cancers 13 02776 g004
Cancers 13 02776 g005 550
Figure 5. Kaplan–Meier 5-year (A) overall survival (OS) and (B) progression-free survival (PFS) analysis based on targeted versus traditional therapy use in solid pan-caner metastatic patients. Patients who received targeted therapies as the first-line treatment had improved but not statistically significant improvement in OS compared to patients who did not (median survival 45 months compared to 37 months for patients who received traditional therapies, p = 0.072). However, PFS was significantly improved in patients receiving targeted therapies (median PFS 9.1 months [CI: 8.1–10.1 months] compared to 10.9 months [CI: 9.7–11.9 months] in patients receiving traditional therapies, p = 0.0029). CI: confidence interval.
Figure 5. Kaplan–Meier 5-year (A) overall survival (OS) and (B) progression-free survival (PFS) analysis based on targeted versus traditional therapy use in solid pan-caner metastatic patients. Patients who received targeted therapies as the first-line treatment had improved but not statistically significant improvement in OS compared to patients who did not (median survival 45 months compared to 37 months for patients who received traditional therapies, p = 0.072). However, PFS was significantly improved in patients receiving targeted therapies (median PFS 9.1 months [CI: 8.1–10.1 months] compared to 10.9 months [CI: 9.7–11.9 months] in patients receiving traditional therapies, p = 0.0029). CI: confidence interval.
Cancers 13 02776 g005
Table
Table 1. Baseline characteristics of 957 patients with seven metastatic solid cancers. IQR: interquartile range; MSI: microsatellite instability.
Table 1. Baseline characteristics of 957 patients with seven metastatic solid cancers. IQR: interquartile range; MSI: microsatellite instability.
Characteristic Pan-Cancer
(n = 957)		Breast (n = 87)Colorectal (n = 322)Gastrointestinal
(n = 123)		Genitourinary (n = 118)Gynecological (n = 59)Head and Neck (n = 19)Thoracic (n = 229)
Age at diagnosis, median (IQR)60 (50–69)50 (42–56)56 (48–66)63 (49–70)67 (60–73)61 (51.5–71)61 (50.5–65.5)64 (54–72)
Race, n (%)
White667 (69.7)50 (57.5)230 (71.4)74 (60.2) 98 (83.1)47 (79.7)17 (89.5)151 (65.9)
Asian194 (20.3)18 (20.7)60 (18.6)35 (28.5)13 (11)8 (13.6)1 (5.3)59 (25.9)
African American47 (4.9)12 (13.8)15 (4.7)4 (3.3)2 (1.7)2 (3.4)1 (5.3)11 (4.8)
Other18 (1.9)4 (4.6)7 (2.1)7 (5.7)3 (2.5)2 (3.7)0 (0)4 (1.7)
Unknown22 (2.3)3 (3.5)10 (3.1)3 (2.4)2 (1.6)0 (0)0(0)4 (1.7)
Ethnicity, n (%)
Hispanic or Latino193 (20.2)24 (27.6)74 (23)37 (30.1)10 (8.5)16 (27.1)3 (15.8)29 (12.7)
Not Hispanic or Latino752 (78.6)61(70.1)244 (75.8)83 (67.5)108 (91.5)43 (72.9)16 (84.2)197 (86)
Unknown or not disclosed12 (1.2)2(2.3)4 (1.2)3 (2.4)0 (0)0 (0)0 (0)3 (1.3)
Sex, n (%)
Female497 (51.9)87 (100)138 (42.9)54 (43.9)27 (22.9)59 (100)8 (42.1)124 (54.1)
Male460 (48.1)0 (0)184 (57.1)69 (56.1)91 (77.1)0 (0)11 (57.9)105 (45.9)
Tumor Burden, n (%)
Low204 (21.2)19 (21.8)70 (21.7)20 (16.316 (13.1)21 (35.6)10 (52.6)48 (21.1)
Intermediate115 (12)6 (6.9)37 (11.5)17 (13.8)8 (6.6)10 (16.9)3 (15.8)34 (14.9)
High17 (1.8)0 (0)2 (0.6)1 (0.8)1 (0.8)1 (1.7)2 (10.5)10 (4.4)
Unknown624 (65)62 (71.2)221 (66.1)85 (69.1)97 (79.5)27 (45.8)4 (21.1)136 (59.6)
Stage at initial diagnosis, n (%)
I41 (4.3)16 (18.4)9 (2.8)0 (0)0 (0)5 (8.5)1 (5.3)10 (4.4)
II81 (8.5)32 (36.8)29 (9)2 (1.6)0 (0)4 (6.8)3 (15.8)11 (4.8)
III134 (14)23 (26.4)62 (19.3)3 (2.4)4 (3.9)23 (39)0 (0)19 (8.3)
IV701 (73.3)16 (18.4)222 (68.9)118 (95.9)114 (96.6)27 (45.8)15 (78.9)189 (82.5)
MSI status, n (%)
High8 (0.8)0 (0)2 (0.6)2 (1.6)0 (0)1 (1.7)1 (5.3)2 (0.9)
Stable351 (36.7)26 (29.9)116 (36)39 (31.7)25 (21.2)32 (54.8)15 (78.9)98 (42.8)
Unknown598 (62.5)60 (70.1)204 (63.4)82 (66.7)93 (78.8)26 (44.1)3 (15.8)129 (56.3)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Roosan, M.R.; Mambetsariev, I.; Pharaon, R.; Fricke, J.; Baroz, A.R.; Chao, J.; Chen, C.; Nasser, M.W.; Chirravuri-Venkata, R.; Jain, M.; et al. Evaluation of Somatic Mutations in Solid Metastatic Pan-Cancer Patients. Cancers 2021, 13, 2776. https://doi.org/10.3390/cancers13112776

AMA Style

Roosan MR, Mambetsariev I, Pharaon R, Fricke J, Baroz AR, Chao J, Chen C, Nasser MW, Chirravuri-Venkata R, Jain M, et al. Evaluation of Somatic Mutations in Solid Metastatic Pan-Cancer Patients. Cancers. 2021; 13(11):2776. https://doi.org/10.3390/cancers13112776

Chicago/Turabian Style

Roosan, Moom R., Isa Mambetsariev, Rebecca Pharaon, Jeremy Fricke, Angel R. Baroz, Joseph Chao, Chen Chen, Mohd W. Nasser, Ramakanth Chirravuri-Venkata, Maneesh Jain, and et al. 2021. "Evaluation of Somatic Mutations in Solid Metastatic Pan-Cancer Patients" Cancers 13, no. 11: 2776. https://doi.org/10.3390/cancers13112776

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop