Search

Empirical determination of the limit of detection. Defined oncogenic driver and background SNVs (a), oncogene and tumor suppressor indels (b), and oncogenic fusions (c) were titrated to specified VAFs, and multiple replicates were measured to determine the 95% LOD. (d) Cell lines comprising defined gene amplifications were titrated to bracket the expected LOD, and the 95% LOD was determined from statistical detection models as per the CLSI Classical Approach.

Alterations detectable by Guardant360 gene panel

Supplemental Methods, Supplemental References, Supplemental Figure Legends, Figure S1-13, Table S1-5 and S14

Supplementary legend

Digital Sequencing analytical performance metrics summary table.

cfDNA variants associated with resistance to on-label targeted therapies across six cancer types

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.

Supplemental Table 1

Purpose:PARP inhibitors (PARPi) are efficacious in multiple cancers harboring germline (and possibly somatic) BRCA1/2 mutations. Acquired reversions can restore BRCA1/2 function, causing resistance to PARPi and/or platinum-based chemotherapy. The optimal method of identifying patients with germline, somatic, and/or reversion mutations in BRCA1/2 has not been established. Next-generation sequencing (NGS) of cell-free DNA (cfDNA) provides a platform to identify these three types of BRCA1/2 mutations.Experimental Design:Patients with advanced breast, ovarian, prostate, or pancreatic cancer were tested using a clinically validated 73-gene cfDNA assay that evaluates single-nucleotide variants and insertion–deletion mutations (indels) in BRCA1/2, and distinguishes somatic/reversion from germline mutations with high accuracy.Results:Among 828 patients, one or more deleterious BRCA1/2 mutations were detected in 60 (7.2%) patients, including germline (n = 42) and somatic (n = 18) mutations. Common coexisting mutations included TP53 (61.6%), MYC (30%), PIK3CA (26.6%), BRAF (15%), and ESR1 (11.5%). Polyclonal reversion mutations (median, 5) were detected in 9 of 42 (21.4%) germline BRCA1/2-mutant patients, the majority (77.7%) of whom had prior PARPi exposure (median duration, 10 months). Serial cfDNA demonstrated emergence of reversion BRCA mutations under therapeutic pressure from initial PARPi exposure, which contributed to subsequent resistance to PARPi and platinum therapy.Conclusions:cfDNA NGS identified high rates of therapeutically relevant mutations without foreknowledge of germline or tissue-based testing results, including deleterious somatic BRCA1/2 mutations missed by germline testing and reversion mutations that can have important treatment implications. Further research is needed to confirm clinical utility of these findings to guide precision medicine approaches for patients with advanced malignancies.

Supplementary Figure Legends

Supplemental Figure 1: Case of advanced NSCLC with high level EGFR T790M mutation on plasma next-generation sequencing (NGS). Supplemental Figure 2: Distribution of variant AFs between 25% and 75% across 105 cases of plasma NGS. Supplemental Figure 3: Plasma NGS cases positive for EGFR T790M submitted for blinded validation.

Purpose: Plasma cell-free DNA (cfDNA) analysis is increasingly used clinically for cancer genotyping, but may lead to incidental identification of germline-risk alleles. We studied EGFR T790M mutations in non–small cell lung cancer (NSCLC) toward the aim of discriminating germline and cancer-derived variants within cfDNA.Experimental Design: Patients with EGFR-mutant NSCLC, some with known germline EGFR T790M, underwent plasma genotyping. Separately, deidentified genomic data and buffy coat specimens from a clinical plasma next-generation sequencing (NGS) laboratory were reviewed and tested.Results: In patients with germline T790M mutations, the T790M allelic fraction (AF) in cfDNA approximates 50%, higher than that of EGFR driver mutations. Review of plasma NGS results reveals three groups of variants: a low-AF tumor group, a heterozygous group (∼50% AF), and a homozygous group (∼100% AF). As the EGFR driver mutation AF increases, the distribution of the heterozygous group changes, suggesting increased copy number variation from increased tumor content. Excluding cases with high copy number variation, mutations can be differentiated into somatic variants and incidentally identified germline variants. We then developed a bioinformatic algorithm to distinguish germline and somatic mutations; blinded validation in 21 cases confirmed a 100% positive predictive value for predicting germline T790M. Querying a database of 31,414 patients with plasma NGS, we identified 48 with germline T790M, 43 with nonsquamous NSCLC (P < 0.0001).Conclusions: With appropriate bioinformatics, plasma genotyping can accurately predict the presence of incidentally detected germline risk alleles. This finding in patients indicates a need for genetic counseling and confirmatory germline testing. Clin Cancer Res; 23(23); 7351–9. ©2017 AACR.

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.

Mutual exclusivity statistics for cfDNA alterations in lung adenocarcinoma, breast cancer, and colorectal cancer

Purpose:Molecular mechanisms of acquired resistance to MET tyrosine kinase inhibitors (TKI) are poorly understood. We aimed to characterize the genomic mechanisms of resistance to type I and type II MET TKIs and their impact on sequential MET TKI therapy outcomes in patients with metastatic MET exon 14–mutant NSCLC.Experimental Design:Genomic alterations occurring at the time of progression on MET TKIs were studied using plasma and tissue next-generation sequencing (NGS).Results:A total of 20 patients had tissue or plasma available for analysis at the time of acquired resistance to a MET TKI. Genomic alterations known or suspected to be mechanisms of resistance were detected in 15 patients (75%). On-target acquired mechanisms of resistance, including single and polyclonal MET kinase domain mutations in codons H1094, G1163, L1195, D1228, Y1230, and high levels of amplification of the MET exon 14–mutant allele, were observed in 7 patients (35%). A number of off-target mechanisms of resistance were detected in 9 patients (45%), including KRAS mutations and amplifications in KRAS, EGFR, HER3, and BRAF; one case displayed both on- and off-target mechanisms of resistance. In 2 patients with on-target resistant mutations, switching between type I and type II MET TKIs resulted in second partial responses.Conclusions:On-target secondary mutations and activation of bypass signaling drive resistance to MET TKIs. A deeper understanding of these molecular mechanisms can support the development of sequential or combinatorial therapeutic strategies to overcome resistance.

Purpose:Molecular mechanisms of acquired resistance to MET tyrosine kinase inhibitors (TKI) are poorly understood. We aimed to characterize the genomic mechanisms of resistance to type I and type II MET TKIs and their impact on sequential MET TKI therapy outcomes in patients with metastatic MET exon 14–mutant NSCLC.Experimental Design:Genomic alterations occurring at the time of progression on MET TKIs were studied using plasma and tissue next-generation sequencing (NGS).Results:A total of 20 patients had tissue or plasma available for analysis at the time of acquired resistance to a MET TKI. Genomic alterations known or suspected to be mechanisms of resistance were detected in 15 patients (75%). On-target acquired mechanisms of resistance, including single and polyclonal MET kinase domain mutations in codons H1094, G1163, L1195, D1228, Y1230, and high levels of amplification of the MET exon 14–mutant allele, were observed in 7 patients (35%). A number of off-target mechanisms of resistance were detected in 9 patients (45%), including KRAS mutations and amplifications in KRAS, EGFR, HER3, and BRAF; one case displayed both on- and off-target mechanisms of resistance. In 2 patients with on-target resistant mutations, switching between type I and type II MET TKIs resulted in second partial responses.Conclusions:On-target secondary mutations and activation of bypass signaling drive resistance to MET TKIs. A deeper understanding of these molecular mechanisms can support the development of sequential or combinatorial therapeutic strategies to overcome resistance.

