A Novel Next-Generation Sequencing Approach to Detecting Microsatellite Instability and Pan-Tumor Characterization of 1000 Microsatellite Instability-High Cases in 67,000 Patient Samples.

Microsatellite instability (MSI) is an important biomarker for predicting response to immune checkpoint inhibitor therapy, as emphasized by the recent checkpoint inhibitor approval for MSI-high (MSI-H) solid tumors. Herein, we describe and validate a novel method for determining MSI status from a next-generation sequencing comprehensive genomic profiling assay using formalin-fixed, paraffin-embedded samples. This method is 97% (65/67) concordant with current standards, PCR and immunohistochemistry. We further apply this method to >67,000 patient tumor samples to identify genes and pathways that are enriched in MSI-stable or MSI-H tumor groups. Data show that although rare in tumors other than colorectal and endometrial carcinomas, MSI-H samples are present in many tumor types. Furthermore, the large sample set revealed that MSI-H tumors selectively share alterations in genes across multiple common pathways, including WNT, phosphatidylinositol 3-kinase, and NOTCH. Last, MSI is sufficient, but not necessary, for a tumor to have elevated tumor mutation burden. Therefore, MSI can be determined from comprehensive genomic profiling with high accuracy, allowing for efficient MSI-H detection across all tumor types, especially those in which routine use of immunohistochemistry or PCR-based assays would be impractical because of a rare incidence of MSI. MSI-H tumors are enriched in alterations in specific signaling pathways, providing a rationale for investigating directed immune checkpoint inhibitor therapies in combination with pathway-targeted therapies.

BRCA Reversion Mutations in Circulating Tumor DNA Predict Primary and Acquired Resistance to the PARP Inhibitor Rucaparib in High-Grade Ovarian Carcinoma.

A key resistance mechanism to platinum-based chemotherapies and PARP inhibitors in BRCA-mutant cancers is the acquisition of BRCA reversion mutations that restore protein function. To estimate the prevalence of BRCA reversion mutations in high-grade ovarian carcinoma (HGOC), we performed targeted next-generation sequencing of circulating cell-free DNA (cfDNA) extracted from pretreatment and postprogression plasma in patients with deleterious germline or somatic BRCA mutations treated with the PARP inhibitor rucaparib. BRCA reversion mutations were identified in pretreatment cfDNA from 18% (2/11) of platinum-refractory and 13% (5/38) of platinum-resistant cancers, compared with 2% (1/48) of platinum-sensitive cancers (P = 0.049). Patients without BRCA reversion mutations detected in pretreatment cfDNA had significantly longer rucaparib progression-free survival than those with reversion mutations (median, 9.0 vs. 1.8 months; HR, 0.12; P < 0.0001). To study acquired resistance, we sequenced 78 postprogression cfDNA, identifying eight additional patients with BRCA reversion mutations not found in pretreatment cfDNA. SIGNIFICANCE: BRCA reversion mutations are detected in cfDNA from platinum-resistant or platinum-refractory HGOC and are associated with decreased clinical benefit from rucaparib treatment. Sequencing of cfDNA can detect multiple BRCA reversion mutations, highlighting the ability to capture multiclonal heterogeneity.This article is highlighted in the In This Issue feature, p. 151.

A computational approach to distinguish somatic vs. germline origin of genomic alterations from deep sequencing of cancer specimens without a matched normal.

