Next Article in Journal
Cisplatin Resistance in Testicular Germ Cell Tumors: Current Challenges from Various Perspectives
Next Article in Special Issue
Primary Driver Mutations in GTF2I Specific to the Development of Thymomas
Previous Article in Journal
Hereditary Gastric and Breast Cancer Syndromes Related to CDH1 Germline Mutation: A Multidisciplinary Clinical Review
Previous Article in Special Issue
Impact and Diagnostic Gaps of Comprehensive Genomic Profiling in Real-World Clinical Practice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Whole Exome Sequencing of Multi-Regional Biopsies from Metastatic Lesions to Evaluate Actionable Truncal Mutations Using a Single-Pass Percutaneous Technique

1
Department of Haematology-Oncology, National University Cancer Institute, Singapore 119074, Singapore
2
Cancer Science Institute of Singapore, National University of Singapore, Singapore 117599, Singapore
3
Department of Radiology, National University Hospital, Singapore 119074, Singapore
4
Department of Pathology, National University of Singapore, Singapore 119077, Singapore
*
Author to whom correspondence should be addressed.
These authors contributed equally to this paper.
Cancers 2020, 12(6), 1599; https://doi.org/10.3390/cancers12061599
Submission received: 17 April 2020 / Revised: 27 May 2020 / Accepted: 7 June 2020 / Published: 17 June 2020
(This article belongs to the Special Issue Precision Medicine in Solid Tumors)

Abstract

:
We investigate the feasibility of obtaining multiple spatially-separated biopsies from a single lesion to explore intratumor heterogeneity and identify actionable truncal mutations using whole exome sequencing (WES). A single-pass radiologically-guided percutaneous technique was used to obtain four spatially-separated biopsies from a single metastatic lesion. WES was performed to identify putative truncal variants (PTVs), defined as a non-synonymous somatic (NSS) variant present in all four spatially separated biopsies. Actionable truncal mutations—filtered using the FoundationOne panel—were defined as clinically relevant PTVs. Mutational landscapes of each biopsy and their association with patient outcomes were assessed. WES on 50 biopsied samples from 13 patients across six cancer types were analyzed. Actionable truncal mutations were identified in 9/13 patients; 31.1 ± 5.12 more unique NSS variants were detected with every additional multi- region tumor biopsy (MRTB) analyzed. The number of PTVs dropped by 16.1 ± 17.9 with every additional MRTB, with the decrease most pronounced (36.8 ± 19.7) when two MRTB were analyzed compared to one. MRTB most reliably predicted PTV compared to in silico analysis of allele frequencies and cancer cell fraction based on one biopsy sample. Three patients treated with actionable truncal mutation-directed therapy derived clinical benefit. Multi-regional sampling for genomics analysis is feasible and informative to help prioritize precision-therapy strategies.

1. Introduction

Intratumor heterogeneity is a key challenge in precision cancer therapy, contributing to treatment resistance, therapeutic failure and poor prognosis [1,2]. With the growing use and reducing cost of next generation sequencing, the full extent of the complexity and genomic diversity within tumors are becoming more apparent [1,3,4]. Genotype-directed targeted therapies are becoming the standard of care, and as tumor molecular profiling becomes more widely used in routine practice, physicians will require the necessary tools to translate genomic information into clinically actionable results. The consequences of intratumor heterogeneity, such as resistance to drug therapy [3,4,5] leading to disease recurrence and death, is at least partially the result of limitations in the ability to define the clonal frequency of driver events for prioritization of drug targeting in tumors. Furthermore, it has been demonstrated that high levels of ITH results in poorer survival outcomes across a wide range of cancer types [6,7]. To mitigate this challenge, a more comprehensive view of the mutational diversity of each tumor lesion is required.
The mutational diversity attributed to ITH limits our ability to resolve the full spectrum of cancer pathway aberrations through a single biopsy of the tumor lesion and may under/overestimate driver alterations [4,8,9]. Therefore, multi-region tumor biopsies (MRTBs) are highly beneficial to attenuate the challenge of estimating the prevalence of oncogenic clonal driver mutations. Targeting clonal driver (truncal) mutations would potentially be more effective than targeting subclonal (branch) mutations in a tumor [10,11]. Yap et al. proposed the targeting of genetic alterations located on the trunk of an individual’s phylogenetic tree as a more effective clinical strategy [10] as truncal mutations are more likely to represent the core driver mutations within the tumor [10]. In view of the importance of identifying truncal mutations, in silico approaches such as the ABSOLUTE algorithm [12] have been developed to predict truncal variants from a single biopsy sample. However, their ability to identify actionable truncal mutations that would be clinically relevant is hitherto unknown. Similarly, various gene panels have been utilized for diagnostic purposes but the minimum number of MRTB samples needed to address issues associated with ITH remains unknown. This study—conducted across six major cancer types—aims to outline: (a) the safety and significance of MRTB to help navigate the complexities of ITH, (b) the minimum number of MRTB samples required when different gene panels were used for clinical assessment, and (c) the feasibility and clinical efficacy of the approach for identifying clinically actionable truncal mutations (i.e., mutations present in all MRTB obtained from a single tumor lesion) and their outcomes when targeted for strategic therapeutic intervention.

2. Results

A cohort of 15 patients with metastatic colorectal carcinoma (CRC; n = 1), non-small cell lung cancer (NSCLC; n = 6), ovarian carcinoma (OV; n = 3), breast carcinoma (BC, n = 1), uterine carcinoma (UC, n = 2), hepatocellular carcinoma (HCC; n = 1), or cervical cancer (CC, n = 1) were recruited to the study. A single-pass radiologically-guided percutaneous biopsy technique was used to obtain MRTBs from a dominantly-progressing metastatic lesion in each patient with core biopsies taken at least 2 mm apart within the same metastatic lesion. Two patients (one with NSCLC and another with cervical cancer) were excluded from analysis as all the MRTB samples collected from them failed quality control (QC). One patient (UC) (P11; Figure 1) had two (out of four) biopsy samples that failed QC which were subsequently excluded from the analysis. Similarly, one patient’s (P01) germline sample (i.e., buccal swab) failed QC and was replaced with whole blood sample. All 15 patients tolerated the procedure well with no significant adverse events, except for one patient (P09) who developed a moderately-sized right sided pneumothorax requiring observation overnight and serial imaging to ensure spontaneous resolution of the pneumothorax.
The investigation pipeline adopted for analyzing whole-exome sequencing (WES) data that passed QC is illustrated in Figure 1. All processed samples had a DNA concentration greater than 4 ng/μL. Average sequencing depth, Q30 percentage and uniformity of coverage obtained were 128X ± 29.2, 87.5% ± 4.35, and 91.2% ± 1.94 respectively.

2.1. Tumour Variant Load

The mutational landscape of each patient was examined to evaluate the extent of ITH across different cancer types. The non-synonymous somatic mutational load (nssML)—defined as the total number of non-synonymous somatic (NSS) variants present—was scrutinized for each biopsy sample (Figure 2A). Results indicate that patient P04 has the highest average non-synonymous somatic mutational load (608.0 ± 41.7), while patient P05 with NSCLC has the highest diversity (i.e., difference in non-synonymous somatic mutational load) among the four MRTB samples analyzed (104.3 ± 49.3). The median diversity across all patients was 7.46 (range: 0.957 to 49.3). Friedman test of difference among the different MRTB samples indicated no statistically significant difference between the number of non-synonymous somatic variants present in each MRTB sample (p = 0.691).

2.2. Intratumor Heterogeneity

