Open Access

The impact of tumor profiling approaches and genomic data strategies for cancer precision medicine

  • Andrea Garofalo1, 2,
  • Lynette Sholl3,
  • Brendan Reardon1, 2,
  • Amaro Taylor-Weiner2,
  • Ali Amin-Mansour2,
  • Diana Miao1, 2,
  • David Liu1, 2,
  • Nelly Oliver1,
  • Laura MacConaill1, 3,
  • Matthew Ducar3,
  • Vanesa Rojas-Rudilla3,
  • Marios Giannakis1, 2,
  • Arezou Ghazani1,
  • Stacy Gray1,
  • Pasi Janne1,
  • Judy Garber1,
  • Steve Joffe4,
  • Neal Lindeman3,
  • Nikhil Wagle1, 2, 5,
  • Levi A. Garraway1, 2, 5Email author and
  • Eliezer M. Van Allen1, 2, 5Email author
Contributed equally
Genome Medicine20168:79

https://doi.org/10.1186/s13073-016-0333-9

Received: 9 March 2016

Accepted: 8 July 2016

Published: 26 July 2016

Abstract

Background

The diversity of clinical tumor profiling approaches (small panels to whole exomes with matched or unmatched germline analysis) may engender uncertainty about their benefits and liabilities, particularly in light of reported germline false positives in tumor-only profiling and use of global mutational and/or neoantigen data. The goal of this study was to determine the impact of genomic analysis strategies on error rates and data interpretation across contexts and ancestries.

Methods

We modeled common tumor profiling modalities—large (n = 300 genes), medium (n = 48 genes), and small (n = 15 genes) panels—using clinical whole exomes (WES) from 157 patients with lung or colon adenocarcinoma. We created a tumor-only analysis algorithm to assess germline false positive rates, the impact of patient ancestry on tumor-only results, and neoantigen detection.

Results

After optimizing a germline filtering strategy, the germline false positive rate with tumor-only large panel sequencing was 14 % (144/1012 variants). For patients whose tumor-only results underwent molecular pathologist review (n = 91), 50/54 (93 %) false positives were correctly interpreted as uncertain variants. Increased germline false positives were observed in tumor-only sequencing of non-European compared with European ancestry patients (p < 0.001; Fisher’s exact) when basic germline filtering approaches were used; however, the ExAC database (60,706 germline exomes) mitigated this disparity (p = 0.53). Matched and unmatched large panel mutational load correlated with WES mutational load (r2 = 0.99 and 0.93, respectively; p < 0.001). Neoantigen load also correlated (r2 = 0.80; p < 0.001), though WES identified a broader spectrum of neoantigens. Small panels did not predict mutational or neoantigen load.

Conclusions

Large tumor-only targeted panels are sufficient for most somatic variant identification and mutational load prediction if paired with expanded germline analysis strategies and molecular pathologist review. Paired germline sequencing reduced overall false positive mutation calls and WES provided the most neoantigens. Without patient-matched germline data, large germline databases are needed to minimize false positive mutation calling and mitigate ethnic disparities.

Keywords

Genomics Precision medicine Disparities Immuno-oncology Neoantigens Panel testing

Background

The mapping of the human genome, together with the advent of massively parallel sequencing, has accelerated discovery of driver genetic alterations in cancer and the development of drugs to target or otherwise exploit these events [1]. Multiple tumor profiling approaches that leverage these advances have entered the clinic. Such assays often consist of targeted sequencing panels that query a subset of typically 200–500 genes implicated in cancer biology or clinical management [28]. Alternatively, panels that emphasize rapid turnaround time by profiling smaller gene sets (n = 15–48 genes) have also emerged [9, 10]. On the other end of the spectrum, clinical whole-exome sequencing (WES; n ~ 20,000 genes) of matched tumor and germline samples has been studied through prospective sequencing efforts [1113]. However, the benefits and limitations of these different sequencing strategies remain incompletely understood.

Understanding the differences in genomic results between different tumor profiling approaches will become increasingly important as the cancer genome is leveraged to stratify patients for new therapeutic strategies. For example, unlike targeted therapies linked to specific genetic lesions (e.g., epidermal growth factor receptor mutations and inhibitors), immune targeting strategies, such as checkpoint blockade or personalized cancer vaccines, may require large-scale ascertainment of mutational and neoantigen loads and individual mutation-associated neoantigens for personalized cancer vaccine development [1418]. One effort demonstrated the ability of two large gene panels (315 or 573 genes) to predict mutational load for immunotherapy response in pilot patient cohorts [19], and another effort demonstrated the ability of one large gene panel (341 genes) to predict DNA mismatch repair protein deficient tumors through mutational load [20], although a systematic characterization of different tumor profiling strategies for both mutation load and personal neoantigen identification should inform their relative utilities for stratifying patients in emerging cancer precision medicine frameworks.

Moreover, although sequencing of paired normal blood or tissue samples is standard practice for research-oriented WES applications, many targeted panel approaches do not include matched normal samples [2, 3, 9, 21, 22]. Together with the limited ancestral diversity in many existing germline databases, this absence of paired normals has raised concerns for the potential of increased false positive somatic mutation calls that are actually germline [23, 24].

To investigate these issues, we analyzed clinical sequencing data from 157 patients with advanced lung and colon adenocarcinoma to ascertain the relative merits of distinct tumor profiling approaches.

Methods

Patients and tumor specimens

