The vast majority of disease-associated single nucleotide polymorphisms (SNP) identified from genome-wide association studies (GWAS) are localized in non-coding regions. A significant fraction of these variants impact transcription factors binding to enhancer elements and alter gene expression. To functionally interrogate the activity of such variants we developed snpSTARRseq, a high-throughput experimental method that can interrogate the functional impact of hundreds to thousands of non-coding variants on enhancer activity. snpSTARRseq dramatically improves signal-to-noise by utilizing a novel sequencing and bioinformatic approach that increases both insert size and the number of variants tested per loci. Using this strategy, we interrogated known prostate cancer (PCa) risk-associated loci and demonstrated that 35% of them harbor SNPs that significantly altered enhancer activity. Combining these results with chromosomal looping data we could identify interacting genes and provide a mechanism of action for 20 PCa GWAS risk regions. When benchmarked to orthogonal methods, snpSTARRseq showed a strong correlation with in vivo experimental allelic-imbalance studies whereas there was no correlation with predictive in silico approaches. Overall, snpSTARRseq provides an integrated experimental and computational framework to functionally test non-coding genetic variants.

Androgen receptor (AR) drives prostate cancer (PCa) development and progression. AR chromatin binding profiles are highly plastic and form recurrent programmatic changes that differentiate disease stages, subtypes and patient outcomes. While prior studies focused on concordance between patient subgroups, inter-tumor heterogeneity of AR enhancer selectivity remains unexplored. Here we report high levels of AR chromatin binding heterogeneity in human primary prostate tumors, that overlap with heterogeneity observed in healthy prostate epithelium. Such heterogeneity has functional consequences, as somatic mutations converge on commonly-shared AR sites in primary over metastatic tissues. In contrast, less-frequently shared AR sites associate strongly with AR-driven gene expression, while such heterogeneous AR enhancer usage also distinguishes patients outcome. These findings indicate that epigenetic heterogeneity in primary disease is directly informative for risk of biochemical relapse. Cumulatively, our results illustrate a high level of AR enhancer heterogeneity in primary PCa driving differential expression and clinical impact.

Many genetic variants affect disease risk by altering context-dependent gene regulation. Such variants are difficult to study mechanistically using current methods that link genetic variation to steady-state gene expression levels, such as expression quantitative trait loci (eQTLs). To address this challenge, we developed the cistrome-wide association study (CWAS), a framework for identifying genotypic and allele-specific effects on chromatin that are also associated with disease. In prostate cancer, CWAS identified regulatory elements and androgen receptor-binding sites that explained the association at 52 of 98 known prostate cancer risk loci and discovered 17 additional risk loci. CWAS implicated key developmental transcription factors in prostate cancer risk that are overlooked by eQTL-based approaches due to context-dependent gene regulation. We experimentally validated associations and demonstrated the extensibility of CWAS to additional epigenomic datasets and phenotypes, including response to prostate cancer treatment. CWAS is a powerful and biologically interpretable paradigm for studying variants that influence traits by affecting transcriptional regulation.

BACKGROUNDAndrogen receptor signaling inhibitors (ARSIs) have improved outcomes for patients with metastatic castration-resistant prostate cancer (mCRPC), but their clinical benefit is limited by treatment resistance.METHODSTo investigate the mechanisms of ARSI resistance, we analyzed the whole-genome (n = 45) and transcriptome (n = 31) sequencing data generated from paired metastatic biopsies obtained before initiation of first-line ARSI therapy for mCRPC and after radiographic disease progression. We investigated the effects of genetic and pharmacologic modulation of SSTR1 in 22Rv1 cells, a representative mCRPC cell line.RESULTSWe confirmed the predominant role of tumor genetic alterations converging on augmenting androgen receptor (AR) signaling and the increased transcriptional heterogeneity and lineage plasticity during the emergence of ARSI resistance. We further identified amplifications involving a putative enhancer downstream of the AR and transcriptional downregulation of SSTR1, encoding somatostatin receptor 1, in ARSI-resistant tumors. We found that patients with SSTR1-low mCRPC tumors derived less benefit from subsequent ARSI therapy in a retrospective cohort. We showed that SSTR1 was antiproliferative in 22Rv1 cells and that the FDA-approved drug pasireotide suppressed 22Rv1 cell proliferation.CONCLUSIONOur findings expand the knowledge of ARSI resistance and point out actionable next steps, exemplified by potentially targeting SSTR1, to improve patient outcomes.FUNDINGNational Cancer Institute (NCI), NIH; Prostate Cancer Foundation; Conquer Cancer, American Society of Clinical Oncology Foundation; UCSF Benioff Initiative for Prostate Cancer Research; Netherlands Cancer Institute.

In prostate cancer, androgen receptor (AR)-targeting agents are very effective in various disease stages. However, therapy resistance inevitably occurs, and little is known about how tumor cells adapt to bypass AR suppression. Here, we performed integrative multiomics analyses on tissues isolated before and after 3 months of AR-targeting enzalutamide monotherapy from patients with high-risk prostate cancer enrolled in a neoadjuvant clinical trial. Transcriptomic analyses demonstrated that AR inhibition drove tumors toward a neuroendocrine-like disease state. Additionally, epigenomic profiling revealed massive enzalutamide-induced reprogramming of pioneer factor FOXA1 from inactive chromatin sites toward active cis-regulatory elements that dictate prosurvival signals. Notably, treatment-induced FOXA1 sites were enriched for the circadian clock component ARNTL. Posttreatment ARNTL levels were associated with patients' clinical outcomes, and ARNTL knockout strongly decreased prostate cancer cell growth. Our data highlight a remarkable cistromic plasticity of FOXA1 following AR-targeted therapy and revealed an acquired dependency on the circadian regulator ARNTL, a novel candidate therapeutic target.

Significance

Understanding how prostate cancers adapt to AR-targeted interventions is critical for identifying novel drug targets to improve the clinical management of treatment-resistant disease. Our study revealed an enzalutamide-induced epigenomic plasticity toward prosurvival signaling and uncovered the circadian regulator ARNTL as an acquired vulnerability after AR inhibition, presenting a novel lead for therapeutic development. See related commentary by Zhang et al., p. 2017. This article is highlighted in the In This Issue feature, p. 2007.

While mutations affecting protein-coding regions have been examined across many cancers, structural variants at the genome-wide level are still poorly defined. Through integrative deep whole-genome and -transcriptome analysis of 101 castration-resistant prostate cancer metastases (109X tumor/38X normal coverage), we identified structural variants altering critical regulators of tumorigenesis and progression not detectable by exome approaches. Notably, we observed amplification of an intergenic enhancer region 624 kb upstream of the androgen receptor (AR) in 81% of patients, correlating with increased AR expression. Tandem duplication hotspots also occur near MYC, in lncRNAs associated with post-translational MYC regulation. Classes of structural variations were linked to distinct DNA repair deficiencies, suggesting their etiology, including associations of CDK12 mutation with tandem duplications, TP53 inactivation with inverted rearrangements and chromothripsis, and BRCA2 inactivation with deletions. Together, these observations provide a comprehensive view of how structural variations affect critical regulators in metastatic prostate cancer.