To investigate the extent of the intratumoral heterogeneity, the amount of truncal (i.e., ubiquitous non-synonymous somatic variants that occur in all MRTB samples analyzed) and branch (i.e., non-synonymous somatic variants that do not occur in all MRTB samples) variants were analyzed (Figure 2B). Phylogenetic trees were also constructed to illustrate this phenomenon graphically (Figure 2C, Figure S1). As demonstrated, two patients (P01 and P05) did not have any truncal variants while patients P03 and P13 only had two and one truncal variant(s), respectively. On average, 24.1% ± 20.7, 14.7% ± 13.7 and 61.2% ± 20.6 of non-synonymous somatic variants were truncal, branch and private mutations, respectively. A high level of intratumoral heterogeneity (75.5% ± 34.6) across different tumors was observed, with private mutations dominating the mutational landscape (p < 0.05). When copy number alterations (CNAs) were interrogated, a moderate degree of diversity (branch amplification: 54.3% ± 34.7, p = 0.083; branch deletion: 59.4% ± 34.8, p = 0.050) was observed (Figure 2D, Figure S2).
Statistically significant somatic cancer driver mutations (ssCDMs) were juxtaposed with non-synonymous somatic variants identified in our study cohort. Results indicate that detectable somatic cancer driver mutations were more likely to be truncal variants (68.2%; Figure 2E, Figures S3 and S4); however, the difference was not statistically significant (p = 0.177). Truncal somatic cancer driver mutations across this study cohort include AKT1, ATM, BCOR, CHD4, KRAS, MAP3K1, and PIK3CA; conversely, branch somatic cancer driver mutations include ERBB2, FOXA1, and PPM1D. In our cohort, EGFR mutations were confined to lung cancers, with three out of four (75%) NSCLC patients found to harbor a truncal variant in at least one reportable mutation in EGFR [13].

2.3. Statistical Saturation Analysis

The relationship between truncal variants and the number of MRTB samples analyzed was examined. A unique variant in this case refers to a distinct non-synonymous somatic variant that appears in at least one of the MRTB samples analyzed simultaneously. In general, with an increasing number of MRTB samples, a monotonically increasing trend in the number of unique variants and correspondingly decreasing number of truncal variants can be observed (Figure S5).
At the exome level (i.e., WES NSS gene panel), on average 31.1 ± 5.12 more unique non-synonymous somatic variants were detected with every additional MRTB sample analyzed. Conversely, using the number of MRTB samples analyzed simultaneously as the baseline reference to determine putative truncal variants (PTVs), a monotonically decreasing number of PTVs can be observed with an increasing number of MRTB samples. The number of PTVs dropped by 16.1 ± 17.9 on average with every additional MRTB sample, with the decrease most pronounced (36.8 ± 19.7) when two MRTB samples were analyzed compared to just one. Similar trends were observed from filtered variants when the WES data was mapped to genes matching four cancer gene panels—namely COSMIC Cancer Gene Census (CGC), Ion AmpliSeq™ Cancer Hotspot Panel v2 (Life Technologies, Carlsbad, CA, USA), TruSight® Cancer panel (Illumina Inc, San Diego, CA, USA) and FoundationOne™ cancer gene panel (Foundation Medicine, Cambridge, MA, USA) (Figure S5, Tables S1 and S2). However, the change in the number of unique/truncal mutations with increasing MRTB samples was less pronounced (<2 variants on average).
To quantitatively corroborate the minimum number of MRTB samples required across different gene panels, statistical saturation analysis was conducted. Results (Figure S6) indicate that every additional MRTB sample analyzed increases the ability to detect unique variants when WES non-synonymous somatic, TruSight® and FoundationOne™ gene panels were used; as for CGC and AmpliSeq™ gene panels, at least two or three MRTB samples (depending on the panel used) were required, respectively, before changes in the number of unique variants became statistically not significant. To identify PTVs, results (Figure 3A) suggest that at least two MRTB samples were required for the CGC, AmpliSeq™ and TruSight® gene panels, and three samples were required for the FoundationOne gene panel; the WES NSS gene panel, on other hand, required four or more MRTB samples based on our analysis. In addition, the positive predictive value (PPV) was determined to evaluate the extent to which truncal variants (defined using four MRTB samples as the baseline reference) can be identified among all variants found in a set of less than four MRTB samples. Results (Figure 3B) indicate that four (or more) and two MRTB samples are needed for WES NSS and FoundationOne™ gene panels, respectively, while CGC, AmpliSeq™ and TruSight® cancer gene panels only require a single biopsy sample. The greatest significant increase in PPV based on WES NSS and FoundationOne panels were from one analyzed sample to two samples (Figure 3B).

2.4. Prediction of Truncal Mutations

To evaluate the ability to identify truncal variants (defined based on four MRTB samples) using a single biopsy sample, two metrics were used as classification thresholds, namely allele frequency (AF) and cancer cell fraction (CCF). Computation of CCF values—the proportion of cancer cells within which the variant is present—for patient P07 (OV) was unable to be performed due to inadequate information related to somatic copy number alteration. Hence, patient P07 was excluded for the purpose of this analysis.
First, the threshold value (for each respective metric) that produces the best average prediction accuracy across all patients was examined. Results, as illustrated in Figure S7, suggest that AF generally outperformed CCF across different patients. Average prediction accuracy improved between 2.7% and 15.3% (across different gene panels) when AF was used as the classification threshold. However, statistical significance of difference was achieved for the WES NSS gene panel only (p = 0.021).
Next, the threshold value (for each respective metric) that produces the best average prediction accuracy across patients with the same cancer type was scrutinized. Results, as demonstrated in Figure 3C, show that AF outperformed CCF by 15.6% to 30.4% across the different gene panels, with the FoundationOne™ cancer gene panel having the largest difference. Statistical significance of difference was achieved for the WES NSS gene panel only (p = 0.031).

2.5. Clinical Therapeutic Intervention

Three (23.1%) patients (P06, P10, and P11) received an inhibitor targeting an actionable truncal mutation based on molecular profiling while another six (46.2%) patients were treated with non-actionable truncal mutation-directed therapy either because they did not have any actionable truncal mutations or there was no available therapy to target the actionable truncal mutation at our center (Table 1). An illustrative example of the patients’ mutational profile can be found in Figures S8 and S9.
Using each patient as his/her own control as a strategy to attenuate confounding factors resulting from the diverse patient population and tumor types, we assessed the clinical efficacy of actionable truncal mutation-directed therapy by comparing progression free survival (PFS) on actionable truncal mutation-directed therapy (PFS-actionable truncal mutation-directed therapy) or non-actionable truncal mutation-directed therapy with the PFS for the most recent prior therapy (PFS-A) in each of these patients [14]. Two NSCLC patients harboring an Epidermal growth factor receptor (EGFR)_T790M mutation were treated with a single agent EGFR_T790M specific tyrosine kinase inhibitor [15] with differing clinical outcomes. Patient P06 had a truncal EGFR_T790M mutation while patient P05 had an EGFR_T790M mutation as a branch mutation. Patient P06 had a partial response and was still on active treatment at last review with a PFS of >25 months (Figure S10a), while patient P05 developed worsening neuro-cognitive defects resulting in cessation of treatment after two months. The PFS ratio for patients P05 and P06 was 0.06 and 10.2, respectively.
Patient P10 with breast cancer and a phosphoinositide-3-kinase catalytic alpha polypeptide (PIK3CA)_H1047R truncal mutation was enrolled into a highly selective PI3Kα/β inhibitor phase 1 (dose escalation) trial [16], but progressed shortly after with a PFS of 1.9 months and PFS ratio of 0.95 despite deriving symptomatic benefit while on the trial. Lastly, patient P11 with uterine carcinoma harbored a truncal RAC-alpha serine/threonine-protein kinase (AKT1) E17K mutation in two out of two of her MRTB cores analyzed (only two of four cores had DNA of sufficient quality for analysis in her case). She had significant sacral bone pain from bone metastasis and received a pan-AKT inhibitor. Strikingly, her pain significantly improved and subsequent scans revealed a 21% reduction in the sum of target lesions with a PFS of 6.1 months (Figure S10b). When compared to the PFS from her most recent physician’s choice therapy, a PFS ratio of 1.5 was observed.
It is noteworthy that the CCF metric performed relatively well in predicting the truncal status of the variants targeted. The median PFS-actionable truncal mutation-directed therapy for the small number of patients treated with actionable truncal mutation-directed therapy was 6.1 months, with a median PFS ratio of 1.5. These findings do suggest that the truncal status of tumors influences response and, if validated, could potentially be used in personalized cancer treatment to help prioritize therapeutic strategies.

