Citations

Abstract Background: Tumor-infiltrating lymphocytes (TIL) have prognostic and potentially predictive significance in the (neo)adjuvant treatment of high-risk breast cancer. However, quantitative TIL measurement is not routinely performed. It is unclear why some tumors attract large quantities of TIL while others do not. We sought to confirm the association between TIL and pathologic complete response rate (pCR) and to further use next generation sequencing (NGS) to identify genes and gene pathways associated with the presence/absence of TIL. Methods: We studied 140 women with high risk stage I-III breast cancer, enrolled in the Breast Cancer Genome Guided Therapy Study (BEAUTY), obtaining serial biopsies for DNA/RNA sequencing and MRI imaging to assess response to neoadjuvant chemotherapy (NAC) with taxane (+/- trastuzumab+/-pertuzumab for HER2+ disease) followed by AC or (F)EC. Diagnostic pre-NAC core needle biopsies and surgical resection specimens post-NAC were available from 110 patients. Stromal TIL were semi-quantitated on a scale of 1-4 (with 1: ≤10/hpf, 2: subtle infiltrate >10/hpf, 3: moderate infiltrate readily visible at low power magnification, 4: dense infiltrate with innumerable lymphocytes). For this analysis, low TIL was defined as scores of 1-2 vs. high defined as 3-4. Using pre-NAC biopsies, RNAseq was performed using the Illumina HiSeq2000 and the Mayo Analysis Pipeline for RNAseq (MAP-Rseq) for quality control, sequence alignment, and gene counts. The quantity of TIL was associated with transcripts across the transcriptome after conditional quantile normalization. Differentially expressed genes were obtained using EdgeR analysis, using a false discovery rate of 0.05, and pathways were evaluated using GAGE methods. Results: The pCR and residual cancer burden (RCB)-0/I rates by stromal TIL status within each molecular subtype are presented in the table. A diverse spectrum of 1344 genes with differential expression between tumors with high vs. low stromal TIL was identified. The genes with >2.0-fold change (FC) and p<1e-09 included S100A7 (4.49 FC), LCN2 (2.48 FC), and ART3 (2.82 FC) (genes known to be involved in immune regulation), as well as TDRD1 (2.71 FC) (a gene related to ERG [ETS-related gene] expression). In addition, the "regulation of actin cytoskeleton" pathway was upregulated in tumors with high TIL, while the "Hedgehog signaling" and "Wnt signaling" pathways were downregulated. Molecular SubtypeStromal TILspCR rate n (%)RCB-0/I rateLuminal AHigh------Luminal ALow0/9 (0%)0/9 (0%)Luminal BHigh1/9 (11.1%)1/8 (12.5%)Luminal BLow3/24 (12.5%)6/23 (26.1%)ER+/HER2+High3/9 (33.3%)4/9 (44.4%)ER+/HER2+Low1/6 (16.7%)1/6 (16.7%)ER-/HER2+High8/9 (88.9%)7/7 (100%)ER-/HER2+Low4/8 (50.0%)6/8 (75.0%)Triple NegativeHigh10/19 (52.6%)13/19 (68.4%)Triple NegativeLow7/14 (50.0%)9/13 (69.2%) Conclusions: We identified genes and gene pathways associated with high TIL expression in breast tumors prior to NAC that provide insight into the interactions between TIL and tumors. TIL can be easily semi-quantitated on H&E and along with these novel biomarkers, may contribute to the personalization of breast cancer therapy. Citation Format: Moyer AM, Boughey JC, Kalari KR, Suman VJ, McLaughlin SA, Moreno-Aspitia A, Northfelt DW, Gray RJ, Sinnwell JP, Carlson EE, Dockter TJ, Jones KN, Felten SJ, Conners AL, Wieben ED, Ingle JN, Wang L, Weinshilboum RM, Visscher DW, Goetz MP. Differential mRNA expression patterns in breast tumors with high vs. low quantity of stromal tumor–Infiltrating lymphocytes. [abstract]. In: Proceedings of the Thirty-Eighth Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2015 Dec 8-12; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2016;76(4 Suppl):Abstract nr P4-04-05.