Purpose: PARP inhibitors (PARPi) are efficacious in multiple cancers harboring germline (and possibly somatic) BRCA1/2 mutations. Acquired reversions can restore BRCA1/2 function, causing resistance to PARPi and/or platinum-based chemotherapy. The optimal method of identifying patients with germline, somatic, and/or reversion mutations in BRCA1/2 has not been established. Next-generation sequencing (NGS) of cell-free DNA (cfDNA) provides a platform to identify these three types of BRCA1/2 mutations. Experimental Design: Patients with advanced breast, ovarian, prostate, or pancreatic cancer were tested using a clinically validated 73-gene cfDNA assay that evaluates single-nucleotide variants and insertion–deletion mutations (indels) in BRCA1/2, and distinguishes somatic/reversion from germline mutations with high accuracy. Results: Among 828 patients, one or more deleterious BRCA1/2 mutations were detected in 60 (7.2%) patients, including germline (n = 42) and somatic (n = 18) mutations. Common coexisting mutations included TP53 (61.6%), MYC (30%), PIK3CA (26.6%), BRAF (15%), and ESR1 (11.5%). Polyclonal reversion mutations (median, 5) were detected in 9 of 42 (21.4%) germline BRCA1/2-mutant patients, the majority (77.7%) of whom had prior PARPi exposure (median duration, 10 months). Serial cfDNA demonstrated emergence of reversion BRCA mutations under therapeutic pressure from initial PARPi exposure, which contributed to subsequent resistance to PARPi and platinum therapy. Conclusions: cfDNA NGS identified high rates of therapeutically relevant mutations without foreknowledge of germline or tissue-based testing results, including deleterious somatic BRCA1/2 mutations missed by germline testing and reversion mutations that can have important treatment implications. Further research is needed to confirm clinical utility of these findings to guide precision medicine approaches for patients with advanced malignancies.

Purpose:To analytically and clinically validate microsatellite instability (MSI) detection using cell-free DNA (cfDNA) sequencing.Experimental Design:Pan-cancer MSI detection using Guardant360 was analytically validated according to established guidelines and clinically validated using 1,145 cfDNA samples for which tissue MSI status based on standard-of-care tissue testing was available. The landscape of cfDNA-based MSI across solid tumor types was investigated in a cohort of 28,459 clinical plasma samples. Clinical outcomes for 16 patients with cfDNA MSI-H gastric cancer treated with immunotherapy were evaluated.Results:cfDNA MSI evaluation was shown to have high specificity, precision, and sensitivity, with a limit of detection of 0.1% tumor content. In evaluable patients, cfDNA testing accurately detected 87% (71/82) of tissue MSI-H and 99.5% of tissue microsatellite stable (863/867) for an overall accuracy of 98.4% (934/949) and a positive predictive value of 95% (71/75). Concordance of cfDNA MSI with tissue PCR and next-generation sequencing was significantly higher than IHC. Prevalence of cfDNA MSI for major cancer types was consistent with those reported for tissue. Finally, robust clinical activity of immunotherapy treatment was seen in patients with advanced gastric cancer positive for MSI by cfDNA, with 63% (10/16) of patients achieving complete or partial remission with sustained clinical benefit.Conclusions:cfDNA-based MSI detection using Guardant360 is highly concordant with tissue-based testing, enabling highly accurate detection of MSI status concurrent with comprehensive genomic profiling and expanding access to immunotherapy for patients with advanced cancer for whom current testing practices are inadequate.See related commentary by Wang and Ajani, p. 6887

Purpose: RET is an emerging oncogenic target showing promise in phase I/II clinical trials. An understudied aspect of RET-driven cancers is the extent to which co-occurring genomic alterations exist and how they may impact prognosis or therapeutic response. Experimental Design: Somatic activating RET alterations were identified among 32,989 consecutive patients with metastatic solid tumors tested with a clinical cell-free circulating tumor DNA (cfDNA) assay. This comprehensive next-generation sequencing (NGS) assay evaluates single-nucleotide variants, and select indels, fusions, and copy number gains in 68–73 clinically relevant cancer genes. Results: A total of 176 somatic activating RET alterations were detected in 170 patients (143 fusions and 33 missense mutations). Patients had non–small cell lung (NSCLC, n = 125), colorectal (n = 15), breast (n = 8), thyroid (n = 8), or other (n = 14) cancers. Alterations in other oncogenic signaling pathway genes were frequently identified in RET-positive samples and varied by specific RET fusion gene partner. RET fusions involving partners other than KIF5B were enriched for alterations in MAPK pathway genes and other bona fide oncogenic drivers of NSCLC, particularly EGFR. Molecular and clinical data revealed that these variants emerged later in the genomic evolution of the tumor as mechanisms of resistance to EGFR tyrosine kinase inhibitors. Conclusions: In the largest cancer cohort with somatic activating RET alterations, we describe novel co-occurrences of oncogenic signaling pathway aberrations. We find that KIF5B-RET fusions are highly specific for NSCLC. In our study, only non-KIF5B-RET fusions contributed to anti-EGFR therapy resistance. Knowledge of specific RET fusion gene partner may have clinical significance.

PARP1/2 inhibitors are effective against BRCA2-deficient tumors. The PARP inhibitor (PARPi) olaparib received FDA breakthrough designation for treatment of metastatic castration-resistant prostate cancers (CRPC) carrying mutations in BRCA1/2 or ATM genes. Emergent resistance to PARPi has been associated with tumor-specific BRCA2 mutations that revert the normal open reading frame rescuing homologous recombination. We describe a case of metastatic CRPC with germline BRCA2 mutation with acquired resistance to olaparib related to biallelic BRCA2 reversion mutations of both the germline and somatic loss of function alleles detected by circulating tumor DNA testing. We also summarize a retrospective analysis of 1,534 prostate cancer cases with ctDNA analysis showing a 1.6% incidence of germline BRCA2 mutations. Within the germline BRCA2-positive cases exposed to platinum chemotherapy or PARP inhibition, the prevalence of reversion mutations was 40%. This report documents the frequency of reversion mutations in a large cohort of prostate cancer patients carrying of BRCA mutations. It also shows the potential utility of ctDNA analyses for early detection of reversion mutation driving tumor resistance.

Cell-free DNA (cfDNA) next-generation sequencing has the potential to capture tumor heterogeneity and genomic evolution under treatment pressure in a non-invasive manner. Here, we report the detection of EGFR L792 mutations, a non-covalent mechanism of osimertinib resistance, using Guardant360 cfDNA testing in a patient with metastatic EGFR-mutant non-small cell lung cancer (NSCLC) whose disease progressed on osimertinib. We subsequently analyzed a large cohort of over 1800 additional patient samples harboring an EGFR T790M mutation and identified a concomitant L792 mutation in a total of 22 (1.2%) cases. In vitro functional assays demonstrated that the EGFR L858R/T790M/L792F/H mutations conferred intermediate-level resistance to osimertinib. Further understanding of potential acquired resistance mechanisms to targeted therapy may help inform treatment strategy in EGFR-mutant NSCLC.

