Extensive androgen receptor enhancer heterogeneity in primary prostate cancers underlies transcriptional diversity and metastatic potential

Estrogen receptor 1 chromatin profiling in human breast tumors reveals high inter-patient heterogeneity with enrichment of risk SNPs and enhancer activity at most-conserved regions

Estrogen Receptor 1 (ESR1; also known as ERα, encoded by ESR1 gene) is the main driver and prime drug target in luminal breast cancer. ESR1 chromatin binding is extensively studied in cell lines and a limited number of human tumors, using consensi of peaks shared among samples. However, little is known about inter-tumor heterogeneity of ESR1 chromatin action, along with its biological implications. Here, we use a large set of ESR1 ChIP-seq data from 70 ESR1+ breast cancers to explore inter-patient heterogeneity in ESR1 DNA binding to reveal a striking inter-tumor heterogeneity of ESR1 action. Of note, commonly shared ESR1 sites show the highest estrogen-driven enhancer activity and are most engaged in long-range chromatin interactions. In addition, the most commonly shared ESR1-occupied enhancers are enriched for breast cancer risk SNP loci. We experimentally confirm SNVs to impact chromatin binding potential for ESR1 and its pioneer factor FOXA1. Finally, in the TCGA breast cancer cohort, we can confirm these variations to associate with differences in expression for the target gene. Cumulatively, we reveal a natural hierarchy of ESR1-chromatin interactions in breast cancers within a highly heterogeneous inter-tumor ESR1 landscape, with the most common shared regions being most active and affected by germline functional risk SNPs for breast cancer development.

Canonical AREs are tumor suppressive regulatory elements in the prostate

The androgen receptor (AR) is the central determinant of prostate tissue identity and differentiation, controlling normal, growth-suppressive prostate-specific gene expression 1 . It is also a key driver of prostate tumorigenesis, becoming "hijacked" to drive oncogenic transcription 2-5 . However, the regulatory elements determining the execution of the growth suppressive AR transcriptional program, and whether this can be reactivated in prostate cancer (PCa) cells remains unclear. Canonical androgen response element (ARE) motifs are the classic DNA binding element for AR 6 . Here, we used a genome-wide strategy to modulate regulatory elements containing AREs to define distinct AR transcriptional programs. We find that activation of these AREs is specifically associated with differentiation and growth suppressive transcription, and this can be reactivated to cause death in AR + PCa cells. In contrast, repression of AREs is well tolerated by PCa cells, but deleterious to normal prostate cells. Finally, gene expression signatures driven by ARE activity are associated with improved prognosis and luminal phenotypes in human PCa patients. This study demonstrates that canonical AREs are responsible for a normal, growth-suppressive, lineage-specific transcriptional program, that this can be reengaged in PCa cells for potential therapeutic benefit, and genes controlled by this mechanism are clinically relevant in human PCa patients.

Increased chromatin accessibility drives transition to androgen receptor splice variant dependence in castration-resistant prostate cancer

Androgen receptor (AR) splice variants, of which ARv7 is the most common, are increased in prostate cancer (PC) that develops resistance to androgen signaling inhibitor drugs, but the extent to which these variants drive AR activity, and whether they have novel functions or dependencies, remain to be determined. We generated a subline of VCaP PC cells (VCaP16) that is resistant to the AR inhibitor enzalutamide (ENZ) and found that AR activity was independent of the full-length AR (ARfl), despite its continued high-level expression, and was instead driven by ARv7. The ARv7 cistrome and transcriptome in VCaP16 cells mirrored that of the ARfl in VCaP cells, although ARv7 chromatin binding was weaker, and strong ARv7 binding sites correlated with higher affinity ARfl binding sites across multiple models and clinical samples. Notably, although ARv7 expression in VCaP cells increased rapidly in response to ENZ, there was a long lag before it gained chromatin binding and transcriptional activity. This lag was associated with an increase in chromatin accessibility, with the AR and nuclear factor I (NFI) motifs being most enriched at these more accessible sites. Moreover, the transcriptional effects of combined NFIB and NFIX knockdown versus ARv7 knockdown were highly correlated. These findings indicate that ARv7 can drive the AR program, but that its activity is dependent on adaptations that increase chromatin accessibility to enhance its intrinsically weak chromatin binding.

