Citations

Abstract Background Genomic intratumor heterogeneity is an increasingly recognized phenomenon in various solid tumors and is now extensively being investigated by next-generation sequencing (NGS) approaches. However, most methods are not designed to resolve the complexity of mixed cell populations. Therefore important information may be lost, eventually leading to assumptions that may not reflect the true nature of the disease. We have therefore developed the method of Multiparameter Ploidy Profiling (MPP). We applied this technology on tumor specimens from malignant melanoma and from non-small cell lung cancer (NSCLC). Methods MPP involves isolation of nuclei from tissues, multiparameter flow-sorting of tumor populations by DNA content and specific tumor markers, and profiling genomes using high resolution array comparative genomic hybridization (aCGH) and NGS. As some tumors are diploid by flow cytometry, specific lineage markers, e.g. SOX10 (melanoma) and TTF1 (NSCLC) were used to separate these from normal cells. Fluorescence-in-situ hybridization (FISH) and immunohistochemistry (IHC) were performed on sorted nuclei and on whole tissue sections to validate results. Results MPP allowed us to sort different tumor populations within a single tumor specimen, while excluding non-neoplastic contamination in downstream analysis with aCGH and NGS. We separated three different tumor populations from a single primary NSCLC with a TTF1 multiparameter sort. We detected significant differences in copy numbers, e.g. a CDKN2A homozygous deletion only present in the diploid tumor population as well as differences in allelic frequencies of gene mutations (e.g. EGFR exon 20). Further sequencing of the diploid TTF1 negative population proved to be concordant to genomic DNA obtained from a tumor free lymph node. Finally, we used the population-specific genomic data in order to infer the clonal relationship of all tumor populations and to postulate the evolution of these tumor populations within single tumor specimens. Conclusions The power of MPP consists of sorting pure tumor populations and of deciphering the level of intratumor heterogeneity based on ploidy. Most importantly MPP enables the separation of diploid tumor populations from diploid, benign stromal cells that can serve as reference for genomic profiling studies, where germline DNA is otherwise not available. To our knowledge, diploid tumor populations were for first time purifyed from normal cells in a multiparameter sort and genomically characterized. We also demonstrated that the analysis of diploid tumor populations within a single biopsy is crucial for the understanding of the clonal composition and evolution of a given tumor. MPP on multiple biopsies from individual patients is necessary to deeply investigate the clonal evolution and the behavior of a diploid tumor population and their potential impact on metastasis and therapy response. Citation Format: Thomas Lorber, Sabrina Rau, Valeria Perrina, Michael Barrett, Christian Ruiz, Lukas Bubendorf. Multiparameter ploidy profiling: a powerful tool to investigate the genomics of diploid tumor populations. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4924. doi:10.1158/1538-7445.AM2015-4924