3. Discussion

ITH represents a significant challenge to precision medicine and contributes to drug resistance. Several studies employing multi-region tumor sampling from post-surgical samples have greatly increased our understanding of tumor evolution and highlighted the importance of tumor sampling from spatially distinct areas in order to avoid erroneous interpretation of genomic data from single sampling bias [8,9,10]. In clinical practice, however, a systematic regional analysis of resected tumor specimens is unfeasible in the majority of patients with metastatic or recurrent cancer who may only have limited accessible intracorporeal tissue for sampling/biopsy. Hence, high quality patient samples across six major cancer types were analyzed to address certain exigent issues related to ITH and devise a potential novel solution to tackling the complexities of tumor heterogeneity when confronted with the reality of treatment decision-making based on limited access to tumor tissue.
Results in our small cohort demonstrate a high degree of ITH (>65% branch mutations) across the majority of patients, with private mutations dominating the mutational landscape. Clearly, this indicates that ITH is a ubiquitous issue that would confound the ability to identify bona fide truncal variants. Statistical saturation analysis demonstrates that for small targeted cancer gene panels like CGC, AmpliSeq™ and TruSight®, a minimum of two MRTB samples are required to identify PTVs; for a larger cancer gene panel like FoundationOne™, at least three MRTB samples are needed. The determination of the minimum number of MRTB samples required is highly valuable as it enables clinicians to find the equilibrium between cost and accuracy (of identifying bona fide truncal variants), and allows the choice of which cancer gene panel to use with the amount of tumor tissue available.
Examination of the nssML of each patient indicates that individual intratumor biopsy samples comprise a similar amount of NSS variants (Figure 2A) while the aggregated nssML shows that every additional MRTB sample would offer a statistically significant increment in the total number of NSS variants (Figure S6). Correlation analysis indicates a strong correlation between the average number of NSS variants among individual intratumor biopsy samples and aggregated NSS variants across all four MRTB samples (Pearson’s rho = 0.975, p < 0.001). This suggests minimal intratumoral variation in the ML based on our series and that mutational burden is less likely to be impaired by sampling bias.
Given the coveted utopia of making informed clinical decisions based on a single biopsy sample, in silico methods for predicting truncal variants are of particular interest. AF and CCF are two favored metrics commonly used. Empirical experiments indicate that when AF was used as the threshold to classify truncal variants, it achieved comparable, if not better, accuracy compared to CCF; although both approaches were less compelling for some cancer types. Of note, different cancer types (at the whole-exome level) favor different threshold values for segregating truncal from branch variants, suggesting that each tumor type exhibits distinct biological characteristics that require dedicated data analytics.
Improved clinical outcomes were observed in two out of three patients whose truncal mutations were selectively targeted. Remarkably, all patients in our series treated with actionable truncal mutation-directed therapy derived symptomatic benefit with improvement in their performance status. The small patient numbers across a diverse spectrum of tumors limits our ability to draw significant conclusions within each tumor type, but nonetheless demonstrates its applicability in a variety of tumor types and preliminary evidence of clinical benefit when used for therapeutic prioritization in selected patients. It is noteworthy that patient P05—who had an EGFR_T790M branch mutation that was targeted—did not respond well to the treatment (Table 1). This reaffirms the hypothesis that increased therapeutic efficacy can be achieved by targeting truncal mutations within a tumor, and that targeting branch mutations may result in only partial treatment efficacy and/or accelerated growth in non-targeted subpopulations [17]. Undeniably, the cost per patient of this approach is high; it has been estimated at USD$5000 per patient for the acquisition of biopsy samples and profiling of four biopsy core samples as well as a germline control, but this is likely to be mitigated in the future as next generation sequencing technologies become more widely used and cost of sequencing gradually decreases. Crucially, the data provided by multi-region sequencing of a tumor could have important implications for the prioritization of druggable targets in the clinical setting. To the best of our knowledge, this is the first study assessing the feasibility and utility of obtaining tissue biopsies from multiple spatially separated regions from a single metastatic site percutaneously. A limitation of this study is the small sample size. Nevertheless, it provides adequate resolution into the complexity and management of ITH. In addition, the ITH analysis performed in this study is based on the construction of phylogenetic trees with the implicit assumption that a tumor sample can be meaningfully summarized as the collection of mutations observed in that sample, or that only a single or dominant clone exists per sample that carries all mutations, which could lead to biased inferences.
In our study, we only analyzed single nucleotide variants (SNV) which may potentially underestimate the frequency of clinically actionable mutations and the mutational load of the tumor. We focused solely on SNV mainly because they make up the majority of pathogenic variants relevant in solid tumor malignancies (59.39%) compared to other genomic alterations such as indels, structural variants and copy number loss [18]. Indeed, the majority of annotated variants in oncogenic and actionable target databases such as OncoKB and cancer hotspots consist of predominantly SNV. In addition, the majority of approved targeted inhibitors available for solid cancers currently are also mainly directed at aberrations associated with SNV. As our results relied entirely on WES analysis, and so one of the limitations of our study is the dependence on the size of the panel testing. As the number of variants being considered increases, so does the required number of samples. Our study also used fresh frozen tissue for WES analysis, which resulted in 13 of the 81 samples collected failing quality assurance due to degradation of DNA. As we continue to expand our taxonomy of tumors and seek to enhance the applicability of this approach in clinical practice, it would be ideal to optimize this approach for the clinical grade analysis of formalin-fixed, paraffin-embedded tumor samples, to enable histological and immunohistochemical analyses of samples to be performed in parallel with genomic analysis in the future. It has been suggested that liquid biopsies based on genomic analyses of circulating cell-free tumor DNA (ctDNA) and circulating tumor cells (CTC) may obviate the need for tumor biopsies [19]. Liquid biopsy platforms offer the potential for real-time sampling and resampling of tumor material for monitoring of therapeutic efficacy [20] and early detection of resistance subclones [21]. Furthermore, the MRTB approach we have used in this study will not be feasible in patients with inaccessible lesions. However, the inability to characterize liquid biopsies histologically limits the extent of biomarker analyses, particularly where tumor microenvironmental features (e.g., programme cell death-1/programme cell death ligand-1(PD1/PDL1) protein expression), immune cell infiltrates and stromal content are concerned. Further studies comparing the clinical utility of multi-spatial or multi-lesional biopsy approaches with that of liquid biopsies in monitoring the emergence of resistance and therapeutic efficacy are eagerly awaited

4. Materials and Methods

4.1. Patients and Specimens Collection

All patients were recruited and treated at the National University Cancer Institute (NCIS), Singapore, between December 2014 and May 2016. WES and data analytics were performed at the National University of Singapore (NUS), Singapore. All procedures were conducted in accordance with the approved protocols and written informed consent was provided by the patients (File S1). Eligible patients were at least 21 years of age, had a histological or cytological diagnosis of advanced or metastatic solid malignancy and had recurrent disease for which tissue biopsy was indicated as part of routine clinical practice. This study was approved by the National Health Group Domain Specific Review Board IRB number 2014/00665
Each patient had 5 tumor biopsy samples obtained from one metastatic lesion using a single-pass radiologically-guided percutaneous biopsy technique. This technique involves the insertion of a coaxial needle together with its trocar into a lesion. The trocar is subsequently removed to allow for the introduction of a biopsy device that is composed of a needle with a 1.5 cm throw to facilitate multiple passes along acute angles from a single lesion via a single percutaneous access. Each biopsy sample was obtained at least 2 mm apart. One biopsy core was sent to the histopathology lab as part of routine clinical management while the remaining four tumor biopsy samples were analyzed using WES. Four patients had biopsies obtained from the lung (P3, P4, P5 and P6), four had peritoneal nodes biopsied (P1, P7, P8 and P9), three patients underwent a liver biopsy (P10, P11 and P12) and one patient each had a bone (P13) and supraclavicular lymph node (P2) biopsied. Germline samples were collected from each patient in the form of a buccal swab. If the germline sample failed quantitative QC, whole blood would be used as replacement. Samples with DNA concentration <4 ng/μL would be deemed to have failed QC. Radiological images were obtained as part of clinical care.