Breast cancer risk SNPs converge on estrogen receptor binding sites commonly shared between breast tumors to locally alter estrogen signalling output

Estrogen Receptor alpha (ERα) is the main driver and prime drug target in luminal breast. ERα chromatin binding is extensively studied in cell lines and a limited number of human tumors, using consensi of peaks shared among samples. However, little is known about inter-tumor heterogeneity of ERα chromatin action, along with its biological implications. Here, we use a large set of ERα ChIP-seq data from 70 ERα+ breast cancers to explore inter-patient heterogeneity in ERα DNA binding, to reveal a striking inter-tumor heterogeneity of ERα action. Interestingly, commonly-shared ERα sites showed the highest estrogen-driven enhancer activity and were most-engaged in long-range chromatin interactions. In addition, the most-commonly shared ERα-occupied enhancers were enriched for breast cancer risk SNP loci. We experimentally confirm SNVs to impact chromatin binding potential for ERα and its pioneer factor FOXA1. Finally, in the TCGA breast cancer cohort, we could confirm these variations to associate with differences in expression for the target gene. Cumulatively, we reveal a natural hierarchy of ERα-chromatin interactions in breast cancers within a highly heterogeneous inter-tumor ERα landscape, with the most-common shared regions being most active and affected by germline functional risk SNPs for breast cancer development.

Functional implications and therapeutic targeting of androgen response elements in prostate cancer

The androgen receptor (AR) plays an essential role in the growth and progression of prostate cancer (CaP). Ligand-activated AR inside the nucleus binds to the androgen response element (ARE) of the target genes in dimeric form and recruits transcriptional machinery to facilitate gene transcription. Pharmacological compounds that inhibit the AR action either bind to the ligand binding domain (LBD) or interfere with the interactions of AR with other co-regulatory proteins, slowing the progression of the disease. However, the emergence of resistance to conventional treatment makes clinical management of CaP difficult. Resistance has been associated with activation of androgen/AR axis that restores AR transcriptional activity. Activated AR signaling in resistance cases can be mediated by several mechanisms including AR amplification, gain-of-function AR mutations, androgen receptor variant (ARVs), intracrine androgen production, and overexpression of AR coactivators. Importantly, in castration resistant prostate cancer, ARVs lacking the LBD become constitutively active and promote hormone-independent development, underlining the need to concentrate on the other domain or the AR-DNA interface for the identification of novel actionable targets. In this review, we highlight the plasticity of AR-DNA binding and explain how fine-tuning AR's cooperative interactions with DNA translate into developing an alternative strategy to antagonize AR activity.

Unravelling the Role of Iroquois Homeobox 4 and its Interplay with Androgen Receptor in Prostate Cancer

Genome-wide association studies have linked Iroquois-Homeobox 4 (IRX4) as a robust expression quantitative-trait locus associated with prostate cancer (PCa) risk. However, the intricate mechanism and regulatory factors governing IRX4 expression in PCa remain poorly understood. Here, we unveil enrichment of androgen-responsive gene signatures in metastatic prostate tumors exhibiting heightened IRX4 expression. Furthermore, we uncover a novel interaction between IRX4 and the androgen receptor (AR) co-factor, FOXA1, suggesting that IRX4 modulates PCa cell behavior through AR cistrome alteration. Remarkably, we identified a distinctive short insertion-deletion polymorphism (INDEL), upstream of the IRX4 gene that differentially regulates IRX4 expression through the disruption of AR binding. This INDEL emerges as the most significant PCa risk-associated variant within the 5p15 locus, in a genetic analysis involving 82,591 PCa cases and 61,213 controls and was associated with PCa survival in patients undergoing androgen-deprivation therapy. These studies suggest the potential of this INDEL as a prognostic biomarker for androgen therapy in PCa and IRX4 as a potential therapeutic target in combination with current clinical management.

Lineage Plasticity and Stemness Phenotypes in Prostate Cancer: Harnessing the Power of Integrated "Omics" Approaches to Explore Measurable Metrics