Purpose: Plasma cell-free DNA (cfDNA) analysis is increasingly used clinically for cancer genotyping, but may lead to incidental identification of germline-risk alleles. We studied EGFR T790M mutations in non–small cell lung cancer (NSCLC) toward the aim of discriminating germline and cancer-derived variants within cfDNA. Experimental Design: Patients with EGFR-mutant NSCLC, some with known germline EGFR T790M, underwent plasma genotyping. Separately, deidentified genomic data and buffy coat specimens from a clinical plasma next-generation sequencing (NGS) laboratory were reviewed and tested. Results: In patients with germline T790M mutations, the T790M allelic fraction (AF) in cfDNA approximates 50%, higher than that of EGFR driver mutations. Review of plasma NGS results reveals three groups of variants: a low-AF tumor group, a heterozygous group (∼50% AF), and a homozygous group (∼100% AF). As the EGFR driver mutation AF increases, the distribution of the heterozygous group changes, suggesting increased copy number variation from increased tumor content. Excluding cases with high copy number variation, mutations can be differentiated into somatic variants and incidentally identified germline variants. We then developed a bioinformatic algorithm to distinguish germline and somatic mutations; blinded validation in 21 cases confirmed a 100% positive predictive value for predicting germline T790M. Querying a database of 31,414 patients with plasma NGS, we identified 48 with germline T790M, 43 with nonsquamous NSCLC (P < 0.0001). Conclusions: With appropriate bioinformatics, plasma genotyping can accurately predict the presence of incidentally detected germline risk alleles. This finding in patients indicates a need for genetic counseling and confirmatory germline testing. Clin Cancer Res; 23(23); 7351–9. ©2017 AACR.