Abstract Background: Next generation sequencing (NGS) has revealed that the genetic profiles of individual tumors are highly heterogeneous due to their subclonal composition. Tumor heterogeneity has profound clinical implications affecting differences in treatment response and therapeutic resistance. Methods: An estrogen receptor positive, Her2 normal breast cancer patient underwent definitive surgical treatment with modified radical mastectomy. Multiple samples were obtained from the primary cancer and metastatic lymph nodes. Flow cytometry based methods were used to isolate and classify distinct cell subpopulations according to their ploidy, and samples were processed for whole exome sequencing (WES). Comparative Genomic Hybridization (CGH) methods were used to identify somatic unique copy-number alterations (CNAs). WES-Single Nucleotide Variations (SNVs) were used to infer subclonal phylogenetic relationships using Maximum Parsimony methods and annotated with variation-cluster-barcodes (vc-barcodes) generated using a variation Bayesian mixture model (VBMM) approach that focuses on copy-number neutral sequence variations for subclone identification (SciClone). Our comparative phylogenetic-VBMM clustering method was used to identify CNA and/or SNV drivers that are key in differentiating subclonal populations at multiple tumor sites. Results: CNA-neutral WES-SNVs from 25 subclonal populations (isolated from 4 distinct sites) clustered into 6 major groups and were used to generate vc-barcodes. A phylogenetic reconstruction of subclonal architecture annotated with vc-barcode information expedited identification of key variations underlying the differentiation of subclones at distinct tumor sites. For example, we identified an amplification event involving the Anaplastic Lymphoma Kinase (ALK) gene that was more frequent in primary biopsy subclones (12/17, ∼71%) in comparison to metastatic subclones (2/6, ∼34%). As the ALK amplification is lost, we find that most of the metastatic subclones also acquire a predicted damaging mutation in this oncogene (4/6, ∼67%). In addition, we identified a potential driver metastatic cell lineage that carries the ALK amplification in the absence of the nonsynonymous mutation. Conclusions: We developed a hybrid clustering methodology and used it to reconstruct the subclonal architecture in primary and metastatic tumors from a single breast cancer patient. Our methods have identified several CNAs and/or SNVs underlying subclonal differentiation, as well as potential biomarkers of disease progression and indicators of emerging resistance. Application of these methods to larger cohorts and types of tumors should be conducted to ascertain more precise estimates of the predictive accuracy of our processes. Citation Format: Mia D. Champion, Princy Francis, Barbara A. Pockaj, Michael T. Barrett. Hybrid clustering methodologies to distinguish CNAs and/or SNVs that drive subclonal differentiation in samples from a breast cancer patient primary tumor and metastatic lymph nodes. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4861. doi:10.1158/1538-7445.AM2015-4861

Abstract Introduction: Overexpression and amplification of PD-L1 has been reported in a subset of cancers, including lymphoma, and gastric cancer. We have recently described amplification of chromosome 9p24 targeting PD-L1 in a subset of triple negative breast cancers (TNBC), as well as 1-3% of colorectal carcinomas and glioblastomas. The clinical impact of these amplifications is not known. Methods: Tumor samples from 41 subjects with TNBC were evaluated for copy number aberrations. We performed DNA-based flow cytometry to sort nuclei from diploid, tetraploid, and aneuploid populations of tumor cells. The genomes of each sorted tumor cell population were then interrogated with oligonucleotide CGH arrays. Copy number intervals were identified and their genomic boundaries mapped using a step gram algorithm. Clinical data was available on 36 of 45 (80%) patients. Results: A recurring top ranked high level (log2ratio >1) amplicon that spanned 9p24.1 was detected in 12/41 (29.3%) of TNBCs. In contrast this amplicon was absent in ER+ (n = 8) and Her2+ (N = 15) breast tumors. The shortest region of overlap (SRO) spanned 777 kb and always included the PD-1 ligands PD-L1, PD-L2, and Janus kinase 2 (JAK2) loci, termed the “PDJ” amplicon. Additional genes in this amplicon include RCL1, MIR101-2, INSL6, INSL4, RLN2, RLN1, and C9orf46. The height of the PDJ amplicon included mean log2ratios of >4, consistent with genomic drivers such as Her2 and Myc. Patients with PD-L1 amplification (n = 8) were noted to have larger tumors (mean 3.9 vs. 1.9 cm, p = 0.04) and a higher incidence of lymph node metastases (75% vs. 26%, p = 0.01). Lymphocytic infiltration was only noted in 4 patients, all without PDJ amplification. Almost all patients received chemotherapy after definitive surgical therapy. The recurrence rate at 5 years was 62.5% in the PDJ amplified patients, and 25% in the unamplified patients (p = 0.02). Overall survival at 5 years was 37.5% in the PDJ amplified patients, compared with 78.2% in the unamplified patients (p = 0.02). Conclusion: We identified a focal amplification of chromosome 9p24.1 involving PD-L1, PD-L2, and JAK2 in a significant proportion of TNBCs. The high level of amplification and inverse correlation with survival suggests that the PDJ amplicon is a genetic driver event in these cancers, and warrants confirmation in larger studies. Citation Format: Karen S. Anderson, Christine L. Klassen, Amylou C. Dueck, Ann E. McCullough, Donald W. Northfelt, Mariacarla Andreozzi, Barbara A. Pockaj, Michael T. Barrett. High-level amplification of chromosome 9p24 targeting PD-L1 and JAK2 correlates with worse DFS and OS in triple negative breast cancer. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 2861. doi:10.1158/1538-7445.AM2015-2861