Prostate cancer (PCa), the most frequent and second most lethal cancer type in men in developed countries, is a highly heterogeneous disease. PCa heterogeneity, therapy resistance, stemness, and lethal progression have been attributed to lineage plasticity, which refers to the ability of neoplastic cells to undergo phenotypic changes under microenvironmental pressures by switching between developmental cell states. What remains to be elucidated is how to identify measurements of lineage plasticity, how to implement them to inform preclinical and clinical research, and, further, how to classify patients and inform therapeutic strategies in the clinic. Recent research has highlighted the crucial role of next-generation sequencing technologies in identifying potential biomarkers associated with lineage plasticity. Here, we review the genomic, transcriptomic, and epigenetic events that have been described in PCa and highlight those with significance for lineage plasticity. We further focus on their relevance in PCa research and their benefits in PCa patient classification. Finally, we explore ways in which bioinformatic analyses can be used to determine lineage plasticity based on large omics analyses and algorithms that can shed light on upstream and downstream events. Most importantly, an integrated multiomics approach may soon allow for the identification of a lineage plasticity signature, which would revolutionize the molecular classification of PCa patients.

African American Prostate Cancer Displays Quantitatively Distinct Vitamin D Receptor Cistrome-transcriptome Relationships Regulated by BAZ1A

African American (AA) prostate cancer associates with vitamin D3 deficiency, but vitamin D receptor (VDR) genomic actions have not been investigated in this context. We undertook VDR proteogenomic analyses in European American (EA) and AA prostate cell lines and four clinical cohorts. Rapid immunoprecipitation mass spectrometry of endogenous protein (RIME) analyses revealed that nonmalignant AA RC43N prostate cells displayed the greatest dynamic protein content in the VDR complex. Likewise, in AA cells, Assay for Transposase-Accessible Chromatin using sequencing established greater 1α,25(OH)2D3-regulated chromatin accessibility, chromatin immunoprecipitation sequencing revealed significant enhancer-enriched VDR cistrome, and RNA sequencing identified the largest 1α,25(OH)2D3-dependent transcriptome. These VDR functions were significantly corrupted in the isogenic AA RC43T prostate cancer cells, and significantly distinct from EA cell models. We identified reduced expression of the chromatin remodeler, BAZ1A, in three AA prostate cancer cohorts as well as RC43T compared with RC43N. Restored BAZ1A expression significantly increased 1α,25(OH)2D3-regulated VDR-dependent gene expression in RC43T, but not HPr1AR or LNCaP cells. The clinical impact of VDR cistrome-transcriptome relationships were tested in three different clinical prostate cancer cohorts. Strikingly, only in AA patients with prostate cancer, the genes bound by VDR and/or associated with 1α,25(OH)2D3-dependent open chromatin (i) predicted progression from high-grade prostatic intraepithelial neoplasia to prostate cancer; (ii) responded to vitamin D3 supplementation in prostate cancer tumors; (iii) differentially responded to 25(OH)D3 serum levels. Finally, partial correlation analyses established that BAZ1A and components of the VDR complex identified by RIME significantly strengthened the correlation between VDR and target genes in AA prostate cancer only. Therefore, VDR transcriptional control is most potent in AA prostate cells and distorted through a BAZ1A-dependent control of VDR function.

Significance

Our study identified that genomic ancestry drives the VDR complex composition, genomic distribution, and transcriptional function, and is disrupted by BAZ1A and illustrates a novel driver for AA prostate cancer.

Let’s Go 3D! New Generation of Models for Evaluating Drug Response and Resistance in Prostate Cancer

Prostate cancer (PC) is the third most frequently diagnosed cancer worldwide and the second most frequent in men. Several risk factors can contribute to the development of PC, and those include age, family history, and specific genetic mutations. So far, drug testing in PC, as well as in cancer research in general, has been performed on 2D cell cultures. This is mainly because of the vast benefits these models provide, including simplicity and cost effectiveness. However, it is now known that these models are exposed to much higher stiffness; lose physiological extracellular matrix on artificial plastic surfaces; and show changes in differentiation, polarization, and cell–cell communication. This leads to the loss of crucial cellular signaling pathways and changes in cell responses to stimuli when compared to in vivo conditions. Here, we emphasize the importance of a diverse collection of 3D PC models and their benefits over 2D models in drug discovery and screening from the studies done so far, outlining their benefits and limitations. We highlight the differences between the diverse types of 3D models, with the focus on tumor–stroma interactions, cell populations, and extracellular matrix composition, and we summarize various standard and novel therapies tested on 3D models of PC for the purpose of raising awareness of the possibilities for a personalized approach in PC therapy.