Molecular response in newly diagnosed chronic-phase chronic myeloid leukemia: prediction modeling and pathway analysis.

Tyrosine kinase inhibitor therapy revolutionized chronic myeloid leukemia treatment and showed how targeted therapy and molecular monitoring could be used to substantially improve survival outcomes. We used chronic myeloid leukemia as a model to understand a critical question: why do some patients have an excellent response to therapy, while others have a poor response? We studied gene expression in whole blood samples from 112 patients from a large phase III randomized trial (clinicaltrials gov. Identifier: NCT00471497), dichotomizing cases into good responders (BCR::ABL1 ≤10% on the International Scale by 3 and 6 months and ≤0.1% by 12 months) and poor responders (failure to meet these criteria). Predictive models based on gene expression demonstrated the best performance (area under the curve =0.76, standard deviation =0.07). All of the top 20 pathways overexpressed in good responders involved immune regulation, a finding validated in an independent data set. This study emphasizes the importance of pretreatment adaptive immune response in treatment efficacy and suggests biological pathways that can be targeted to improve response.

Contribution and clinical relevance of germline variation to the cancer transcriptome.

BACKGROUND: Somatic alterations in the cancer genome, some of which are associated with changes in gene expression, have been characterized in multiple studies across diverse cancer types. However, less is known about germline variants that influence tumor biology by shaping the cancer transcriptome. METHODS: We performed expression quantitative trait loci (eQTL) analyses using multi-dimensional data from The Cancer Genome Atlas to explore the role of germline variation in mediating the cancer transcriptome. After accounting for associations between somatic alterations and gene expression, we determined the contribution of inherited variants to the cancer transcriptome relative to that of somatic variants. Finally, we performed an interaction analysis using estimates of tumor cellularity to identify cell type-restricted eQTLs. RESULTS: The proportion of genes with at least one eQTL varied between cancer types, ranging between 0.8% in melanoma to 28.5% in thyroid cancer and was correlated more strongly with intratumor heterogeneity than with somatic alteration rates. Although contributions to variance in gene expression was low for most genes, some eQTLs accounted for more than 30% of expression of proximal genes. We identified cell type-restricted eQTLs in genes known to be cancer drivers including LPP and EZH2 that were associated with disease-specific mortality in TCGA but not associated with disease risk in published GWAS. Together, our results highlight the need to consider germline variation in interpreting cancer biology beyond risk prediction.

Adjuvant dabrafenib plus trametinib versus placebo in patients with resected, BRAFV600-mutant, stage III melanoma (COMBI-AD): exploratory biomarker analyses from a randomised, phase 3 trial.