A key constraint in genomic testing in oncology is that matched normal specimens are not commonly obtained in clinical practice. Thus, while well-characterized genomic alterations do not require normal tissue for interpretation, a significant number of alterations will be unknown in whether they are germline or somatic, in the absence of a matched normal control. We introduce SGZ (somatic-germline-zygosity), a computational method for predicting somatic vs. germline origin and homozygous vs. heterozygous or sub-clonal state of variants identified from deep massively parallel sequencing (MPS) of cancer specimens. The method does not require a patient matched normal control, enabling broad application in clinical research. SGZ predicts the somatic vs. germline status of each alteration identified by modeling the alteration's allele frequency (AF), taking into account the tumor content, tumor ploidy, and the local copy number. Accuracy of the prediction depends on the depth of sequencing and copy number model fit, which are achieved in our clinical assay by sequencing to high depth (>500x) using MPS, covering 394 cancer-related genes and over 3,500 genome-wide single nucleotide polymorphisms (SNPs). Calls are made using a statistic based on read depth and local variability of SNP AF. To validate the method, we first evaluated performance on samples from 30 lung and colon cancer patients, where we sequenced tumors and matched normal tissue. We examined predictions for 17 somatic hotspot mutations and 20 common germline SNPs in 20,182 clinical cancer specimens. To assess the impact of stromal admixture, we examined three cell lines, which were titrated with their matched normal to six levels (10-75%). Overall, predictions were made in 85% of cases, with 95-99% of variants predicted correctly, a significantly superior performance compared to a basic approach based on AF alone. We then applied the SGZ method to the COSMIC database of known somatic variants in cancer and found >50 that are in fact more likely to be germline.

First-in-human trial of multikinase VEGF inhibitor regorafenib and anti-EGFR antibody cetuximab in advanced cancer patients.

BACKGROUND: The combination of multikinase VEGF inhibitor regorafenib and anti-EGFR antibody cetuximab overcomes intrinsic and acquired resistance in both EGFR-sensitive and EGFR-resistant preclinical models of colorectal cancer (CRC). METHODS: Utilizing a standard 3+3 design, a phase I study was designed to determine safety, maximum tolerated dose (MTD), and dose-limiting toxicities (DLTs) of the regorafenib plus cetuximab combination among patients with advanced cancer including CRC. Comprehensive genomic profiling was performed on the exceptional responder. RESULTS: Among the 27 patients enrolled the median age was 54 years. None of 19 patients treated at dose level 1 (cetuximab i.v. 200 mg/m2 followed by 150 mg/m2 weekly + regorafenib 80 mg daily) experienced a DLT, and 2 of 5 patients treated at dose level 2 (cetuximab i.v. 200 mg/m2 followed by 150 mg/m2 weekly + regorafenib 120 mg daily) experienced a DLT (grade 3 thrombocytopenia [n = 1] and grade 3 intra-abdominal bleed [n = 1]). Most common adverse events were grade 1 or 2 rash (20 patients). Of 24 evaluable patients, 11 (46%) patients had clinical benefit (stable disease > 6 cycles or partial response [PR]) (CRC n = 8, one patient each with head and neck cancer, carcinoma of unknown primary, and glioblastoma). A CRC patient, who progressed on anti-EGFR and regorafenib, achieved a PR (46% decrease per RECIST v1.1) lasting 15 months. Genomic profiling of an exceptional responder with response for over 27 cycles revealed hypermutated genotype with microsatellite instability (MSI). CONCLUSION: Regorafenib 80 mg daily plus cetuximab 200 mg/m2 loading dose, followed by 150 mg/m2 every week is the MTD/recommended phase II dose. The combination demonstrated early signals of activity in wild-type CRC, including 1 exceptional responder with MSI high. TRIAL REGISTRATION: clinicaltrials.gov NCT02095054FUNDING. The University of Texas MD Anderson Cancer Center is supported by the NIH Cancer Center Support Grant CA016672. This work was supported in part by the Cancer Prevention Research Institute of Texas grant RP110584 and National Center for Advancing Translational Sciences grant UL1 TR000371 (Center for Clinical and Translational Sciences).

Genomic Profiling of a Large Set of Diverse Pediatric Cancers Identifies Known and Novel Mutations across Tumor Spectra.