Abstract Background Cancer genomic investigations have identified recurrent genomic aberrations critical for cancer initiation, progression, and metastases. However, these investigations are typically performed in isolation, and the effects of treatment on the clonal selection of tumor cells are mostly unknown. We hypothesized that molecular profiling of residual tumors after neoadjuvant chemotherapy (NAC) would identify new drug targets/pathways in patients at high risk for disease recurrence. To better identify clonal populations of resistant breast cancer cells, we utilized DNA content-based flow sorting of nuclei to identify and isolate clonal populations for aCGH and next generation sequencing (NGS) both before and after NAC. Methods The Breast Cancer Genome Guided Therapy Study (BEAUTY) (NCT 02022202) is a prospective study of patients with high-risk breast cancer treated with neoadjuvant 12 weekly paclitaxel (T) +/- trastuzumab followed by 4 cycles of anthracycline based chemotherapy. Tumor tissue from baseline, residual disease from surgery, distant metastases, and patient derived xenografts (PDX) are obtained for cell sorting by DNA ploidy, aCGH, RNA and exome sequencing. Results: 140 patients have been enrolled, 104 have completed surgery and 30 unique PDX have been established corresponding to 26 patients prior to chemotherapy and 4 from residual disease at surgery. Baseline exome and RNA sequencing is complete in 140. Currently, genomic analyses of flow sorted matched baseline, surgical, PDX, and distant disease samples are available in 6 patients. Substantial genomic variation was observed in the surgical sample compared to the primary tumor including gain of oncogenic drivers (EGFR) and loss of negative regulators (ATG5) (Table). The PDX recapitulated these events with excellent fidelity compared to the corresponding human tumor. In patients with TNBC, RNA seq obained from matched samples demonstrated changes in immune related pathways. Evaluation of drug targets/pathways identified in the resistant tumors are ongoing using the PDX and sequencing of the remaining matched baseline/surgical disease will be reported. Clonal changes occurring over time in patients with residual disease or disease recurrence after NACTumor SubtypeResidual Cancer BurdenDisease StatusClonal Aberration Changes (baseline and post NAC)TNBC (AR subtype)3Contralateral lymph node recurrence at 150 days2p25.2 - p25.1 amplicon lost at recurrence; 6p21.32 -p21.31 amplicon lost at recurrenceTNBC (BL 1)3Bone/liver/lymph node recurrence at 135 days5q11.2 amplicon gained at surgery; 12p13.33 - p13.2 amplicon gained at surgeryTNBC (BL 1)3Disease-free at day 1506q21 amplicon lost at sugeryTNBC (BL 2)0Brain recurrence at 390 daysChr 2 chromothripsis at surgeryLuminal B3Progression during chemotherapy7p12.1 - p11.2 amplicon gained at surgeryLuminal HER23Disease-free at day 3609q33.2 amplicon lost at surgery; 15q13.1 - q13.3 amplicon lost at surgery; 15q22.2 - q24.1 amplicon lost at surgery Conclusions: We observed substantial evolutionary changes in residual breast tumors remaining after NAC. Our findings suggest that a comprehensive assessment of the mutational landscape that has evolved during NAC can inform drug development in high risk breast cancer patients. Citation Format: Matthew P Goetz, Michael T Barrett, Krishna R Kalari, Vera J Suman, Sarah A McLaughlin, Alvaro Moreno-Aspitia, Ann M Moyer, Donald W Northfelt, Richard J Gray, Jason Sinnwell, Douglas Mahoney, Poulami Barman, Peter Vedell, Xiaojia Tang, Kevin Thompson, Travis Dockter, Katie Jones, Sara J Felten, Amy Conners, Jeanette Eckel-Passow, Hughes Sicotte, Steven N Hart, Jia Yu, Daniel W Visscher, Eric D Wieben, Cloann Schultz, Minetta C Liu, James N Ingle, Liewei Wang, Richard W Weinshilboum, Judy C Boughey. Impact of neoadjuvant chemotherapy on the clonal composition of breast cancer [abstract]. In: Proceedings of the Thirty-Seventh Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2014 Dec 9-13; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2015;75(9 Suppl):Abstract nr PD3-3.

