Kalari, Krishna R.

Abstract BACKGROUND: Based on prospective clinical trials, there is no evidence for differences in efficacy between the 3 aromatase inhibitors (AIs) anastrozole, exemestane, and letrozole. The purpose of this study was to identify germline genetic variants associated with response to AIs and to help identify novel mechanisms associated with drug disease efficacy. METHODS: A genome-wide association study (GWAS) was performed for 624 patients (Steroids 2015;99:32-38) to identify SNPs associated with estrogen level change in women with estrogen receptor (ER) positive breast cancer treated with anastrozole. Replication of associated SNPs was performed in a GWAS from the MA.27 trial that compared adjuvant anastrozole and exemestane treatment of post-menopausal women with ER+ breast cancer. Functional studies were subsequently performed to determine SNP effects and underlying mechanisms. RESULTS: Our initial GWAS identified SNPs within CSMD1 that were associated with changes in estrogen levels during anastrozole therapy. An additional SNP in CSMD1 was also associated with breast cancer events in CCTG MA.27. Functionally, we showed that CSMD1 regulates CYP19 expression in a SNP-, and in an anastrozole- dependent fashion. These phenomena were not observed for either letrozole or exemestane. In MA.27, an anastrozole- specific effect was also seen with the minor allele having a protective effect on time to distant metastasis (HR=0.49, p=0.00259), but this was not the case for exemestane (HR=0.71, p=0.111). Our in vitro functional studies indicated that overexpression of CSMD1 sensitized anastrozole or letrozole resistant cells to anastrozole but not to the other two AIs. The SNP in CSMD1 that was associated with increased CSMD1 and CYP19 expression levels increased anastrozole sensitivity, but not letrozole or exemestane in lymphoblastoid cell lines (LCLs) homozygous for either WT or variant CSMD1 SNP genotypes. Based on these observations, we explored whether anastrozole has additional mechanisms beyond its function as a CYP19 inhibitor. Utilizing an estrogen response element (ERE) luciferase reporter assay in a CYP19 CRISPR knockout breast cancer T47D cell line and a surface plasmon resonance (SPR) assay, we found that anastrozole can also function as an ERα agonist, and can bind to, and result in, proteasome dependent ERα degradation, especially in the presence of E2. Treatment of these CYP19 CRISPR knockout cells with anastrozole in the presence of increasing concentrations of E2 results in greater sensitivity compared with anastrozole alone, while the addition of E2, as expected, does not improve letrozole or exemestane sensitivity. These same observations were also seen in letrozole and anastrazole resistant cells. CONCLUSIONS: Our findings suggest that anastrozole might be more effective than letrozole or exemestane in patients with the CSMD1 SNP. Furthermore, anastrozole can function as an ERα agonist, binding to ERα and resulting in its degradation, especially in the presence of E2. These findings should help to make it possible to develop precision endocrine therapies for women who are candidates for AIs. Citation Format: Cairns J, Ingle J, Dudenkov T, Kalari K, Buzdar A, Kubo M, Robson M, Ellis M, Goss P, Shepherd L, Goetz M, Weinshilboum R, Wang L. CSMD1 SNPs selectively affect anastrozole response in postmenopausal breast cancer patients [abstract]. In: Proceedings of the 2016 San Antonio Breast Cancer Symposium; 2016 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2017;77(4 Suppl):Abstract nr PD1-04.

