Spectrum of driver mutations and clinical impact of circulating tumor DNA analysis in non-small cell lung cancer: Analysis of over 8000 cases.

BACKGROUND: Circulating cell-free tumor DNA (ctDNA)-based mutation profiling, if sufficiently sensitive and comprehensive, can efficiently identify genomic targets in advanced lung adenocarcinoma. Therefore, the authors investigated the accuracy and clinical utility of a commercially available digital next-generation sequencing platform in a large series of patients with non-small cell lung cancer (NSCLC). METHODS: Plasma-based comprehensive genomic profiling results from 8388 consecutively tested patients with advanced NSCLC were analyzed. Driver and resistance mutations were examined with regard to their distribution, frequency, co-occurrence, and mutual exclusivity. RESULTS: Somatic alterations were detected in 86% of samples. The median variant allele fraction was 0.43% (range, 0.03%-97.62%). Activating alterations in actionable oncogenes were identified in 48% of patients, including EGFR (26.4%), MET (6.1%), and BRAF (2.8%) alterations and fusions (ALK, RET, and ROS1) in 2.3%. Treatment-induced resistance mutations were common in this cohort, including driver-dependent and driver-independent alterations. In the subset of patients who had progressive disease during EGFR therapy, 64% had known or putative resistance alterations detected in plasma. Subset analysis revealed that ctDNA increased the identification of driver mutations by 65% over standard-of-care, tissue-based testing at diagnosis. A pooled data analysis on this plasma-based assay demonstrated that targeted therapy response rates were equivalent to those reported from tissue analysis. CONCLUSIONS: Comprehensive ctDNA analysis detected the presence of therapeutically targetable driver and resistance mutations at the frequencies and distributions predicted for the study population. These findings add support for comprehensive ctDNA testing in patients who are incompletely tested at the time of diagnosis and as a primary option at the time of progression on targeted therapies.

Co-evolution of transcriptional silencing proteins and the DNA elements specifying their assembly.

Co-evolution of transcriptional regulatory proteins and their sites of action has been often hypothesized but rarely demonstrated. Here we provide experimental evidence of such co-evolution in yeast silent chromatin, a finding that emerged from studies of hybrids formed between two closely related Saccharomyces species. A unidirectional silencing incompatibility between S. cerevisiae and S. bayanus led to a key discovery: asymmetrical complementation of divergent orthologs of the silent chromatin component Sir4. In S. cerevisiae/S. bayanus interspecies hybrids, ChIP-Seq analysis revealed a restriction against S. cerevisiae Sir4 associating with most S. bayanus silenced regions; in contrast, S. bayanus Sir4 associated with S. cerevisiae silenced loci to an even greater degree than did S. cerevisiae's own Sir4. Functional changes in silencer sequences paralleled changes in Sir4 sequence and a reduction in Sir1 family members in S. cerevisiae. Critically, species-specific silencing of the S. bayanus HMR locus could be reconstituted in S. cerevisiae by co-transfer of the S. bayanus Sir4 and Kos3 (the ancestral relative of Sir1) proteins. As Sir1/Kos3 and Sir4 bind conserved silencer-binding proteins, but not specific DNA sequences, these rapidly evolving proteins served to interpret differences in the two species' silencers presumably involving emergent features created by the regulatory proteins that bind sequences within silencers. The results presented here, and in particular the high resolution ChIP-Seq localization of the Sir4 protein, provided unanticipated insights into the mechanism of silent chromatin assembly in yeast.

Systematic discovery of neoepitope-HLA pairs for neoantigens shared among patients and tumor types.

The broad application of precision cancer immunotherapies is limited by the number of validated neoepitopes that are common among patients or tumor types. To expand the known repertoire of shared neoantigen-human leukocyte antigen (HLA) complexes, we developed a high-throughput platform that coupled an in vitro peptide-HLA binding assay with engineered cellular models expressing individual HLA alleles in combination with a concatenated transgene harboring 47 common cancer neoantigens. From more than 24,000 possible neoepitope-HLA combinations, biochemical and computational assessment yielded 844 unique candidates, of which 86 were verified after immunoprecipitation mass spectrometry analyses of engineered, monoallelic cell lines. To evaluate the potential for immunogenicity, we identified T cell receptors that recognized select neoepitope-HLA pairs and elicited a response after introduction into human T cells. These cellular systems and our data on therapeutically relevant neoepitopes in their HLA contexts will aid researchers studying antigen processing as well as neoepitope targeting therapies.