BACKGROUND: Adjuvant dabrafenib plus trametinib reduced the risk of relapse versus placebo in patients with resected, BRAFV600-mutant, stage III melanoma in the phase 3 COMBI-AD trial. This prespecified exploratory biomarker analysis aimed to evaluate potential prognostic or predictive factors and mechanisms of resistance to adjuvant targeted therapy. METHODS: COMBI-AD is a randomised, double-blind, placebo-controlled, phase 3 trial comparing dabrafenib 150 mg orally twice daily plus trametinib 2 mg orally once daily versus two matched placebos. Study participants were at least 18 years of age and underwent complete resection of stage IIIA (lymph node metastases >1 mm), IIIB, or IIIC cutaneous melanoma, per American Joint Committee on Cancer 7th edition criteria, with a BRAFV600E or BRAFV600K mutation. Patients were randomly assigned (1:1) to the two treatment groups by an interactive voice response system, stratified by mutation type and disease stage. Patients, physicians, and the investigators who analysed data were masked to treatment allocation. The primary outcome was relapse-free survival, defined as the time from randomisation to disease recurrence or death from any cause. Biomarker assessment was a prespecified exploratory outcome of the trial. We assessed intrinsic tumour genomic features by use of next-generation DNA sequencing and characteristics of the tumour microenvironment by use of a NanoString RNA assay, which might provide prognostic and predictive information. This trial is registered with ClinicalTrials.gov, number NCT01682083, and is ongoing but no longer recruiting participants. FINDINGS: Between Jan 31, 2013, and Dec 11, 2014, 870 patients were enrolled in the trial. Median follow-up at data cutoff (April 30, 2018) was 44 months (IQR 38-49) in the dabrafenib plus trametinib group and 42 months (21-49) in the placebo group. Intrinsic tumour genomic features were assessed in 368 patients (DNA sequencing set) and tumour microenvironment characteristics were assessed in 507 patients (NanoString biomarker set). MAPK pathway genomic alterations at baseline did not affect treatment benefit or clinical outcome. An IFNγ gene expression signature higher than the median was prognostic for prolonged relapse-free survival in both treatment groups. Tumour mutational burden was independently prognostic for relapse-free survival in the placebo group (high TMB, top third; hazard ratio [HR] 0·56, 95% CI 0·37-0·85, p=0·0056), but not in the dabrafenib plus trametinib group (0·83, 95% CI 0·53-1·32, p=0·44). Patients with tumour mutational burden in the lower two terciles seem to derive a substantial long-term relapse-free survival benefit from targeted therapy (HR [versus placebo] 0·49, 95% CI 0·35-0·68, p<0·0001). However, patients with high tumour mutational burden seem to have a less pronounced benefit with targeted therapy (HR [versus placebo] 0·75, 95% CI 0·44-1·26, p=0·27), especially if they had an IFNγ signature lower than the median (HR 0·88 [95% CI 0·40-1·93], p=0·74). INTERPRETATION: Tumour mutational burden alone or in combination with IFNγ gene expression signature or other markers for an adaptive immune response might be of relevance for identifying patients with stage III melanoma who might derive clinical benefit from targeted therapy. Further validation in prospective clinical trials is warranted. FUNDING: Novartis Pharmaceuticals.

Interchromosomal core duplicons drive both evolutionary instability and disease susceptibility of the Chromosome 8p23.1 region.

Recurrent rearrangements of Chromosome 8p23.1 are associated with congenital heart defects and developmental delay. The complexity of this region has led to inconsistencies in the current reference assembly, confounding studies of genetic variation. Using comparative sequence-based approaches, we generated a high-quality 6.3-Mbp alternate reference assembly of an inverted Chromosome 8p23.1 haplotype. Comparison with nonhuman primates reveals a 746-kbp duplicative transposition and two separate inversion events that arose in the last million years of human evolution. The breakpoints associated with these rearrangements map to an ape-specific interchromosomal core duplicon that clusters at sites of evolutionary inversion (P = 7.8 × 10-5). Refinement of microdeletion breakpoints identifies a subgroup of patients that map to the same interchromosomal core involved in the evolutionary formation of the duplication blocks. Our results define a higher-order genomic instability element that has shaped the structure of specific chromosomes during primate evolution contributing to rearrangements associated with inversion and disease.

Human copy number variation and complex genetic disease.

Copy number variants (CNVs) play an important role in human disease and population diversity. Advancements in technology have allowed for the analysis of CNVs in thousands of individuals with disease in addition to thousands of controls. These studies have identified rare CNVs associated with neuropsychiatric diseases such as autism, schizophrenia, and intellectual disability. In addition, copy number polymorphisms (CNPs) are present at higher frequencies in the population, show high diversity in copy number, sequence, and structure, and have been associated with multiple phenotypes, primarily related to immune or environmental response. However, the landscape of copy number variation still remains largely unexplored, especially for smaller CNVs and those embedded within complex regions of the human genome. An integrated approach including characterization of single nucleotide variants and CNVs in a large number of individuals with disease and normal genomes holds the promise of thoroughly elucidating the genetic basis of human disease and diversity.

Properties and rates of germline mutations in humans.