Abstract Background: Breast tumors exhibit intratumor heterogeneity resulting in targeted therapy resistance and other challenges in disease management. To address the sources of heterogeneity, we performed a unique, in-depth analysis of clonal architecture in primary chemoradiation-naïve breast cancers. We combined DNA content-based flow cytometry and ploidy analysis with aCGH and next-generation sequencing (NGS) in multiple biopsies from the tumors and involved lymph nodes (LNs). Material and methods: We used DNA content-based flow sorting to isolate nuclei from distinct populations of diploid and aneuploid tumor cells in surgical tumor samples from two chemoradiation-naïve patients. Each sorted tumor cell population was interrogated with aCGH and exome NGS. In Patient #1, we used 12 fresh frozen sections morphologically mapped from within a HER2+, ER+, PR- primary invasive ductal carcinoma (IDC) of histological grade 3 with LN involvement and 2-3 sections from 2 out of 5 LNs. In Patient #2, 11 morphologically mapped fresh frozen sections were analyzed from a grade 2, ER+, PR+, HER2+, BRCA2 mutant LN- IDC. In parallel, matching samples were processed for IHC assays. Results: We identified multiple co-existing aneuploid populations within the biopsies. The 17 primary and LN biopsies from Patient #1 fell into 6 distinct ploidy groups albeit with aberrant but homogenous aCGH profiles, characterized by SARC amplification and homozygous deletions in ROBO1 and ROBO2. In contrast a dominant ploidy was identified throughout Patient #2 but with heterogeneous aCGH profiles. Mutation profiles obtained through exome sequencing further confirmed that ploidy was the main driver in Patient #1 whereas copy number aberrations played the key role in Patient #2 with the BRCA2 mutation (R3129X). A dendrogram based on exome variant calls of the aneuploid populations in Patient #1 strongly correlated with ploidy group and further revealed the specific clonal population characterized by a 5N ploidy and homozygous mutations in TP53 and PIK3CA as the progenitor to the ploidies present in the distant LNs. Strikngly, both patients had no HER2 amplification or mutation across their clonal populations, contradicting the initial IHC staining in a single core biopsy. Conclusions Rather than inferring the presence of distinct tumor cell populations, our novel flow-sorting based approach of first identifying the clonal populations and then interrogating their genomes, provides an objective method of exploring the sources and clinical significance of tumor heterogeneity. Our approach of clonal analysis has broad implications in the study of tumor heterogeneity and the identification of drivers in breast and other solid tumors that can advance more effective treatment and clinical management of patients with this disease. Citation Format: Princy Francis, Lucretia M Alvarez, Elizabeth Lenkiewicz, Mia Champion, Lisa Evers, Karen L Anderson, Ann E McCullough, Michael T Barrett, Barbara Pockaj. Deep clonal profiling identifies distinct mechanisms of heterogeneity and evolution in breast cancer [abstract]. In: Proceedings of the Thirty-Seventh Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2014 Dec 9-13; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2015;75(9 Suppl):Abstract nr P2-03-07.