Pediatric cancers are generally characterized by low mutational burden and few recurrently mutated genes. Recent studies suggest that genomic alterations may help guide treatment decisions and clinical trial selection. Here, we describe genomic profiles from 1,215 pediatric tumors representing sarcomas, extracranial embryonal tumors, brain tumors, hematologic malignancies, carcinomas, and gonadal tumors. Comparable published datasets identified similar frequencies of clinically relevant alterations, validating this dataset as biologically relevant. We identified novel ALK fusions in a neuroblastoma (BEND5-ALK) and an astrocytoma (PPP1CB-ALK), novel BRAF fusions in an astrocytoma (BCAS1-BRAF) and a ganglioglioma (TMEM106B-BRAF), and a novel PAX3-GLI2 fusion in a rhabdomyosarcoma. Previously characterized ALK, NTRK1, and PAX3 fusions were observed in unexpected malignancies, challenging the "disease-specific" alterations paradigm. Finally, we identified recurrent variants of unknown significance in MLL3 and PRSS1 predicted to have functional impact. Data from these 1,215 tumors are publicly available for discovery and validation. Cancer Res; 77(2); 509-19. ©2017 AACR.

Targeted Next Generation Sequencing Identifies Markers of Response to PD-1 Blockade.

Therapeutic antibodies blocking programmed death-1 and its ligand (PD-1/PD-L1) induce durable responses in a substantial fraction of melanoma patients. We sought to determine whether the number and/or type of mutations identified using a next-generation sequencing (NGS) panel available in the clinic was correlated with response to anti-PD-1 in melanoma. Using archival melanoma samples from anti-PD-1/PD-L1-treated patients, we performed hybrid capture-based NGS on 236-315 genes and T-cell receptor (TCR) sequencing on initial and validation cohorts from two centers. Patients who responded to anti-PD-1/PD-L1 had higher mutational loads in an initial cohort (median, 45.6 vs. 3.9 mutations/MB; P = 0.003) and a validation cohort (37.1 vs. 12.8 mutations/MB; P = 0.002) compared with nonresponders. Response rate, progression-free survival, and overall survival were superior in the high, compared with intermediate and low, mutation load groups. Melanomas with NF1 mutations harbored high mutational loads (median, 62.7 mutations/MB) and high response rates (74%), whereas BRAF/NRAS/NF1 wild-type melanomas had a lower mutational load. In these archival samples, TCR clonality did not predict response. Mutation numbers in the 315 genes in the NGS platform strongly correlated with those detected by whole-exome sequencing in The Cancer Genome Atlas samples, but was not associated with survival. In conclusion, mutational load, as determined by an NGS platform available in the clinic, effectively stratified patients by likelihood of response. This approach may provide a clinically feasible predictor of response to anti-PD-1/PD-L1. Cancer Immunol Res; 4(11); 959-67. ©2016 AACR.

Broad Detection of Alterations Predicted to Confer Lack of Benefit From EGFR Antibodies or Sensitivity to Targeted Therapy in Advanced Colorectal Cancer.