All genetic variation arises via new mutations; therefore, determining the rate and biases for different classes of mutation is essential for understanding the genetics of human disease and evolution. Decades of mutation rate analyses have focused on a relatively small number of loci because of technical limitations. However, advances in sequencing technology have allowed for empirical assessments of genome-wide rates of mutation. Recent studies have shown that 76% of new mutations originate in the paternal lineage and provide unequivocal evidence for an increase in mutation with paternal age. Although most analyses have focused on single nucleotide variants (SNVs), studies have begun to provide insight into the mutation rate for other classes of variation, including copy number variants (CNVs), microsatellites, and mobile element insertions (MEIs). Here, we review the genome-wide analyses for the mutation rate of several types of variants and suggest areas for future research.

Refinement and discovery of new hotspots of copy-number variation associated with autism spectrum disorder.

Rare copy-number variants (CNVs) have been implicated in autism and intellectual disability. These variants are large and affect many genes but lack clear specificity toward autism as opposed to developmental-delay phenotypes. We exploited the repeat architecture of the genome to target segmental duplication-mediated rearrangement hotspots (n = 120, median size 1.78 Mbp, range 240 kbp to 13 Mbp) and smaller hotspots flanked by repetitive sequence (n = 1,247, median size 79 kbp, range 3-96 kbp) in 2,588 autistic individuals from simplex and multiplex families and in 580 controls. Our analysis identified several recurrent large hotspot events, including association with 1q21 duplications, which are more likely to be identified in individuals with autism than in those with developmental delay (p = 0.01; OR = 2.7). Within larger hotspots, we also identified smaller atypical CNVs that implicated CHD1L and ACACA for the 1q21 and 17q12 deletions, respectively. Our analysis, however, suggested no overall increase in the burden of smaller hotspots in autistic individuals as compared to controls. By focusing on gene-disruptive events, we identified recurrent CNVs, including DPP10, PLCB1, TRPM1, NRXN1, FHIT, and HYDIN, that are enriched in autism. We found that as the size of deletions increases, nonverbal IQ significantly decreases, but there is no impact on autism severity; and as the size of duplications increases, autism severity significantly increases but nonverbal IQ is not affected. The absence of an increased burden of smaller CNVs in individuals with autism and the failure of most large hotspots to refine to single genes is consistent with a model where imbalance of multiple genes contributes to a disease state.

Genome-wide association study and admixture mapping reveal new loci associated with total IgE levels in Latinos.

BACKGROUND: IgE is a key mediator of allergic inflammation, and its levels are frequently increased in patients with allergic disorders. OBJECTIVE: We sought to identify genetic variants associated with IgE levels in Latinos. METHODS: We performed a genome-wide association study and admixture mapping of total IgE levels in 3334 Latinos from the Genes-environments & Admixture in Latino Americans (GALA II) study. Replication was evaluated in 454 Latinos, 1564 European Americans, and 3187 African Americans from independent studies. RESULTS: We confirmed associations of 6 genes identified by means of previous genome-wide association studies and identified a novel genome-wide significant association of a polymorphism in the zinc finger protein 365 gene (ZNF365) with total IgE levels (rs200076616, P = 2.3 × 10(-8)). We next identified 4 admixture mapping peaks (6p21.32-p22.1, 13p22-31, 14q23.2, and 22q13.1) at which local African, European, and/or Native American ancestry was significantly associated with IgE levels. The most significant peak was 6p21.32-p22.1, where Native American ancestry was associated with lower IgE levels (P = 4.95 × 10(-8)). All but 22q13.1 were replicated in an independent sample of Latinos, and 2 of the peaks were replicated in African Americans (6p21.32-p22.1 and 14q23.2). Fine mapping of 6p21.32-p22.1 identified 6 genome-wide significant single nucleotide polymorphisms in Latinos, 2 of which replicated in European Americans. Another single nucleotide polymorphism was peak-wide significant within 14q23.2 in African Americans (rs1741099, P = 3.7 × 10(-6)) and replicated in non-African American samples (P = .011). CONCLUSION: We confirmed genetic associations at 6 genes and identified novel associations within ZNF365, HLA-DQA1, and 14q23.2. Our results highlight the importance of studying diverse multiethnic populations to uncover novel loci associated with total IgE levels.

