An Atlas of Accessible Chromatin in Advanced Prostate Cancer Reveals the Epigenetic Evolution during Tumor Progression.

Metastatic castration-resistant prostate cancer (mCRPC) is a lethal disease that resists therapy targeting androgen signaling, the primary driver of prostate cancer. mCRPC resists androgen receptor (AR) inhibitors by amplifying AR signaling or by evolving into therapy-resistant subtypes that do not depend on AR. Elucidation of the epigenetic underpinnings of these subtypes could provide important insights into the drivers of therapy resistance. In this study, we produced chromatin accessibility maps linked to the binding of lineage-specific transcription factors (TF) by performing ATAC sequencing on 70 mCRPC tissue biopsies integrated with transcriptome and whole genome sequencing. mCRPC had a distinct global chromatin accessibility profile linked to AR function. Analysis of TF occupancy across accessible chromatin revealed 203 TFs associated with mCRPC subtypes. Notably, ZNF263 was identified as a putative prostate cancer TF with a significant impact on gene activity in the double-negative (AR- neuroendocrine-) subtype, potentially activating MYC targets. Overall, this analysis of chromatin accessibility in mCRPC provides valuable insights into epigenetic changes that occur during progression to mCRPC.

Integrated analyses highlight interactions between the three-dimensional genome and DNA, RNA and epigenomic alterations in metastatic prostate cancer.

The impact of variations in the three-dimensional structure of the genome has been recognized, but solid cancer tissue studies are limited. Here, we performed integrated deep Hi-C sequencing with matched whole-genome sequencing, whole-genome bisulfite sequencing, 5-hydroxymethylcytosine (5hmC) sequencing and RNA sequencing across a cohort of 80 biopsy samples from patients with metastatic castration-resistant prostate cancer. Dramatic differences were present in gene expression, 5-methylcytosine/5hmC methylation and in structural variation versus mutation rate between A and B (open and closed) chromatin compartments. A subset of tumors exhibited depleted regional chromatin contacts at the AR locus, linked to extrachromosomal circular DNA (ecDNA) and worse response to AR signaling inhibitors. We also identified topological subtypes associated with stark differences in methylation structure, gene expression and prognosis. Our data suggested that DNA interactions may predispose to structural variant formation, exemplified by the recurrent TMPRSS2-ERG fusion. This comprehensive integrated sequencing effort represents a unique clinical tumor resource.

A prostate cancer gastrointestinal transcriptional phenotype may be associated with diminished response to AR-targeted therapy.

Background: Prostate cancer is a heterogenous disease, but once it becomes metastatic it eventually becomes treatment resistant. One mechanism of resistance to AR-targeting therapy is lineage plasticity, where the tumor undergoes a transformation to an AR-indifferent phenotype, most studied in the context of neuroendocrine prostate cancer (NEPC). However, activation of additional de- or trans-differentiation programs, including a gastrointestinal (GI) gene expression program, has been suggested as an alternative method of resistance. In this study, we explored the previously identified GI prostate cancer phenotype (PCa-GI) in a large cohort of metastatic castration-resistant prostate cancer (mCRPC) patient biopsy samples. Methods: We analyzed a dataset of 634 mCRPC samples with batch effect corrected gene expression data from the West Coast Dream Team (WCDT), the East Coast Dream Team (ECDT), the Fred Hutchinson Cancer Research Center (FHCRC) and the Weill Cornell Medical center (WCM). Survival data was available from the WCDT and ECDT cohorts. We calculated a gene expression GI score using the sum of z-scores of genes from a published set of PCa-GI-defining genes (N=38). Survival analysis was performed using the Kaplan-Meier method and Cox proportional hazards regression with endpoint overall survival from time of biopsy to death of any cause. Results: We found that the PCa-GI score had a bimodal distribution, identifying a distinct set of tumors with an activated GI expression pattern. Approximately 35% of samples were classified as PCa-GI high, which was concordant with prior reports. Liver metastases had the highest median score but after excluding liver samples, 29% of the remaining samples were still classified as PCa-GI high, suggesting a distinct phenotype not exclusive to liver metastases. No correlation was observed between GI score and proliferation, AR signaling, or NEPC scores. Furthermore, the PCa-GI score was not associated with genomic alterations in AR, FOXA1, RB1, TP53 or PTEN. However, tumors with MYC amplifications showed significantly higher GI scores (p=0.0001). Patients with PCa-GI tumors had a shorter survival (HR=1.5 [1.1-2.1], p=0.02), but this result was not significant after adjusting for the liver as metastatic site (HR=1.2 [0.82-1.7], p=0.35). Patients with PCa-GI low samples had a better outcome after androgen receptor signaling inhibitors (ASI, abiraterone or enzalutamide) than other therapies (HR=0.37 [0.22-0.61], p=0.0001) while the benefit of ASI was smaller and non-significant for PCa-GI high samples (HR=0.55 [0.29-1.1], p=0.07). A differential pathway analysis identified FOXA2 signaling to be upregulated PCa-GI high tumors (FDR = 3.7 × 10-13). Conclusions: The PCa-GI phenotype is prevalent in clinical mCRPC samples and may represent a distinct biological entity. PCa-GI tumors may respond less to ASI and could offer a strategy to study novel therapeutic targets.