Abstract Genetic risks in breast cancer remain only partly understood. Here, we report the results of a genome-wide association study of germline DNA from 4,658 women, including 252 women experiencing a breast cancer recurrence, who were entered on the MA.27 adjuvant trial comparing the aromatase inhibitors (AI) anastrozole and exemestane. Single-nucleotide polymorphisms (SNP) of top significance were identified in the gene encoding MIR2052HG, a long noncoding RNA of unknown function. Heterozygous or homozygous individuals for variant alleles exhibited a ∼40% or ∼63% decrease, respectively, in the hazard of breast cancer recurrence relative to homozygous wild-type individuals. Functional genomic studies in lymphoblastoid cell lines and ERα-positive breast cancer cell lines showed that expression from MIR2052HG and the ESR1 gene encoding estrogen receptor-α (ERα) was induced by estrogen and AI in a SNP-dependent manner. Variant SNP genotypes exhibited increased ERα binding to estrogen response elements, relative to wild-type genotypes, a pattern that was reversed by AI treatment. Further, variant SNPs were associated with lower expression of MIR2052HG and ERα. RNAi-mediated silencing of MIR2052HG in breast cancer cell lines decreased ERα expression, cell proliferation, and anchorage-independent colony formation. Mechanistic investigations revealed that MIR2052HG sustained ERα levels both by promoting AKT/FOXO3-mediated ESR1 transcription and by limiting ubiquitin-mediated, proteasome-dependent degradation of ERα. Taken together, our results define MIR2052HS as a functionally polymorphic gene that affects risks of breast cancer recurrence in women treated with AI. More broadly, our results offer a pharmacogenomic basis to understand differences in the response of breast cancer patients to AI therapy. Cancer Res; 76(23); 7012–23. ©2016 AACR.

Abstract INTRODUCTION: The tumor microenvironment consists of a complicated interplay between a heterogeneous mixture of tumor cells and the surrounding tissue and there is a burgeoning interest in the role of the tissue microbiota within this microenvironment. Our group has investigated the breast tumor microbiome using the unaligned (non-human) reads from 796 RNASequencing data obtained from The Cancer Genome Atlas (TCGA) repository of Breast Cancer. METHODS: Unmapped sequences for 720 TCGA tumor samples and 76 paired normal samples were analyzed with Kraken v0.10.5 using the RefSeq 16S and standard genomes. To minimize false discoveries Kraken results were filtered to species that have bacterial genomic and 16S read support. Cluster analysis of the T-distributed Stochastic Neighbor Embedding (tsne) was performed and consensus voting of 26 maximum scoring metrics identified 3 distinct clusters. Differential abundance analysis was performed with the edgeR package. Covariate association analysis was performed for tumor clinical characteristics, as well as sample heterogeneity characteristics. RESULTS: There were three microbiome clusters we identified using TCGA breast samples. The average silhouette width of the third cluster was deemed statistically insignificant (0.071), suggesting further stratification of this cluster is necessary. The remaining two TCGA clusters (467) shared an average silhouette width of 0.39, using Bray's dissimilarity index. Differential abundance analysis of the two clusters after Bonferroni correction identified 450 significant species. The biomes were represented by 63 triple negative, 93 HER2+, 258 Luminal, and 50 tumor adjacent normals. Concordant cluster assignment was observed in 93.7 +/- 7.3% of the normal-tumor pairings, averaged across clinical subtypes. We observed clear delineation of biome populations at the order taxonomic level, which was not associated with the clinical subtype or the tumor/tumor adjacent pairings. We found enrichment of Burkholderiales in the first biome and enrichment of Bacillales and Lactobacillales within the second biome. The second biome appears to be deplete of representatives from the Asian population (p value 0.005) and survival associations were observed to be significantly reduced for HER2 patients associated with Biome 1 (p value 0.025). Geneset enrichment analysis of the expression signatures indicated drug metabolism and chemical carcinogenesis pathways to be associated with biome 2. CONCLUSIONS: In our study, we observed two significantly different microbiomes among TCGA breast cancer samples that were not associated with the clinical subtypes. Current efforts include validating these biomes, using 16S sequencing on matching fresh frozen paraffin embedded samples. Additionally, we are evaluating the influence of the biomes in independent adjuvant and neo-adjuvant studies with next generation sequencing data. Citation Format: Kevin J. Thompson, Jason Sinnwell, Xiaojia Tang, James N. Ingle, Matthew P. Goetz, Krishna R. Kalari. Distinct microbiome populations within breast cancer microenvironments. [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 3298.

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