Population-genetic properties of differentiated human copy-number polymorphisms.

Copy-number variants (CNVs) can reach appreciable frequencies in the human population, and recent discoveries have shown that several of these copy-number polymorphisms (CNPs) are associated with human diseases, including lupus, psoriasis, Crohn disease, and obesity. Despite new advances, significant biases remain in terms of CNP discovery and genotyping. We developed a method based on single-channel intensity data and benchmarked against copy numbers determined from sequencing read depth to successfully obtain CNP genotypes for 1495 CNPs from 487 human DNA samples of diverse ethnic backgrounds. This microarray contained CNPs in segmental duplication-rich regions and insertions of sequences not represented in the reference genome assembly or on standard SNP microarray platforms. We observe that CNPs in segmental duplications are more likely to be population differentiated than CNPs in unique regions (p = 0.015) and that biallelic CNPs show greater stratification when compared to frequency-matched SNPs (p = 0.0026). Although biallelic CNPs show a strong correlation of copy number with flanking SNP genotypes, the majority of multicopy CNPs do not (40% with r > 0.8). We selected a subset of CNPs for further characterization in 1876 additional samples from 62 populations; this revealed striking population-differentiated structural variants in genes of clinical significance such as OCLN, a tight junction protein involved in hepatitis C viral entry. Our microarray design allows these variants to be rapidly tested for disease association and our results suggest that CNPs (especially those that cannot be imputed from SNP genotypes) might have contributed disproportionately to human diversity and selection.

Demonstrating stratification in a European American population.

Population stratification occurs in case-control association studies when allele frequencies differ between cases and controls because of ancestry. Stratification may lead to false positive associations, although this issue remains controversial. Empirical studies have found little evidence of stratification in European-derived populations, but potentially significant levels of stratification could not be ruled out. We studied a European American panel discordant for height, a heritable trait that varies widely across Europe. Genotyping 178 SNPs and applying standard analytical methods yielded no evidence of stratification. But a SNP in the gene LCT that varies widely in frequency across Europe was strongly associated with height (P < 10(-6)). This apparent association was largely or completely due to stratification; rematching individuals on the basis of European ancestry greatly reduced the apparent association, and no association was observed in Polish or Scandinavian individuals. The failure of standard methods to detect this stratification indicates that new methods may be required.

Next-generation sequencing of baseline genetic mutations and outcomes of eltrombopag and azacitidine therapy in patients with myelodysplastic syndromes and thrombocytopenia: Data from the SUPPORT clinical trial.

Eltrombopag has been previously shown to be effective in reversing azacitidine-mediated thrombocytopenia. This was further investigated in the SUPPORT trial, a phase III study assessing the efficacy/safety of eltrombopag plus azacitidine in patients with intermediate- to high-risk myelodysplastic syndromes and thrombocytopenia. The results did not support a clinical benefit for the addition of eltrombopag to azacitidine. We investigated if the somatic mutational profiles in the patient cohort were associated with treatment outcomes. Based on the available data, we observed no imbalance in the mutational profiles between treatment arms or a clear association between identified somatic mutations and clinical outcomes.

Clonal hematopoiesis detection in patients with cancer using cell-free DNA sequencing.