4.2. Whole-Exome Sequencing

Extraction of genomic DNA from tumor samples was carried out using the Qiagen Allprep DNA/RNA Micro Kit (Qiagen, Hilden, Germany). The MasterAmp Buccal Swab Kit (Epicenter, Madison, WI, USA) was used to extract DNA from buccal swab samples while the Qiagen EZ1 DNA Blood 350 µL Kit (Qiagen, Hilden, Germany) was used to process whole blood samples. An Illumina NextSeq 500 Sequencing System (Illumina, San Diego, CA, USA) was utilized to perform 150 base pair paired-end WES. All experiments were conducted in accordance to manufacturer guidelines at the Cancer Science Institute of Singapore, NUS. Full patient data can be found in the National Centre for Biotechnology Information (NCBI) Sequence Read Archive (SRA) under accession number SRP137039

4.3. Sequencing Reads Alignment and Somatic Variant Detection

Sequence reads were aligned to the hg19 reference genome using the Burrows–Wheeler Aligner (BWA) v0.7.7 [22] and realignment to the hg19 reference genome was performed by using the Genome Analysis ToolKit (GATK) v3.3.0 [23]. Variant calling was performed with duplication removal and base recalibration prior to variant calling using MuTect somatic variant caller v1.1.7 [24] and annotated using Oncotator v1.8.0.0 [25]. Only NSS variants were filtered out for analysis independent of CCF or AF. All sequencing data have been made available in the NCBI Sequence Read Archive (SRA) under accession number SRP137039.

4.4. Copy Number Alterations

Somatic CNAs were detected using several algorithms. Succinctly, the computation of raw copy number calls and the adjustment of GC content of the raw copy number calls were performed using VarScan2 v2.3.9 [26]. Re-centering and segmentation of the adjusted copy number calls were conducted using DNAcopy v1.44.0 [27]. Sample purity was assessed using the ABSOLUTE algorithm (Appendix A, Table A1).

4.5. Cancer Gene Panels

Five gene panels—WES NSS, COSMIC Cancer Gene Census (CGC), the Ion AmpliSeq™ Cancer Hotspot Panel v2, the TruSight® Cancer panel and the FoundationOne™ cancer gene panel—were examined. The WES NSS gene panel comprises of all protein-coding genes (with NSS property) in the genome, while the CGC gene panel consists of all statistically significant cancer-specific genes curated from the COSMIC Cancer Gene Census (CGC) [28,29] and The Cancer Genome Atlas (TCGA) [30,31]. The complete list of CGC interrogated genes for CRC (n = 29), NSCLC (n = 30), OV (n = 10), BC (n = 52), EC (n = 58), and HCC (n = 26) is available in Figures S3 and S4. Commercial cancer gene panels like AmpliSeq™, TruSight® and FoundationOne™ are comprised of 50, 94, and 315 cancer-related genes, respectively; their interrogated gene lists are available at ThermoFisher [32], Illumina [33], and Foundation Medicine [34], respectively.
All NSS variants identified in the WES NSS gene panel were subsequently juxtaposed with individual targeted cancer gene panels (i.e., CGC, AmpliSeq™, TruSight® and FoundationOne™) and variants from mismatched genes were winnowed out. The resulting list of variants for each cancer gene panel was used for subsequent downstream analysis.

4.6. Construction of Phylogenetic Trees

The construction of phylogenetic trees was carried out based on a binary table that represents the presence or absence of variants across all MRTB samples. Using the PHYLogeny Inference Package v3.695 (PHYLIP) [35] and matched germline information as the outgroup root, discrete character parsimony was used to generate the topology of the phylogenetic trees. Based on the computed mutation counts, the length of the trunk, shared and private branches were drawn accordingly.

4.7. Statistical Saturation Analysis

The computation of the average number of unique variants present when ‘k’ number of MRTB samples were analyzed concurrently was performed based on the following Formula (1):
A v e r a g e   u n i q u e   v a r i a n t s k = 1 n i = 1 n x i ,
where ‘n’ denotes the number of permutation combinations available when ‘k’ number of MRTB samples were selected from 4 MRTB samples, and ‘xi’ refers to the number of variants that are present in at least one of the MRTB samples examined in combination set ‘i’. Similarly, the average number of PTVs was calculated based on the formula above but with ‘xi’ defined as the number of variants present in all MRTB samples scrutinized in combination set ‘i’.
PPV, on other hand, was computed based on the following Formula (2):
A v e r a g e   P P V k = 1 n i = 1 n C P C V i , k ,
where ‘n’ denotes the number of permutation combinations available when ‘k’ number of MRTB samples were selected from 4 MRTB samples, ‘C’ represents the number of variants that occur in all 4 MRTB samples, and ‘PTVi,k’ (i.e., putative truncal variant) refers to the number of variants present in all ‘k’ number of MRTB samples examined simultaneously in combination set ‘i’.
All statistical hypothesis tests were conducted using the Wilcoxon signed rank test unless otherwise stated. Statistical significance is considered when the p-value is less than 0.05.

4.8. Cancer Cell Fraction and Allele Frequency

The CCF value was estimated using the ABSOLUTE algorithm [13] for each somatic single nucleotide variant (SNV) site based on its AF, CNAs, ploidy and purity of the tumor tissue analyzed. Based on the computed CCF values, a range of thresholds—from 0.90 to 1.00 incremented at a step size of 0.01—were used to classify variants into either truncal or branch.
AF was calculated by dividing the number of alternative sequence read counts with the total number of (alternative and reference) sequence read counts. Likewise, the clonality of variants was determined by comparing the variant’s AF with a range of thresholds (from 0.01 to 0.55 incremented at a step size of 0.01). The respective thresholds used for both the cancer cell fraction and allele frequency analysis according to the individual panels are shown in Table 2 below:

4.9. Prediction of Truncal Mutations

To stratify variants into either truncal or branch based on a single biopsy sample, different (AF and CCF) threshold values were investigated. For each threshold value examined, the following formula was employed to compute the average classification accuracy, by which 4 MRTB samples were used as the baseline reference for defining bona fide truncal variants (3).
A v e r a g e   A c c u r a c y k = 1 4 i = 1 4 x i , k y i ,
where ‘k’ denotes the examined threshold value, ‘xi,k’ represents the number of correct classification made for biopsy sample ‘i’ using threshold value ‘k’, and ‘yi’ refers to the total number of variants assessed for biopsy sample ‘i’.

4.10. Assessing Targeted Therapy Outcomes

Treatment decisions were made using molecular profiling results from one core biopsy, as reported in the Intergrated Molecular Analysis of Cancer (IMAC) study [14] while the remaining four core biopsies were analyzed for this study. Using each patient as his/her own control as a strategy to attenuate confounding factors resulting from the diverse patient population and tumor types, we assessed the clinical efficacy of actionable truncal mutation-directed therapy by comparing the PFS for each patient who received actionable truncal mutation-directed therapy (PFS-actionable truncal mutation-directed therapy) with the PFS for the therapy immediately before actionable truncal mutation-directed therapy (PFS-A) [36]. If the PFS of PFS-actionable truncal mutation-directed therapy/PFS-A ratio was ≥1.3, then the molecular profiling-selected actionable truncal mutation-directed therapy was defined as having benefit for the patient compared to the physician’s choice chemotherapy. The PFS ratio was defined as PFS actionable   truncal   mutation directed   therapy PFS   for   the   therapy   immediately   before   actionable   truncal   mutation directed   therapy   ( PFS A ) and was used to evaluate the efficiency of the therapeutic intervention [36].

5. Conclusions