Origins and Timing of Emerging Lesions in Advanced Renal Cell Carcinoma.

Renal cell carcinoma (RCC) with venous tumor thrombus (VTT) arising from the primary tumor occurs in approximately 10% of cases and is thought to represent more advanced disease. The intravascular nature of VTT suggests that it may serve as a source for hematogenous metastases. RCC with VTT and distant metastasis provides unique opportunities to examine the origins and emergence timing of these distinct tumor lesions, and to identify molecular correlates with disease state. We performed multi-region exome and RNA-sequencing analysis of 16 patients with RCC with VTT, with eight patients also having sequenced metastasis, to identify genomic alterations, biological pathways, and evolutionary processes contributing to VTT and metastasis, and to ask whether metastasis arises directly from or independent of VTT. No specific genomic alterations were associated with VTT. Hallmark copy-number alterations (deletions of 14q, 8p, and 4q) were associated with metastasis and disease recurrence, and secondary driver alterations tended to accumulate in metastatic lineages. Mismatch repair mutational signatures co-occurred across most tumors, suggesting a role for intracellular DNA damage in RCC. Robust phylogenetic timing analysis indicated that metastasis typically emerged before VTT, rather than deriving from it, with the earliest metastases predicted to emerge years before diagnosis. As a result, VTT in metastatic cases frequently derived from a metastatic lineage. Relative to the primary tumor, VTT upregulated immediate-early genes and transcriptional targets of the TNFα/NF-κB pathway, whereas metastases upregulated MTOR and transcriptional targets downstream of mTORC1 activation. IMPLICATIONS: These results suggest that VTT and metastasis formation occur independently, VTT presence alone does not necessarily imply more advanced disease with inevitably poor prognosis.

Single-cell analysis reveals clonally expanded tumor-associated CD57+ CD8 T cells are enriched in the periphery of patients with metastatic urothelial cancer responding to PD-L1 blockade.

BACKGROUND: A growing body of evidence suggests that T-cell responses against neoantigens are critical regulators of response to immune checkpoint blockade. We previously showed that circulating neoantigen-specific CD8 T cells in patients with lung cancer responding to anti-Programmed death-ligand 1 (PD-L1) (atezolizumab) exhibit a unique phenotype with high expression of CD57, CD244, and KLRG1. Here, we extended our analysis on neoantigen-specific CD8 T cells to patients with metastatic urothelial cancer (mUC) and further profiled total CD8 T cells to identify blood-based predictive biomarkers of response to atezolizumab. METHODS: We identified tumor neoantigens from 20 patients with mUC and profiled their peripheral CD8 T cells using highly multiplexed combinatorial tetramer staining. Another set of patients with mUC treated with atezolizumab (n=30) or chemotherapy (n=40) were selected to profile peripheral CD8 T cells by mass cytometry. Using single-cell transcriptional analysis (single-cell RNA sequencing (scRNA-seq)), together with CITE-seq (cellular indexing of transcriptomes and epitopes by sequencing) and paired T-cell receptor (TCR) sequencing, we further characterized peripheral CD8 T cells in a subset of patients (n=16). RESULTS: High frequency of CD57 was observed in neoantigen-specific CD8 T cells in patients with mUC responding to atezolizumab. Extending these findings to bulk CD8 T cells, we found higher frequency of CD57 expressing CD8 T cells before treatment in patients responding to atezolizumab (n=20, p<0.01) but not to chemotherapy. These findings were corroborated in a validation cohort (n=30, p<0.01) and notably were independent of known biomarkers of response. scRNA-seq analysis identified a clonally expanded cluster enriched within CD57+ CD8 T cells in responding patients characterized by higher expression of genes associated with activation, cytotoxicity, and tissue-resident memory markers. Furthermore, compared with CD57- CD8 T cells, TCRs of CD57+ CD8 T cells showed increased overlap with the TCR repertoire of tumor-infiltrating T cells. CONCLUSIONS: Collectively, we show high frequencies of CD57 among neoantigen-specific and bulk CD8 T cells in patients responding to atezolizumab. The TCR repertoire overlap between peripheral CD57+ CD8 T cells and tumor-infiltrating lymphocytes suggest that accumulation of peripheral CD57+ CD8 T cells is reflective of an ongoing antitumor T-cell response. Our findings provide evidence and rationale for using circulating CD8 T cells expressing CD57 as a readily accessible blood-based biomarker for selecting patients with mUC for atezolizumab therapy.