In the context of cancer, clonal hematopoiesis of indeterminate potential (CHIP) is associated with the development of therapy-related myeloid neoplasms and shorter overall survival. Cell-free DNA (cfDNA) sequencing is becoming widely adopted for genomic screening of patients with cancer but has not been used extensively to determine CHIP status because of a requirement for matched blood and tumor sequencing. We present an accurate classification approach to determine the CH status from cfDNA sequencing alone, applying our model to 4324 oncology clinical cfDNA samples. Using this method, we determined that 30.3% of patients in this cohort have evidence of CH, and the incidence of CH varies by tumor type. Matched RNA sequencing data show evidence of increased inflammation, especially neutrophil activation, within the tumors and tumor microenvironments of patients with CH. In addition, patients with CH had evidence of neutrophil activation systemically, pointing to a potential mechanism of action for the worse outcomes associated with CH status. Neutrophil activation may be one of many mechanisms, however, because patients with estrogen receptor-positive breast cancer harboring TET2 frameshift mutations had worse outcomes but similar neutrophil frequencies to patients without CH. Together, these data show the feasibility of detecting CH through cfDNA sequencing alone and an application of this method, demonstrating increased inflammation in patients with CH both systemically and in the tumor microenvironment.

Initial Evidence for the Efficacy of Naporafenib in Combination With Trametinib in NRAS-Mutant Melanoma: Results From the Expansion Arm of a Phase Ib, Open-Label Study.

PURPOSE: No approved targeted therapy for the treatment of patients with neuroblastoma RAS viral (v-ras) oncogene homolog (NRAS)-mutant melanoma is currently available. PATIENTS AND METHODS: In this phase Ib escalation/expansion study (ClinicalTrials.gov identifier: NCT02974725), the safety, tolerability, and preliminary antitumor activity of naporafenib (LXH254), a BRAF/CRAF protein kinases inhibitor, were explored in combination with trametinib in patients with advanced/metastatic KRAS- or BRAF-mutant non-small-cell lung cancer (escalation arm) or NRAS-mutant melanoma (escalation and expansion arms). RESULTS: Thirty-six and 30 patients were enrolled in escalation and expansion, respectively. During escalation, six patients reported grade ≥3 dose-limiting toxicities, including dermatitis acneiform (n = 2), maculopapular rash (n = 2), increased lipase (n = 1), and Stevens-Johnson syndrome (n = 1). The recommended doses for expansion were naporafenib 200 mg twice a day plus trametinib 1 mg once daily and naporafenib 400 mg twice a day plus trametinib 0.5 mg once daily. During expansion, all 30 patients experienced a treatment-related adverse event, the most common being rash (80%, n = 24), blood creatine phosphokinase increased, diarrhea, and nausea (30%, n = 9 each). In expansion, the objective response rate, median duration of response, and median progression-free survival were 46.7% (95% CI, 21.3 to 73.4; 7 of 15 patients), 3.75 (95% CI, 1.97 to not estimable [NE]) months, and 5.52 months, respectively, in patients treated with naporafenib 200 mg twice a day plus trametinib 1 mg once daily, and 13.3% (95% CI, 1.7 to 40.5; 2 of 15 patients), 3.75 (95% CI, 2.04 to NE) months, and 4.21 months, respectively, in patients treated with naporafenib 400 mg twice a day plus trametinib 0.5 mg once daily. CONCLUSION: Naporafenib plus trametinib showed promising preliminary antitumor activity in patients with NRAS-mutant melanoma. Prophylactic strategies aimed to lower the incidence of skin-related events are under investigation.

Combined PD-1, BRAF and MEK inhibition in BRAFV600E colorectal cancer: a phase 2 trial.

