Viktor A. Adalsteinsson

Whole exome sequencing of circulating tumor cells provides a window into metastatic prostate cancer

Comprehensive analyses of cancer genomes promise to inform prognoses and precise cancer treatments. A major barrier, however, is inaccessibility of metastatic tissue. A potential solution is to characterize circulating tumor cells (CTCs), but this requires overcoming the challenges of isolating rare cells and sequencing low-input material. Here we report an integrated process to isolate, qualify and sequence whole exomes of CTCs with high fidelity using a census-based sequencing strategy. Power calculations suggest that mapping of >99.995% of the standard exome is possible in CTCs. We validated our process in two patients with prostate cancer, including one for whom we sequenced CTCs, a lymph node metastasis and nine cores of the primary tumor. Fifty-one of 73 CTC mutations (70%) were present in matched tissue. Moreover, we identified 10 early trunk and 56 metastatic trunk mutations in the non-CTC tumor samples and found 90% and 73% of these mutations, respectively, in CTC exomes. This study establishes a foundation for CTC genomics in the clinic.

EGFR variant heterogeneity in glioblastoma resolved through single-nucleus sequencing

UNLABELLED Glioblastomas (GBM) with EGFR amplification represent approximately 50% of newly diagnosed cases, and recent studies have revealed frequent coexistence of multiple EGFR aberrations within the same tumor, which has implications for mutation cooperation and treatment resistance. However, bulk tumor sequencing studies cannot resolve the patterns of how the multiple EGFR aberrations coexist with other mutations within single tumor cells. Here, we applied a population-based single-cell whole-genome sequencing methodology to characterize genomic heterogeneity in EGFR-amplified glioblastomas. Our analysis effectively identified clonal events, including a novel translocation of a super enhancer to the TERT promoter, as well as subclonal LOH and multiple EGFR mutational variants within tumors. Correlating the EGFR mutations onto the cellular hierarchy revealed that EGFR truncation variants (EGFRvII and EGFR carboxyl-terminal deletions) identified in the bulk tumor segregate into nonoverlapping subclonal populations. In vitro and in vivo functional studies show that EGFRvII is oncogenic and sensitive to EGFR inhibitors currently in clinical trials. Thus, the association between diverse activating mutations in EGFR and other subclonal mutations within a single tumor supports an intrinsic mechanism for proliferative and clonal diversification with broad implications in resistance to treatment. SIGNIFICANCE We developed a novel single-cell sequencing methodology capable of identifying unique, nonoverlapping subclonal alterations from archived frozen clinical specimens. Using GBM as an example, we validated our method to successfully define tumor cell subpopulations containing distinct genetic and treatment resistance profiles and potentially mutually cooperative combinations of alterations in EGFR and other genes.

Resistance to checkpoint blockade therapy through inactivation of antigen presentation

Treatment with immune checkpoint blockade (CPB) therapies often leads to prolonged responses in patients with metastatic melanoma, but the common mechanisms of primary and acquired resistance to these agents remain incompletely characterized and have yet to be validated in large cohorts. By analyzing longitudinal tumor biopsies from 17 metastatic melanoma patients treated with CPB therapies, we observed point mutations, deletions or loss of heterozygosity (LOH) in beta-2-microglobulin (B2M), an essential component of MHC class I antigen presentation, in 29.4% of patients with progressing disease. In two independent cohorts of melanoma patients treated with anti-CTLA4 and anti-PD1, respectively, we find that B2M LOH is enriched threefold in non-responders (~30%) compared to responders (~10%) and associated with poorer overall survival. Loss of both copies of B2M is found only in non-responders. B2M loss is likely a common mechanism of resistance to therapies targeting CTLA4 or PD1.Resistance to immune-checkpoint blockade often occurs in treated patients. Here, the authors demonstrate that B2M loss is a mechanism of primary and acquired resistance to therapies targeting CTLA4 or PD-1 in melanoma patients.

Scalable whole-exome sequencing of cell-free DNA reveals high concordance with metastatic tumors

Whole-exome sequencing of cell-free DNA (cfDNA) could enable comprehensive profiling of tumors from blood but the genome-wide concordance between cfDNA and tumor biopsies is uncertain. Here we report ichorCNA, software that quantifies tumor content in cfDNA from 0.1× coverage whole-genome sequencing data without prior knowledge of tumor mutations. We apply ichorCNA to 1439 blood samples from 520 patients with metastatic prostate or breast cancers. In the earliest tested sample for each patient, 34% of patients have ≥10% tumor-derived cfDNA, sufficient for standard coverage whole-exome sequencing. Using whole-exome sequencing, we validate the concordance of clonal somatic mutations (88%), copy number alterations (80%), mutational signatures, and neoantigens between cfDNA and matched tumor biopsies from 41 patients with ≥10% cfDNA tumor content. In summary, we provide methods to identify patients eligible for comprehensive cfDNA profiling, revealing its applicability to many patients, and demonstrate high concordance of cfDNA and metastatic tumor whole-exome sequencing.Identifying the mutational landscape of tumours from cell-free DNA in the blood could help diagnostics in cancer. Here, the authors present ichorCNA, software that quantifies tumour content in cell free DNA, and they demonstrate that cell-free DNA whole-exome sequencing is concordant with metastatic tumour whole-exome sequencing.