All patients consented to an institutional review board-approved protocol that allows comprehensive genetic analysis of tumor and germline samples (Dana-Farber Cancer Institute #12-078). Ancestry status was self-reported. Samples were selected from pathology archives by a board-certified anatomic pathologist based on sample size, tumor purity, and timing relative to date of study enrollment and analyzed by the Center for Advanced Molecular Diagnostics (CAMD) at Brigham and Women’s Hospital (BWH), a Clinical Laboratory Improvements Amendments (CLIA)-certified laboratory. Tumor content was estimated by an anatomic pathologist from corresponding stained slides and only samples with at least 20 % malignant cells were analyzed. DNA was isolated with a commercial kit (QIAamp DNA Mini Kit, Qiagen, Valencia, CA, USA) following the manufacturer’s instructions. DNA was quantified (PicoGreen, ThermoFisher Scientific, Waltham, MA, USA) and samples with at least 50 ng/μL of DNA proceeded to library preparation.

Whole exome sequencing

WES from formalin-fixed, paraffin embedded (FFPE) samples was performed as described previously [12]. Whole-exome capture libraries were constructed from tumor and normal DNA after sample shearing, end repair, phosphorylation, and ligation to barcoded sequencing adaptors. DNA was then subjected to solution-phase hybrid capture using Agilent baits. The samples were multiplexed and sequenced using Illumina HiSeq technology. All BAM files were deposited in dbGap phs001075.v1.p1.

Genomic analysis

Sequence data processing and quality control

WES data were processed using established analytical pipelines at the Broad Institute [12]. A BAM file was produced using the Picard pipeline (http://broadinstitute.github.io/picard/), which aligns tumor and normal sequences to the hg19 human genome build from raw Illumina reads using the BWA aligner (version 0.5.9-tpx [0.5.9 with an internal patch to support threading]). BAM files were uploaded into the Firehose pipeline (http://www.broadinstitute.org/cancer/cga/Firehose), which manages input and output files to be executed by GenePattern [25]. Quality control modules within Firehose were applied to all sequencing data for comparison of the origin for tumor and normal genotypes and to assess fingerprinting concordance. Cross-contamination of samples was estimated using ContEst [26]; those with >5 % contamination were excluded from subsequent analysis.

Somatic alterations and downsampling

MuTect (version 1.1.6) [27] was applied to identify somatic single-nucleotide variants. Artifacts introduced by DNA oxidation during sequencing or from FFPE were computationally removed using a filter-based method [28]. Annotation of identified variants was done using Oncotator (version 1.2.7.0) [29]. Representative large (n = 300) [3], medium (n = 48), and small (n = 15) [9] gene sets were defined through review of literature (Additional file 1: Table S1; Additional file 2: Table S2; Additional file 3: Table S3). The aggregate mutation data across the whole exome were collected in two files, one containing all tumor-germline matched calls (Additional file 4: Table S4) and a second containing tumor-only calls (Additional file 5: Table S5). “Downsampling” was performed on the aggregate somatic mutation alteration data files to derive subsets of WES data for gene sets represented by the large, medium, and small gene lists defined by the respective panels. For example, to model the 300 gene panel, the set of mutations from the entire WES data was restricted to only consider events in those 300 genes. The analyses were performed using the R statistical software.

Tumor-only and germline analysis

Tumor-normal paired mutation data were taken to be the set of all true somatic mutations for each patient. To ensure only high-confidence mutation calls were considered, only mutations with an allelic fraction ≥5 % and sequencing or FFPE artifact filtering strategies described above were considered. Tumor-only mutation calling was performed by using MuTect and pairing the tumor WES with an FFPE germline whole exome from another patient to reduce false positive calls introduced by artifacts from the sequencing process, as described previously [23]. Variants were removed if they were present in combinations of dbSNP (build 134) [30] and 1000 Genomes (phase 1, version 3) [31] using the Oncotator annotation algorithm [29], along with the ExAC (version 0.3) [32] databases. Mutations were rescued if listed as somatic in COSMIC (version 74) at least one, three, five, or ten times, for increasing stringency [33]. Positive predictive values for each filter were calculated by dividing the number of true somatic mutations in the post-filtering mutation data by the total number of unfiltered mutations. Sensitivity was calculated by dividing the number of true somatic mutations in the tumor-only post-filtering data by the total number of known somatic mutations in the paired mutation data. To obtain the set of known germline variants, the GATK HaplotypeCaller (version 3.1.1) [34, 35] was applied to germline sequence BAMs to identify germline single-nucleotide polymorphisms (SNPs) using WES data from each patient. Unfiltered germline variants in tumor-only targeted panel data after application of various database filters were identified by comparison with exome germline SNP data.

Orthogonal large panel molecular pathology review

A subset of cases (n = 91) underwent separate testing with an academic lab large gene panel (“OncoPanel”), followed by molecular pathology review [3]. Specifically, after variants were identified by computational approaches, an individual molecular pathologist reviewed each variant and assigned a tier based on clinical actionability and to determine whether the variant was likely somatic or germline and whether there were any clinical actions for the variant. The four tiers in this system are:
  1. Tier 1

    The alteration has well-established published evidence confirming clinical utility in this tumor type in at least one of the following contexts: predicting response to treatment with a US Food and Drug Administration (FDA)-approved therapy; assessing prognosis; establishing a definitive diagnosis; or conferring an inherited increased risk of cancer to this patient and family.

     
  2. Tier 2

    The alteration may have clinical utility in at least one of the following contexts: selection of an investigational therapy in clinical trials for this cancer type; limited evidence of prognostic association; supportive of a specific diagnosis; proven association of response to treatment with an FDA-approved therapy in a different type of cancer; or similar to a different mutation with a proven association with response to treatment with an FDA-approved therapy in this type of cancer.

     
  3. Tier 3

    The alteration is of uncertain clinical utility but may have a role as suggested by at least one of the following: demonstration of association with response to treatment in this cancer type in preclinical studies (e.g., in vitro studies or animal models); alteration in a biochemical pathway that has other known, therapeutically targetable alterations; alteration in a highly conserved region of the protein predicted, in silico, to alter protein function; or selection of an investigational therapy for a different cancer type.

     
  4. Tier 4

    The alteration is novel or its significance has not been studied in cancer.

     

Mutation rates and neoantigens

The mutation rate for each sample was calculated by dividing the number of bona fide mutation calls post-filtering by the total genomic territory sequenced (in megabases). Germline exome data from each patient were used to genotype human leukocyte antigen (HLA) loci with POLYSOLVER [36]. Patient HLA genotypes and matched exome mutation data were used as inputs for NetMHCPan [37] to generate predicted binding affinities of somatic mutations linked to specific MHC class I molecules. Predicted mutation-associated neoantigens were defined as all 9- and 10-amino-acid peptides resulting from tumor-specific mutations with predicted HLA binding affinities <500 nM using NetMHCpan (v2.4). Downsampled panel data were queried for exome neoantigens to determine the fraction of putative neoantigens observed in WES that were recoverable from panel mutation calls.

Statistical tests

The sample size was based on available material, and thus there was not an a priori power calculation. Comparisons of germline false positive rates between ancestry groups was performed with two-sided Fisher’s Exact test. Pearson correlation tests were performed for mutational and neoantigen load comparisons in the three panel settings.

Results

Reducing false-positive germline variants in tumor-only analysis

The study included 157 patients, who underwent clinical sequencing: 75 with colorectal adenocarcinoma and 82 with lung adenocarcinoma. Clinical tumor and germline WES [12] produced a mean coverage of 154× and 133× in tumor and normal DNA, respectively (Additional file 6: Table S6). The combination of sequencing depth and tumor purity (“Methods”) enabled mutation detection in these cases. Large (n = 300) [3], medium (n = 48), and small (n = 15) [9] gene panel data were produced by creating subsets of the whole exome mutation data (“downsampling”) to simultaneously model the different gene sets captured in multiple representative academic and commercial efforts. Matched tumor–germline sequencing revealed a median of 75 (interquartile range of first and third quartiles [IQR] = 55–134), 4 (IQR = 3–6), 2 (IQR = 1–3), and 1 (IQR = 1–2) mutations per patient for WES, large, medium, and small panels, respectively. Unmatched tumor-only sequencing produced a median of 445 (IQR = 404–531), 10 (IQR = 8–13), 3 (IQR = 2–5), and 2 (IQR = 1–2) mutations per patient. Thus, the proportion of putative somatic variants was increased in tumor-only sequencing data under all conditions.

For all targeted panel options, both sensitivity and positive predictive value (PPV) could be optimized by using an analytical pipeline that consisted of an unmatched germline sample, the largest publically available germline WES database (ExAC) [32], and recovery of somatic mutations with COSMIC at the highest frequency threshold (n ≥ 10 events) to recover mutational hotspots [38] (Fig. 1a, b; “Methods”; Additional file 7: Table S7). For large tumor-only targeted panels, this filtering led to 14 % (144/1012 variants) of putative somatic mutations that were actually germline false positive variants (Fig 1c). For unmatched WES, the germline false positive rate with this approach was even higher (18 %; 5282/29,738 variants).
Fig. 1

Germline false positives in tumor-only clinical sequencing. Sensitivity and positive predictive value (PPV) curves for multiple germline filtering strategies identifies optimal approaches for unmatched large targeted panel testing (a) and whole-exome sequencing (b). For 91 patients, germline exome data were used to identify false positives post-filtering. Subsequent molecular pathologist review of variants was performed on individual cases to further classify putative germline variants. With molecular pathology review, 50/54 false positive variants were correctly classified as unknown (“tier 4”), with the remaining variants classified as having uncertain (n = 3; “tier 3”) or potential (n = 1; “tier 2”) clinical utility (c, d). Please see “Methods” for detailed descriptions of the four-tier classification schema

A subset of these cases (n = 91) underwent orthogonal molecular pathologist review of variants (“Methods”). The addition of pathologist review after in silico analysis resulted in 50/54 (93 %) false positives interpreted as unknown variants (tier 4) that may be germline false positives rather than somatic alterations (Fig. 1c, d). Of the four remaining germline false positive results, three (RB1 K874Q, MSH6 S405C, BRCA2 splice site mutation at the junction of exons 2 and 3 (g.32890665G > A)) were classified as uncertain clinical utility and only one, a known pathogenic variant associated with hereditary cancer syndromes (TP53 R248Q), was classified as having potential clinical utility based on negative prognostic implications (Fig. 1d).

The impact of ancestry on germline false positives in tumor-only analysis

The use of germline databases is a critical component for removing false positive germline calls in tumor-only panel sequencing (“Methods”); however, the representation of non-European ancestry in these databases is incomplete [39]. Therefore, we next sought to measure the variation in false positive rates in populations with different self-reported ancestries. When the dbSNP database was used as a filter in large panel tumor profiling analysis [30], germline false positives were significantly increased in tumor-only sequencing of non-European compared with European ancestry patients (odds ratio [OR] = 2.52, p < 0.001; Fisher’s exact; Fig. 2a). While the median number of false positives among the two populations was the same (n = 2), 32 % (6/19) of non-white patients had more than five false positives, whereas 5 % (7/132) of white patients had more than five germline false positives.. A similar relationship was observed with the use of 1000 Genomes (OR = 1.83, p < 0.001; Fisher’s exact; Fig. 2b). However, use of ExAC [32], a public database of 60,706 germline exomes that represents an order-of-magnitude increase in germline variant data compared with other databases, mitigated this disparity (OR = 1.19, p = 0.53; Fig. 2c). Therefore, tumor-only mutational profiles require the use of germline databases with sufficiently broad representation to minimize the elevations in false positives that might otherwise be seen in patients with diverse ethnic backgrounds.
Fig. 2

Ancestry and germline false positives using different analysis strategies. a The use of dbSNP as the primary germline filtration strategy results in a significant increase in false positives among non-white patients (p < 0.001). b A similar increase was observed with the use of 1000 Genomes (p < 0.001). c With larger germline databases such as ExAC, this disparity is mitigated

Impact of gene panel size on mutational load and neoantigen prediction

Given the potential utility for immuno-oncology patient stratification, we next investigated the extent to which various tumor profiling platforms might inform genome-wide properties. First, we explored whether mutational loads measured directly (using paired WES data) correlated with mutational loads inferred from targeted panel data (e.g., <2 % of the genomic territory covered by WES). To test this, we divided the number of mutations observed in the panel by the genomic territory covered (in megabases) by that panel (“Methods”). Consistent with previous reports [19], we found that large-panel mutational loads correlated strongly with WES-based mutational load regardless of whether tumor-only or paired data were used (r2 = 0.99 for matched and 0.93 for unmatched; p < 0.001), with median nonsynonymous mutation rates of 2.3/Mb (IQR = 1.7–4.2) and 5.5/Mb (IQR = 4.1–8.2) in WES and panels, respectively (Fig. 3a, b). The ability of medium gene panels to predict the WES mutational load was somewhat reduced (r2 = 0.84 and 0.71, respectively). Small panels were poor predictors of overall mutational load (0.4 < r2 ≤ 0.6 for all conditions). When analyzed separately, the lung and colon cancer subsets achieved comparable results across the three settings (Additional file 8: Figure S1; Additional file 9: Figure S2). Thus, large matched or unmatched panels successfully recapitulated the WES mutational load.
Fig. 3

Mutational load predictions with different panel tests. Comparison of mutational load predictions using WES or either matched (a) or unmatched (b) large panel tests (n = 300 genes) demonstrates both can reliably predict the mutational load. The linear regression line is shown in black with 95 % confidence bands shaded in grey. The identity line (dashed) is shown for comparison. With medium sized panels (n = 48 genes), this ability decreases in both the matched and unmatched setting and is not possible with small (n = 15) gene panels

Given the potential importance of identifying patient-specific neoantigens—novel protein sequences absent from the normal human genome that arise from somatic mutations [40]—for immuno-oncology applications [1618, 41], we sought to identify neoantigens derived from the different panels and WES. To do this, we integrated patient human leukocyte antigen (HLA) typing [42] with the set of all potential neoantigens to identify those with predicted high affinities to the patient’s MHC class I alleles (“Methods”). The median neoantigen load across all patients as determined by WES data was 38 neoantigens per exome, while the median number of those called by large panel targeted sequencing was one neoantigen per panel (Fig. 4a). The median number of neoantigens called in both small and medium panels was zero. No correlation was discovered between small/medium panel neoantigen calls and exome neoantigen calls (r2 = 0.24 and r2 = 0.62, respectively); however, we did observe a correlation between large panel neoantigen and exome neoantigens (r2 = 0.81) (Fig. 4b–d). In the WES data, 5511 neoantigens with binding affinities of <500 nM were identified across the 157 patient samples. Of these, 229 (4.1 %) were observed through matched targeted panel sequencing data (Fig. 4a). Thus, while large panels were able to recapitulate mutation and predicted neoantigen loads, most potentially immunogenic neoantigens occurred in genes that were not represented in these cancer panels and would not be observable for cancer vaccine strategies.
Fig. 4

Neoantigen predictions in panels. a The proportion of neoantigens called in large panel targeted sequencing data demonstrates an inability to identify as a broad spectrum of neoantigens compared to WES. b Nonetheless, there is a linear relationship between large panel neoantigens recovered from exome and germline-matched large panel data. This linear relationship no longer holds when considering neoantigen data from medium (c) and small (d) targeted panels

Discussion

As precision medicine efforts proliferate in clinical oncology settings, a spectrum of tumor gene profiling strategies—from individual variant testing to comprehensive WES—are being utilized. We sought to assess the relative merits of each by analyzing genomic data from 157 lung and colon adenocarcinoma cases, followed by in silico modeling of different forms of targeted panel testing to determine their analytical strengths and weaknesses.

One priority was to determine best practices for germline filtering in the setting of tumor-only profiling. While our results identified a consistent (albeit lower) rate of germline false positive findings when leveraging large, publically available germline variant data sets, we found that the addition of molecular pathologist review was highly effective in reducing false positive errors germane to unmatched sequencing. This observation may be informative in centers where assembling a molecular tumor board may not be practically feasible but individual molecular pathologists can act as reviewers.

Even so, the higher false positive germline variants in non-white ancestries highlights the limitations of utilizing germline genomic databases for such filtering, as these cohorts may not represent the clinical population being tested. Tumor-only analytical pipelines that do not anticipate diverse ancestry could unwittingly produce a higher rate of germline false positives in some ancestral backgrounds. Expansion of germline databases to represent the diversity of patients tested is necessary to mitigate this source of false positives, and this strategy may contribute to the improvement of precision medicine health disparities resulting from analytical features of the human genome. Indeed, as clinical genomic profiling becomes increasingly expansive technologically, with whole-genome and whole-transcriptome sequencing being performed in clinical settings, the need to capture a diverse set of patients is especially relevant.

Furthermore, because of the emerging therapeutic avenues associated with tumor neoantigens and mutational load, we sought to analyze the ability of targeted panels to identify patients who might benefit from such treatments. Indeed, targeted large panel mutation rate approximated the WES mutation rate well for most samples, whether the targeted panel was matched or unmatched. This suggests that large targeted panels may be useful for flagging patients with exceptionally high mutation rates for specific clinical investigations. Yet, targeted panel sequencing results failed to recapitulate the neoantigen load estimated from WES data, which may be as relevant as mutational load data when combined with immunohistochemistry markers (i.e., PD-L1 staining) to stratify patients for immunotherapies [43].

In addition, since only a small fraction of total neoantigens fell within the genomic regions covered by the targeted panel (average tumor sample had less than three neoantigens called in the targeted sequencing data), it is likely impossible to stratify patients based on relative neoantigen loads from targeted sequencing data alone. Also, the targeted panel does not call most of the patient-specific neoantigens themselves, which may become increasingly relevant as personalized cancer vaccine strategies requiring knowledge of specific neoantigens expand across many clinical settings [4447].

The main limitation of this study is that, due to the rapidly expanding diversity of panel-based sequencing approaches offered in commercial and academic labs, it was not possible to directly and comprehensively compare outputs of all available approaches with these clinical samples. As a result, certain components of the workflow could not be examined in this context, such as the impact of higher sequencing depth on variant detection sensitivity [8] or the differences in germline false positive results from different lab analysis processes. Indeed, since not all labs report specific details about analysis methods, such as how exactly germline variants are filtered, we could not confirm whether the approach outlined in this effort is consistent across vendors. This highlights the importance of encouraging transparency in analytical efforts given how widely variable results may seem depending on which approach is used. Since germline variants may also have immediate clinical implications for assigning cancer risk [48] and therapeutic strategies [49], distinguishing somatic and germline events is especially relevant in this context. Furthermore, this study highlights certain benefits of WES, although WES compared with panel testing has additional costs (i.e., financial, interpretive) beyond analytical.

Conclusions

Broadly, our work highlights the relative advantages and disadvantages of WES and targeted panel sequencing for clinical precision oncology. Targeted panel sequencing maintains an advantage over WES for variant identification in a small set of known clinically informative cancer genes and utilization of germline enhances somatic mutation identification. Additionally, prior studies have demonstrated that targeted panels enable more sequencing depth compared with WES. Since the WES obtained for this study was considered sufficient to enable the subsequent analyses, this technical component of targeted panels did not require further exploration to enable the studies described herein.

Even so, the breadth and adaptability of WES may ultimately offer advantages over targeted panels for certain immunotherapy regimens. As treatment paradigms shift and require detailed assessments of global genomic changes for immunotherapy purposes, as well as deep clonal architecture of tumors feasibly enabled through deeper targeted sequencing, a combination of these strategies may prove most effective for genomic analysis in the clinic. When paired with up-to-date bioinformatics and database filtering, along with molecular pathology assessment, this strategy may inform wider analytical standardization for genomic analysis.

Declarations

Acknowledgements

We thank the patients for participating in this study. We also thank Carrie Sougnez for project management support, the Broad Genomics Platform for sequencing contributions, and the Broad Cancer Genome Analysis group for analysis methodology and feedback. Drs. Van Allen and Garraway had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.

Funding

This work was supported by NIH U01HG006492. This work was also supported by K08 CA188615-02 (E.M.V.) and BroadNext10 (E.M.V.). The funding organizations were not responsible for design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.

Availability of data and materials

All BAM files are deposited in dbGap phs001075.v1.p1.

Authors’ contributions

AG, BR, ATW, AAM, DM, DL, LAG, and EMV performed genomic analysis of whole exome data and analyses throughout this effort. LS, LM, MD, VRR, and NL performed analysis of clinical panel tests. NO, AZ, SG, PJ, JG, SJ, NL, NW, LAG, and EMV coordinated sample acquisition, accrual, and project management of effort. MG, DM, AG, BR, and EMV performed mutational and neoantigen load comparison analysis. AG, LAG, and EMV wrote the manuscript. All authors read and approved the final manuscript.

Competing interests

Drs. Garraway and Wagle are equity holders in and consultants to Foundation Medicine. Dr. Garraway is a consultant to Novartis, Millenium/Takeda, and Boehringer Ingelheim and a recipient of a grant from Novartis. Dr. Van Allen is a consultant to Roche Ventana, Takeda Pharmaceuticals, and Third Rock Ventures. No other competing interests.

Ethics approval and consent to participate

All patients consented to an institutional review board-approved protocol that allows comprehensive genetic analysis of tumor and germline samples (Dana-Farber Cancer Institute #12-078). This study conforms to the Declaration of Helsinki.

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

Authors’ Affiliations

(1)
Department of Medical Oncology, Dana-Farber Cancer Institute, Harvard Medical School
(2)
Broad Institute of MIT and Harvard
(3)
Department of Pathology, Brigham and Women’s Hospital
(4)
Department of Medical Ethics and Health Policy, University of Pennsylvania
(5)
Center for Cancer Precision Medicine, Dana-Farber Cancer Institute

References

  1. Garraway LA. Genomics-driven oncology: framework for an emerging paradigm. J Clin Oncol. 2013. doi:10.1200/JCO.2012.46.8934.Epub2013/04/17.Google Scholar
  2. Frampton GM, Fichtenholtz A, Otto GA, Wang K, Downing SR, He J, et al. Development and validation of a clinical cancer genomic profiling test based on massively parallel DNA sequencing. Nat Biotechnol. 2013;31(11):1023–31. doi:10.1038/nbt.2696.View ArticlePubMedGoogle Scholar
  3. Cryan JB, Haidar S, Ramkissoon LA, Bi WL, Knoff DS, Schultz N, et al. Clinical multiplexed exome sequencing distinguishes adult oligodendroglial neoplasms from astrocytic and mixed lineage gliomas. Oncotarget. 2014;5(18):8083–92. doi:10.18632/oncotarget.2342. PubMed PMID: 25257301; PubMed Central PMCID: PMC4226668.View ArticlePubMedPubMed CentralGoogle Scholar
  4. Pritchard CC, Salipante SJ, Koehler K, Smith C, Scroggins S, Wood B, et al. Validation and implementation of targeted capture and sequencing for the detection of actionable mutation, copy number variation, and gene rearrangement in clinical cancer specimens. J Mol Diagn. 2014;16(1):56–67. doi:10.1016/j.jmoldx.2013.08.004. PubMed PMID: 24189654, PubMed Central PMCID: PMC3873496.View ArticlePubMedPubMed CentralGoogle Scholar
  5. Cheng DT, Mitchell TN, Zehir A, Shah RH, Benayed R, Syed A, et al. Memorial Sloan Kettering-Integrated Mutation Profiling of Actionable Cancer Targets (MSK-IMPACT): a hybridization capture-based next-generation sequencing clinical assay for solid tumor molecular oncology. J Mol Diagn. 2015;17(3):251–64. doi:10.1016/j.jmoldx.2014.12.006.View ArticlePubMedGoogle Scholar
  6. Wagle N, Berger MF, Davis MJ, Blumenstiel B, Defelice M, Pochanard P, et al. High-throughput detection of actionable genomic alterations in clinical tumor samples by targeted, massively parallel sequencing. Cancer Discov. 2012;2(1):82–93. doi:10.1158/2159-8290.CD-11-0184. PubMed PMID: 22585170, PubMed Central PMCID: PMC3353152.View ArticlePubMedGoogle Scholar
  7. Sohal DPS, BI Rini, Khorana AA, Dreicer R, Abraham J, Procop GW, Saunthararajah Y, Pennell NA, Stevenson JP, Pelley R, Estfan B, Shepard D, Funchain P, Elson P, Adelstein DJ, Brian J. Bolwell prospective clinical study of precision oncology in solid tumors. JNCI J Natl Cancer Inst. 2016;108(3):djv332. doi:10.1093/jnci/djv332.
  8. Chen K, Meric-Bernstam F, Zhao H, Zhang Q, Ezzeddine N, Tang LY, et al. Clinical actionability enhanced through deep targeted sequencing of solid tumors. Clin Chem. 2015;61(3):544–53. doi:10.1373/clinchem.2014.231100. PubMed PMID: 25626406, PubMed Central PMCID: PMC4511273.View ArticlePubMedPubMed CentralGoogle Scholar
  9. Fisher KE, Zhang L, Wang J, Smith GH, Newman S, Schneider TM, et al. Clinical validation and implementation of a targeted next-generation sequencing assay to detect somatic variants in non-small cell lung, melanoma, and gastrointestinal malignancies. J Mol Diagn. 2016;18(2):299–315. doi:10.1016/j.jmoldx.2015.11.006.View ArticlePubMedGoogle Scholar
  10. Cottrell CE, Al-Kateb H, Bredemeyer AJ, Duncavage EJ, Spencer DH, Abel HJ, et al. Validation of a next-generation sequencing assay for clinical molecular oncology. J Mol Diagn. 2014;16(1):89–105. doi:10.1016/j.jmoldx.2013.10.002.View ArticlePubMedGoogle Scholar
  11. Roychowdhury S, Iyer MK, Robinson DR, Lonigro RJ, Wu YM, Cao X, et al. Personalized oncology through integrative high-throughput sequencing: a pilot study. Sci Transl Med. 2011;3(111):111ra21. doi:10.1126/scitranslmed.3003161. Epub 2011/12/03. PubMed PMID: 22133722; PubMed Central PMCID: PMC3476478.
  12. Van Allen EM, Wagle N, Stojanov P, Perrin DL, Cibulskis K, Marlow S, et al. Whole-exome sequencing and clinical interpretation of formalin-fixed, paraffin-embedded tumor samples to guide precision cancer medicine. Nat Med. 2014;20(6):682–8. doi:10.1038/nm.3559. PubMed PMID: 24836576, PubMed Central PMCID: PMC4048335.View ArticlePubMedPubMed CentralGoogle Scholar
  13. Beltran H, Eng K, Mosquera JM, Sigaras A, Romanel A, Rennert H, et al. Whole-exome sequencing of metastatic cancer and biomarkers of treatment response. JAMA Oncol. 2015;1(4):466–74. doi:10.1001/jamaoncol.2015.1313. PubMed PMID: 26181256, PubMed Central PMCID: PMC4505739.View ArticlePubMedPubMed CentralGoogle Scholar
  14. Schumacher TN, Kesmir C, van Buuren MM. Biomarkers in cancer immunotherapy. Cancer Cell. 2015;27(1):12–4. doi:10.1016/j.ccell.2014.12.004.View ArticlePubMedGoogle Scholar
  15. Le DT, Uram JN, Wang H, Bartlett BR, Kemberling H, Eyring AD, et al. PD-1 Blockade in tumors with mismatch-repair deficiency. N Engl J Med. 2015. doi:10.1056/NEJMoa1500596.Google Scholar
  16. Snyder A, Makarov V, Merghoub T, Yuan J, Zaretsky JM, Desrichard A, et al. Genetic basis for clinical response to CTLA-4 blockade in melanoma. N Engl J Med. 2014;371(23):2189–99. doi:10.1056/NEJMoa1406498.View ArticlePubMedPubMed CentralGoogle Scholar
  17. Rizvi NA, Hellmann MD, Snyder A, Kvistborg P, Makarov V, Havel JJ, et al. Cancer immunology. Mutational landscape determines sensitivity to PD-1 blockade in non-small cell lung cancer. Science. 2015;348(6230):124–8. doi:10.1126/science.aaa1348.View ArticlePubMedGoogle Scholar
  18. Van Allen EM, Miao D, Schilling B, Shukla SA, Blank C, Zimmer L, et al. Genomic correlates of response to CTLA-4 blockade in metastatic melanoma. Science. 2015;350(6257):207–11. doi:10.1126/science.aad0095.View ArticlePubMedGoogle Scholar
  19. Campesato LF, Barroso-Sousa R, Jimenez L, Correa BR, Sabbaga J, Hoff PM, et al. Comprehensive cancer-gene panels can be used to estimate mutational load and predict clinical benefit to PD-1 blockade in clinical practice. Oncotarget. 2015;6(33):34221–7. doi:10.18632/oncotarget.5950.Google Scholar
  20. Stadler ZK, Battaglin F, Middha S, Hechtman JF, Tran C, Cercek A, et al. Reliable detection of mismatch repair deficiency in colorectal cancers using mutational load in next-generation sequencing panels. J Clin Oncol. 2016. doi:10.1200/JCO.2015.65.1067.Google Scholar
  21. Meric-Bernstam F, Brusco L, Shaw K, Horombe C, Kopetz S, Davies MA, et al. Feasibility of large-scale genomic testing to facilitate enrollment onto genomically matched clinical trials. J Clin Oncol. 2015;33(25):2753–62. doi:10.1200/JCO.2014.60.4165. PubMed PMID: 26014291, PubMed Central PMCID: PMC4550690.View ArticlePubMedGoogle Scholar
  22. MacConaill LE, Garcia E, Shivdasani P, Ducar M, Adusumilli R, Breneiser M, et al. Prospective enterprise-level molecular genotyping of a cohort of cancer patients. J Mol Diagn. 2014;16(6):660–72. doi:10.1016/j.jmoldx.2014.06.004. PubMed PMID: 25157968, PubMed Central PMCID: PMC4210463.View ArticlePubMedPubMed CentralGoogle Scholar
  23. Jones S, Anagnostou V, Lytle K, Parpart-Li S, Nesselbush M, Riley DR, et al. Personalized genomic analyses for cancer mutation discovery and interpretation. Sci Transl Med. 2015;7(283):283ra53. doi:10.1126/scitranslmed.aaa7161. PubMed PMID: 25877891; PubMed Central PMCID: PMC4442685.
  24. Schrader KA, Cheng DT, Joseph V, Prasad M, Walsh M, Zehir A, et al. Germline variants in targeted tumor sequencing using matched normal DNA. JAMA Oncol. 2016;2(1):104–11. doi:10.1001/jamaoncol.2015.5208.View ArticlePubMedGoogle Scholar
  25. Reich M, Liefeld T, Gould J, Lerner J, Tamayo P, Mesirov JP. GenePattern 2.0. Nat Genet. 2006;38(5):500–1. doi:10.1038/ng0506-500. Epub 2006/04/28.View ArticlePubMedGoogle Scholar
  26. Cibulskis K, McKenna A, Fennell T, Banks E, DePristo M, Getz G. ContEst: estimating cross-contamination of human samples in next-generation sequencing data. Bioinformatics. 2011;27(18):2601–2. doi:10.1093/bioinformatics/btr446. PubMed PMID: 21803805, PubMed Central PMCID: PMC3167057, Epub 2011/08/02.PubMedPubMed CentralGoogle Scholar
  27. Cibulskis K, Lawrence MS, Carter SL, Sivachenko A, Jaffe D, Sougnez C, et al. Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples. Nat Biotechnol. 2013. doi:10.1038/nbt.2514.Epub2013/02/12.PubMedPubMed CentralGoogle Scholar
  28. Costello M, Pugh TJ, Fennell TJ, Stewart C, Lichtenstein L, Meldrim JC, et al. Discovery and characterization of artifactual mutations in deep coverage targeted capture sequencing data due to oxidative DNA damage during sample preparation. Nucleic Acids Res. 2013;41(6):e67. doi:10.1093/nar/gks1443. PubMed PMID: 23303777, PubMed Central PMCID: PMC3616734, Epub 2013/01/11.View ArticlePubMedPubMed CentralGoogle Scholar
  29. Ramos AH, Lichtenstein L, Gupta M, Lawrence MS, Pugh TJ, Saksena G, et al. Oncotator: cancer variant annotation tool. Hum Mutat. 2015;36(4):E2423–9. doi:10.1002/humu.22771.View ArticlePubMedGoogle Scholar
  30. Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, et al. dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 2001;29(1):308–11. PubMed PMID: 11125122, PubMed Central PMCID: PMC29783.View ArticlePubMedPubMed CentralGoogle Scholar
  31. Genomes Project C, Abecasis GR, Auton A, Brooks LD, DePristo MA, Durbin RM, et al. An integrated map of genetic variation from 1,092 human genomes. Nature. 2012;491(7422):56–65. doi:10.1038/nature11632. PubMed PMID: 23128226, PubMed Central PMCID: PMC3498066.View ArticleGoogle Scholar
  32. ExAC. Large-scale aggregation of human genomic data. 2016. http://exac.broadinstitute.org.
  33. Forbes SA, Bindal N, Bamford S, Cole C, Kok CY, Beare D, et al. COSMIC: mining complete cancer genomes in the Catalogue of Somatic Mutations in Cancer. Nucleic Acids Res. 2011;39(Database issue):D945–50. doi:10.1093/nar/gkq929. PubMed PMID: 20952405, PubMed Central PMCID: PMC3013785.View ArticlePubMedGoogle Scholar
  34. DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011;43(5):491–8. doi:10.1038/ng.806. PubMed PMID: 21478889, PubMed Central PMCID: PMC3083463, Epub 2011/04/12.View ArticlePubMedPubMed CentralGoogle Scholar
  35. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303. doi:10.1101/gr.107524.110. PubMed PMID: 20644199, PubMed Central PMCID: PMC2928508, Epub 2010/07/21.View ArticlePubMedPubMed CentralGoogle Scholar
  36. Shukla SA, Rooney MS, Rajasagi M, Tiao G, Dixon PM, Lawrence MS, et al. Comprehensive analysis of cancer-associated somatic mutations in class I HLA genes. Nat Biotechnol. 2015;33(11):1152–8. doi:10.1038/nbt.3344. PubMed PMID: 26372948, PubMed Central PMCID: PMC4747795.View ArticlePubMedPubMed CentralGoogle Scholar
  37. Hoof I, Peters B, Sidney J, Pedersen LE, Sette A, Lund O, et al. NetMHCpan, a method for MHC class I binding prediction beyond humans. Immunogenetics. 2009;61(1):1–13. doi:10.1007/s00251-008-0341-z. PubMed PMID: 19002680, PubMed Central PMCID: PMC3319061.View ArticlePubMedGoogle Scholar
  38. Forbes SA, Beare D, Gunasekaran P, Leung K, Bindal N, Boutselakis H, et al. COSMIC: exploring the world’s knowledge of somatic mutations in human cancer. Nucleic Acids Res. 2015;43(Database issue):D805–11. doi:10.1093/nar/gku1075. PubMed PMID: 25355519, PubMed Central PMCID: PMC4383913.View ArticlePubMedGoogle Scholar
  39. Bustamante CD, Burchard EG, De la Vega FM. Genomics for the world. Nature. 2011;475(7355):163–5. doi:10.1038/475163a. PubMed PMID: 21753830, PubMed Central PMCID: PMC3708540.View ArticlePubMedPubMed CentralGoogle Scholar
  40. Schumacher TN, Schreiber RD. Neoantigens in cancer immunotherapy. Science. 2015;348(6230):69–74. doi:10.1126/science.aaa4971.View ArticlePubMedGoogle Scholar
  41. Carreno BM, Magrini V, Becker-Hapak M, Kaabinejadian S, Hundal J, Petti AA, et al. Cancer immunotherapy. A dendritic cell vaccine increases the breadth and diversity of melanoma neoantigen-specific T cells. Science. 2015;348(6236):803–8. doi:10.1126/science.aaa3828.View ArticlePubMedPubMed CentralGoogle Scholar
  42. Shukla SA, Rooney MS, Rajasagi M, Tiao G, Dixon PM, Lawrence M, et al. Sensitive detection of somatic mutations in class I HLA genes reveals enrichment for functional mutations in cancer. Nat Biotechnol. 2015; in press.Google Scholar
  43. Tumeh PC, Harview CL, Yearley JH, Shintaku IP, Taylor EJ, Robert L, et al. PD-1 blockade induces responses by inhibiting adaptive immune resistance. Nature. 2014;515(7528):568–71. doi:10.1038/nature13954. PubMed PMID: 25428505, PubMed Central PMCID: PMC4246418.View ArticlePubMedPubMed CentralGoogle Scholar
  44. Rajasagi M, Shukla SA, Fritsch EF, Keskin DB, DeLuca D, Carmona E, et al. Systematic identification of personal tumor-specific neoantigens in chronic lymphocytic leukemia. Blood. 2014;124(3):453–62. doi:10.1182/blood-2014-04-567933. PubMed PMID: 24891321, PubMed Central PMCID: PMC4102716.View ArticlePubMedPubMed CentralGoogle Scholar
  45. Gubin MM, Artyomov MN, Mardis ER, Schreiber RD. Tumor neoantigens: building a framework for personalized cancer immunotherapy. J Clin Invest. 2015;125(9):3413–21. doi:10.1172/JCI80008. PubMed PMID: 26258412, PubMed Central PMCID: PMC4588307.View ArticlePubMedGoogle Scholar
  46. Yoshimura K, Minami T, Nozawa M, Kimura T, Egawa S, Fujimoto H, et al. A phase 2 randomized controlled trial of personalized peptide vaccine immunotherapy with low-dose dexamethasone versus dexamethasone alone in chemotherapy-naive castration-resistant prostate cancer. Eur Urol. 2016. doi:10.1016/j.eururo.2015.12.050.PubMedGoogle Scholar
  47. Noguchi M, Matsumoto K, Uemura H, Arai G, Eto M, Naito S, et al. An open-label, randomized phase II trial of personalized peptide vaccination in patients with bladder cancer that progressed after platinum-based chemotherapy. Clin Cancer Res. 2016;22(1):54–60. doi:10.1158/1078-0432.CCR-15-1265.View ArticlePubMedGoogle Scholar
  48. Lu C, Xie M, Wendl MC, Wang J, McLellan MD, Leiserson MD, et al. Patterns and functional implications of rare germline variants across 12 cancer types. Nat Commun. 2015;6:10086. doi:10.1038/ncomms10086. PubMed PMID: 26689913, PubMed Central PMCID: PMC4703835.View ArticlePubMedPubMed CentralGoogle Scholar
  49. Mateo J, Carreira S, Sandhu S, Miranda S, Mossop H, Perez-Lopez R, et al. DNA-repair defects and olaparib in metastatic prostate cancer. N Engl J Med. 2015;373(18):1697–708. doi:10.1056/NEJMoa1506859.View ArticlePubMedGoogle Scholar

Copyright

© The Author(s). 2016

Advertisement