In conclusion, this study has demonstrated: (i) the importance of performing multiple biopsies despite extant in silico prediction methods, (ii) the minimum number of MRTB samples required to alleviate challenges related to ITH is dependent on the tested hypothesis and the examined gene panel, but that at least two biopsies should be submitted for analysis to achieve a PPV of >90% identifying AT mutations, and (iii) the feasibility and clinical efficacy of adopting the proposed approach for strategic therapeutic intervention. Further validation of this approach for identifying and targeting AT mutations in larger cohorts will be required to fully assess its potential value as a precision medicine strategy to circumvent the challenges of intratumoral heterogeneity in cancer therapy.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6694/12/6/1599/s1, Figure S1: Phylogenetic trees constructed from NSS mutations, Figure S2: Heatmap visualization illustrating the presence/absence of copy number alterations in relation to each region, Figure S3: ssCDMs for CRC, NSCLC, and OV, Figure S4: ssCDMs for BC, UC, and HCC, Figure S5: Mutational characteristics across different number of MRTB samples and gene panels, Figure S6: Boxplot illustrating the average number of unique variants across different gene panels and number of MRTB samples, Figure S7: Best average prediction accuracy of PTVs across different patients, Figure S8: Illustrative example of patients P01-P06 mutational profile, Figure S9: Illustrative example of patients P07-P13 mutational profile, Figure S10: Pre - treatment and 8- week post treatment images from, Table S1: Average number of unique variants detected across different number of MRTB samples and gene panels, Table S2: Average number of PTVs detected across different number of MRTB samples and gene panels, File S1: Study Protocol.

Author Contributions

Conceptualization, D.S.P.T., R.S. and A.G.; methodology, A.G., R.S., M.F.M.O. and D.S.P.T.; software, D.T., V.H., T.Z.T. and S.E.G.; validation, D.T., T.Z.T. and S.E.G.; formal analysis, D.T., V.H., T.Z.T., S.E.G., A.G., R.S. (Richie Soong), M.F.M.O. and D.S.P.T.; investigation, V.H., D.T., T.Z.T., S.E.G., B.W., A.G., R.S. (Richie Soong), M.F.M.O., D.S.P.T., D.L., B.P., A.W., R.S. (Ross Soo), B.C.G., S.C.L., C.E.C., W.P.Y. and S.O. resources, A.G., R.S. (Richie Soong) and D.S.P.T.; data curation, V.H., D.T., T.Z.T., S.E.G., R.S. (Richie Soong), M.F.M.O. and D.S.P.T.; writing— V.H., D.T., S.E.G., R.S. (Richie Soong) and D.S.P.T.; writing—review and editing V.H., D.T., T.Z.T. and D.S.P.T., visualization D.S.P.T., R.S. (Richie Soong) and A.G. supervision, D.S.P.T., R.S. (Richie Soong) and A.G.; project administration V.H., D.T., T.Z.T., S.E.G., R.S. (Richie Soong), M.F.M.O. and D.S.P.T.; funding acquisition, D.S.P.T., R.S. (Richie Soong) and A.G. All authors have read and agree to the published version of the manuscript.

Funding

This study was supported by the National Medical Research Council Singapore, National Research Foundation Singapore and the Singapore Ministry of Education under its Research Centres of Excellence initiative and a Yong Loo Lin fellowship grant as well as the Singapore Ministry of Health’s National Medical Research Council Clinician Scientist Award (NMRC/CSA-INV/0016/2017) to DSPT.

Conflicts of Interest

D. Tan consultancy fees from Astra Zeneca, Roche, MSD, Merck Serono, Tessa Therapeutics, Eisai and Genmab. Reasearch funding from AstraZeneca, Bayer and Karyopharm and V. Heong consult on the advisory board for Astra Zeneca and Pfizer; the sponsor had no involvement in the design and conduct of the study; collection, management, analysis and interpretation of the data; preparation, review, or approval of the manuscript; or the decision to submit the manuscript for publication.

Appendix A