Genomic Heterogeneity as a Barrier to Precision Medicine in Gastroesophageal Adenocarcinoma

Gastroesophageal adenocarcinoma (GEA) is a lethal disease where targeted therapies, even when guided by genomic biomarkers, have had limited efficacy. A potential reason for the failure of such therapies is that genomic profiling results could commonly differ between the primary and metastatic tumors. To evaluate genomic heterogeneity, we sequenced paired primary GEA and synchronous metastatic lesions across multiple cohorts, finding extensive differences in genomic alterations, including discrepancies in potentially clinically relevant alterations. Multiregion sequencing showed significant discrepancy within the primary tumor (PT) and between the PT and disseminated disease, with oncogene amplification profiles commonly discordant. In addition, a pilot analysis of cell-free DNA (cfDNA) sequencing demonstrated the feasibility of detecting genomic amplifications not detected in PT sampling. Lastly, we profiled paired primary tumors, metastatic tumors, and cfDNA from patients enrolled in the personalized antibodies for GEA (PANGEA) trial of targeted therapies in GEA and found that genomic biomarkers were recurrently discrepant between the PT and untreated metastases. Divergent primary and metastatic tissue profiling led to treatment reassignment in 32% (9/28) of patients. In discordant primary and metastatic lesions, we found 87.5% concordance for targetable alterations in metastatic tissue and cfDNA, suggesting the potential for cfDNA profiling to enhance selection of therapy.Significance: We demonstrate frequent baseline heterogeneity in targetable genomic alterations in GEA, indicating that current tissue sampling practices for biomarker testing do not effectively guide precision medicine in this disease and that routine profiling of metastatic lesions and/or cfDNA should be systematically evaluated. Cancer Discov; 8(1); 37-48. ©2017 AACR.See related commentary by Sundar and Tan, p. 14See related article by Janjigian et al., p. 49This article is highlighted in the In This Issue feature, p. 1.

Structural Alterations Driving Castration-Resistant Prostate Cancer Revealed by Linked-Read Genome Sequencing

Nearly all prostate cancer deaths are from metastatic castration-resistant prostate cancer (mCRPC), but there have been few whole-genome sequencing (WGS) studies of this disease state. We performed linked-read WGS on 23 mCRPC biopsy specimens and analyzed cell-free DNA sequencing data from 86 patients with mCRPC. In addition to frequent rearrangements affecting known prostate cancer genes, we observed complex rearrangements of the AR locus in most cases. Unexpectedly, these rearrangements include highly recurrent tandem duplications involving an upstream enhancer of AR in 70%-87% of cases compared with <2% of primary prostate cancers. A subset of cases displayed AR or MYC enhancer duplication in the context of a genome-wide tandem duplicator phenotype associated with CDK12 inactivation. Our findings highlight the complex genomic structure of mCRPC, nominate alterations that may inform prostate cancer treatment, and suggest that additional recurrent events in the non-coding mCRPC genome remain to be discovered.

Tumor cells are dislodged into the pulmonary vein during lobectomy.

OBJECTIVE Intraoperative tumor shedding may facilitate tumor dissemination. In earlier studies, shed tumor cells were defined primarily by cytomorphological examination, and normal epithelial cells could not always be distinguished from tumor cells. We sought to accurately identify tumor cells using single-cell sequencing and determine whether these cells were mobilized into the circulation during pulmonary lobectomy. METHODS Forty-two blood samples collected from the tumor-draining pulmonary vein at the end of lobectomy procedures were analyzed. Arrays of nanowells were used to enumerate and retrieve single EpCAM(+) cells. Targeted sequencing of 10 to 15 cells and nested polymerase chain reaction of single cells detected somatic mutations in shed epithelial cells consistent with patient-matched tumor but not normal tissue. RESULTS The mean number of EpCAM(+) cells in video-assisted thoracoscopy (VATS) lobectomy (no wedge) specimens (n = 16) was 165 (median, 115; range, 0-509) but sampling cells from 3 patients indicated that only 0% to 38% of the EpCAM(+) cells were tumor cells. The mean number of EpCAM(+) cells in VATS lobectomy (wedge) specimens (n = 12) was 1128 (median, 197; range, 47-9406) and all of the EpCAM(+) cells were normal epithelial cells in 2 patients sampled. The mean number of EpCAM(+) cells in thoracotomy specimens (n = 14) was 238 (median, 22; range, 9-2920) and 0% to 50% of total EpCAM(+) cells were tumor cells based on 4 patients sampled. CONCLUSIONS Surgery mobilizes tumor cells into the pulmonary vein, along with many normal epithelial cells. EpCAM alone cannot differentiate between normal and tumor cells. On the other hand, single-cell genetic approaches with patient-matched normal and tumor tissues can accurately quantify the number of shed tumor cells.