ARID1A mutations confer intrinsic and acquired resistance to cetuximab treatment in colorectal cancer.

Most colorectal (CRC) tumors are dependent on EGFR/KRAS/BRAF/MAPK signaling activation. ARID1A is an epigenetic regulator mutated in approximately 5% of non-hypermutated CRC tumors. Here we show that anti-EGFR but not anti-VEGF treatment enriches for emerging ARID1A mutations in CRC patients. In addition, we find that patients with ARID1A mutations, at baseline, are associated with worse outcome when treated with cetuximab- but not bevacizumab-containing therapies; thus, this suggests that ARID1A mutations may provide both an acquired and intrinsic mechanism of resistance to anti-EGFR therapies. We find that, ARID1A and EGFR-pathway genetic alterations are mutually exclusive across lung and colorectal cancers, further supporting a functional connection between these pathways. Our results not only suggest that ARID1A could be potentially used as a predictive biomarker for cetuximab treatment decisions but also provide a rationale for exploring therapeutic MAPK inhibition in an unexpected but genetically defined segment of CRC patients.

Indication-specific tumor evolution and its impact on neoantigen targeting and biomarkers for individualized cancer immunotherapies.

BACKGROUND: Individualized neoantigen-specific immunotherapy (iNeST) requires robustly expressed clonal neoantigens for efficacy, but tumor mutational heterogeneity, loss of neoantigen expression, and variable tissue sampling present challenges. It is assumed that clonal neoantigens are preferred targets for immunotherapy, but the distributions of clonal neoantigens are not well characterized across cancer types. METHODS: We combined multiregion sequencing (MR-seq) analysis of five untreated, synchronously sampled metastatic solid tumors with re-analysis of published MR-seq data from 103 patients in order to characterize their globally clonal neoantigen content and factors that would impact neoantigen targeting. RESULTS: Branching evolution in colorectal cancer and renal cell carcinoma led to fewer clonal neoantigens and to clade-specific neoantigens (those shared across a subset of tumor regions but not fully clonal), with the latter not being readily distinguishable in single tumor samples. In colorectal, renal, and bladder cancer, most tumors had few globally clonal neoantigens. Prioritizing mutations with higher purity-adjusted and ploidy-adjusted variant allele frequency enriched for globally clonal neoantigens (those found in all tumor regions), whereas estimated cancer cell fraction derived from clustering-based tools, surprisingly, did not. Neoantigen quality was associated with loss of neoantigen expression in the bladder cancer case, and HLA-allele loss was observed in the renal and non-small cell lung cancer cases. CONCLUSIONS: We show that tumor type, multilesion sampling, neoantigen expression, and HLA allele retention are important factors for iNeST targeting and patient selection, and may also be important factors to consider in the development of biomarker strategies.

The Landscape of Actionable Genomic Alterations in Cell-Free Circulating Tumor DNA from 21,807 Advanced Cancer Patients.