Table A1. Assessment of tumor purity for each sample using ABSOLUTE algorithm. Abbrev: VAF, variant allele frequency; BRCA, breast carcinoma; CRC, colorectal carcinoma; OV, ovarian carcinoma; UCEC, uterine carcinoma; LIHC, hepatocellular carcinoma; NaN, not computable.
Table A1. Assessment of tumor purity for each sample using ABSOLUTE algorithm. Abbrev: VAF, variant allele frequency; BRCA, breast carcinoma; CRC, colorectal carcinoma; OV, ovarian carcinoma; UCEC, uterine carcinoma; LIHC, hepatocellular carcinoma; NaN, not computable.
SamplePatientPatient labelDiseaseABSOLUTE_call Status (Called = Clonal)ABSOLUTE_purityABSOLUTE_ploidyABSOLUTECancer DNA FactionABSOLUTECoverage for 80% PowerVAF Range-All VariantsVAF Range-Putative Truncal VariantsVAF Range-Private VariantsVAF Range-Branch VariantsMclust.wCN.Cluster#Mclust.noCN.Cluster#VAF_ROC_AUC_TruncalVsNonTruncalVAF_YoudenThreshold_TruncalVsNonTruncalVAF_YoudenThreshold_TruncalVsNonTruncal_AccuracyMclust.cn.threshold_clusterWHighestVAFMclust.cn.threshold_TruncalVsNonTruncal_AccuracyRemark
RE3F1RE3P01CRCcalled0.255.860.49640.14038NaN0.1355319150.244622NaNNaNNaN0.241285714NaNNo Truncal Mutation detected
RE3F2RE3P01CRCcalled0.211.940.2520.134576087NaN0.0992539680.15296694212NaNNaNNaN0.1086NaNNo Truncal Mutation detected
RE3F3RE3P01CRCcalled0.284.020.44490.146770833NaN0.1438241760.200422NaNNaNNaN0.348666667NaNNo Truncal Mutation detected
RE3F4RE3P01CRChigh non-clonal12.06180.383224409NaN0.303893130.46771544722NaNNaNNaN0.447333333NaNNo Truncal Mutation detected
RE4F1RE4P02LUNGcalled0.393.830.55360.1835919280.1914909090.1639444440.1565220.6751306170.07662.277951930.41966666750.49111808---
RE4F2RE4P02LUNGcalled0.53.810.65290.2050121950.2268484850.1451666670.177076923220.6884773660.09168.249158250.213557.35129068---
RE4F3RE4P02LUNGcalled0.225.430.43680.1769419090.1922666670.1230454550.17203125130.6805821370.08968.125996810.2057555.2830941---
RE4F4RE4P02LUNGcalled0.423.850.58340.1800340430.1949393940.14381250.145815789220.7145454550.09667.705627710.599549.89177489---
RE6F1RE6P03LUNGcalled0.24.030.33660.1064318180.1560.0854827590.1505020.6593750.364500.17368.75Only 2 Truncal Mutations detected
RE6F2RE6P03LUNGcalled0.161.960.16660.1067941180.121250.1020370370.130333333030.5750.5500.134563.33333333Only 2 Truncal Mutations detected
RE6F3RE6P03LUNGcalled0.232.120.24470.1055853660.2360.08390.124020.6925675680.52462.50.32733333373.64864865Only 2 Truncal Mutations detected
RE6F4RE6P03LUNGcalled0.161.070.09610.1064814810.218250.0796666670.191375120.75250.52500.07260.5Only 2 Truncal Mutations detected
RE8F1RE8P04LUNGcalled0.324.210.49460.2220098210.2253070270.2108363640.204115108230.5895632470.05250.407472560.29691891953.80586048---
RE8F2RE8P04LUNGhigh non-clonal0.584.340.75280.291516330.3161794590.2018333330.234288194230.6555375570.08455.154487640.56573333354.26535877---
RE8F3RE8P04LUNGhigh entropy0.464.390.65350.2551517790.2752659460.1567164180.210709559230.6580034660.08353.600349140.35998305157.10528998---
RE8F4RE8P04LUNGcalled0.273.330.38470.1921841650.1983135140.1136774190.187166667340.6191405540.04950.581597540.33236363651.06664348---
RE9F1RE9P05LUNGcalled0.223.870.36580.1520.4080.1179387760.15736309513NaNNaNNaN0.117214286NaNNo Truncal Mutation detected
RE9F2RE9P05LUNGcalled0.315.630.56550.2299216420.4143333330.1413882350.26865555623NaNNaNNaN0.416NaNNo Truncal Mutation detected
RE9F3RE9P05LUNGcalled0.253.940.4530.124543860.4070.1112653060.085222NaNNaNNaN0.08NaNNo Truncal Mutation detected
RE9F4RE9P05LUNGcalled0.34.850.51510.1561373390.4083333330.1120222220.16277837812NaNNaNNaN0.1255NaNNo Truncal Mutation detected
RE10F1RE10P06LUNGcalled0.356.420.63550.2335894740.4202424240.1060769230.2807240.93350.223589.710.431571.11---
RE10F2RE10P06LUNGcalled0.243.980.38560.1545392160.2298787880.1177619050.126333333220.83047870.16780.570.19314285772.99---
RE10F3RE10P06LUNGcalled0.265.30.48600.2037619050.4021212120.1030645160.1735220.9621212120.31188.131313130.33671428682.07070707---
RE10F4RE10P06LUNGcalled0.354.340.54430.2038478260.3762727270.082160.247666667220.9422187980.27391.884951210.33285714376.91319979---
RE5F1RE5P07OVnon-aneuploidNaNNaNNaNNaN0.3415153660.4567923810.1197829460.288920635430.9530099390.27893.777481090.54755.3600356---
RE5F2RE5P07OVnon-aneuploidNaNNaNNaNNaN0.3518592810.3854971430.1760891090.354071429230.7675058280.18876.351648350.36418421167.12687313---
RE5F3RE5P07OVnon-aneuploidNaNNaNNaNNaN0.3724596890.427596190.1384380950.315649351130.8586446890.25882.340659340.37880769279.35531136---
RE5F4RE5P07OVnon-aneuploidNaNNaNNaNNaN0.3912122570.4315447620.174281250.3839375260.818432540.2878.180555560.42409677468.52380952---
RE11F1RE11P08OVcalled0.273.60.4480.1826776320.2147070710.1093809520.174272727230.8259958070.09380.055269680.45853.09700781---
RE11F2RE11P08OVcalled0.253.020.34480.220969880.2732121210.1191041670.206105263430.808306950.12279.059249210.38357.59837178---
RE11F3RE11P08OVcalled0.255.060.46590.2191648940.2796262630.1353731340.202272727130.8060379070.12379.048916130.1517577.1422086---
RE11F4RE11P08OVcalled0.224.740.4650.2253370790.2940505050.10380.220416667130.8374248820.15280.392532920.265567.40825981---
RE13F1RE13P09OVcalled0.475.470.71410.1966666670.2425781250.1261363640.215220.7831662740.13878.198702830.37033333355.76356132---
RE13F2RE13P09OVcalled0.325.190.55500.2546319440.3173593750.1031515150.275574468130.7411132810.22971.8750.25941176566.40625---
RE13F3RE13P09OVcalled0.315.30.55520.2673333330.344781250.1128809520.299893617130.7388518260.17272.059339890.22209090964.9315309---
RE13F4RE13P09OVcalled0.335.210.57490.2958759690.3555156250.1727619050.267886364120.7052884620.14569.194711540.22958333364.42307692---
RE7F1RE7P10BRCAcalled0.822.30.84110.317900990.448260870.17731250.425285714120.8328063240.28982.747035570.32579.48616601---
RE7F2RE7P10BRCAcalled0.782.320.81120.30820.4233478260.1318108110.483857143120.8631422920.30486.561264820.222583.05335968---
RE7F3RE7P10BRCAcalled0.295.980.55590.3119325840.4336304350.111593750.385818182120.8612234580.389.762386250.26986.27401416---
RE7F4RE7P10BRCAcalled0.223.610.33590.191711340.2379782610.1436818180.189571429120.8056265980.13283.013640240.11366666781.15942029---
RE12F1RE12P11UCEChigh non-clonal0.632.240.66140.221503448* 0.3232266666666670.112514286NaN220.9031428570.16488.285714290.3107573.85714286---
RE12F4RE12P11UCECcalled0.664.450.81260.267685185* 0.40720.14737931NaN330.8639846740.19482.597701150.916552.09195402---
RE14F1RE14P12UCECcalled0.24.620.37680.2174022990.276520.1100967740.279166667430.8851351350.12186.837837840.51852.2972973---
RE14F2RE14P12UCECcalled0.265.790.51620.2543218390.304180.13120.303083333230.8175675680.17882.837837840.522551.94594595---
RE14F3RE14P12UCECcalled0.234.90.43620.2950.41680.1018484850.3215130.8821739130.17485.043478260.459560.65217391---
RE14F4RE14P12UCECcalled0.225.480.44680.2955212770.4230.1108571430.305444444130.8886363640.19788.909090910.49563.72727273---
RE15F1RE15P13LIHCcalled0.434.050.61350.2661760.350.2255967740.304672131130.626016260.75500.252746.54471545Only 1 Truncal Mutation detected
RE15F2RE15P13LIHCcalled0.231.960.23460.1792962960.24750.1768421050.164333333040.720.667500.42866666773Only 1 Truncal Mutation detected
RE15F3RE15P13LIHCcalled0.232.040.23470.1660759490.36350.150.163171875000.948051948150NaNNaNOnly 1 Truncal Mutation detected
RE15F4RE15P13LIHCcalled0.23.910.32660.11890.280.0848888890.1479000.9464285710.66750NaNNaNOnly 1 Truncal Mutation detected
* data analysed based on two sample biopsies instead of four.

