Kohli, Manish

Abstract Purpose: Intratumoral androgen synthesis (IAS) is a key mechanism promoting androgen receptor (AR) reactivation and antiandrogen resistance in castration-resistant prostate cancer (CRPC). However, signaling pathways driving aberrant IAS remain poorly understood. Experimental Design: The effect of components of the AKT-RUNX2-osteocalcin (OCN)–GPRC6A–CREB signaling axis on expression of steroidogenesis genes CYP11A1 and CYP17A1 and testosterone level were examined in PTEN-null human prostate cancer cell lines. Pten knockout mice were used to examine the effect of Runx2 heterozygous deletion or abiraterone acetate (ABA), a prodrug of the CYP17A1 inhibitor abiraterone on Cyp11a1 and Cyp17a1 expression, testosterone level and tumor microenvironment (TME) remodeling in vivo. Results: We uncovered that activation of the AKT–RUNX2–OCN–GPRC6A–CREB signaling axis induced expression of CYP11A1 and CYP17A1 and testosterone production in PTEN-null prostate cancer cell lines in culture. Deletion of Runx2 in Pten homozygous knockout prostate tumors decreased Cyp11a1 and Cyp17a1 expression, testosterone level, and tumor growth in castrated mice. ABA treatment also inhibited testosterone synthesis and alleviated Pten loss-induced tumorigenesis in vivo. Pten deletion induced TME remodeling, but Runx2 heterozygous deletion or ABA treatment reversed the effect of Pten loss by decreasing expression of the collagenase Mmp9. Conclusions: Abnormal RUNX2 activation plays a pivotal role in PTEN loss-induced IAS and TME remodeling, suggesting that the identified signaling cascade represents a viable target for effective treatment of PTEN-null prostate cancer, including CRPC. Clin Cancer Res; 24(4); 834–46. ©2017 AACR.

The Memorial Sloan Kettering Cancer Center (MSKCC) prognostic score is based on clinical parameters. We analyzed whole blood mRNA expression in metastatic clear cell renal cell carcinoma (mCCRCC) patients and compared it to the MSKCC score for predicting overall survival. In a discovery set of 19 patients with mRCC, we performed whole transcriptome RNA sequencing and selected eighteen candidate genes for further evaluation based on associations with overall survival and statistical significance. In an independent validation of set of 47 patients with mCCRCC, transcript expression of the 18 candidate genes were quantified using a customized NanoString probeset. Cox regression multivariate analysis confirmed that two of the candidate genes were significantly associated with overall survival. Higher expression of BAG1 [hazard ratio (HR) of 0.14, p < 0.0001, 95% confidence interval (CI) 0.04–0.36] and NOP56 (HR 0.13, p < 0.0001, 95% CI 0.05–0.34) were associated with better prognosis. A prognostic model incorporating expression of BAG1 and NOP56 into the MSKCC score improved prognostication significantly over a model using the MSKCC prognostic score only (p < 0.0001). Prognostic value of using whole blood mRNA gene profiling in mCCRCC is feasible and should be prospectively confirmed in larger studies.

The testis‐specific Y‐encoded‐like protein (TSPYL) gene family includes TSPYL1 to TSPYL6. We previously reported that TSPYL5 regulates cytochrome P450 (CYP) 19A1 expression. Here we show that TSPYLs, especially TSPYL 1, 2, and 4, can regulate the expression of many CYP genes, including CYP17A1, a key enzyme in androgen biosynthesis, and CYP3A4, an enzyme that catalyzes the metabolism of abiraterone, a CYP17 inhibitor. Furthermore, a common TSPYL1 single nucleotide polymorphism (SNP), rs3828743 (G/A) (Pro62Ser), abolishes TSPYL1's ability to suppress CYP3A4 expression, resulting in reduced abiraterone concentrations and increased cell proliferation. Data from a prospective clinical trial of 87 metastatic castration‐resistant prostate cancer patients treated with abiraterone acetate/prednisone showed that the variant SNP genotype (A) was significantly associated with worse response and progression‐free survival. In summary, TSPYL genes are novel CYP gene transcription regulators, and genetic alteration within these genes significantly influences response to drug therapy through transcriptional regulation of CYP450 genes.