Purpose: Cell-free DNA (cfDNA) sequencing provides a noninvasive method for obtaining actionable genomic information to guide personalized cancer treatment, but the presence of multiple alterations in circulation related to treatment and tumor heterogeneity complicate the interpretation of the observed variants.Experimental Design: We describe the somatic mutation landscape of 70 cancer genes from cfDNA deep-sequencing analysis of 21,807 patients with treated, late-stage cancers across >50 cancer types. To facilitate interpretation of the genomic complexity of circulating tumor DNA in advanced, treated cancer patients, we developed methods to identify cfDNA copy-number driver alterations and cfDNA clonality.Results: Patterns and prevalence of cfDNA alterations in major driver genes for non-small cell lung, breast, and colorectal cancer largely recapitulated those from tumor tissue sequencing compendia (The Cancer Genome Atlas and COSMIC; r = 0.90-0.99), with the principal differences in alteration prevalence being due to patient treatment. This highly sensitive cfDNA sequencing assay revealed numerous subclonal tumor-derived alterations, expected as a result of clonal evolution, but leading to an apparent departure from mutual exclusivity in treatment-naïve tumors. Upon applying novel cfDNA clonality and copy-number driver identification methods, robust mutual exclusivity was observed among predicted truncal driver cfDNA alterations (FDR = 5 × 10-7 for EGFR and ERBB2), in effect distinguishing tumor-initiating alterations from secondary alterations. Treatment-associated resistance, including both novel alterations and parallel evolution, was common in the cfDNA cohort and was enriched in patients with targetable driver alterations (>18.6% patients).Conclusions: Together, these retrospective analyses of a large cfDNA sequencing data set reveal subclonal structures and emerging resistance in advanced solid tumors. Clin Cancer Res; 24(15); 3528-38. ©2018 AACR.

Validation of a Plasma-Based Comprehensive Cancer Genotyping Assay Utilizing Orthogonal Tissue- and Plasma-Based Methodologies.

Purpose: To analytically and clinically validate a circulating cell-free tumor DNA sequencing test for comprehensive tumor genotyping and demonstrate its clinical feasibility.Experimental Design: Analytic validation was conducted according to established principles and guidelines. Blood-to-blood clinical validation comprised blinded external comparison with clinical droplet digital PCR across 222 consecutive biomarker-positive clinical samples. Blood-to-tissue clinical validation comprised comparison of digital sequencing calls to those documented in the medical record of 543 consecutive lung cancer patients. Clinical experience was reported from 10,593 consecutive clinical samples.Results: Digital sequencing technology enabled variant detection down to 0.02% to 0.04% allelic fraction/2.12 copies with ≤0.3%/2.24-2.76 copies 95% limits of detection while maintaining high specificity [prevalence-adjusted positive predictive values (PPV) >98%]. Clinical validation using orthogonal plasma- and tissue-based clinical genotyping across >750 patients demonstrated high accuracy and specificity [positive percent agreement (PPAs) and negative percent agreement (NPAs) >99% and PPVs 92%-100%]. Clinical use in 10,593 advanced adult solid tumor patients demonstrated high feasibility (>99.6% technical success rate) and clinical sensitivity (85.9%), with high potential actionability (16.7% with FDA-approved on-label treatment options; 72.0% with treatment or trial recommendations), particularly in non-small cell lung cancer, where 34.5% of patient samples comprised a directly targetable standard-of-care biomarker.Conclusions: High concordance with orthogonal clinical plasma- and tissue-based genotyping methods supports the clinical accuracy of digital sequencing across all four types of targetable genomic alterations. Digital sequencing's clinical applicability is further supported by high rates of technical success and biomarker target discovery. Clin Cancer Res; 24(15); 3539-49. ©2018 AACR.

Cell-Free DNA Next-Generation Sequencing in Pancreatobiliary Carcinomas.

UNLABELLED: Patients with pancreatic and biliary carcinomas lack personalized treatment options, in part because biopsies are often inadequate for molecular characterization. Cell-free DNA (cfDNA) sequencing may enable a precision oncology approach in this setting. We attempted to prospectively analyze 54 genes in tumor and cfDNA for 26 patients. Tumor sequencing failed in 9 patients (35%). In the remaining 17, 90.3% (95% confidence interval, 73.1%-97.5%) of mutations detected in tumor biopsies were also detected in cfDNA. The diagnostic accuracy of cfDNA sequencing was 97.7%, with 92.3% average sensitivity and 100% specificity across five informative genes. Changes in cfDNA correlated well with tumor marker dynamics in serial sampling (r = 0.93). We demonstrate that cfDNA sequencing is feasible, accurate, and sensitive in identifying tumor-derived mutations without prior knowledge of tumor genotype or the abundance of circulating tumor DNA. cfDNA sequencing should be considered in pancreatobiliary cancer trials where tissue sampling is unsafe, infeasible, or otherwise unsuccessful. SIGNIFICANCE: Precision medicine efforts in biliary and pancreatic cancers have been frustrated by difficulties in obtaining adequate tumor tissue for next-generation sequencing. cfDNA sequencing reliably and accurately detects tumor-derived mutations, paving the way for precision oncology approaches in these deadly diseases.