While BRAF inhibitor combinations with EGFR and/or MEK inhibitors have improved clinical efficacy in BRAFV600E colorectal cancer (CRC), response rates remain low and lack durability. Preclinical data suggest that BRAF/MAPK pathway inhibition may augment the tumor immune response. We performed a proof-of-concept single-arm phase 2 clinical trial of combined PD-1, BRAF and MEK inhibition with sparatlizumab (PDR001), dabrafenib and trametinib in 37 patients with BRAFV600E CRC. The primary end point was overall response rate, and the secondary end points were progression-free survival, disease control rate, duration of response and overall survival. The study met its primary end point with a confirmed response rate (24.3% in all patients; 25% in microsatellite stable patients) and durability that were favorable relative to historical controls of BRAF-targeted combinations alone. Single-cell RNA sequencing of 23 paired pretreatment and day 15 on-treatment tumor biopsies revealed greater induction of tumor cell-intrinsic immune programs and more complete MAPK inhibition in patients with better clinical outcome. Immune program induction in matched patient-derived organoids correlated with the degree of MAPK inhibition. These data suggest a potential tumor cell-intrinsic mechanism of cooperativity between MAPK inhibition and immune response, warranting further clinical evaluation of optimized targeted and immune combinations in CRC. ClinicalTrials.gov registration: NCT03668431.

TET2-Driven Clonal Hematopoiesis and Response to Canakinumab: An Exploratory Analysis of the CANTOS Randomized Clinical Trial.

Importance: Clonal hematopoiesis of indeterminate potential (CHIP) is associated with increased risk of atherosclerotic cardiovascular disease, and mouse experiments suggest that CHIP related to Tet2 loss of function in myeloid cells accelerates atherosclerosis via augmented interleukin (IL) 1ß signaling. Objective: To assess whether individuals with CHIP have greater cardiovascular event reduction in response to IL-1ß neutralization in the Canankinumab Anti-inflammatory Thrombosis Outcomes Trial (CANTOS). Design, Setting, and Participants: This randomized clinical trial took place from April 2011 to June 2017 at more than 1000 clinical sites in 39 countries. Targeted deep sequencing of genes previously associated with CHIP in a subset of trial participants using genomic DNA prepared from baseline peripheral blood samples were analyzed. All participants had prior myocardial infarction and elevated high-sensitivity C-reactive protein level above 0.20 mg/dL. Analysis took place between June 2017 and December 2021. Interventions: Canakinumab, an anti-IL-1ß antibody, given at doses of 50, 150, and 300 mg once every 3 months. Main Outcomes and Measures: Major adverse cardiovascular events (MACE). Results: A total of 338 patients (8.6%) were identified in this subset with evidence for clonal hematopoiesis. As expected, the incidence of CHIP increased with age; the mean (SD) age of patients with CHIP was 66.3 (9.2) years and 61.5 (9.6) years in patients without CHIP. Unlike other populations that were not preselected for elevated C-reactive protein, in the CANTOS population variants in TET2 were more common than DNMT3A (119 variants in 103 patients vs 86 variants in 85 patients). Placebo-treated patients with CHIP showed a nonsignificant increase in the rate of MACE compared with patients without CHIP using a Cox proportional hazard model (hazard ratio, 1.32 [95% CI, 0.86-2.04]; P = .21). Exploratory analyses of placebo-treated patients with a somatic variant in either TET2 or DNMT3A (n = 58) showed an equivocal risk for MACE (hazard ratio, 1.65 [95% CI, 0.97-2.80]; P = .06). Patients with CHIP due to somatic variants in TET2 also had reduced risk for MACE while taking canakinumab (hazard ratio, 0.38 [95% CI, 0.15-0.96]) with equivocal difference compared with others (P for interaction = .14). Conclusions and Relevance: These results are consistent with observations of increased risk for cardiovascular events in patients with CHIP and raise the possibility that those with TET2 variants may respond better to canakinumab than those without CHIP. Future studies are required to further substantiate this hypothesis. Trial Registration: ClinicalTrials.gov Identifier: NCT01327846.

Cell-free DNA captures tumor heterogeneity and driver alterations in rapid autopsies with pre-treated metastatic cancer.

