Wang, Liewei

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 Purpose: Tamoxifen and raloxifene, are primary prevention strategies for women at high risk of breast cancer (Visvanathan, 2013). Recent advances in genetic studies of breast cancer risk have identified common susceptibility loci that explain 14% of familial risk for breast cancer in the general population (Michailidou, 2013). However, it is not known if these loci are risk factors for breast cancer among high-risk women treated with SERMs for breast cancer prevention. We hypothesized that the large risk reduction associated with SERMs, coupled with the fact that several breast cancer loci correlate with family history, may limit the contribution of these common genetic loci to breast cancer in this high risk population. We present the first report to evaluate 75 established breast cancer susceptibility loci, in the context of a polygenic risk score (PRS), as a risk factor for breast cancer among high risk women taking raloxifene and tamoxifen for breast cancer prevention. Methods: We conducted a matched case-control study of 594 cases (i.e., participants who developed breast cancer while on SERM therapy) and 1,171 matched controls selected from the 33,000 participants enrolled in the NSABP P-1 and P-2 breast cancer prevention trials. Genotypes of 75 single nucleotide polymorphisms (SNPs) were available from a genome-wide association study conducted at the RIKEN Center for Genomic Medicine. We formed a quantitative PRS from reported per-SNP odds ratios (OR) for the 75 susceptibility loci. Conditional logistic regression was used to examine the PRS as a risk factor for breast cancer and to assess whether the PRS and breast cancer association differed by treatment type, family history, or other clinical characteristics. Analyses also examined associations of PRS with invasive vs. in situ cancer and ER-positive vs. ER-negative cancer. Results: The PRS ranged from 3.98 to 7.74, and a one unit change in PRS was associated with a 42% increase in breast cancer (OR=1.42; 95% CI: 1.18-1.70; P = 0.0002). There was evidence of a stronger association of PRS with breast cancer among women with no first-degree family history (OR=1.62, 95% CI: 1.18-2.21) compared to those with a positive family history (OR=1.32, 95% CI: 1.06-1.66) (Pintx<0.05). The PRS also appeared a stronger risk factor for ER-positive (OR=1.59, 95% CI: 1.25-2.02, P < 0.0002) vs. ER-negative (OR=1.05, 95% CI: 0.68-1.62, P=0.84) breast cancer, although differences did not reach statistical significance (Pintx=0.10). PRS and breast cancer associations were similar across tamoxifen and raloxifene treatments, age at trial entry, 5-year predicted Gail model risk, hysterectomy status, BMI, presence of atypical hyperplasia and invasive vs. in situ cancer. Conclusion: A polygenic risk score composed of 75 loci was a risk factor for ER-positive breast cancer, especially in the absence of a first-degree family history of breast cancer. Further, the PRS associations with breast cancer were similar for women taking tamoxifen or raloxifene for prevention. These data suggest that common genetic variation adds information on risk of ER-positive breast cancer in a high-risk population receiving SERMs. Citation Format: Celine M Vachon, Daniel J Schaid, James N Ingle, Matthew P Goetz, Donald L Wickerham, Michiaki Kubo, Erin E Carlson, Soonmyung Paik, Norman Wolmark, Yusuke Nakamura, Liewei Wang, Richard M Weinshilboum, Fergus J Couch. The contribution of common genetic variation to breast cancer risk among women receiving tamoxifen or raloxifene within the National Surgical Adjuvant Breast and Bowel Project (NSABP) P-1 and P-2 trials [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 P6-10-03.

Abstract INTRODUCTION Patient derived xenografts (PDX) may better reflect individual patient (pt) tumor biology; however, the feasibility of collecting PDX from percutaneous tumor biopsies (PTB) in the neoadjuvant setting is unknown. Furthermore, drug response phenotypes observed in PDX have not been prospectively compared to the corresponding pt clinical outcomes. METHODS The Breast Cancer Genome Guided Therapy Study (BEAUTY) is a prospective Mayo study of pts with high-risk breast cancer treated with neoadjuvant weekly paclitaxel (T) +/- trastuzumab followed by anthracycline based chemotherapy. PTB (at baseline) and residual surgical tissue (after all chemotherapy) are obtained for next generation sequencing (NGS) and PDX. Tumor biopsies (1-2 cores from 14 gauge needle) were implanted with Matrigel <1 hour of collection in the flanks of NOD-SCID or NSG mice. Low dose estradiol was supplemented in the drinking water. Primary outgrowth rate was defined as PDX tumor volume >50 mm3. Take rate was defined as development of at least 1 stably transplantable xenograft line/pt. To determine whether clinical T response assessed by MRI corresponded with in vivo T response, pretreatment PDX from 5 pts were injected into NOD-SCID mice (20 mice per pt PDX) and when tumors reached 100-200mm3, mice were randomized to no treatment vs T (20 mg/kg, ip. every 3-4 days). Two of these 5 patients had a MRI response defined as >30% decrease in longest lesion. RESULTS Pretreatment PTB from 81 unique pts were implanted in 251 mice (2-4 mice/pt). PDX outgrowth rates were 33.3% (27/81 pts) and 22 stable PDX were established (overall take rate 27.2%). Take rates were as follows: triple negative breast cancer (46%; 13/28); HER2 (27%; 6/22), Luminal B (13%; 3/22), and luminal A (0%; 0/9). Residual surgical tumor (after all treatment) from 17 pts was injected into 85 mice (average 5 mice/pt) and the initial outgrowth rate was 23% (4/17) with 3 stably transplantable lines established. PDX, derived from pretreatment PTB of 5 pts (2 responders and 3 non-responders), were assessed for in vivo T response. The size of the T treated group was significantly smaller than the no treatment group for the PDX derived from the 2 clinical responders, with complete disappearance of tumor by 18 days. In contrast, the PDX derived from the 3 clinical non-responders had no evidence for T response. CONCLUSIONS We have demonstrated the feasibility of using PTB to establish PDX in a prospective neoadjuvant clinical study and have demonstrated similar T drug response phenotypes in in the PDX as seen in the corresponding pt. These data suggest that PDX generated prospectively may be useful for biomarker validation and the development and individualization of new drug therapy. Funded by the Mayo Clinic Center for Individualized Medicine and the Mayo Clinic Cancer Center Citation Format: Jia Yu, Ping Yin, Bowen Gao, Jason P. Sinnwell, Ann M. Moyer, Daniel W. Visscher, Amy L. Conners, Travis J. Dockter, Krishna R. Kalari, Xiaojia Tang, Kevin J. Thompson, Hugues Sicotte, Douglas W. Mahoney, Steven N. Hart, Peter T. Vedell, Poulami Barman, Katie N. Jones, Sarah A. McLaughlin, John A. Copland, Alvaro Moreno Aspitia, Donald W. Northfelt, Richard J. Gray, Vera J. Suman, Jeanette E. Eckel Passow, Eric D. Wieben, James N. Ingle, Zhenkun Lou, Gianrico Farrugia, Richard Weinshilboum, Matthew P. Goetz, Judy C. Boughey, Liewei Wang. Feasibility of using percutaneous tumor biopsies from a prospective neoadjuvant breast cancer study to develop patient derived xenografts and assess in vivo chemotherapy sensitivity. [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 1195. doi:10.1158/1538-7445.AM2014-1195

Abstract Rapid development of next generation sequencing technology has enabled the identification of genomic alterations from short sequencing reads. There are a number of software pipelines available for calling single nucleotide variants from genomic DNA but, no comprehensive pipelines to identify, annotate and prioritize expressed SNVs (eSNVs) from non-directional paired-end RNA-Seq data. We have developed the eSNV-Detect, a novel computational system, which utilizes data from multiple aligners to call, even at low read depths, and rank variants from RNA-Seq. Multi-platform comparisons with the eSNV-Detect variant candidates were performed. The method was first applied to RNA-Seq from a lymphoblastoid cell-line, achieving 99.7% precision and 91.0% sensitivity in the expressed SNPs for the matching HumanOmni2.5 BeadChip data. Comparison of RNA-Seq eSNV candidates from 25 ER+ breast tumors from The Cancer Genome Atlas (TCGA) project with whole exome coding data showed 90.6–96.8% precision and 91.6–95.7% sensitivity. Contrasting single-cell mRNA-Seq variants with matching traditional multicellular RNA-Seq data for the MD-MB231 breast cancer cell-line delineated variant heterogeneity among the single-cells. Further, Sanger sequencing validation was performed for an ER+ breast tumor with paired normal adjacent tissue validating 29 out of 31 candidate eSNVs. The source code and user manuals of the eSNV-Detect pipeline for Sun Grid Engine and virtual machine are available at http://bioinformaticstools.mayo.edu/research/esnv-detect/.

Abstract BACKGROUND Metformin, a widely used anti-diabetic drug, is being considered as a highly promising agent for treatment and prevention of many types of cancer, including breast cancer. Although the exact mechanism of action of metformin is unknown, it is thought that metformin is an AMPK activator. To better individualize metformin therapy and explore additional underlying molecular mechanism associated with metformin action, we conducted a metformin pharmacogenomic study using 266 genomic data-enriched lymphoblastoid cell lines (LCLs). METHODS Genome-wide microarray data were generated for all LCLs, including 1.3 million SNPs, 485K DNA methylation probes and 54K mRNA expression probe sets. Metformin cytotoxicity assay was performed using MTS assay. Genome-wide association (GWA) analyses and integrated analyses using all the genomic data and metformin cytotoxicity (IC50 value) were performed to identify top candidate genes for functional validation in breast cancer cell lines, MDA-MB 231 and Hs578T, using siRNA screening approach. Further mechanistic study was performed to investigate the role of candidates in regulation of AMPK activity. RESULTS GWA analyses in LCLs identified 198 mRNA expression probe sets and 210 DNA methylation probes associated with metformin IC50 with p-value <10-4, respectively. Analysis of 1.3 million SNPs found 12 loci (a region containing at least 1 SNP of p-value <10-5 and 1 SNP of p-value <10-3 within 50kb) that were associated with metformin IC50. Integrated SNP loci-mRNA expression-IC50 analysis indicated that the SNP loci on chromosome 16 was associated with metformin IC50 through a trans-regulation of expression of 7 genes with p-value <10-4. Integrated methylation-mRNA expression-IC50 analysis showed that 15 DNA methylation probes in 6 genes were associated with metformin IC50 through both cis- and trans-regulation of expression of 48 genes with p-value <10-4, including 4 methylation probes in 3 of those 6 genes that were cis-correlated with its own gene expression. Ingenuity pathway analysis of top 62 candidate genes identified four major pathways and 55 genes within those pathways were selected for functional validation in breast cancer cell lines. The knockdown experiment showed that down regulation of 25 genes significantly altered metformin sensitivity in MDA-MB231, and 14 of them also showed the same effect in Hs578T. Preliminary mechanistic experiment indicated that knockdown of 7 genes significantly reduced AMPK activity in MDA-MB231 and 3 of them had same effect in Hs578T. CONCLUSIONS GWAS using a genomic data-enriched LCL model system, together with functional validation in breast cancer cell lines, could help us to identify novel genetic and epigenetic biomarkers involved in metformin response and help us to better understand the mechanisms of metformin in cancer treatment. Citation Format: Nifang Niu, Xianglin Tan, Brooke L. Fridley, Daniel J. Schaid, Ryan P. Abo, Anthony Batzler, Erin E. Carlson, Gregory Jenkins, Sebastian Moran, Holger Andrea Heyn, Manel Esteller Badosa, Liewei Wang. Metformin pharmacogenomics: A genome-wide associate study to identify genetic and epigenetic biomarkers involved in metformin response. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 2271. doi:10.1158/1538-7445.AM2013-2271

Abstract Definitive radiotherapy combined with chemotherapy is associated with enhanced locoregional control and clinical outcomes in non-small cell lung cancer (NSCLC) patients who are not candidates for surgery. However, radiation-induced pneumonitis and esophagitis are common dose-limiting inflammatory conditions that occur in these patients. These events can compromise the effect of the overall treatment and reduce quality of life in these patients. It is not clear what factors contribute to differences in the development of toxicities among patients with same clinical or dosimetric parameters. In this study, we conducted a large scale pathway-based analysis to identify inflammation-related genetic variants associated with radiation-induced pneumonitis or esophagitis. In the initial discovery phase, 11,930 inflammation-related genetic variants were genotyped in 243 stage I-III NSCLC patients who were treated with definitive radiotherapy at MD Anderson Cancer Center. A fast-track validation was performed for the top candidate SNPs in an additional 212 NSCLC cases following the same criteria. A total of 14 SNPs and 12 SNPs were validated for pneumonitis and esophagitis, respectively. Meta-analysis was used to summarize the results from two phases. The most significant SNPs for esophagitis and pneumonitis mapped to chromosome 18p (OR=9.05, 95%CI=2.68-30.59, P-meta=3.90×10−4) and 14q (OR=4.75, 95%CI=2.17-10.41, P-meta =9.61×10−5), respectively. To explore the potential underlying mechanisms, we analyzed radiation sensitivity and eQTL relationships of the identified SNPs and their host genes in 277 ethnically defined human lymphoblastoid cell lines. Among the 24 host genes for the 26 identified SNPs, three genes (DDX58, GZMB, CDC2) showed significant association with radiation responses in these cell lines (P<0.05). Moreover, a significant eQTL relationship was identified between CDC2 expression and the genotype of the identified CDC2 SNP (P=0.022). Interestingly, this CDC2 SNP is predicted to create a new miRNA binding site. Luciferase reporter assays to test the influence of this SNP on miRNA regulation is currently ongoing. Our results provide evidence that inflammation-related genetic variants contribute to the development of radiation-induced pneumonitis and esophagitis. The results also suggest that gene expression changes through altered miRNA binding may be an underlying mechanism for these events. Citation Format: Xia Pu, Liewei Wang, Michelle AT Hildebrandt, Yuanqing Ye, Joe Y. Chang, Charles Lu, Heath D. Skinner, Nifang Niu, Gregory D. Jenkins, Ritsuko Komaki, John D. Minna, Jack A. Roth, Richard M. Weinshilboum, Xifeng Wu. Inflammation-related genetic variants predict radiation-induced toxicity following definitive radiotherapy for non-small cell lung cancer. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 3626. doi:10.1158/1538-7445.AM2013-3626

Abstract The FKBP5 gene encodes the FKBP51 immunophilin, a negative regulator of AKT. Variability in level of FKBP5 expression is a major factor contributing to variation in response to several classes of chemotherapeutic agents including gemcitabine, first line treatment for pancreatic cancer. Genetic variation of FKBP5 could influence its expression, function, and, ultimately, treatment response of pancreatic cancer. We set out to comprehensively study the role of genetic variation in FKBP5 identified by Next Generation DNA resequencing in response to gemcitabine treatment of pancreatic cancer by utilizing both tumor and germline DNA samples from 43 pancreatic cancer patients, including 19 paired tumor-normal samples. Genotype-phenotype association studies with overall survival (in tumor only) as well as with FKBP5 gene expression were performed using the same samples in which Next Generation DNA resequencing had been performed, followed by functional genomic studies of top selected single nucleotide polymorphisms (SNPs). In-depth resequencing identified 404 FKBP5 SNPs in both normal and tumor DNA patient samples. SNPs with the smallest p-values for association with FKBP5 survival or expression were subjected to functional genomic studies. Electromobility shift assay (EMSA) showed that the rs73748206 “T” allele altered DNA-protein binding patterns, consistent with a significantly increased reporter gene activity. The effect of rs73748206 was also confirmed by mRNA expression studies using lymphoblastoid cell lines (LCL) derived from the same patients from whom DNA used for resequencing had been obtained. This comprehensive FKBP5 resequencing study provides insight into the role of inheritance in variation in response to gemcitabine therapy of pancreatic cancer. Citation Format: Katarzyna A. Ellsworth, Bruce W. Eckloff, Liang Li, Irene Moon, Brooke L. Fridley, Greg D. Jenkins, Erin Carlson, Abra Brisbin, Ryan Abo, William Bamlet, Gloria Petersen, Eric D. Wieben, Liewei Wang. Contribution of FKBP5 genetic variation to gemcitabine treatment and survival in pancreatic cancer. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 2209. doi:10.1158/1538-7445.AM2013-2209