Abstract Background: Tumorigenesis involves both genetic and epigenetic changes. Epigenetic alterations are reversible and are promising cancer therapeutic targets. Decitabine (5-aza-2'-deoxycytidine), a DNA methyltransferase inhibitor, is FDA approved for hematological malignancies. However, the effect of decitabine in breast cancer is not completely understood. Previous reports indicated that one decitabine mechanism involves regulation of protein levels for DNMT1, the major DNA methyltransferase that methylates hemimethylated CpG di-nucleotides in DNA. However, the E3 ligase involved in this process has not been identified. Whether decitabine also regulates DNMT3A and 3B in a similar fashion remains unclear. Therefore, our goals were to 1) understand mechanisms underlying decitabine action, 2) test the antitumor activity of decitabine in breast cancer models and 3) identify biomarkers associated with response to decitabine. Methods and Results: Western blots of breast cancer cell lines showed that DNMT1, DNMT3A, and DNMT3B protein levels decreased following decitabine treatment without a reduction in mRNA levels. Bioinformatic analysis of DNA methyltransferase sequences revealed a potential TRAF6 binding motif, and the interaction with TRAF6 (TNF receptor-associated factor 6) was confirmed by IP. TRAF6 functions as an E3 ligase. To determine whether TRAF6 might be the E3 ligase responsible for the degradation of DNMTs after decitabine treatment, we knocked down TRAF6 by RNA interference or knocked out the TRAF6 gene by CRISPR/Cas9. Down regulation of TRAF6 attenuated DNMT ubiquitination and increased DNMT protein levels, suggesting that TRAF6 might mediate proteasome-dependent degradation of all three DNMTs. This was further confirmed by reconstituting the knockout cells with WT and a TRAF6-C70A mutant, followed by assessing DNMT protein levels. Global DNA methylation was also increased after TRAF6 depletion and was confirmed in TRAF6 knock out cells in which DNMT levels were unaffected by decitabine. Cell cytotoxicity and colony forming assays showed that TRAF6 knockout cells were resistant to decitabine, suggesting that a major decitabine mechanism of action is through the regulation of TRAF6 which, in turn, degrades DNMTs, leading to decreased global methylation. Finally, decitabine significantly induced TRAF6 at both mRNA and protein levels, a process that might create positive feedback leading to increased degradation of DNMT proteins upon decitabine treatment. Based on these results, we further hypothesized that levels of the three DNMTs might influence decitabine response. Using 18 breast cancer patient derived xenograft (PDX) models, we found a wide range of DNMT protein levels regardless of ER/HER2 status. DNMT levels in the PDX models were directly associated with sensitivity to decitabine treatment, confirming our hypothesis. Conclusion: Our data showed that decitabine might be an effective agent for treating breast cancer and revealed a novel mechanism underlying decitabine treatment. Baseline DNMT protein levels may serve as a biomarker for predicting decitabine drug response. Citation Format: Yu J, Qin B, Boughey JC, Moyer AM, Visscher DW, Sinnwell JP, Yin P, Thompson KJ, Docter TJ, Kalari KR, Suman VJ, Wieben ED, Felten SJ, Conners AL, Jones KN, McLaughlin SA, Copland JA III, Moreno Aspitia A, Northfelt DW, Gray RJ, Ingle JN, Lou Z, Weinshilboum R, Goetz MP, Wang L. Regulation of DNA methyltransferases via TRAF6 determines breast cancer response to decitabine. [abstract]. In: Proceedings of the Thirty-Eighth Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2015 Dec 8-12; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2016;76(4 Suppl):Abstract nr P3-07-51.

Abstract Background: Both HR deficiency and BRCA mutation status predict response to platinum-based therapy and BRCA mutation status predicts docetaxel resistance. However, the association of either biomarker with response to the individual elements of either AC or taxanes (T) is unknown since T is commonly given concomitantly with or after anthracyclines (A). We evaluated the association of HRD and BRCA mutation status with response to neoadjuvant weekly T followed by AC or (F)EC in high-risk breast cancer. Methods: We studied 140 high risk Stage I-III breast cancer patients (pts), enrolled in the breast cancer genome guided therapy study (BEAUTY), obtaining biopsies for DNA/RNA sequencing and MRI imaging to assess response to neoadjuvant weekly T (+trastuzumab+/-pertuzumab for HER2+ disease) followed by AC or (F)EC. Germline BRCA status and HR status of tumor samples (Myriad laboratories) were obtained. HR deficient tumor was defined as HRD score ≥42 or BRCA mutation. MRI response by changes in tumor size after 12 weeks of T was classified by WHO criteria. pCR was defined as ypT0/Tis ypN0. Both MRI response after T and pCR (after T and AC) were examined in terms of germline BRCA mutation (gBRCAmut vs. gBRCAwt) and tumor HR deficiency. Results: Of 140 pts enrolled, 8 withdrew consent and 2 carboplatin treated pts were excluded. Germline data were available for 124/130 pts. 12 patients had BRCA deleterious germline mutations (4 BRCA1, 8 BRCA2). MRI partial (PR)/complete response (CR) rate to T was 47.3% (95% CI: 37.8-57.0%) in the BRCAwt group and 66.7% (95% CI: 34.9-90.1%) in the BRCAmut group. No MRI CR's were observed in BRCA1 mut pts. In contrast, pCR rate was 50% in the 12 gBRCAmut pts (95% CI: 21.1-78.9%) and 31.3% in the 112 gBRCAwt pts (95% CI: 22.8-40.7%). HR deficiency status has thus far been determined for 74 pts: 26 pts have HD deficient tumors: 18 TNBC, 5 Luminal B, 2 ER-/HER2+; and 1 ER+/HER2+. Determination of HR deficiency is ongoing and will be reported for the full cohort in terms of 12 week MRI response to T and pCR to T+AC. HR deficientMolecular Subtypeyes (%)no (%)TBD (%)Luminal A0/112/11 (18.2)9/11 (81.8)Luminal B5/37 (13.5)13/37 (35.1)19/37 (51.3)Luminal NOS0/21/2 (50)1/2 (50)ER+/Her2+1/17 (5.8)14/17 (82.4)2/17 (11.8)ER-/Her2+2/20 (10)11/20 (55)7/20 (35)Triple Negative18/43 (41.9)6/43 (18.6)17/43 (39.5)germline BRCA statusMRI partial response after T (%)MRI complete response after T (%)pCR after T&AC (%)BRCA11/4 (25)0/42/4 (50)BRCA25/8 (62.5)2/8 (25)4/8 (50)BRCAwt35/112 (31.3)18/112 (16.1)35/112 (31.3) Conclusion: In the setting of neoadjuvant weekly T followed by AC, pCR rates were non-significantly higher in pts with BRCA1 mutations. While we observed no overall association between BRCA mutation status and response rates to taxanes; nearly all MRI responses to taxanes (partial and complete) were observed in the BRCA2 group. Prospective studies are needed to validate these findings and to determine whether BRCA status can be used to select therapy. HR deficiency is uncommon in luminal A and HER2+, frequent in TNBC, and the association of HRD with both MRI response to taxanes and pCR will be reported at the meeting. Citation Format: Boughey JC, Kalari KR, Suman VJ, McLaughlin SA, Moreno Aspitia A, Moyer AM, Northfelt DW, Gray RJ, Vedell PT, Tang X, Dockter TJ, Jones KN, Felten SJ, Conners AL, Hart SN, Visscher DW, Wieben ED, Ingle JN, Hartman A-R, Timms K, Elkin E, Jones J, Wang L, Weinshilboum RW, Goetz MP. Role of germline BRCA status and tumor homologous recombination (HR) deficiency in response to neoadjuvant weekly paclitaxel followed by anthracycline-based chemotherapy. [abstract]. In: Proceedings of the Thirty-Eighth Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2015 Dec 8-12; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2016;76(4 Suppl):Abstract nr P3-07-29.

Abstract The addition of trastuzumab to chemotherapy extends survival among patients with HER2+ breast cancer. Prior work showed that trastuzumab and chemotherapy augments HER2 extracellular domain (ECD)-specific antibodies. The current study investigated whether combination therapy induced immune responses beyond HER2-ECD and, importantly, whether those immune responses were associated with survival. Pretreatment and posttreatment sera were obtained from 48 women with metastatic HER2+ breast cancer on NCCTG (now Alliance for Clinical Trials in Oncology) studies, N0337 and N983252. IgG to HER2 intracellular domain (ICD), HER2-ECD, p53, IGFBP2, CEA, and tetanus toxoid were examined. Sera from 25 age-matched controls and 26 surgically resected HER2+ patients were also examined. Prior to therapy, some patients with metastatic disease had elevated antibodies to IGFBP2, p53, HER2-ICD, HER2-ECD, and CEA, but not to tetanus toxin, relative to controls and surgically resected patients. Treatment augmented antibody responses to HER2-ICD in 69% of metastatic patients, which was highly associated with improved progression-free survival (PFS; HR = 0.5, P = 0.0042) and overall survival (OS; HR = 0.7, P = 0.038). Augmented antibody responses to HER2-ICD also correlated (P = 0.03) with increased antibody responses to CEA, IGFBP2, and p53, indicating that treatment induces epitope spreading. Paradoxically, patients who already had high preexisting immunity to HER2-ICD did not respond to therapy with increased antibodies to HER2-ICD and demonstrated poorer PFS (HR = 1.6, P < 0.0001) and OS (HR = 1.4, P = 0.0006). Overall, the findings further demonstrate the importance of the adaptive immune system in the efficacy of trastuzumab-containing regimens. Cancer Res; 76(13); 3702–10. ©2016 AACR.

Abstract Motivation: Single-cell RNA-sequencing technology allows detection of gene expression at the single-cell level. One typical feature of the data is a bimodality in the cellular distribution even for highly expressed genes, primarily caused by a proportion of non-expressing cells. The standard and the over-dispersed gamma-Poisson models that are commonly used in bulk-cell RNA-sequencing are not able to capture this property. Results: We introduce a beta-Poisson mixture model that can capture the bimodality of the single-cell gene expression distribution. We further integrate the model into the generalized linear model framework in order to perform differential expression analyses. The whole analytical procedure is called BPSC. The results from several real single-cell RNA-seq datasets indicate that ∼90% of the transcripts are well characterized by the beta-Poisson model; the model-fit from BPSC is better than the fit of the standard gamma-Poisson model in > 80% of the transcripts. Moreover, in differential expression analyses of simulated and real datasets, BPSC performs well against edgeR, a conventional method widely used in bulk-cell RNA-sequencing data, and against scde and MAST, two recent methods specifically designed for single-cell RNA-seq data. Availability and Implementation: An R package BPSC for model fitting and differential expression analyses of single-cell RNA-seq data is available under GPL-3 license at https://github.com/nghiavtr/BPSC . Contact: yudi.pawitan@ki.se or mattias.rantalainen@ki.se Supplementary information: Supplementary data are available at Bioinformatics online.

Abstract We previously reported, on the basis of a genome-wide association study for aromatase inhibitor-induced musculoskeletal symptoms, that single-nucleotide polymorphisms (SNPs) near the T-cell leukemia/lymphoma 1A (TCL1A) gene were associated with aromatase inhibitor-induced musculoskeletal pain and with estradiol (E2)-induced TCL1A expression. Furthermore, variation in TCL1A expression influenced the downstream expression of proinflammatory cytokines and cytokine receptors. Specifically, the top hit genome-wide association study SNP, rs11849538, created a functional estrogen response element (ERE) that displayed estrogen receptor (ER) binding and increased E2 induction of TCL1A expression only for the variant SNP genotype. In the present study, we pursued mechanisms underlying the E2-SNP-dependent regulation of TCL1A expression and, in parallel, our subsequent observations that SNPs at a distance from EREs can regulate ERα binding and that ER antagonists can reverse phenotypes associated with those SNPs. Specifically, we performed a series of functional genomic studies using a large panel of lymphoblastoid cell lines with dense genomic data that demonstrated that TCL1A SNPs at a distance from EREs can modulate ERα binding and expression of TCL1A as well as the expression of downstream immune mediators. Furthermore, 4-hydroxytamoxifen or fulvestrant could reverse these SNP-genotype effects. Similar results were found for SNPs in the IL17A cytokine and CCR6 chemokine receptor genes. These observations greatly expand our previous results and support the existence of a novel molecular mechanism that contributes to the complex interplay between estrogens and immune systems. They also raise the possibility of the pharmacological manipulation of the expression of proinflammatory cytokines and chemokines in a SNP genotype-dependent fashion.

Abstract Breast cancer is the most common invasive cancer in women worldwide and treatment related toxicity has high inter-patient variability. Our goal was to identify additional genetic predictors of toxicity by completing a pathway based analysis taking advantage of the genotyping results obtained through a GWAS of the same patients. For our analysis, we have utilized GWA data from the SUCCESS-A trial in which 3,754 breast cancer patients were treated with adjuvant chemotherapy. We have included SNPs associated with 112 genes involved in drug disposition and 375 genes in selected cancer pathways. We identified 35 SNPs associated with neutropenia (P-value <1×10−4). These patients initially received three cycles consisting of 5-fluorouracil (5-FU), cyclophosphamide, and epirubicin. Currently we are characterizing the function of three SNPs located within an intron of the hyaluronan mediated motility receptor (HMMR) gene that are associated with the neutropenia phenotype (P-value <1×10−4). These SNPs are highly associated with the expression of a gene encoding a B-cell survival factor, TNFSF13B, (P-value <1×10−5). Our in vitro studies using lymphoblastoid cell lines (LCLs) selected based on the genotypes of HMMR gene showed that heterozygous and homozygous variant cells (n = 6) have an increased sensitivity to 5-fluorouracil, but not cyclophosphamide or epirubicin when compared with homozygous wild-type cells (n = 6). LCLs exposed to the previously mentioned chemotherapeutic agents resulted in an upregulation of TNFSF13B gene expression and downregulation of HMMR gene expression. We have utilized a pathway based approach to analyzing GWA data to identify several SNPs in the HMMR gene highly associated with the expression of one of the top signals in the cytopenia GWAS. Our goal is to further characterize the role of HMMR gene variants in chemotherapy response by assessing whether combinations of these agents in vitro (including 5-FU) also have a genotype dependent response to cytotoxicity. Citation Format: Behzad Bidadi, Krishna R. Kalari, Matthias Rubner, Alexander Hein, Matthias W. Beckmann, Peter A. Fasching, Brigitte Rack, Wolfgang Janni, Richard Weinshilboum, Liewei Wang. Utilizing a pathway based analysis of genome wide association data to identify biomarkers of toxicity in breast cancer patients. [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 5488. doi:10.1158/1538-7445.AM2015-5488

Abstract Introduction: Triple negative breast cancers (TNBC) are characterized as lacking estrogen receptor (ER), progesterone receptor (PR) and human epidermal growth factor receptor 2 (HER2) expression and TNBC patients have higher rates of recurrence and death compared with other breast cancer subtypes. For TNBC patients who fail standard chemotherapy, there are a lack of novel drug therapies, given the absence of well-defined molecular targets. Recently, a microarray meta-analysis identified 7 triple negative subtypes, including the validation of the luminal androgen receptor (LAR) positive subtype [Lehmann, 2011]. However, microarray technology is dependent on probe-target specificity and the 7 subtypes have yet to be validated using RNA sequencing data, and the presences of recurrent genomic alterations in the 7 subtypes are unknown. METHODS: We obtained 1106 breast cancer RNA-Seq bam files from The Cancer Genome Atlas (TCGA) and aligned with Tophat v1.3. The PAM50 intrinsic gene signature was used to extract a cohort of 128 TNBC samples. Consensus clustering of genes, greater than 75th percentile of variance, was performed using Kmeans clustering in Spearman's correlation space. A nearest centroid prediction model was developed from genes differentially expressed among the clusters. Eighty independent TNBC RNA sequencing samples were obtained (British Colombia; BC) [Shah, 2012] which were calibrated to our TNBC conditional quantile normalized cohort and sub-typed by our model. RESULTS: Using RNA-Seq gene expression count data, we identified 5 clusters, all of which were stable, including the LAR cluster. Signaling pathway impact analysis (SPIA) implicated cytokine-cytokine receptor interaction, leukocyte transendothelial migration, and regulation of actin cytoskeleton pathways commonly altered in the non-LAR TNBC subtypes. In contrast, cell cycle, ECM-receptor interaction, endocrine regulated calcium reabsorption, and insulin signaling pathways were altered in the LAR versus non-LAR subtypes. Neuroactive ligand-receptor interactions were observed to be altered commonly between all sub-types. We then applied our model to the Shah, et al cohort. In this cohort, the LAR subtype was consistent with Shah's classification of ‘other’ TNBC and contained no basal samples by PAM50 intrinsic modeling. Analysis of sub-type specific mutation data from the BC cohort demonstrates an increased mutational load in ECM-related proteins, particularly the myosins, along with increased TP53 clonality in the non-LAR subtypes. CONCLUSIONS: Using TCGA RNASeq data, we have confirmed the presence of 5 major TNBC subtypes, including the LAR; which was negligible in basal composition by PAM50 intrinsic modeling. SPIA pathway analysis indicates a core set of pathways demonstrating altered expression across the TNBC sub-types and the identification of molecular targets within each subtype is ongoing. Citation Format: Kevin J. Thompson, Xiaojia Tang, Zhifu Sun, Jason P. Sinnwell, Hugues Sicotte, Douglas W. Mahoney, Steven Hart, Peter T. Vedell, Poulami Barman, Jeanette E. Eckel Passow, Eric D. Wieben, James N. Ingle, Judy C. Boughey, Liewei Wang, Richard Weinshilboum, Krishna R. Kalari, Matthew P. Goetz. Molecular classification of triple negative breast cancer via RNA-sequencing data. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 5592. doi:10.1158/1538-7445.AM2014-5592

Abstract INTRODUCTION Next generation sequencing (NGS) of patients has significantly changed our ways to study cancer genomics as it provides precise estimates of gene expression, fusion transcripts, expressed single nucleotide variants (eSNVs), splice variants and copy number variants. The Breast Cancer Genome Guided Therapy (BEAUTY) is an ongoing clinical study in which RNA sequencing (RNAseq) and whole exome sequencing (WES) are performed prior to, during and after neoadjuvant chemotherapy. Here we report the use of these sequencing technologies to investigate gene expression levels and mutational profiles in a triple negative breast cancer (TNBC) patient enrolled in BEAUTY whose disease did not respond to neoadjuvant paclitaxel and anthracycline/cycphosphamide. METHODS Computational approaches were used to integrate WES and RNAseq data obtained before therapy (V1T), after 12 weekly paclitaxel treatments (V2T) and after anthracycline-based regimen at surgical resection (V3T) to study a single patient with persistent TNBC disease after neoadjuvant chemotherapy. RESULTS Using RNA-Seq data, we identified an inter-chromosomal fusion transcript between chromosome 20 and 22 (GNAS-TTC38) that was highly expressed at V1T, V2T and V3T. We also identified intra-chromosomal fusion transcripts that were expressed at two time points, such as fusion transcript (KANSL1-ARL17A) on chromosome 17 for V1T and V2T and fusion transcript (RBM12B-LINC00535) on chromosome 8 for V2T and V3T time points. Several gene expression changes were also observed. Gene expression analysis of V1T, V2T and V3T tumors was performed. Differential temporal gene expression profiles of 9884 genes that were significant and varying at different time points were obtained for pathway analysis. Pathway analysis of 9884 genes identified up regulation and down regulation of several transcription factors with a fold change of 2x or more. When compared to blood, DNA tumor and RNA-Seq data, we identified 81 common somatic eSNVs that were expressed in both V1T and V2T time points and we are in the process of investigating V3T data. We found alterations of key transporter domains (CD225, Coatamer_beta_C, DUF2435, Dynamitin, EI24, GLTP, LMF1, Porin_3, V-ATPase_C) in our V1T and V2T SNV data. Similar to gene expression analysis, we are in the process of obtaining the list of mutations at various time points to identify driver and passenger mutation candidate genes for this specific TNBC patient. CONCLUSIONS Our initial time-series analysis of eSNV, fusion transcripts and gene expression data demonstrate that intensive analysis for individual patients is feasible. Further investigation of drug transporters and transcription regulators may help develop personalized treatment strategies for patients with disease resistant to current regimens. Citation Format: Krishna R. Kalari, Xiaojia Tang, Kevin J. Thompson, Douglas W. Mahoney, Poulami Barman, Jason P. Sinnwell, Hugues Sicotte, Peter Vedell, Steven N. Hart, Travis J. Dockter, Katie N. Jones, Amy L. Conners, Ann M. Moyer, Daniel W. Visscher, Jia Yu, Bowen Gao, Sarah A. McLaughlin, John A. Copland, Alvaro Moreno-Aspitia, Donald W. Northfelt, Richard J. Gray, Vera J. Suman, Jeanette E. Eckel Passow, Jean-Pierre A. Kocher, Eric D. Wieben, Gianrico Farrugia, Cloann G. Schultz, James N. Ingle, Richard Weinshilboum, Matthew P. Goetz, Liewei Wang, Judy C. Boughey. Analysis of sequencing data to identify potential drug targets for an individual newly diagnosed with basal breast cancer who failed to respond to current standard neoadjuvant chemotherapy. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 4185. doi:10.1158/1538-7445.AM2014-4185

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.