In patients with metastatic cancer, spatial heterogeneity of somatic alterations may lead to incomplete assessment of a cancer's mutational profile when analyzing a single tumor biopsy. In this study, we perform sequencing of cell-free DNA (cfDNA) and distinct metastatic tissue samples from ten rapid autopsy cases with pre-treated metastatic cancer. We show that levels of heterogeneity in genetic biomarkers vary between patients but that gene expression signatures representative of the tumor microenvironment are more consistent. Across nine patients with plasma samples available, we are able to detect 62/62 truncal and 47/121 non-truncal point mutations in cfDNA. We observe that mutation clonality in cfDNA is correlated with the number of metastatic lesions in which the mutation is detected and use this result to derive a clonality threshold to classify truncal and non-truncal driver alterations with reasonable specificity. In contrast, mutation truncality is more often incorrectly assigned when studying single tissue samples. Our results demonstrate the utility of a single cfDNA sample relative to that of single tissue samples when treating patients with metastatic cancer.

Combined PD-1, BRAF and MEK inhibition in advanced BRAF-mutant melanoma: safety run-in and biomarker cohorts of COMBI-i.

Immune and targeted therapies achieve long-term survival in metastatic melanoma; however, new treatment strategies are needed to improve patients' outcomes1,2. We report on the efficacy, safety and biomarker analysis from the single-arm safety run-in (part 1; n = 9) and biomarker (part 2; n = 27) cohorts of the randomized, placebo-controlled, phase 3 COMBI-i trial (NCT02967692) of the anti-PD-1 antibody spartalizumab, in combination with the BRAF inhibitor dabrafenib and MEK inhibitor trametinib. Patients (n = 36) had previously untreated BRAF V600-mutant unresectable or metastatic melanoma. In part 1, the recommended phase 3 regimen was identified based on the incidence of dose-limiting toxicities (DLTs; primary endpoint): 400 mg of spartalizumab every 4 weeks plus 150 mg of dabrafenib twice daily plus 2 mg of trametinib once daily. Part 2 characterized changes in PD-L1 levels and CD8+ cells following treatment (primary endpoint), and analyzed additional biomarkers. Assessments of efficacy and safety were key secondary endpoints (median follow-up, 24.3 months). Spartalizumab plus dabrafenib and trametinib led to an objective response rate (ORR) of 78%, including 44% complete responses (CRs). Grade ≥3 treatment-related adverse events (TRAEs) were experienced by 72% of patients. All patients had temporary dose modifications, and 17% permanently discontinued all three study drugs due to TRAEs. Early progression-free survival (PFS) events were associated with low tumor mutational burden/T cell-inflamed gene expression signature (GES) or high immunosuppressive tumor microenvironment (TME) GES levels at baseline; an immunosuppressive TME may also preclude CR. Overall, the efficacy, safety and on-treatment biomarker modulations associated with spartalizumab plus dabrafenib and trametinib are promising, and biomarkers that may predict long-term benefit were identified.

BRAF-Mutant Transcriptional Subtypes Predict Outcome of Combined BRAF, MEK, and EGFR Blockade with Dabrafenib, Trametinib, and Panitumumab in Patients with Colorectal Cancer.

PURPOSE: The influence of the transcriptional and immunologic context of mutations on therapeutic outcomes with targeted therapy in cancer has not been well defined. BRAF V600E-mutant (BM) colorectal cancer comprises two main transcriptional subtypes, BM1 and BM2. We sought to determine the impact of BM subtype, as well as distinct biological features of those subtypes, on response to BRAF/MEK/EGFR inhibition in patients with colorectal cancer. PATIENTS AND METHODS: Paired fresh tumor biopsies were acquired at baseline and on day 15 of treatment from all consenting patients with BM colorectal cancer enrolled in a phase II clinical trial of dabrafenib, trametinib, and panitumumab. For each sample, BM subtype, cell cycle, and immune gene signature expression were determined using RNA-sequencing (RNA-seq), and a Cox proportional hazards model was applied to determine association with progression-free survival (PFS). RESULTS: Confirmed response rates, median PFS, and median overall survival (OS) were higher in BM1 subtype patients compared with BM2 subtype patients. Evaluation of immune contexture identified greater immune reactivity in BM1, whereas cell-cycle signatures were more highly expressed in BM2. A multivariate model of PFS incorporating BM subtype plus immune and cell-cycle signatures revealed that BM subtype encompasses the majority of the effect. CONCLUSIONS: BM subtype is significantly associated with the outcome of combination dabrafenib, trametinib, and panitumumab therapy and may serve as a standalone predictive biomarker beyond mutational status. Our findings support a more nuanced approach to targeted therapeutic decisions that incorporates assessment of transcriptional context.