References

  1. Jamal-Hanjani, M.; Wilson, G.A. Tracking the Evolution of Non–Small-Cell Lung Cancer. N. Engl. J. Med. 2017, 376, 2109–2121. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Greaves, M. Evolutionary determinants of cancer. Cancer Discov. 2015, 5, 806–820. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Burrell, R.A.; McGranahan, N. The causes and consequences of genetic heterogeneity in cancer evolution. Nature 2013, 501, 338–345. [Google Scholar] [CrossRef] [PubMed]
  4. Gerlinger, M.; Rowan, A.J. Intratumor Heterogeneity and Branched Evolution Revealed by Multiregion Sequencing. N. Engl. J. Med. 2012, 366, 883–892. [Google Scholar] [CrossRef] [Green Version]
  5. Dagogo-Jack, I.; Shaw, A.T. Tumour heterogeneity and resistance to cancer therapies. Nat. Rev. Clin. Oncol. 2018, 15, 81–94. [Google Scholar] [CrossRef]
  6. Turner, N.C.; Reis-Filho, J.S. Genetic heterogeneity and cancer drug resistance. Lancet Oncol. 2012, 13, 178–185. [Google Scholar] [CrossRef]
  7. Morris, L.G.T.; Riaz, N. Pan-cancer analysis of intratumor heterogeneity as a prognostic determinant of survival. Oncotarget 2016, 7, 10051–10063. [Google Scholar] [CrossRef] [Green Version]
  8. McGranahan, N.; Swanton, C. Biological and therapeutic impact of intratumor heterogeneity in cancer evolution. Cancer Cell 2015, 27, 15–26. [Google Scholar] [CrossRef] [Green Version]
  9. Gerlinger, M.; Horswell, S. Genomic architecture and evolution of clear cell renal cell carcinomas defined by multiregion sequencing. Nat. Genet. 2014, 46, 225–233. [Google Scholar] [CrossRef]
  10. Yap, T.A.; Gerlinger, M. Intratumor heterogeneity: Seeing the wood for the trees. Sci. Transl. Med. 2012, 4, 127–137. [Google Scholar] [CrossRef] [Green Version]
  11. Schmitt, M.W.; Loeb, L.A. The influence of subclonal resistance mutations on targeted cancer therapy. Nat. Rev. Clin. Oncol. 2016, 13, 335–347. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Carter, S.L.; Cibulskis, K. Absolute quantification of somatic DNA alterations in human cancer. Nat. Biotechnol. 2012, 30, 413–421. [Google Scholar] [CrossRef] [PubMed]
  13. Shi, Y.; Au, J.S.-K. A prospective, molecular epidemiology study of EGFR mutations in Asian patients with advanced non-small-cell lung cancer of adenocarcinoma histology (PIONEER). J. Thorac. Oncol. 2014, 9, 154–162. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Heong, V.; Syn, N.L. Value of a molecular screening program to support clinical trial enrollment in Asian cancer patients: The Integrated Molecular Analysis of Cancer (IMAC) Study. Int. J. Cancer 2018, 142, 1890–1900. [Google Scholar] [CrossRef] [Green Version]
  15. Wang, S.; Cang, S. Third-generation inhibitors targeting EGFR T790M mutation in advanced non-small cell lung cancer. J. Hematol. Oncol. 2016, 9, 34. [Google Scholar] [CrossRef] [Green Version]
  16. Janku, F.; Wheler, J.J. PIK3CA mutation H1047R is associated with response to PI3K/AKT/mTOR signaling pathway inhibitors in early-phase clinical trials. Cancer Res. 2013, 73, 276–284. [Google Scholar] [CrossRef] [Green Version]
  17. Lohr, J.G.; Stojanov, P. Widespread genetic heterogeneity in multiple myeloma: Implications for targeted therapy. Cancer Cell 2014, 25, 91–101. [Google Scholar] [CrossRef] [Green Version]
  18. Bailey, M.H.; Tokheim, C. Comprehensive Characterization of Cancer Driver Genes and Mutations. Cell 2018, 173, 371–385.e18. [Google Scholar] [CrossRef] [Green Version]
  19. Jovelet, C.; Ileana, E. Circulating Cell-Free Tumor DNA Analysis of 50 Genes by Next-Generation Sequencing in the Prospective MOSCATO Trial. Clin. Cancer Res. 2016, 22, 2960–2968. [Google Scholar] [CrossRef] [Green Version]
  20. Marchetti, A.; Palma, J.F. Early Prediction of Response to Tyrosine Kinase Inhibitors by Quantification of EGFR Mutations in Plasma of NSCLC Patients. J. Thorac. Oncol. 2015, 10, 1437–1443. [Google Scholar] [CrossRef] [Green Version]
  21. Sorensen, B.S.; Wu, L. Monitoring of epidermal growth factor receptor tyrosine kinase inhibitor-sensitizing and resistance mutations in the plasma DNA of patients with advanced non-small cell lung cancer during treatment with erlotinib. Cancer 2014, 120, 3896–3901.20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [Green Version]
  23. McKenna, A.; Hanna, M. The genome analysis toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010, 20, 1297–1303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Cibulskis, K.; Lawrence, M.S. Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples. Nat. Biotechnol. 2013, 31, 213–219. [Google Scholar] [CrossRef] [PubMed]
  25. Ramos, A.H.; Lichtenstein, L. Oncotator: Cancer variant annotation tool. Hum. Mutat. 2015, 36, E2423–E2429. [Google Scholar] [CrossRef]
  26. Koboldt, D.C.; Zhang, Q. VarScan 2: Somatic mutation and copy number alteration discovery in cancer by exome sequencing. Genome Res. 2012, 22, 568–576. [Google Scholar]
  27. Seshan, V.E.; Olshen, A. DNAcopy: DNA copy number data analysis. R package version 1.62.0. Available online: https://bioconductor.org/biocLite.R (accessed on 14 August 2017).
  28. Futreal, P.A.; Coin, L. A census of human cancer genes. Nat. Rev. Cancer 2004, 4, 177–183. [Google Scholar] [CrossRef]
  29. Wellcome Sanger Institute. Cancer Gene Census. 2011. Available online: http://www.sanger.ac.uk/genetics/CGP/Census/ (accessed on 10 June 2018).
  30. Lawrence, M.S.; Stojanov, P. Discovery and saturation analysis of cancer genes across 21 tumour types. Nature 2014, 505, 495–501. [Google Scholar] [CrossRef] [Green Version]
  31. Ally, A.; Balasundaram, M. Comprehensive and Integrative Genomic Characterization of Hepatocellular Carcinoma. Cell 2017, 169, 1327–1341.e23. [Google Scholar] [CrossRef]
  32. Thermo Fisher. Ion AmpliSeq Cancer Hotspot Panel v2. Available online: https://www.thermofisher.com/order/catalog/product/4475346 (accessed on 28 November 2018).
  33. Illumina. TruSight® Cancer Panel. Available online: https://www.illumina.com/products/by-type/clinical-research-products/trusight-cancer.html (accessed on 26 July 2018).
  34. Foundation Medicine. FoundationOneTM Cancer Gene Panel. Available online: https://www.foundationmedicine.com/genomic-testing/foundation-one (accessed on 11 August 2018).
  35. Felsenstein, J. PHYLIP (Phylogeny Inference Package) Version 3.6; Distributed by the Author; Department of Genome Sciences, University of Washington: Seattle, WA, USA, 2005; Available online: http://evolution.genetics.washington.edu/phylip.html (accessed on 15 August 2017).
  36. Von Hoff, D.D.; Stephenson, J.J. Pilot study using molecular profiling of patients’ tumors to find potential targets and select treatments for their refractory cancers. J. Clin. Oncol. 2010, 28, 4877–4883. [Google Scholar] [CrossRef]