Analytical and Clinical Validation of a Digital Sequencing Panel for Quantitative, Highly Accurate Evaluation of Cell-Free Circulating Tumor DNA.

Next-generation sequencing of cell-free circulating solid tumor DNA addresses two challenges in contemporary cancer care. First this method of massively parallel and deep sequencing enables assessment of a comprehensive panel of genomic targets from a single sample, and second, it obviates the need for repeat invasive tissue biopsies. Digital Sequencing™ is a novel method for high-quality sequencing of circulating tumor DNA simultaneously across a comprehensive panel of over 50 cancer-related genes with a simple blood test. Here we report the analytic and clinical validation of the gene panel. Analytic sensitivity down to 0.1% mutant allele fraction is demonstrated via serial dilution studies of known samples. Near-perfect analytic specificity (> 99.9999%) enables complete coverage of many genes without the false positives typically seen with traditional sequencing assays at mutant allele frequencies or fractions below 5%. We compared digital sequencing of plasma-derived cell-free DNA to tissue-based sequencing on 165 consecutive matched samples from five outside centers in patients with stage III-IV solid tumor cancers. Clinical sensitivity of plasma-derived NGS was 85.0%, comparable to 80.7% sensitivity for tissue. The assay success rate on 1,000 consecutive samples in clinical practice was 99.8%. Digital sequencing of plasma-derived DNA is indicated in advanced cancer patients to prevent repeated invasive biopsies when the initial biopsy is inadequate, unobtainable for genomic testing, or uninformative, or when the patient's cancer has progressed despite treatment. Its clinical utility is derived from reduction in the costs, complications and delays associated with invasive tissue biopsies for genomic testing.

Prospective blinded study of somatic mutation detection in cell-free DNA utilizing a targeted 54-gene next generation sequencing panel in metastatic solid tumor patients.

Sequencing of the mutant allele fraction of circulating cell-free DNA (cfDNA) derived from tumors is increasingly utilized to detect actionable genomic alterations in cancer. We conducted a prospective blinded study of a comprehensive cfDNA sequencing panel with 54 cancer genes. To evaluate the concordance between cfDNA and tumor DNA (tDNA), sequencing results were compared between cfDNA from plasma and genomic tumor DNA (tDNA). Utilizing next generation digital sequencing technology (DST), we profiled approximately 78,000 bases encoding 512 complete exons in the targeted genes in cfDNA from plasma. Seventy-five patients were prospectively enrolled between February 2013 and March 2014, including 61 metastatic cancer patients and 14 clinical stage II CRC patients with matched plasma and tissue samples. Using the 54-gene panel, we detected at least one somatic mutation in 44 of 61 tDNA (72.1%) and 29 of 44 (65.9%) cfDNA. The overall concordance rate of cfDNA to tDNA was 85.9%, when all detected mutations were considered. We collected serial cfDNAs during cetuximab-based treatment in 2 metastatic KRAS wild-type CRC patients, one with acquired resistance and one with primary resistance. We demonstrate newly emerged KRAS mutation in cfDNA 1.5 months before radiologic progression. Another patient had a newly emerged PIK3CA H1047R mutation on cfDNA analysis at progression during cetuximab/irinotecan chemotherapy with gradual increase in allele frequency from 0.8 to 2.1%. This blinded, prospective study of a cfDNA sequencing showed high concordance to tDNA suggesting that the DST approach may be used as a non-invasive biopsy-free alternative to conventional sequencing using tumor biopsy.