Tracking the Evolution of Resistance to ALK Tyrosine Kinase Inhibitors through Longitudinal Analysis of Circulating Tumor DNA.

PURPOSE: ALK rearrangements predict for sensitivity to ALK tyrosine kinase inhibitors (TKIs). However, responses to ALK TKIs are generally short-lived. Serial molecular analysis is an informative strategy for identifying genetic mediators of resistance. Although multiple studies support the clinical benefits of repeat tissue sampling, the clinical utility of longitudinal circulating tumor DNA analysis has not been established in ALK-positive lung cancer. METHODS: Using a 566-gene hybrid-capture next-generation sequencing (NGS) assay, we performed longitudinal analysis of plasma specimens from 22 ALK-positive patients with acquired resistance to ALK TKIs to track the evolution of resistance during treatment. To determine tissue-plasma concordance, we compared plasma findings to results of repeat biopsies. RESULTS: At progression, we detected an ALK fusion in plasma from 19 (86%) of 22 patients, and identified ALK resistance mutations in plasma specimens from 11 (50%) patients. There was 100% agreement between tissue- and plasma-detected ALK fusions. Among 16 cases where contemporaneous plasma and tissue specimens were available, we observed 100% concordance between ALK mutation calls. ALK mutations emerged and disappeared during treatment with sequential ALK TKIs, suggesting that plasma mutation profiles were dependent on the specific TKI administered. ALK G1202R, the most frequent plasma mutation detected after progression on a second-generation TKI, was consistently suppressed during treatment with lorlatinib. CONCLUSIONS: Plasma genotyping by NGS is an effective method for detecting ALK fusions and ALK mutations in patients progressing on ALK TKIs. The correlation between plasma ALK mutations and response to distinct ALK TKIs highlights the potential for plasma analysis to guide selection of ALK-directed therapies.

PureCN: copy number calling and SNV classification using targeted short read sequencing.

BACKGROUND: Matched sequencing of both tumor and normal tissue is routinely used to classify variants of uncertain significance (VUS) into somatic vs. germline. However, assays used in molecular diagnostics focus on known somatic alterations in cancer genes and often only sequence tumors. Therefore, an algorithm that reliably classifies variants would be helpful for retrospective exploratory analyses. Contamination of tumor samples with normal cells results in differences in expected allelic fractions of germline and somatic variants, which can be exploited to accurately infer genotypes after adjusting for local copy number. However, existing algorithms for determining tumor purity, ploidy and copy number are not designed for unmatched short read sequencing data. RESULTS: We describe a methodology and corresponding open source software for estimating tumor purity, copy number, loss of heterozygosity (LOH), and contamination, and for classification of single nucleotide variants (SNVs) by somatic status and clonality. This R package, PureCN, is optimized for targeted short read sequencing data, integrates well with standard somatic variant detection pipelines, and has support for matched and unmatched tumor samples. Accuracy is demonstrated on simulated data and on real whole exome sequencing data. CONCLUSIONS: Our algorithm provides accurate estimates of tumor purity and ploidy, even if matched normal samples are not available. This in turn allows accurate classification of SNVs. The software is provided as open source (Artistic License 2.0) R/Bioconductor package PureCN (http://bioconductor.org/packages/PureCN/).