Abstract Prostate cancer is the most frequently diagnosed male cancer. For prostate cancer that has progressed to an advanced or metastatic stage, androgen deprivation therapy (ADT) is the standard of care. ADT inhibits activity of the androgen receptor (AR), a master regulator transcription factor in normal and cancerous prostate cells. The major limitation of ADT is the development of castration-resistant prostate cancer (CRPC), which is almost invariably due to transcriptional re-activation of the AR. One mechanism of AR transcriptional re-activation is expression of AR-V7, a truncated, constitutively active AR variant (AR-V) arising from alternative AR pre-mRNA splicing. Noteworthy, AR-V7 is being developed as a predictive biomarker of primary resistance to androgen receptor (AR)-targeted therapies in CRPC. Multiple additional AR-V species are expressed in clinical CRPC, but the extent to which these may be co-expressed with AR-V7 or predict resistance is not known. Here we utilized long read sequencing to identify and quantify AR isoforms expressed in CRPC. To unambiguously characterize all AR isoforms, we prepared Iso-Seq™ libraries via 3’ rapid amplification of cDNA ends (RACE) with RNA isolated from prostate cancer cell lines and xenograft tissues using a forward primer anchored in AR exon 1. 3’RACE reactions were subjected to single molecule, real-time (SMRT®) long-read sequencing with a Pacific Biosciences RSII System. Our work identified AR-V9 as a truncated isoform that is frequently co-expressed with AR-V7 in CRPC. Mechanistically, our work re-annotated AR-V7 and AR-V9 mRNAs, showing these two species shared a common 3’ terminal exon containing separate splice acceptor sites. Taking into account this new information, novel siRNAs and antibodies which could distinguish between AR-V7 and AR-V9 were designed, validated and used to measure the relative expression of these two AR isoforms in CRPC cells with a view to determining the potential of AR-V9 as a predictive biomarker of primary resistance to AR-targeted therapies. Citation Format: Manish Kohli, Yeung Ho, David W. Hillman, Jamie L. Van Etten, Christine Henzler, Rendong Yang, Yingming Li, Elizabeth Tseng, Ting Hon, Tyson A. Clark, Liguo Wang, Kevin Silverstein, Liewei Wang, Scott M. Dehm. SMRT® Sequencing of full-length androgen receptor isoforms in prostate cancer reveals previously hidden drug resistant variants [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 425. doi:10.1158/1538-7445.AM2017-425

Abstract Background: Circulating cell-free DNA (cfDNA) analysis is emerging as a less invasive approach to assessing tumor genomic alterations in cancer patients. Although high concordance has been reported between tumor tissue NGS and cfDNA in studies investigating specific genetic alterations, the fidelity of cfDNA to tumor tissue DNA in the global genomic scale is largely unknown. In a correlative study of a prospective clinical trial (NCT# 01953640) conducted in mCRPC stage patients treated with abiraterone acetate/prednisone (A/P), we evaluated correlation of genomic CNA in tumor DNA obtained from biopsy of metastatic lesions and matched plasma cfDNA. Methods: mCRPC patients (pts) underwent two image-guided core biopsies of metastases prior to initiation of A/P (visit 1, V1) and 12 weeks after treatment (visit 2, V2). At both time points plasma was obtained at the time of the core biopsies and cfDNA was extracted using DNA Blood Mini Kit (Qiagen, Valencia, CA). High coverage (for tumor tissue) and low coverage (for cfDNA) whole genome sequencing reads were first mapped to the human genome hg19. Read counts (RC) from the mapped sequence files were then binned into 1Mb windows. The RC ratio in each genomic bin was calculated by comparing tumor tissue DNA to lymphocyte gDNA derived from the same patient, and was further log2 transformed, corrected for GC content, and normalized by CGHnormaliter. The fully normalized log2 ratios data was subjected to segmentation using DNAcopy algorithm. Results: Between 05/2013 and 09/2015, 92 patients (pts) were enrolled of which tissue and plasma NGS data both visits was available for 18 pts. The correlation of CNAs between tumor tissue and its cfDNA counterpart ranges from 0.013 to 0.83 for V1 samples, and -0.05 to 0.9 for V2 samples. The decreased correlation in some pairs of samples is largely due to low tumor content and heterogeneity in the cfDNA. Although there is a wide range of correlation, commonly shared CNAs were identified in multiple chromosomal regions, including loss in 8p, gain in 8q and chromosome 5 and X. On the other hand, several genomic regions show inconsistent CNAs between tumor tissue and cfDNA among 15 pairs of samples in which cfDNA have large tumor content. These include: loss in 6q, 1p, 2p, 9q, 11q, 13p, 15p, 20p and gain in 8q and 10p in cfDNA but not in tumor tissue; marked loss in chromosome 22 and gain in p arm of X chromosome in tumor tissue but much less evident in cfDNA. Conclusion: High concordance of CNAs between cfDNA and tumor tissue DNA can be achieved given sufficient tumor content of cfDNA. However, more CNAs can be identified in cfDNA than its tumor tissue counterpart. Our result suggests that cfDNA NGS is a useful tool to investigate clonal evolution associated with cancer progression. Citation Format: Chiang-Ching Huang, Meijun Du, Liguo Wang, Yen-chen Lin, Hua Huang, Liewei Wang, Liang Wang, Manish Kohli. Comparison of copy number aberrations (CNAs) between plasma cell free DNA (cfDNA) and tissue DNA in metastatic castrate resistant prostate cancer (mCRPC) [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 2733. doi:10.1158/1538-7445.AM2017-2733

Abstract Purpose: Androgen receptor (AR) variant AR-V7 is a ligand-independent transcription factor that promotes prostate cancer resistance to AR-targeted therapies. Accordingly, efforts are under way to develop strategies for monitoring and inhibiting AR-V7 in castration-resistant prostate cancer (CRPC). The purpose of this study was to understand whether other AR variants may be coexpressed with AR-V7 and promote resistance to AR-targeted therapies. Experimental Design: We utilized complementary short- and long-read sequencing of intact AR mRNA isoforms to characterize AR expression in CRPC models. Coexpression of AR-V7 and AR-V9 mRNA in CRPC metastases and circulating tumor cells was assessed by RNA-seq and RT-PCR, respectively. Expression of AR-V9 protein in CRPC models was evaluated with polyclonal antisera. Multivariate analysis was performed to test whether AR variant mRNA expression in metastatic tissues was associated with a 12-week progression-free survival endpoint in a prospective clinical trial of 78 CRPC-stage patients initiating therapy with the androgen synthesis inhibitor, abiraterone acetate. Results: AR-V9 was frequently coexpressed with AR-V7. Both AR variant species were found to share a common 3′ terminal cryptic exon, which rendered AR-V9 susceptible to experimental manipulations that were previously thought to target AR-V7 uniquely. AR-V9 promoted ligand-independent growth of prostate cancer cells. High AR-V9 mRNA expression in CRPC metastases was predictive of primary resistance to abiraterone acetate (HR = 4.0; 95% confidence interval, 1.31–12.2; P = 0.02). Conclusions: AR-V9 may be an important component of therapeutic resistance in CRPC. Clin Cancer Res; 23(16); 4704–15. ©2017 AACR.

500 Background: The Memorial Sloan Kettering Cancer Center (MSKCC) prognostic score is based on clinical parameters and has several limitations. We analyzed whole blood mRNA expression in mCCRCC patients (pts) and compared it to MSKCC for predicting overall survival (OS). Methods: A prospective collection of whole blood in PAXgene RNA tubes from mCCRCC pts was performed between 2011 and 2015 with uniform processing. A discovery set (DS) of 19 randomly selected pt samples was used to generate whole transcriptome RNA sequencing and signatures describing differential gene expression patterns associated with OS. We calculated the fragment per kilobase of transcript per million (FPKM) of each transcript and gene to identify correlations between the differentially expressed genes and MSKCC score. Candidate genes were selected based on statistical significance (p<0.05) and associations with OS. In an independent validation set (VS) of 47 pts with mCCRCC, we used the nCounter Analysis System (NanoString Technologies, Seattle, WA, USA) to quantify gene transcript expression selected from the DS. Cox regression and Kaplan-Meier analysis were used for association of candidate mRNAs in the VS with OS (time of metastatic recurrence to death or last follow-up). We compared our gene signature to the MSKCC prognostic score for poor versus good survival. Results: In the DS, at a median follow-up of 4.13 years (yrs) (range: 1.3-4.5), 11 patients had died. We identified eighteen genes in the DS with statistically significant differential gene expression patterns and ability to discriminate patient OS. In the VS at a median follow-up of 4.9 yrs (range: 0.1-11.3), twenty-five patients had died. We confirmed two of the eighteen selected genes that had strong prognostic significance: BAG1 [hazard ratio (HR) 0.31, p=0.007] and NOP56 [HR 0.28, p=0.004]. Comparing OS between good and poor outcome groups for these 2 scores, our two-gene panel had higher prognostic ability (HR 20.07, p <0.0001) than the MSKCC prognostic score (HR 3.19, p 0.04). Conclusions: Prognostic value of using whole blood mRNA gene profiling in mCCRCC is feasible and should be prospectively confirmed in larger studies.

181 Background: To evaluate whether metformin (Met) a widely-used, nontoxic oral antidiabetic drug with putative anticancer properties leads to improvements in prostate cancer (PC) outcomes in the CHAARTED trial. Methods: In the CHAARTED database where metformin use at baseline was recorded prospectively, we identified patients with metastatic PC who underwent either ADT alone or ADT and docetaxel (D) chemotherapy. Cox proportional hazards models were used to determine the effect of Metformin on outcomes. Results: A total of 788 patients (median age, 63 y) had complete data after randomization. Comparison of ADT+D+Met (n = 39) to ADT+D (n = 357) and ADT+Met (n = 29) to ADT alone (n = 363) revealed similar clinicopathologic characteristics. Cause of death was PC in 13(81%) of ADT+D+Met, 72(85%) ADT+D, 9(82%) ADT+Met and 105(84%) ADT alone groups. See table for PC outcomes and overall survival by metformin use. Cox regression analysis for overall survival stratified by stratification factors at randomization demonstrates Met use was associated with a trend for worse overall survival (HR 1.47 95%CI: [0.95,2.26], p = 0.08) with adjustment for treatment arm and prior local therapy. In contrast, ADT+D use (HR 0.62; 95%CI: [0.47,0.81]) and prior local therapy with surgery or radiation (HR 0.56; 95% CI: [0.38, 0.82]) were associated with improved survival. Conclusions: In this study, baseline metformin did not improve PC outcomes. Partial support and drug supply by Sanofi. Clinical trial information: NCT00309985. [Table: see text]

175 Background: Genome and transcriptome aberrations associated with primary resistance to abiraterone acetate/prednisone (AA/P) in mCRPC are not known. In a prospective trial (NCT#01953640) we assessed whole-exome and RNA-seq based aberrations in metastases of CRPC stage patients (pts) for identifying markers associated with primary resistance to AA/P. Methods: Whole-exome and RNA-seq of biopsies from metastases was performed followed by analyses for association between genomic aberrations & primary resistance. Primary resistance was defined by progression on AA/P after 12-weeks of therapy(non-responders) using PSA, RECIST, bone scan and symptom criteria (per PCWG2). Gene network analysis was performed in genes mutated more frequently in non-responders, and also in genes that were differentially expressed between non-responders and responders and a “risk ratio” (RR) was calculated thereof. Results: Between 6/2013 & 8/2015, 92 pts were enrolled of which 82 had complete whole-exome, RNA-seq and 12-week outcome data available for analysis. Median age was 72.5 yrs (IQR: 68.5-78); median PSA-18 ng/ml (IQR: 8.1- 46.6). At 12-weeks 33/82 had progressed. Using the RR of 2 as threshold, we identified 113 and 292 genes that were more frequently mutated in non-responders & responders respectively. Among the 113 genes, OBSCN, ADAM21, LPHN3, DOCK10 ( P = 0.08, RR= inf) & USP42 ( P = 0.16, RR = 5.71) were the top candidates.Gene network analysis revealed that in non-responders, genes involved in the Wnt/β-catenin pathway were frequently mutated and negative regulators of Wnt pathway ( DKK4, SFRP2 & LRP6) were also frequently deleted. Gene expression analyses revealed the expression levels of Wnt/β-catenin pathway inhibitors were significantly reduced while expression levels of cell cycle regulatory genes were significantly increased in non-responders. Conclusions: In this study we observed Wnt/β-catenin pathway activation to be associated with primary AA/P resistance. This finding offers a potential for the development of predictive biomarkers and modulation of targeted pathways to overcome AA/P resistance. Clinical trial information: NCT# 01953640.