Multiplexed Elimination of Wild-Type DNA and High-Resolution Melting Prior to Targeted Resequencing of Liquid Biopsies

BACKGROUND The use of clinical samples and circulating cell-free DNA (cfDNA) collected from liquid biopsies for diagnostic and prognostic applications in cancer is burgeoning, and improved methods that reduce the influence of excess wild-type (WT) portion of the sample are desirable. Here we present enrichment of mutation-containing sequences using enzymatic degradation of WT DNA. Mutation enrichment is combined with high-resolution melting (HRM) performed in multiplexed closed-tube reactions as a rapid, cost-effective screening tool before targeted resequencing. METHODS We developed a homogeneous, closed-tube approach to use a double-stranded DNA-specific nuclease for degradation of WT DNA at multiple targets simultaneously. The No Denaturation Nuclease-assisted Minor Allele Enrichment with Probe Overlap (ND-NaME-PrO) uses WT oligonucleotides overlapping both strands on putative DNA targets. Under conditions of partial denaturation (DNA breathing), the oligonucleotide probes enhance double-stranded DNA-specific nuclease digestion at the selected targets, with high preference toward WT over mutant DNA. To validate ND-NaME-PrO, we used multiplexed HRM, digital PCR, and MiSeq targeted resequencing of mutated genomic DNA and cfDNA. RESULTS Serial dilution of KRAS mutation-containing DNA shows mutation enrichment by 10- to 120-fold and detection of allelic fractions down to 0.01%. Multiplexed ND-NaME-PrO combined with multiplexed PCR-HRM showed mutation scanning of 10-20 DNA amplicons simultaneously. ND-NaME-PrO applied on cfDNA from clinical samples enables mutation enrichment and HRM scanning over 10 DNA targets. cfDNA mutations were enriched up to approximately 100-fold (average approximately 25-fold) and identified via targeted resequencing. CONCLUSIONS Closed-tube homogeneous ND-NaME-PrO combined with multiplexed HRM is a convenient approach to efficiently enrich for mutations on multiple DNA targets and to enable prescreening before targeted resequencing.

Technological considerations for genome-guided diagnosis and management of cancer

Technological, methodological, and analytical advances continue to improve the resolution of our view into the cancer genome, even as we discover ways to carry out analyses at greater distances from the primary tumor sites. These advances are finally making the integration of cancer genomic profiling into clinical practice feasible. Formalin fixation and paraffin embedding, which has long been the default pathological biopsy medium, is now being supplemented with liquid biopsy as a means to profile the cancer genomes of patients. At each stage of the genomic data generation process—sample collection, preservation, storage, extraction, library construction, sequencing, and variant calling—there are variables that impact the sensitivity and specificity of the analytical result and the clinical utility of the test. These variables include sample degradation, low yields of nucleic acid, and low variant allele fractions (proportions of assayed molecules carrying variant allele(s)). We review here the most common pre-analytical and analytical factors relating to routine cancer patient genome profiling, some solutions to common challenges, and the major sample preparation and sequencing technology choices available today.

Whole-exome sequencing of cell-free DNA and circulating tumor cells in multiple myeloma

Liquid biopsies including circulating tumor cells (CTCs) and cell-free DNA (cfDNA) have enabled minimally invasive characterization of many cancers, but are rarely analyzed together. Understanding the detectability and genomic concordance of CTCs and cfDNA may inform their use in guiding cancer precision medicine. Here, we report the detectability of cfDNA and CTCs in blood samples from 107 and 56 patients with multiple myeloma (MM), respectively. Using ultra-low pass whole-genome sequencing, we find both tumor fractions correlate with disease progression. Applying whole-exome sequencing (WES) to cfDNA, CTCs, and matched tumor biopsies, we find concordance in clonal somatic mutations (~99%) and copy number alterations (~81%) between liquid and tumor biopsies. Importantly, analyzing CTCs and cfDNA together enables cross-validation of mutations, uncovers mutations exclusive to either CTCs or cfDNA, and allows blood-based tumor profiling in a greater fraction of patients. Our study demonstrates the utility of analyzing both CTCs and cfDNA in MM.Circulating tumor cells (CTCs) and cell-free DNA (cfDNA) enables characterization of a patient’s cancer. Here, the authors analyse CTCs, cfDNA, and tumor biopsies from multiple myeloma patients to show these approaches are complementary for mutation detection, together enabling a greater fraction of patient tumors to be profiled.