Evolutionary analysis of heterochromatin protein compatibility by interspecies complementation in Saccharomyces.

The genetic bases for species-specific traits are widely sought, but reliable experimental methods with which to identify functionally divergent genes are lacking. In the Saccharomyces genus, interspecies complementation tests can be used to evaluate functional conservation and divergence of biological pathways or networks. Silent information regulator (SIR) proteins in S. bayanus provide an ideal test case for this approach because they show remarkable divergence in sequence and paralog number from those found in the closely related S. cerevisiae. We identified genes required for silencing in S. bayanus using a genetic screen for silencing-defective mutants. Complementation tests in interspecies hybrids identified an evolutionarily conserved Sir-protein-based silencing machinery, as defined by two interspecies complementation groups (SIR2 and SIR3). However, recessive mutations in S. bayanus SIR4 isolated from this screen could not be complemented by S. cerevisiae SIR4, revealing species-specific functional divergence in the Sir4 protein despite conservation of the overall function of the Sir2/3/4 complex. A cladistic complementation series localized the occurrence of functional changes in SIR4 to the S. cerevisiae and S. paradoxus branches of the Saccharomyces phylogeny. Most of this functional divergence mapped to sequence changes in the Sir4 PAD. Finally, a hemizygosity modifier screen in the interspecies hybrids identified additional genes involved in S. bayanus silencing. Thus, interspecies complementation tests can be used to identify (1) mutations in genetically underexplored organisms, (2) loci that have functionally diverged between species, and (3) evolutionary events of functional consequence within a genus.

The Awesome Power of Yeast Evolutionary Genetics: New Genome Sequences and Strain Resources for the Saccharomyces sensu stricto Genus.

High-quality, well-annotated genome sequences and standardized laboratory strains fuel experimental and evolutionary research. We present improved genome sequences of three species of Saccharomyces sensu stricto yeasts: S. bayanus var. uvarum (CBS 7001), S. kudriavzevii (IFO 1802(T) and ZP 591), and S. mikatae (IFO 1815(T)), and describe their comparison to the genomes of S. cerevisiae and S. paradoxus. The new sequences, derived by assembling millions of short DNA sequence reads together with previously published Sanger shotgun reads, have vastly greater long-range continuity and far fewer gaps than the previously available genome sequences. New gene predictions defined a set of 5261 protein-coding orthologs across the five most commonly studied Saccharomyces yeasts, enabling a re-examination of the tempo and mode of yeast gene evolution and improved inferences of species-specific gains and losses. To facilitate experimental investigations, we generated genetically marked, stable haploid strains for all three of these Saccharomyces species. These nearly complete genome sequences and the collection of genetically marked strains provide a valuable toolset for comparative studies of gene function, metabolism, and evolution, and render Saccharomyces sensu stricto the most experimentally tractable model genus. These resources are freely available and accessible through www.SaccharomycesSensuStricto.org.

Interspecies variation reveals a conserved repressor of alpha-specific genes in Saccharomyces yeasts.

The mating-type determination circuit in Saccharomyces yeast serves as a classic paradigm for the genetic control of cell type in all eukaryotes. Using comparative genetics, we discovered a central and conserved, yet previously undetected, component of this genetic circuit: active repression of alpha-specific genes in a cells. Upon inactivation of the SUM1 gene in Saccharomyces bayanus, a close relative of Saccharomyces cerevisiae, a cells acquired mating characteristics of alpha cells and displayed autocrine activation of their mating response pathway. Sum1 protein bound to the promoters of alpha-specific genes, repressing their transcription. In contrast to the standard model, alpha1 was important but not required for alpha-specific gene activation and mating of alpha cells in the absence of Sum1. Neither Sum1 protein expression, nor its association with target promoters was mating-type-regulated. Thus, the alpha1/Mcm1 coactivators did not overcome repression by occluding Sum1 binding to DNA. Surprisingly, the mating-type regulatory function of Sum1 was conserved in S. cerevisiae. We suggest that a comprehensive understanding of some genetic pathways may be best attained through the expanded phenotypic space provided by study of those pathways in multiple related organisms.