INTRODUCTION: A KRAS mutation represented the first genomic biomarker to predict lack of benefit from anti-epidermal growth factor receptor (EGFR) antibody therapy in advanced colorectal cancer (CRC). Expanded RAS testing has further refined the treatment approach, but understanding of genomic alterations underlying primary and acquired resistance is limited and further study is needed. MATERIALS AND METHODS: We prospectively analyzed 4,422 clinical samples from patients with advanced CRC, using hybrid-capture based comprehensive genomic profiling (CGP) at the request of the individual treating physicians. Comparison with prior molecular testing results, when available, was performed to assess concordance. RESULTS: We identified a RAS/RAF pathway mutation or amplification in 62% of cases, including samples harboring KRAS mutations outside of the codon 12/13 hotspot region in 6.4% of cases. Among cases with KRAS non-codon 12/13 alterations for which prior test results were available, 79 of 90 (88%) were not identified by focused testing. Of 1,644 RAS/RAF wild-type cases analyzed by CGP, 31% harbored a genomic alteration (GA) associated with resistance to anti-EGFR therapy in advanced CRC including mutations in PIK3CA, PTEN, EGFR, and ERBB2. We also identified other targetable GA, including novel kinase fusions, receptor tyrosine kinase amplification, activating point mutations, as well as microsatellite instability. CONCLUSION: Extended genomic profiling reliably detects alterations associated with lack of benefit to anti-EGFR therapy in advanced CRC, while simultaneously identifying alterations potentially important in guiding treatment. The use of CGP during the course of clinical care allows for the refined selection of appropriate targeted therapies and clinical trials, increasing the chance of clinical benefit and avoiding therapeutic futility. IMPLICATIONS FOR PRACTICE: Comprehensive genomic profiling (CGP) detects diverse genomic alterations associated with lack of benefit to anti-epidermal growth factor receptor therapy in advanced colorectal cancer (CRC), as well as targetable alterations in many other genes. This includes detection of a broad spectrum of activating KRAS alterations frequently missed by focused molecular hotspot testing, as well as other RAS/RAF pathway alterations, mutations shown to disrupt antibody binding, RTK activating point mutations, amplifications, and rearrangements, and activating alterations in downstream effectors including PI3K and MEK1. The use of CGP in clinical practice is critical to guide appropriate selection of targeted therapies for patients with advanced CRC.

Clinical Actionability of Comprehensive Genomic Profiling for Management of Rare or Refractory Cancers.

BACKGROUND: The frequency with which targeted tumor sequencing results will lead to implemented change in care is unclear. Prospective assessment of the feasibility and limitations of using genomic sequencing is critically important. METHODS: A prospective clinical study was conducted on 100 patients with diverse-histology, rare, or poor-prognosis cancers to evaluate the clinical actionability of a Clinical Laboratory Improvement Amendments (CLIA)-certified, comprehensive genomic profiling assay (FoundationOne), using formalin-fixed, paraffin-embedded tumors. The primary objectives were to assess utility, feasibility, and limitations of genomic sequencing for genomically guided therapy or other clinical purpose in the setting of a multidisciplinary molecular tumor board. RESULTS: Of the tumors from the 92 patients with sufficient tissue, 88 (96%) had at least one genomic alteration (average 3.6, range 0-10). Commonly altered pathways included p53 (46%), RAS/RAF/MAPK (rat sarcoma; rapidly accelerated fibrosarcoma; mitogen-activated protein kinase) (45%), receptor tyrosine kinases/ligand (44%), PI3K/AKT/mTOR (phosphatidylinositol-4,5-bisphosphate 3-kinase; protein kinase B; mammalian target of rapamycin) (35%), transcription factors/regulators (31%), and cell cycle regulators (30%). Many low frequency but potentially actionable alterations were identified in diverse histologies. Use of comprehensive profiling led to implementable clinical action in 35% of tumors with genomic alterations, including genomically guided therapy, diagnostic modification, and trigger for germline genetic testing. CONCLUSION: Use of targeted next-generation sequencing in the setting of an institutional molecular tumor board led to implementable clinical action in more than one third of patients with rare and poor-prognosis cancers. Major barriers to implementation of genomically guided therapy were clinical status of the patient and drug access. Early and serial sequencing in the clinical course and expanded access to genomically guided early-phase clinical trials and targeted agents may increase actionability. IMPLICATIONS FOR PRACTICE: Identification of key factors that facilitate use of genomic tumor testing results and implementation of genomically guided therapy may lead to enhanced benefit for patients with rare or difficult to treat cancers. Clinical use of a targeted next-generation sequencing assay in the setting of an institutional molecular tumor board led to implementable clinical action in over one third of patients with rare and poor prognosis cancers. The major barriers to implementation of genomically guided therapy were clinical status of the patient and drug access both on trial and off label. Approaches to increase actionability include early and serial sequencing in the clinical course and expanded access to genomically guided early phase clinical trials and targeted agents.