Background: NTRK rearrangements, though rare in common cancers, are clinically actionable targets with two FDA-approved drugs for pan-cancer indications, larotrectinib and entrectinib. Fusion detection from ctDNA provides an opportunity to facilitate screening for this biomarker, though technical challenges such as heterogeneity in fusion partners remain. To address this, we developed an assembly-based de novo fusion algorithm that does not rely on a fixed set of partner genes and applied it to >18,000 clinical samples to detect NTRK1 fusions. Methods: A cohort of 18,867 patients across multiple cancer types (lung adenocarcinoma, breast cancer, and colorectal adenocarcinoma), plus 276 healthy control samples were previously tested with Guardant360(R), a CLIA-validated 74-gene cfDNA NGS-based assay. The median unique molecule coverage was approximately 3,000 molecules sequenced to 15,000x read depth. Samples were reanalyzed in silico using the de novo fusion algorithm: in brief, reads aligned to candidate fusion breakpoints were assembled into de Bruijn graphs. Resulting contigs were aligned to the reference and analytical filters were applied to achieve high specificity. Results: NTRK1 fusions were detected in 0.13% of patients, a similar prevalence to what has been previously described in tissue (Rosen et al., 2020). No fusions were identified in 276 healthy control samples. When examining intergenic fusions, 53% of fusions included previously characterized partner genes (TPM3, TPR, LMNA, and TP53), while the remaining fusions contained novel partners.The majority of NTRK1 fusion partners (76%) were detected only once, consistent with previous studies. KRAS G12D and BRAF V600E were detected in 17% and 13% of NTRK1 fusion-positive patients, respectively, and occurred predominantly in the colorectal cohort; most variants (60%) were subclonal to the NTRK1 fusion, and have been previously shown to be associated with resistance to TRK inhibitors through downstream activation of the MAPK pathway. The on-target resistance variant NTRK1 G595R was detected in 8% of patients positive for NTRK1 fusions, all of which were subclonal to the fusion. EGFR L858R was detected in 13% of NTRK1 fusion-positive patients, all of which were clonal to the fusion and in the lung cohort, in accordance with previous studies showing NTRK1 fusions as a resistance mechanism to TKIs in EGFR-positive NSCLC (Xia et al., 2020). Conclusions: NTRK1 fusions were detected in cfDNA at a similar prevalence to tissue NGS, demonstrating high sensitivity of plasma-based assays to detect these fusions. NTRK1 fusion partners were diverse, with the majority of partner genes observed only once across the cohort. Both clonal and subclonal NTRK1 rearrangements were detected, affirming that this biomarker can emerge as an oncogenic driver or as a mechanism of resistance. Citation Format: Arielle Yablonovitch, Sante Gnerre, Jennifer Yen, Scott Shell, Elena Helman, Stephen Fairclough, Rebecca J. Nagy, Justin Odegaard, Darya Chudova, AmirAli Talasaz. NTRK1 fusion detection from clinical cfDNA NGS using a de novo fusion caller [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 537.

PARP inhibitors (PARPi) are efficacious in multiple cancers harboring germline (and possibly somatic)Patients with advanced breast, ovarian, prostate, or pancreatic cancer were tested using a clinically validated 73-gene cfDNA assay that evaluates single-nucleotide variants and insertion-deletion mutations (indels) inAmong 828 patients, one or more deleteriouscfDNA NGS identified high rates of therapeutically relevant mutations without foreknowledge of germline or tissue-based testing results, including deleterious somatic

Introduction: ESR1 mutations are biological indicators of endocrine resistance observed in 20-30% of patients with metastatic breast cancer treated with aromatase inhibitor (AI) therapy. A number of ESR1 resistance mutations have been characterized, including single nucleotide variants in codons 380, 463, 536, 537, and 538. However, ESR1 resistance mutations have not been characterized using large-scale, real-world population studies. Here we describe the subclonal landscape of ESR1 resistance in a real-world dataset of patients with advanced breast cancer treated with anti-estrogen therapy. Methods: The GuardantINFORM™ clinical-genomic database links circulating tumor DNA (ctDNA) results (Guardant360®, Guardant Health) with de-identified aggregated clinical encounters including diagnosis, treatments, and real-world outcomes. GuardantINFORM was queried for adult patients, with a confirmed diagnosis of primary breast cancer, with at least one claim for endocrine therapy (tamoxifen and/or an AI), and at least one Guardant360 test completed following AI treatment. Data was reviewed retrospectively to determine the landscape of ESR1 alterations identified post AI-treatment. Results: 6,541 patients met inclusion criteria. The patient cohort was predominantly female (99%) with a median age at first anti-estrogen therapy of 59.7 years (range 25 - 91 years). 2,044 patients were positive for at least one ESR1 mutation (31%) on the post AI treatment Guardant360 test. 1,943 patients (95%) had an ESR1 mutation identified that is known to confer resistance to AI therapy while the remainder had ESR1 mutations of unknown clinical significance. The majority of patients (65.8%) had multiple ESR1 mutations one of which was always a canonical ESR1 resistance mutation (Table 1). Preliminary outcome analysis showed no difference in real-world overall survival in those with a canonical ESR1 mutation with or without an additional co-occurring subclonal ESR1 mutation. Additional analyses to understand the impact of co-occurring ESR1 alterations and time to therapy resistance is ongoing. Conclusions: Uniquely well-characterized clinical-genomic data in a proprietary dataset identified that approximately 30% of patients with advanced breast cancer had somatic ESR1 mutations following AI therapy, consistent with previously published data. The majority of patients had multiple subclonal ESR1 resistance mutations following AI treatment, always with a canonical resistance mutation, which did not impact real-world overall survival. Additional work is needed to explore the contribution of these ESR1 subclones to the time to clinical resistance. Table 1: Spectrum of ESR1 variants observed in eligible patients in GuardantINFORMESR1 mutationTotal prevalence in ESR1 positive patientsPrevalence as sole ESR1 variantPrevalence with additional ESR1 variant(s)D538G53.8%13.0%40.8%Y537S38.2%8.8%29.4%Y537N18.4%1.6%16.8%E380Q13.0%2.6%10.4%L536H5.0%0.4%4.6%S463P1.4%0.3%1.1%Other29.5%25.6%3.9% Citation Format: David Hanna, Jamie Hutchins, Junhua Yu, Enrique Marino, Naveen Kumar, Rajesh Kucharlapati, Rajesh Chavali, Gautam Nayak, Aaron Hardin, Stephen Fairclough, Victoria M Raymond, Kathryn Lang. Real-world clinical-genomic data identifies the ESR1 clonal and subclonal circulating tumor DNA (ctDNA) landscape and provides insight into clinical outcomes [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr PS18-15.

Background: Real-world evidence (RWE), specifically the integration of clinical and molecular data, plays an increasingly important role in elucidating the molecular landscape of advanced solid tumors, at diagnosis and at progression, and furthering precision oncology research. However, reliance on tissue genomic testing results within these RWE databases significantly limits our understanding of tumor evolution as tissue profiles are typically limited to a single time point, generally at diagnosis, and biopsies are not routinely completed at progression. Circulating tumor DNA (ctDNA) analysis at diagnosis and progression is an alternative to serial tissue testing. Studying the real-world data derived from ctDNA testing results can provide additional information, specifically around tumor evolution and mechanisms of innate and acquired resistance in response to intervening therapies (targeted and non-targeted). Methods: Anonymized results from patients consecutively tested with Guardant360®, a CLIA-certified, CAP-accredited, NYSDOH-approved ctDNA test indicated for patients with advanced and metastatic solid tumor cancers (Guardant Health, Redwood City, California) were analyzed. Using a secure and HIPAA-compliant methodology, these results were then linked to a de-identified encounters database containing medical and pharmacy claims to provide a longitudinal view of patients' journeys, including diagnosis, treatments, and real-world outcomes. Results: Of 85,874 patients, the disease groups most highly represented were lung, breast, colorectal, prostate, and pancreas cancers. 48,317 patients (56.3%) had a record of treatment with an antineoplastic treatment (i.e. chemo, hormonal, small molecule, or immuno-therapy). Of patients with a record of therapy, 14,608 (30.2%) had a record of therapy following, but none prior to, the Guardant360 test date, 14,493 (30.0%) had a record of therapy prior to, but not following theGuardant360 test date, and 19,216 (39.8%) had a record of therapy both before and after Guardant360 test date. Conclusions: In a cohort of patients who underwent ctDNA testing and for whom treatment information was available, post-therapy molecular results were available in 69.8%. The availability of genomic information following systemic therapy suggest that this dataset is a viable resource to study drug resistance and tumor evolution in the real-world setting. Citation Format: Naveen Kumar, Rajesh Kucharlapati, Aaron Hardin, Victoria M. Raymond, Stephen Fairclough, Kathryn Lang. Development of a clinical-genomic database to study tumor evolution and molecular biomarkers of drug resistance in a real-world setting [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 4421.

Background: PARP inhibition can cause synthetic lethality and increased therapeutic sensitivity in patients with HRR deficiency (HRD), which can be detected through the molecular profiling of HRR genes. Prostate cancer has a high prevalence of HRD (20-25%, Athie 2019). High failure rates for tissue biopsy in metastatic prostate cancer patients (25-75%) (Ross 2005, Spritzer 2013) pose challenges for HRD profiling, underscoring the need for a non-invasive, ctDNA alternative. Copy number loss, a frequent cause of HRD, is further difficult to call due to signal dilution by cell-free leukocytic DNA. We developed a pipeline that detects loss-of-function SNV/Indels, structural rearrangements, and gene deletions to identify HRD on GuardantOMNITM, a 500-gene liquid biopsy panel. We present its performance across >620 prostate cancer GuardantOMNITM samples. Methods: Samples from 627 prostate cancer patients were processed on GuardantOMNITM RUO, with median unique coverage of 4982 molecules sequenced to 20,000x read depth. Somatic and germline SNVs and small indels were called using the Guardant bioinformatics pipeline (Helman 2018). A novel HRD module was developed to annotate pathogenic SNV/Indels and identify structural rearrangements, gene-level homozygous deletions, loss-of-heterozygosity (LOH) and genome-wide LOH, comprising of a novel CNV (Barbacioru 2019) and de-novo fusion caller (Gnerre, submitted). Loss-of-function variants were analyzed in 24 HRD genes. Results: Pathogenic alterations in HRD genes were called in 260/620 (42%) prostate cancer samples with ctDNA detected: 28% of all samples had a pathogenic somatic or germline SNV/Indel, 20% had a homozygous deletion, 3.4% had a rearrangement involving an HRD gene. The majority of SNV/Indels occurred in BRCA2 (31% of all 158 SNV/Indels) and ATM (22%), similar to tissue (Dahwan 2016), but mutations also occurred across an additional 21 genes, including CDK12 (8%), CHEK2 (5%) and NBN (3.8%). Of prostate patients with a germline BRCA1/2 SNV/Indel and sufficient tumor shedding for LOH detection (max MAF>10%), 10/19 (52%) also had LOH, compared to 86% in tissue (Jonsson 2019). Homozygous deletions were enriched in BRCA2 (6.9% of all samples), ATM (4.4%) and CHEK1 (2.3%). Rearrangements, including fusions and multi-exonic deletions, accounted for 6.5% of inactivating HRD mutations detected. In total, 24% of prostate samples had a biallelic inactivation involving an SNV, Indel or deletion. Conclusion: We demonstrate in a prostate cancer cohort that GuardantOMNITM ctDNA profiling calls all classes of mutations contributing to HRD, with relative prevalence of alterations consistent with those in tissue. CfDNA presents a potential alternative for identifying patients who may benefit from PARP or cisplatin/platinum therapies, expanding the prevalence from 28% using small variants to 42% with the complete HRD biomarker set. Citation Format: Jennifer Yen, Sante Gnerre, Catalin Barbacioru, Elena Helman, Ohad Manor, Arielle Yablonovitch, Ravi Vijaya Satya, Leo Liu, Jennifer Saam, Stephen Fairclough, Becky Nagy, Richard Lanman, Darya Chudova, AmirAli Talasaz. Landscape of homologous recombination repair (HRR) mutations in prostate cancer profiled by ctDNA next-generation sequencing [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 729.

3545 Background: FGFR2/ 3 rearrangements are promising therapeutic targets, especially in advanced urothelial cancer (aUC) with FDA-approved erdafitinib. Liquid biopsy is an attractive non-invasive method to identify these fusions, but detection in cfDNA is technically challenging due to low tumor shedding levels, short molecules, and wide variation in gene partners. To address this, we developed an assembly-based fusion detection algorithm to call rearrangements in a de novo fashion without reliance on a fixed partner set and applied it to > 15,000 clinical samples. Methods: A cohort of 15,218 patients with mixed cancer types (including 698 aUC patients, as well as breast, cholangiocarcinoma, colorectal, and gastric), plus 276 healthy control samples were previously tested with Guardant360(R), a clinical 74-gene cfDNA NGS-based assay. The median unique molecule coverage was approximately 3,000 molecules sequenced to 15,000x read depth. Samples were reanalyzed in silico using the novel algorithm: in brief, reads aligned to candidate fusion breakpoints were assembled into de Bruijn graphs. Resulting contigs were aligned to the reference and filters were applied to remove technical artifacts. Results: The majority of FGFR2 (86%) and FGFR3 fusion partners (73%) in the mixed cancer cohort were observed only once, consistent with previous reports (Helsten 2016). FGFR3-TACC3 was the most common fusion, occurring in 72% of FGFR3 fusion-positive patients. In 37% of FGFR2 fusion positive patients, the de novo caller detected partners not previously described. In the aUC cohort, FGFR3 fusions were detected in 3.3% of patients, with 8/10 (80%) partner genes/intergenic regions occurring only once, which is in line with previous reports (Nassar 2018). No fusions were identified in 276 healthy control samples. In the mixed cancer cohort, common mutations co-occurring with FGFR2 fusions were FGFR2 N549K, PIK3CA H1047R, and TP53 R175H (5.6% each); KRAS Q61H was observed in 28% of patients with FGFR3 fusions. Conclusions: FGFR2/ 3 fusion partners detected by a highly specific assembly-based de novo fusion caller were heterogeneous and individually rare, highlighting the importance of a de novo approach. We observed an FGFR3 fusion prevalence in cfDNA from aUC patients that is comparable to previous reports for tissue testing, demonstrating an ability to capture targetable genomic rearrangements with plasma-based NGS in this patient population.

Targeting HER2 is effective in 24% of ERBB2 amplified metastatic colorectal cancer; however, secondary resistance occurs in most of the cases. We studied the evolution of individual metastases during treatment to discover spatially resolved determinants of resistance. Circulating tumor DNA (ctDNA) analysis identified alterations associated with resistance in the majority of refractory patients. ctDNA profiles and lesion-specific radiographic reports revealed organ- or metastasis-private evolutionary patterns. When radiologic assessments documented progressive disease in target lesions, response to HER2 blockade was retained in other metastases. Genomic and functional analyses on samples and cell models from eight metastases of a patient co-recruited to a postmortem study unveiled lesion-specific evolutionary trees and pharmacologic vulnerabilities. Lesion size and contribution of distinct metastases to plasma ctDNA were correlated.

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.

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.

Fixed chromosomal inversions can reduce gene flow and promote speciation in two ways: by suppressing recombination and by carrying locally favoured alleles at multiple loci. However, it is unknown whether favoured mutations slowly accumulate on older inversions or if young inversions spread because they capture pre-existing adaptive quantitative trait loci (QTLs). By genetic mapping, chromosome painting and genome sequencing, we have identified a major inversion controlling ecologically important traits in Boechera stricta. The inversion arose since the last glaciation and subsequently reached local high frequency in a hybrid speciation zone. Furthermore, the inversion shows signs of positive directional selection. To test whether the inversion could have captured existing, linked QTLs, we crossed standard, collinear haplotypes from the hybrid zone and found multiple linked phenology QTLs within the inversion region. These findings provide the first direct evidence that linked, locally adapted QTLs may be captured by young inversions during incipient speciation.Chromosome inversions play an important role in local adaptation and speciation1,2, and selectively important inversions have been identified in many species3,4. Selection due to different environmental factors or stages in the life cycle1 may favour inversions carrying locally adapted alleles at several loci. In addition, established inversions are predicted to accumulate selectively important genetic differences, which may contribute to reproductive isolation during speciation1.Although few studies have identified the actual loci that influence selection on inversions2,4,5, rearrangements may be favoured due to gene alterations near breakpoints6, chromatin changes7 or combinations of advantageous, co-adapted alleles8. Inversions suppress recombination, so locally advantageous alleles may segregate together, causing higher fitness than recombinant haplotypes9. Most evolutionary studies have focused on widespread, older inversions, so we have little knowledge of the evolutionary processes that guide their initial increase in frequency. For example, do inversions drift to higher frequency, and then acquire new advantageous mutations after they are common? Or are multiple linked, advantageous alleles captured in a new inversion, allowing them to spread together? Analysis of younger inversions may elucidate the evolutionary forces controlling the initial spread of chromosome inversions, which therefore influence their role in adaptation and speciation4,8.Related species often differ for chromosome inversions that carry locally favoured alleles at multiple loci10,11. A key distinction among models for the evolution of inversions is whether early frequency increase is due to genetic drift or natural selection. Genetic drift might predominate initially, with subsequent accumulation of advantageous variants12. Alternatively, the Kirkpatrick–Barton model9 argues that linked, locally adapted alleles exist first, and subsequently are captured within a new, selectively favoured inversion13. In this ‘inversion-late’ evolutionary sequence1,5, linked quantitative trait loci (QTLs), similar to the ancestral haplotype that gave rise to the inversion, may still exist in non-inverted genotypes9. Here, we test these predictions of the Kirkpatrick–Barton model. First, we introduce ecologically diverged subspecies of Boechera stricta. Next, we examine a young inversion to infer the selective forces controlling its early increase in frequency. Finally, we cross collinear, standard genotypes from the hybrid zone to ask whether old, linked QTLs can be found within the inversion region.

Background: In addition to BRCA1/2 germline and somatic inactivating mutations, loss of heterozygosity (LOH) and biallellic somatic copy number loss are associated with BRCA1/2 loss of function. Patients with cancers harboring these somatic features may benefit from treatment with PARP inhibitors. ctDNA analysis has proven a viable alternative to tumor tissue genotyping, especially in tissue- or time-limited clinical scenarios. However, ctDNA assessment of LOH and biallelic copy number loss is challenging given that ctDNA is diluted by cell-free leukocytic DNA making it difficult to confidently call these genomic events. We developed a method to identify somatic biallelic copy number loss in ctDNA using targeted sequencing of cell-free (cfDNA) across a wide range of cancer types. Methods: A novel statistical model was developed using coverage profiles and single nucleotide polymorphism (SNP) allelic frequencies estimated from plasma samples to determine the presence of LOH and biallellic copy number loss in BRCA1 or BRCA2. For each gene of interest, the model uses observed coverage and fragment size distribution, together with allele frequency of germline SNPs. The model was applied to plasma samples from 28,000 patients with advanced solid tumors sequenced using a 73-gene next generation sequencing ctDNA panel (Guardant360®, Guardant Health, Redwood City, CA). Results: The model was analytically validated using in-silico simulations in order to assess both the limit of blank and limit of detection. This method shows 95% sensitivity in detecting LOH and bi-allellic copy number loss for samples with a maximum somatic variant allele frequency as low as 9%. Sensitivity is mainly driven by the panel size and the depth of coverage of the gene of interest. The model was then applied to the 28,199 patient cohort. BRCA1 and BRCA2 LOH was observed in 2.4% (134/5568) and 7.4% (415/5568) of classic homologous recombination deficient (HRD) cancers including breast, ovarian, prostate, and pancreas. BRCA1 and BRCA2biallelic copy number loss was observed in 0.3% (19/5568) and 0.6% (31/5568) of this same group of HRD cancers. Discussion: In this cohort of 5,568 patients with classic HRD associated cancers, somatic LOH and biallelic copy number loss was detected in BRCA1 in 2.7% of samples and in BRCA2 in 8.0% of samples, which is aligned with previously reported tissue prevalence. BRCA1/2 somatic LOH and biallelic copy number loss can accurately be detected in ctDNA utilizing likelihood-based scores based on coverage and germline SNPs allele frequencies. The ability to identify this therapeutically targetable genomic alteration through a non-invasive ctDNA assessment has significant clinical implications. Citation Format: Catalin Barbacioru, Victoria M. Raymond, Marcin Sikora, Elena Helman, Sante Gnerre, Stephen Fairclough, Darya Chudova, Richard B. Lanman, AmirAli Talasaz. Cell-free circulating tumor DNA (ctDNA) detects somatic copy number loss in homologous recombination repair genes [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 435.

Background: Late stage therapy selection is advancing beyond the detection of classical driver alterations (SNV, Indel, CNV, Fusion) to the use of aggregate biomarkers comprising combinations of SNVs, Indels, and fusions, that reflect the underlying mutational process, indicate the cellular dysfunction, and suggest the neoantigen signature. GuardantOMNI (OMNI), is a highly sensitive 500-gene cfDNA sequencing test designed to detect primary somatic mutations as well as aggregate biomarkers including TMB and DDR pathway defects. The newest version of OMNI introduces MSI biomarker detection, which is an indication for the first FDA-approved pan-cancer therapy. Here we present the first analytical validation study of OMNI MSI-H detection. Methods: Probabilistic models for MSI detection were designed to interrogate instability across >1300 loci, differentiating PCR and sequencing induced artifacts from true biological signals using molecular barcoding as implemented in the Guardant digital sequencing platform. Cut-offs were selected based on the analysis of the training healthy donor cohort and locked prior to validation. Analytical performance was assessed using pre-characterized cell lines, cancer-free donor-samples, and clinical patient samples. Qualitative and quantitative orthogonal confirmation was provided by next-generation sequencing (NGS) and published information. Results: 350 samples were processed for this study, using both 5ng and 30ng cfDNA input levels, and all samples passed sequencing QC metrics established prior to testing. Cell line-based dilution studies performed in triplicate spanning 5ng and 30ng input range as well as in silico simulation demonstrated 95% limit of detection (LoD) at 0.1% Analysis of 121 healthy donor and 16 MSS advanced cancer patient samples identified no false positives (FP). In 240 patient samples with orthogonal confirmation of MSI status from tissue or panel NGS, GuardantOMNI demonstrated 100% analytic sensitivity (3/3) and 100% analytic specificity (237/237) with a PPV of 1 and and NPV 1. Retrospective analysis of 2000+ clinical samples detected MSI high in 16% of colorectal cancers, 4% of prostate cancers, and 0.6% of non-small cell lung cancers. Conclusions: Here we have demonstrated the performance of MSI-H detection with GuardantOMNI (OMNI), which has the potential to accelerate clinical trial enrollment, research and discovery with a single, non-invasive blood sample. Citation Format: Jing Zhao, Alexander Artyomenko, Carlo Artieri, John Latham, Stephen R. Fairclough, Catalin Barbacioru, Elena Helman, John Strickler, Darya Chudova, Richard Lanman, AmirAli Talasaz. Analytical validation of MSI High detection with GuardantOMNI [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 1675.

The evolution of cadherins, which are essential for metazoan multicellularity and restricted to metazoans and their closest relatives, has special relevance for understanding metazoan origins. To reconstruct the ancestry and evolution of cadherin gene families, we analyzed the genomes of the choanoflagellate Salpingoeca rosetta , the unicellular outgroup of choanoflagellates and metazoans Capsaspora owczarzaki , and a draft genome assembly from the homoscleromorph sponge Oscarella carmela . Our finding of a cadherin gene in C. owczarzaki reveals that cadherins predate the divergence of the C. owczarzaki , choanoflagellate, and metazoan lineages. Data from these analyses also suggest that the last common ancestor of metazoans and choanoflagellates contained representatives of at least three cadherin families, lefftyrin, coherin, and hedgling. Additionally, we find that an O. carmela classical cadherin has predicted structural features that, in bilaterian classical cadherins, facilitate binding to the cytoplasmic protein β-catenin and, thereby, promote cadherin-mediated cell adhesion. In contrast with premetazoan cadherin families (i.e., those conserved between choanoflagellates and metazoans), the later appearance of classical cadherins coincides with metazoan origins.

Background: Cell-free circulating tumor DNA analysis provides a non-invasive method for obtaining actionable genomic information to guide personalized cancer treatment. Deep sequencing of cell-free DNA (cfDNA) can potentially provide insights into tumor heterogeneity across multiple tumor sites in a patient, including emerging treatment-resistant subclones. However, the increased informational complexity of polyclonal cfDNA in circulation poses analysis challenges, particularly in tumors with abundant copy number alterations. To facilitate interpretation of this added complexity, we developed methods to identify cfDNA copy-number driver alterations and cfDNA clonality, Methods: We analyzed a large clinical sequencing database of somatic point mutations and copy number alterations from targeted cfDNA sequencing of 21,807 consecutive patients across >50 cancer types (Guardant Health, CA). We evaluated a minimal cfDNA clonality model that relies on the relationship between variant allele frequency (VAF) in cfDNA and the level of tumor DNA in circulation (ctDNA level), while accounting for copy number alterations. Results: We found that the initial simple model of cfDNA clonality performed well on >90% of samples, given a relatively small targeted genomic region (70 genes, 150 kb). However, normalizing VAF by copy number is subject to error in some samples due to the effect of ctDNA level on variant detection, variable unique molecule coverage across samples, and non-linearity of VAF at high copy number. Therefore, we developed an improved cfDNA clonality model that incorporated these analytical and biological features, which was then trained on a portion of the large cfDNA data set. Our cfDNA clonality model accurately distinguished subclonal resistance from driver alterations in a test set of over 5,000 lung, colorectal, and breast cancer patients. Although numerous subclonal tumor-derived alterations were apparent in the initial test data set, leading to an apparent departure from mutual exclusivity in treatment-naïve tumors, robust mutual exclusivity was observed among cfDNA clonal driver alterations when our cfDNA clonality analysis method was applied. These results suggest our analytical approach can be used to identify treatment-associated emerging resistance alterations in patients from a single blood draw, including parallel evolution of distinct subclonal alterations. Conclusion: Managing cancer will likely depend on identifying emerging treatment-resistant subclones at or in anticipation of progression. Highly accurate deep sequencing of cfDNA, along with comprehensive models of cfDNA clonality, can elucidate subclonal structure of the tumor and identify emerging treatment resistance. Citation Format: Stephen Fairclough, Oliver Zill, Catalin Barbacioru, Justin Odegaard, Richard B. Lanman, AmirAli Talasaz, Darya Chudova. A method for differentiating clonal driver mutations from subclonal emerging resistance mutations in circulating cell-free DNA [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 2190.

12097Background: PARP inhibitors (PARPi) have shown efficacy in breast (BC), ovarian (OC), pancreatic (PaC), and prostate (PrC) cancers harboring deleterious germline or somatic alterations (alts) ...

For more than a century and half choanoflagellates, the closest living relatives of animals, have fascinated evolutionary biologists. By characterizing the similarities and differences between choanoflagellates and animals, biologists have gained perspective on the biology of their last common ancestor, the “Urchoanimal”, as well as the evolutionary foundations of multicellularity and the origin of animals. The best-studied colonial choanoflagellate, Salpingoeca rosetta, forms colonies by cell division and not by cell aggregation. The observation that cytoplasmic bridges connect cells in S. rosetta colonies and other colonial choanoflagellates, as well as cells in sponges, suggests that this mechanism of colony formation may be ancestral within the choanoflagellate lineage and may have been present in the Urchoanimal as well. The comparison of choanoflagellate gene content and gene function with animal gene content and gene function has revealed that many of the basic mechanisms of cell adhesion, signaling, and differentiation that were previously thought to be unique to animals are also present in choanoflagellates, indicating that these genes were present prior to the evolution of animals. These insights refine our understanding of genes that emerged on the stem lineage leading to the last common ancestor of all animals, the “Urmetazoan”. Taken together the data from choanoflagellates have provided deep insights into the biology of the Urchoanimal and the evolution of animal multicellularity.

Fixed chromosomal inversions can reduce gene flow and promote speciation in two ways: by suppressing recombination and by carrying locally favored alleles at multiple loci. However, it is unknown whether favored mutations slowly accumulate on older inversions or if young inversions spread because they capture preexisting adaptive Quantitative Trait Loci (QTLs). By genetic mapping, chromosome painting and genome sequencing we have identified a major inversion controlling ecologically important traits in Boechera stricta. The inversion arose since the last glaciation and subsequently reached local high frequency in a hybrid speciation zone. Furthermore, the inversion shows signs of positive directional selection. To test whether the inversion could have captured existing, linked QTLs, we crossed standard, collinear haplotypes from the hybrid zone and found multiple linked phenology QTLs within the inversion region. These findings provide the first direct evidence that linked, locally adapted QTLs may be captured by young inversions during incipient speciation.

Guardant360 is a cell-free circulating tumor DNA (ctDNA) test that genotypes all guideline-recommended solid tumor somatic genomic treatment targets from a single non-invasive blood draw. The new version, v2.10, was redesigned to enhance sensitivity and specificity across 73 cancer-related genes. It detects all four major variant classes (single nucleotide variants, SNVs, in all 73 genes; indels in 23 genes; gene amplifications, CNAs, in 18 genes; and fusions in 6 genes). Analytical performance was assessed throughout the reportable range via multiple serial dilution studies of orthogonally-characterized contrived and patient samples. Analytical specificity was assessed by calculating the false positive rate in pre-characterized healthy donor sample mixtures serially diluted. Positive predictive value (PPV) was estimated as a function of allelic fraction/copy number from pre-characterized samples and prevalence-adjusted using a cohort of 2,585 consecutive clinical samples. Confirmation was performed using ddPCR. Analytical specificity was 100% for SNVs, fusions, and CNAs and 96% for indels across 25 defined samples. Relative to Guardant360v2.9, v2.10 demonstrated 20-50% increase in fusion molecule recovery. Retrospective in silico analysis of 2,585 consecutive clinical samples demonstrated a 15% increase in actionable fusion detection, a 6%-15% increase in actionable indel detection (excluding newly reportable indels), and a 3%-6% increase in actionable SNV detection. AlterationsReportable Range95% Limit of DetectionAllelic Fraction / Copy NumberAnalytical SensitivityAllelic Fraction / Copy numberPPVSNVs≥0.04%0.25%≥0.25%>99.9%≥0.25%98.7%0.05–0.25%63.8%99.9%≥0.25%98.4%0.05–0.25%67.8% Guardant360 analytical performance characteristics based on standard cfDNA input (30ng). Analytical sensitivity/limit of detection estimates are provided for clinically actionable variants and may vary by sequence context and cfDNA input. PPV is estimated across entire reportable panel space (PPV was 100% for clinically actionable variants). Conclusion: Guardant360 v2.10 comprehensively detects all adult solid tumor guideline-recommended somatic genomic variants with unparalleled sensitivity, accuracy, and specificity. Citation Format: James Vowles, Justin Odegaard, Stefanie Mortimer, Stephen Fairclough, Marcin Sikora, Diana Abdueva, Reza Mokhtari, Arthur Baca, AmirAli Talasaz. Analytical validation of Guardant360 v2.10 [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 5705. doi:10.1158/1538-7445.AM2017-5705

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License, which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited., Filopodia are fine actin-based cellular projections used for both environmental sensing and cell motility, and they are essential organelles for metazoan cells. In this study, we reconstruct the origin of metazoan filopodia and microvilli. We first report on the evolutionary assembly of the filopodial molecular toolkit and show that homologs of many metazoan filopodial components, including fascin and myosin X, were already present in the unicellular or colonial progenitors of metazoans. Furthermore, we find that the actin crosslinking protein fascin localizes to filopodia-like structures and microvilli in the choanoflagellate Salpingoeca rosetta. In addition, homologs of filopodial genes in the holozoan Capsaspora owczarzaki are upregulated in filopodia-bearing cells relative to those that lack them. Therefore, our findings suggest that proteins essential for metazoan filopodia and microvilli are functionally conserved in unicellular and colonial holozoans and that the last common ancestor of metazoans bore a complex and specific filopodial machinery. © 2013 The Author., This work was supported by an Institució Catalana de Recerca i Estudis Avançats contract, a European Research Council Starting Grant (ERC-2007-StG- 206883), a grant (BFU2011-23434) from Ministerio de Economía y Competitividad (MINECO) to I.R.-T., funding from NIH NIGMS GM089977 to N.K., pregraduate Formacion Profesorado Universitario grant from MICINN to A.S-P., a Deutsche Forschungsgemeinschaft (DFG) postdoctoral fellowship to P.B., and the Canadian Research Chair program grant to B.F.L.

Background: Metazoan multicellularity is rooted in mechanisms of cell adhesion, signaling, and differentiation that first evolved in the progenitors of metazoans. To reconstruct the genome composition of metazoan ancestors, we sequenced the genome and transcriptome of the choanoflagellate Salpingoeca rosetta, a close relative of metazoans that forms rosette-shaped colonies of cells. Results: A comparison of the 55 Mb S. rosetta genome with genomes from diverse opisthokonts suggests that the origin of metazoans was preceded by a period of dynamic gene gain and loss. The S. rosetta genome encodes homologs of cell adhesion, neuropeptide, and glycosphingolipid metabolism genes previously found only in metazoans and expands the repertoire of genes inferred to have been present in the progenitors of metazoans and choanoflagellates. Transcriptome analysis revealed that all four S. rosetta septins are upregulated in colonies relative to single cells, suggesting that these conserved cytokinesis proteins may regulate incomplete cytokinesis during colony development. Furthermore, genes shared exclusively by metazoans and choanoflagellates were disproportionately upregulated in colonies and the single cells from which they develop. Conclusions: The S. rosetta genome sequence refines the catalog of metazoan-specific genes while also extending the evolutionary history of certain gene families that are central to metazoan biology. Transcriptome data suggest that conserved cytokinesis genes, including septins, may contribute to S. rosetta colony formation and indicate that the initiation of colony development may preferentially draw upon genes shared with metazoans, while later stages of colony maturation are likely regulated by genes unique to S. rosetta. Background Metazoan multicellularity and development are rooted in basic mechanisms of cell adhesion, signaling, and differentiation that were present in the unicellular and colonial progenitors of metazoans. Reconstructing the evolution of metazoans from their single celled ancestors promises to illuminate one of the major transitions in evolutionary history, while also revealing fundamental mechanisms underlying metazoan cell biology and

Bacterially-produced small molecules exert profound influences on animal health, morphogenesis, and evolution through poorly understood mechanisms. In one of the closest living relatives of animals, the choanoflagellate Salpingoeca rosetta, we find that rosette colony development is induced by the prey bacterium Algoriphagus machipongonensis and its close relatives in the Bacteroidetes phylum. Here we show that a rosette inducing factor (RIF-1) produced by A. machipongonensis belongs to the small class of sulfonolipids, obscure relatives of the better known sphingolipids that play important roles in signal transmission in plants, animals, and fungi. RIF-1 has extraordinary potency (femtomolar, or 10(-15) M) and S. rosetta can respond to it over a broad dynamic range-nine orders of magnitude. This study provides a prototypical example of bacterial sulfonolipids triggering eukaryotic morphogenesis and suggests molecular mechanisms through which bacteria may have contributed to the evolution of animals.DOI:http://dx.doi.org/10.7554/eLife.00013.001.

Analysis of cell-free circulating tumor DNA (ctDNA) by next-generation sequencing (NGS) allows non-invasive real-time profiling of actionable genomic alterations. Liquid biopsy provides an option for disease monitoring in early stage cancer patients post surgical resection, with a potential to aid in adjuvant decision making. However, to be applicable, tests must cover a broad enough genomic footprint to not require a priori knowledge of mutations, have high specificity, and sensitivity higher than conventional methods. NGS is necessary, since inactivating mutations are the most common alteration type in many common cancer types such as colorectal carcinomas (CRC). Here we present a highly efficient and specific NGS assay for detection of ctDNA in early stage cancer patients, capable of detecting single molecule variants across a 12 kb gene panel with an analytical sensitivity of >0.02% for single nucleotide variants (SNVs) and indels. This panel was applied to a clinical study involving 14 early stage (II/III) CRC patients with both pre- and post-op blood draws (up to 7 days post surgery). A subset (6 patients) also had tumor samples collected at the time of the surgical resection of the tumor. Overall, the detection rate of ctDNA in pre-op blood draws was 93%. In the post-op blood draws ctDNA was detectable in 43% of cases. The estimated average minor allele frequency (MAF) is 0.58% (± 0.82%) in pre-op, 0.18% (±0.21%) in post-op, and 40% (±18%) in tumor samples. When tumor tissue was available and used as a reference, the clinical sensitivity, specificity, and accuracy in pre-op blood samples were 83%, 99.995%, and 99.99%, respectively. SNVs with MAF as low as 0.04% were confirmed in tissue data. The clinical specificity of variants detected in post-op blood samples using pre-op samples as the reference is 99.996%. Cohort expansion to 50 patients and follow-up for clinical recurrence in both cohorts is ongoing. In conclusion, we have developed an assay with ultra-high accuracy and specificity, for the detection of ctDNA in early stage CRC patients that is capable of detecting alterations present in the tumor post-surgical resection. This technology allows for a promising non-invasive route for molecular monitoring of residual disease post surgery and for early detection of relapse compared to traditional methodologies. Citation Format: Stefanie A. Mortimer, Katharine Dilger, Stephen Fairclough, Diana Abdueva, Darya Chudova, Ankit Sarin, Jim Leng, Jeeyun Lee, Helmy Eltoukhy, AmirAli Talasaz. Post-surgical resection monitoring in early stage colorectal carcinoma patients using a circulating cell-free DNA assay with ultra-high accuracy and specificity. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 506.

e20643Background: Somatic mutations (SNVs) in ALK are a common mechanism of resistance to crizotinib, and may respond to next-generation ALK inhibitors. Historically, identification of these ALK SN...

e23021Background: In non-small cell lung cancer (NSCLC), targeted treatment with first line tyrosine kinase inhibitors (TKIs) selects for the evolution of resistance at EGFR T790M in 50% of cases. ...

Bacterially-produced small molecules exert profound influences on animal health, morphogenesis, and evolution through poorly understood mechanisms. In one of the closest living relatives of animals, the choanoflagellate Salpingoeca rosetta, we find that rosette colony development is induced by the prey bacterium Algoriphagus machipongonensis and its close relatives in the Bacteroidetes phylum. Here we show that a rosette inducing factor (RIF-1) produced by A. machipongonensis belongs to the small class of sulfonolipids, obscure relatives of the better known sphingolipids that play important roles in signal transmission in plants, animals, and fungi. RIF-1 has extraordinary potency (femtomolar, or 10−15 M) and S. rosetta can respond to it over a broad dynamic range—nine orders of magnitude. This study provides a prototypical example of bacterial sulfonolipids triggering eukaryotic morphogenesis and suggests molecular mechanisms through which bacteria may have contributed to the evolution of animals. DOI: http://dx.doi.org/10.7554/eLife.00013.001, eLife digest All animals, including humans, evolved in a world filled with bacteria. Although bacteria are most familiar as pathogens, some bacteria produce small molecules that are essential for the biology of animals and other eukaryotes, although the details of the ways in which these bacterial molecules are beneficial are not well understood. The choanoflagellates are water-dwelling organisms that use their whip-like flagella to move around, feeding on bacteria. They can exist as one cell or a colony of multiple cells and, perhaps surprisingly, are the closest known living relatives of animals. This means that experiments on these organisms have the potential to improve our understanding of animal development and the transition from egg to embryo to adult. Alegado et al. have explored how the morphology of Salpingoeca rosetta, a colony-forming choanoflagellate, is influenced by its interactions with various species of bacteria. In particular, they find that the development of multicellularity in S. rosetta is triggered by the presence of the bacterium Algoriphagus machipongonensis as well as its close relatives. They also identify the signaling molecule produced by the bacteria to be C32H64NO7S; this lipid molecule is an obscure relative of the sphingolipid molecules that have important roles in signal transmission in animals, plants, and fungi. Moreover, Alegado et al. show that S. rosetta can respond to this molecule – which they call rosette-inducing factor (RIF-1) – over a wide range of concentrations, including concentrations as low as 10−17 M. The work of Alegado et al. suggests that interactions between S. rosetta and Algoriphagus bacteria could be a productive model system for studying the influences of bacteria on animal cell biology, and for investigating the mechanisms of signal delivery and reception. Moreover, the molecular mechanisms revealed by this work leave open the possibility that bacteria might have contributed to the evolution of multicellularity in animals. DOI: http://dx.doi.org/10.7554/eLife.00013.002