Figure 1. Representative workflow of the processing pipeline. (A) Whole-exome sequencing was performed on all germline and MRTB samples obtained from each patient. Bioinformatics analysis was subsequently performed: (B) alignment of sequence reads; (C) somatic variant calling and variant annotation; (D) generation of non-synonymous somatic mutational landscape across all patients; (E) identification of truncal and branch variants present in each patient; (F) curation of statistically significant somatic cancer driver mutations; (G) construction of phylogenetic trees from non-synonymous somatic variants; (H) filtering of genetic variants using AmpliSeq™, TruSight® and FoundationOne™ cancer gene panels; (I) statistical saturation analysis to determine the minimum number of MRTB samples needed (to alleviate challenges associated with ITH) in relation to the gene panel used; (J) copy number alterations analysis; (K) estimation of cancer cell fraction (CCF); (L) prediction of putative truncal variants using two different threshold metrics, namely variant allele frequency and CCF; (M) informed targeted therapies were performed based on patients’ mutational profile that reflects genes from the AmpliSeq™ cancer gene panel. MRTB: multi- region tumor biopsy; ITH: intratumor heterogeneity; TB: the number of MRTB samples resected from the patient; GL: the type of germline sample; BL: whole blood sample; BS: buccal swab sample. CRC: Colorectal cancer; NSCLC: Non-small cell lung cancer; OV: Ovarian Cancer; BC: Breast cancer; UC: Uterine Cancer; HCC: Hepatocellular Carcinoma; P: Patient.
Figure 1. Representative workflow of the processing pipeline. (A) Whole-exome sequencing was performed on all germline and MRTB samples obtained from each patient. Bioinformatics analysis was subsequently performed: (B) alignment of sequence reads; (C) somatic variant calling and variant annotation; (D) generation of non-synonymous somatic mutational landscape across all patients; (E) identification of truncal and branch variants present in each patient; (F) curation of statistically significant somatic cancer driver mutations; (G) construction of phylogenetic trees from non-synonymous somatic variants; (H) filtering of genetic variants using AmpliSeq™, TruSight® and FoundationOne™ cancer gene panels; (I) statistical saturation analysis to determine the minimum number of MRTB samples needed (to alleviate challenges associated with ITH) in relation to the gene panel used; (J) copy number alterations analysis; (K) estimation of cancer cell fraction (CCF); (L) prediction of putative truncal variants using two different threshold metrics, namely variant allele frequency and CCF; (M) informed targeted therapies were performed based on patients’ mutational profile that reflects genes from the AmpliSeq™ cancer gene panel. MRTB: multi- region tumor biopsy; ITH: intratumor heterogeneity; TB: the number of MRTB samples resected from the patient; GL: the type of germline sample; BL: whole blood sample; BS: buccal swab sample. CRC: Colorectal cancer; NSCLC: Non-small cell lung cancer; OV: Ovarian Cancer; BC: Breast cancer; UC: Uterine Cancer; HCC: Hepatocellular Carcinoma; P: Patient.
Cancers 12 01599 g001
Figure 2. Mutational landscape of patients across six cancer types. (A) Boxplot illustrating nssML. A cross (+) represents the mean value of the data. (B) Line chart and stacked bar chart representing the number and proportion of truncal/branch variants, respectively. (C) Representative phylogenetic tree and mutation heatmap for patient P03. Trunk, branch and private branches of the tree signify mutations that occur in all, in some but not all, and only one MRTB sample(s) resected from the patient, respectively. Heatmap demonstrates the presence (green: private; red: branch; blue: trunk) or absence (gray) of NSS mutations in each MRTB sample. Bx denotes an MRTB sample with identification number x. The total number of NSS, truncal (percentage), branch (percentage), and private (percentage) mutations are denoted by ‘n’, ‘C’, ‘S’, and ‘P’, respectively. (D) Heatmap illustrating the presence and absence (gray) of CNAs for patients with OV. Large-scale amplifications and deletions are represented with areas filled with green and blue, respectively. (E) ssCDMs for OC and their associated AF and CCF. CV: clonal (truncal) variant; Y: yes; N: no; AF: allele frequency; CCF: cancer cell fraction.
Figure 2. Mutational landscape of patients across six cancer types. (A) Boxplot illustrating nssML. A cross (+) represents the mean value of the data. (B) Line chart and stacked bar chart representing the number and proportion of truncal/branch variants, respectively. (C) Representative phylogenetic tree and mutation heatmap for patient P03. Trunk, branch and private branches of the tree signify mutations that occur in all, in some but not all, and only one MRTB sample(s) resected from the patient, respectively. Heatmap demonstrates the presence (green: private; red: branch; blue: trunk) or absence (gray) of NSS mutations in each MRTB sample. Bx denotes an MRTB sample with identification number x. The total number of NSS, truncal (percentage), branch (percentage), and private (percentage) mutations are denoted by ‘n’, ‘C’, ‘S’, and ‘P’, respectively. (D) Heatmap illustrating the presence and absence (gray) of CNAs for patients with OV. Large-scale amplifications and deletions are represented with areas filled with green and blue, respectively. (E) ssCDMs for OC and their associated AF and CCF. CV: clonal (truncal) variant; Y: yes; N: no; AF: allele frequency; CCF: cancer cell fraction.
Cancers 12 01599 g002
Figure 3. Number, PPV and in silico prediction accuracy of truncal variants across different gene panels. Five gene panels were scrutinized, namely WES NSS, CGC, AmpliSeq™, TruSight® and FoundationOne™ cancer gene panels. (A) Boxplot illustrating the number of PTVs across different numbers of MRTB samples analyzed concurrently. (B) PPV of PTVs in relation to the number of MRTB samples interrogated simultaneously. (C) Best average prediction accuracy of PTVs across different cancer types. Two types of thresholds were used to classify variants into either truncal or branch, namely AF and CCF. Based on the respective threshold, the best average prediction accuracy achievable (within the defined search domain) among all patients with the same cancer type (across different gene panels) is portrayed above. A single asterisk (*) denotes p < 0.05, double asterisks (**) signify p < 0.01, while triple asterisks (***) indicate p < 0.001. A cross (+) represents the mean value of the data. ‘Not available’ signifies that no variants that are associated with the specific gene panel were found.
Figure 3. Number, PPV and in silico prediction accuracy of truncal variants across different gene panels. Five gene panels were scrutinized, namely WES NSS, CGC, AmpliSeq™, TruSight® and FoundationOne™ cancer gene panels. (A) Boxplot illustrating the number of PTVs across different numbers of MRTB samples analyzed concurrently. (B) PPV of PTVs in relation to the number of MRTB samples interrogated simultaneously. (C) Best average prediction accuracy of PTVs across different cancer types. Two types of thresholds were used to classify variants into either truncal or branch, namely AF and CCF. Based on the respective threshold, the best average prediction accuracy achievable (within the defined search domain) among all patients with the same cancer type (across different gene panels) is portrayed above. A single asterisk (*) denotes p < 0.05, double asterisks (**) signify p < 0.01, while triple asterisks (***) indicate p < 0.001. A cross (+) represents the mean value of the data. ‘Not available’ signifies that no variants that are associated with the specific gene panel were found.
Cancers 12 01599 g003
Table 1. Clinical details of patients who received therapy targeting their actionable truncal mutation.
Table 1. Clinical details of patients who received therapy targeting their actionable truncal mutation.
Cancer TypePatientAgeSexNo. of MRTB Samples with Abnormality of InterestNo. of MRTB Samples that CCF Metric Classified as ClonalTargeted AbnormalityTherapeutic InterventionPFS (Months)PFS RatioRadiological RECIST (v1.1) Response
Initial TherapyActionable Truncal Mutation-Directed Therapy
NSCLCP0674M4/43/4EGFR T790MT790M inhibitor2.525.510.2PR
NSCLCP0543M3/4NAEGFR T790MT790M inhibitor2.13.61.71SD
BCP1041F4/44/4PIK3CA H1047RPI3Kα/β inhibitor21.90.95PD
UCP1146F2/22/2AKT1 E17Kpan-AKT inhibitor46.11.53SD
M: male; F: female; SD: stable disease; PR: partial response; PD: progressive disease (based on RECIST v1.1); PFS: progression free survival. RECIST: Response evaluation criteria in solid tumors; EGFR: Epidermal growth factor receptor; PIK3CA: phosphoinositide-3-kinase catalytic alpha polypeptide; AKT: RAC-alpha serine/threonine-protein kinase.
Table 2. Thresholds used to determine the clonality of the variants.
Table 2. Thresholds used to determine the clonality of the variants.
PanelCCFAF
FoundationOne0.920.13
AmpliSeq0.920.15
TruSight0.960.13
WES10.16

Share and Cite

MDPI and ACS Style

Heong, V.; Tay, D.; Goh, S.E.; Wee, B.; Tan, T.Z.; Soo, R.; Pang, B.; Lim, D.; Gopinathan, A.; Ow, S.; et al. Whole Exome Sequencing of Multi-Regional Biopsies from Metastatic Lesions to Evaluate Actionable Truncal Mutations Using a Single-Pass Percutaneous Technique. Cancers 2020, 12, 1599. https://doi.org/10.3390/cancers12061599

AMA Style

Heong V, Tay D, Goh SE, Wee B, Tan TZ, Soo R, Pang B, Lim D, Gopinathan A, Ow S, et al. Whole Exome Sequencing of Multi-Regional Biopsies from Metastatic Lesions to Evaluate Actionable Truncal Mutations Using a Single-Pass Percutaneous Technique. Cancers. 2020; 12(6):1599. https://doi.org/10.3390/cancers12061599

Chicago/Turabian Style

Heong, Valerie, Darwin Tay, Shane Ee Goh, Bernard Wee, Tuan Zea Tan, Ross Soo, Brendan Pang, Diana Lim, Anil Gopinathan, Samuel Ow, and et al. 2020. "Whole Exome Sequencing of Multi-Regional Biopsies from Metastatic Lesions to Evaluate Actionable Truncal Mutations Using a Single-Pass Percutaneous Technique" Cancers 12, no. 6: 1599. https://doi.org/10.3390/cancers12061599

APA Style

Heong, V., Tay, D., Goh, S. E., Wee, B., Tan, T. Z., Soo, R., Pang, B., Lim, D., Gopinathan, A., Ow, S., Chee, C. E., Goh, B. C., Lee, S. C., Yong, W. P., Wong, A., Omar, M. F. M., Soong, R., & Tan, D. S. (2020). Whole Exome Sequencing of Multi-Regional Biopsies from Metastatic Lesions to Evaluate Actionable Truncal Mutations Using a Single-Pass Percutaneous Technique. Cancers, 12(6), 1599. https://doi.org/10.3390/cancers12061599

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