Comprehensive Genomic Profiling of Clinically Advanced Medullary Thyroid Carcinoma.

OBJECTIVE: The aim of this study was to determine the genomic alterations of cancer-related genes in advanced medullary thyroid carcinoma during the course of clinical care. METHODS: Hybrid-capture-based comprehensive genomic profiling was performed on 34 consecutive medullary thyroid carcinoma cases to identify all four classes of genomic alterations, and outcome for an index patient was collected. RESULTS: RET was mutated in 88% (30/34) of cases, with RET M918T being responsible for 70% (21/30) of the RET alterations. The other RET alterations were RET E632_L633del, C634R, C620R, C618G/R/S, V804M, and RET amplification. Two of the four RET wild-type patients harbored mutations in KRAS or HRAS (1/34 each). The next most frequent genomic alterations were amplifications of CCND1, FGF3, and FGF19 and alterations in CDKN2A (3/34 each). One case with a RET M918T mutation developed acquired resistance to progressively dose-escalated vandetanib. When the mTOR inhibitor everolimus was added to continued vandetanib treatment, the patient achieved a second 25% reduction of tumor volume (RECIST 1.1) for 8 months. CONCLUSIONS: Comprehensive genomic profiling identified the full breadth of RET alterations in metastatic medullary thyroid carcinoma and possible cooperating oncogenic driver alterations. This approach may refine the use of targeted therapy for these patients.

Loss of heterozygosity at the CYP2D6 locus in breast cancer: implications for germline pharmacogenetic studies.

BACKGROUND: Controversy exists regarding the impact of CYP2D6 genotype on tamoxifen responsiveness. We examined loss of heterozygosity (LOH) at the CYP2D6 locus and determined its impact on genotyping error when tumor tissue is used as a DNA source. METHODS: Genomic tumor data from the adjuvant and metastatic settings (The Cancer Genome Atlas [TCGA] and Foundation Medicine [FM]) were analyzed to characterize the impact of CYP2D6 copy number alterations (CNAs) and LOH on Hardy Weinberg equilibrium (HWE). Additionally, we analyzed CYP2D6 *4 genotype from formalin-fixed paraffin-embedded (FFPE) tumor blocks containing nonmalignant tissue and buccal (germline) samples from patients on the North Central Cancer Treatment Group (NCCTG) 89-30-52 tamoxifen trial. All statistical tests were two-sided. RESULTS: In TCGA samples (n =627), the CYP2D6 LOH rate was similar in estrogen receptor (ER)-positive (41.2%) and ER-negative (35.2%) but lower in HER2-positive tumors (15.1%) (P < .001). In FM ER+ samples (n = 290), similar LOH rates were observed (40.8%). In 190 NCCTG samples, the agreement between CYP2D6 genotypes derived from FFPE tumors and FFPE tumors containing nonmalignant tissue was moderate (weighted Kappa = 0.74; 95% CI = 0.63 to 0.84). Comparing CYP2D6 genotypes derived from buccal cells to FFPE tumor DNA, CYP2D6*4 genotype was discordant in six of 31(19.4%). In contrast, there was no disagreement between CYP2D6 genotypes derived from buccal cells with FFPE tumors containing nonmalignant tissue. CONCLUSIONS: LOH at the CYP2D6 locus is common in breast cancer, resulting in potential misclassification of germline CYP2D6 genotypes. Tumor DNA should not be used to determine germline CYP2D6 genotype without sensitive techniques to detect low frequency alleles and quality control procedures appropriate for somatic DNA.

A direct characterization of human mutation based on microsatellites.

Mutations are the raw material of evolution but have been difficult to study directly. We report the largest study of new mutations to date, comprising 2,058 germline changes discovered by analyzing 85,289 Icelanders at 2,477 microsatellites. The paternal-to-maternal mutation rate ratio is 3.3, and the rate in fathers doubles from age 20 to 58, whereas there is no association with age in mothers. Longer microsatellite alleles are more mutagenic and tend to decrease in length, whereas the opposite is seen for shorter alleles. We use these empirical observations to build a model that we apply to individuals for whom we have both genome sequence and microsatellite data, allowing us to estimate key parameters of evolution without calibration to the fossil record. We infer that the sequence mutation rate is 1.4-2.3×10(-8) mutations per base pair per generation (90% credible interval) and that human-chimpanzee speciation occurred 3.7-6.6 million years ago.

Microsatellites are molecular clocks that support accurate inferences about history.

Microsatellite length mutations are often modeled using the generalized stepwise mutation process, which is a type of random walk. If this model is sufficiently accurate, one can estimate the coalescence time between alleles of a locus after a mathematical transformation of the allele lengths. When large-scale microsatellite genotyping first became possible, there was substantial interest in using this approach to make inferences about time and demography, but that interest has waned because it has not been possible to empirically validate the clock by comparing it with data in which the mutation process is well understood. We analyzed data from 783 microsatellite loci in human populations and 292 loci in chimpanzee populations, and compared them with up to one gigabase of aligned sequence data, where the molecular clock based upon nucleotide substitutions is believed to be reliable. We empirically demonstrate a remarkable linearity (r(2) > 0.95) between the microsatellite average square distance statistic and sequence divergence. We demonstrate that microsatellites are accurate molecular clocks for coalescent times of at least 2 million years (My). We apply this insight to confirm that the African populations San, Biaka Pygmy, and Mbuti Pygmy have the deepest coalescent times among populations in the Human Genome Diversity Project. Furthermore, we show that microsatellites support unbiased estimates of population differentiation (F(ST)) that are less subject to ascertainment bias than single nucleotide polymorphism (SNP) F(ST). These results raise the prospect of using microsatellite data sets to determine parameters of population history. When genotyped along with SNPs, microsatellite data can also be used to correct for SNP ascertainment bias.

The cardiac output from blood pressure algorithms trial.

OBJECTIVE: The value of different algorithms that estimate cardiac output (CO) by analysis of a peripheral arterial blood pressure (ABP) waveform has not been definitively identified. In this investigation, we developed a testing data set containing a large number of radial ABP waveform segments and contemporaneous reference CO by thermodilution measurements, collected in an intensive care unit (ICU) patient population during routine clinical operations. We employed this data set to evaluate a set of investigational algorithms, and to establish a public resource for the meaningful comparison of alternative CO-from-ABP algorithms. DESIGN: A retrospective comparative analysis of eight investigational CO-from-ABP algorithms using the Multiparameter Intelligent Monitoring in Intensive Care II database. SETTING: Mixed medical/surgical ICU of a university hospital. PATIENTS: A total of 120 cases. INTERVENTIONS: None. MEASUREMENTS: CO estimated by eight investigational CO-from-ABP algorithms, and CO(TD) as a reference. MAIN RESULTS: All investigational methods were significantly better than mean arterial pressure (MAP) at estimating direction changes in CO(TD). Only the formula proposed by Liljestrand and Zander in 1928 was a significantly better quantitative estimator of CO(TD) compared with MAP (95% limits-of-agreement with CO(TD): -1.76/+1.41 L/min versus -2.20/+1.82 L/min, respectively; p < 0.001, per the Kolmogorov-Smirnov test). The Liljestrand method was even more accurate when applied to the cleanest ABP waveforms. Other investigational algorithms were not significantly superior to MAP as quantitative estimators of CO. CONCLUSIONS: Based on ABP data recorded during routine intensive care unit (ICU) operations, the Liljestrand and Zander method is a better estimator of CO(TD) than MAP alone. Our attempts to fully replicate commercially-available methods were unsuccessful, and these methods could not be evaluated. However, the data set is publicly and freely available, and developers and vendors of CO-from-ABP algorithms are invited to test their methods using these data.