Search

Abstract Androgen receptor (AR)-mediated transcription plays a critical role in development and prostate cancer growth. AR drives gene expression by binding to thousands of cis-regulatory elements (CRE) that loop to hundreds of target promoters. With multiple CREs interacting with a single promoter, it remains unclear how individual AR bound CREs contribute to gene expression. To characterize the involvement of these CREs, we investigate the AR-driven epigenetic and chromosomal chromatin looping changes by generating a kinetic multi-omic dataset comprised of steady-state mRNA, chromatin accessibility, transcription factor binding, histone modifications, chromatin looping, and nascent RNA. Using an integrated regulatory network, we find that AR binding induces sequential changes in the epigenetic features at CREs, independent of gene expression. Further, we show that binding of AR does not result in a substantial rewiring of chromatin loops, but instead increases the contact frequency of pre-existing loops to target promoters. Our results show that gene expression strongly correlates to the changes in contact frequency. We then propose and experimentally validate an unbalanced multi-enhancer model where the impact on gene expression of AR-bound enhancers is heterogeneous, and is proportional to their contact frequency with target gene promoters. Overall, these findings provide insights into AR-mediated gene expression upon acute androgen simulation and develop a mechanistic framework to investigate nuclear receptor mediated perturbations.

Summary: Renal cell carcinoma with sarcomatoid differentiation (sRCC) is associated with poor survival and a heightened response to immune checkpoint inhibitors (ICIs). Two major barriers to improving outcomes for sRCC are the limited understanding of its gene regulatory programs and the low diagnostic yield of tumor biopsies due to spatial heterogeneity. Herein, we characterized the epigenomic landscape of sRCC by profiling 107 epigenomic libraries from tissue and plasma samples from 50 patients with RCC and healthy volunteers. By profiling histone modifications and DNA methylation, we identified highly recurrent epigenomic reprogramming enriched in sRCC. Furthermore, CRISPRa experiments implicated the transcription factor FOSL1 in activating sRCC-associated gene regulatory programs, and FOSL1 expression was associated with the response to ICIs in RCC in two randomized clinical trials. Finally, we established a blood-based diagnostic approach using detectable sRCC epigenomic signatures in patient plasma, providing a framework for discovering epigenomic correlates of tumor histology via liquid biopsy.

Abstract We introduce Promoter-Enhancer-Guided Interaction Networks (PENGUIN), a method for studying protein-protein interaction (PPI) networks within enhancer-promoter interactions. PENGUIN integrates H3K27ac-HiChIP data with tissue-specific PPIs to define enhancer-promoter PPI networks (EPINs). We validated PENGUIN using cancer (LNCaP) and benign (LHSAR) prostate cell lines. Our analysis detected EPIN clusters enriched with the architectural protein CTCF, a regulator of enhancer-promoter interactions. CTCF presence was coupled with the prevalence of prostate cancer (PrCa) single nucleotide polymorphisms (SNPs) within the same EPIN clusters, suggesting functional implications in PrCa. Within the EPINs displaying enrichments in both CTCF and PrCa SNPs, we also show enrichment in oncogenes. We substantiated our identified SNPs through CRISPR/Cas9 knockout and RNAi screens experiments. Here we show that PENGUIN provides insights into the intricate interplay between enhancer-promoter interactions and PPI networks, which are crucial for identifying key genes and potential intervention targets. A dedicated server is available at https://penguin.life.bsc.es/ .

Abstract To date, single-nucleotide polymorphisms (SNPs) have been the most intensively investigated class of polymorphisms in genome wide associations studies (GWAS), however, other classes such as insertion-deletion or multiple nucleotide length polymorphism (MNLPs) may also confer disease risk. Multiple reports have shown that the 5p15.33 prostate cancer risk region is a particularly strong expression quantitative trait locus (eQTL) for Iroquois Homeobox 4 (IRX4) transcripts. Here, we demonstrate using epigenome and genome editing that a biallelic (21 and 47 base pairs (bp)) MNLP is the causal variant regulating IRX4 transcript levels. In LNCaP prostate cancer cells (homozygous for the 21 bp short allele), a single copy knock-in of the 47 bp long allele potently alters the chromatin state, enabling de novo functional binding of the androgen receptor (AR) associated with increased chromatin accessibility, Histone 3 lysine 27 acetylation (H3K27ac), and ~3-fold upregulation of IRX4 expression. We further show that an MNLP is amongst the strongest candidate susceptibility variants at two additional prostate cancer risk loci. We estimated that at least 5% of prostate cancer risk loci could be explained by functional non-SNP causal variants, which may have broader implications for other cancers GWAS. More generally, our results underscore the importance of investigating other classes of inherited variation as causal mediators of human traits.

The epigenomic landscape of renal cell carcinoma (RCC) remains to be explored. Here, integrative epigenomic analysis of primary human RCC samples and RCC GWAS risk SNPs identifies transcription-factor specific subtypes and enrichment of risk variants in allelically-imbalanced peaks.

Epigenetic reprogramming of the androgen receptor (AR) has been identified as an important process driving prostate cancer (PCa) progression. Here, the authors analyze the role of AR chromatin binding heterogeneity in PCa clinical outcomes, metastasis and relapse.

The role of MYC in transcriptional reprogramming in prostate cancer remains poorly characterized. Here, MYC overexpression antagonizes the canonical AR transcriptional program leading to prostate tumor initiation and progression by disrupting transcriptional pause release at AR-regulated genes.

Androgen receptor in prostate cancer (PCa) transcriptionally represses multiple genes including MYC. Here, the authors suggest that increased MYC in response to androgen deprivation contributes to castration-resistant PCa, while decreased MYC may contribute to responses to supraphysiological androgen therapy.

Neuroendocrine carcinomas (NECs) arise from different anatomic sites, but have similar histological and clinical features. Here, the authors show that the epigenetic landscape of a range of NECs converges towards a common epigenetic state, while distinct subtypes occur within neuroendocrine prostate cancer contributing to intratumor heterogeneity in clinical samples.

The molecular processes that lead to neuroendocrine prostate cancer after treating prostate adenocarcinoma (PRAD) are not well understood. Here the authors show that regulation by FOXA1 and changes in the epigenomic profile drive the transition from PRAD to a neuroendocrine phenotype.

The androgen receptor (AR) is a master transcription factor that regulates prostate cancer (PC) development and progression. Inhibition of AR signaling by androgen deprivation is the first-line therapy with initial efficacy for advanced and recurrent PC. Paradoxically, supraphysiological levels of testosterone (SPT) also inhibit PC progression. However, as with any therapy, not all patients show a therapeutic benefit, and responses differ widely in magnitude and duration. In this study, we evaluated whether differences in the AR cistrome before treatment can distinguish between SPT-responding (R) and -nonresponding (NR) tumors. We provide the first preclinical evidence to our knowledge that SPT-R tumors exhibit a distinct AR cistrome when compared with SPT-NR tumors, indicating a differential biological role of the AR. We applied an integrated analysis of ChIP-Seq and RNA-Seq to the pretreatment tumors and identified an SPT-R signature that distinguishes R and NR tumors. Because transcriptomes of SPT-treated clinical specimens are not available, we interrogated available castration-resistant PC (CRPC) transcriptomes and showed that the SPT-R signature is associated with improved survival and has the potential to identify patients who would respond to SPT. These findings provide an opportunity to identify the subset of patients with CRPC who would benefit from SPT therapy.

Metastatic castration-resistant prostate cancers (mCRPCs) are treated with therapies that antagonize the androgen receptor (AR). Nearly all patients develop resistance to AR-targeted therapies (ARTs). Our previous work identified CREB5 as an upregulated target gene in human mCRPC that promoted resistance to all clinically approved ART. The mechanisms by which CREB5 promotes progression of mCRPC or other cancers remains elusive. Integrating ChIP-seq and rapid immunoprecipitation and mass spectroscopy of endogenous proteins, we report that cells overexpressing CREB5 demonstrate extensive reprogramming of nuclear protein–protein interactions in response to the ART agent enzalutamide. Specifically, CREB5 physically interacts with AR, the pioneering actor FOXA1, and other known co-factors of AR and FOXA1 at transcription regulatory elements recently found to be active in mCRPC patients. We identified a subset of CREB5/FOXA1 co-interacting nuclear factors that have critical functions for AR transcription (GRHL2, HOXB13) while others (TBX3, NFIC) regulated cell viability and ART resistance and were amplified or overexpressed in mCRPC. Upon examining the nuclear protein interactions and the impact of CREB5 expression on the mCRPC patient transcriptome, we found that CREB5 was associated with Wnt signaling and epithelial to mesenchymal transitions, implicating these pathways in CREB5/FOXA1-mediated ART resistance. Overall, these observations define the molecular interactions among CREB5, FOXA1, and pathways that promote ART resistance.

The role of non-coding somatic mutations in ovarian cancer is unclear. Here, the authors integrate genomic and epigenomic data from patient samples to show that these mutations frequently converge on the PAX8 transcriptional network.

Summary: Long noncoding RNAs (lncRNAs) have emerged as critical regulators of tumorigenesis, and yet their mechanistic roles remain challenging to characterize. Here, we integrate functional proteomics with lncRNA-interactome profiling to characterize Urothelial Cancer Associated 1 (UCA1), a candidate driver of ovarian cancer development. Reverse phase protein array (RPPA) analysis indicates that UCA1 activates transcription coactivator YAP and its target genes. In vivo RNA antisense purification (iRAP) of UCA1 interacting proteins identified angiomotin (AMOT), a known YAP regulator, as a direct binding partner. Loss-of-function experiments show that AMOT mediates YAP activation by UCA1, as UCA1 enhances the AMOT-YAP interaction to promote YAP dephosphorylation and nuclear translocation. Together, we characterize UCA1 as a lncRNA regulator of Hippo-YAP signaling and highlight the UCA1-AMOT-YAP signaling axis in ovarian cancer development. : Cancer Systems Biology; Omics; Proteomics; Systems Biology Subject Areas: Cancer Systems Biology, Omics, Proteomics, Systems Biology

[This corrects the article DOI: 10.1371/journal.pone.0216997.].

Summary: Fallopian tube secretory epithelial cells (FTSECs) are likely the main precursor cell type of high-grade serous ovarian cancers (HGSOCs), but these tumors may also arise from ovarian surface epithelial cells (OSECs). We profiled global landscapes of gene expression and active chromatin to characterize molecular similarities between OSECs (n = 114), FTSECs (n = 74), and HGSOCs (n = 394). A one-class machine learning algorithm predicts that most HGSOCs derive from FTSECs, with particularly high FTSEC scores in mesenchymal-type HGSOCs (padj < 8 × 10−4). However, a subset of HGSOCs likely derive from OSECs, particularly HGSOCs of the proliferative type (padj < 2 × 10−4), suggesting a dualistic model for HGSOC origins. Super-enhancer (SE) landscapes were also more similar between FTSECs and HGSOCs than between OSECs and HGSOCs (p < 2.2 × 10−16). The SOX18 transcription factor (TF) coincided with a HGSOC-specific SE, and ectopic overexpression of SOX18 in FTSECs caused epithelial-to-mesenchymal transition, indicating that SOX18 plays a role in establishing the mesenchymal signature of fallopian-derived HGSOCs. : Lawrenson et al. profile gene expression and active chromatin in ∼200 ovarian and fallopian epithelial isolates and implement machine learning to demonstrate that most high-grade serous ovarian cancers (HGSOCs) derive from fallopian tube epithelial cells, but a subset may originate from ovarian epithelia. SOX18 induces mesenchymal features to drive early neoplasia in fallopian tube precursors. Keywords: high-grade serous ovarian cancer, ovarian surface epithelial cell, fallopian tube secretory epithelial cell, super enhancers, transcription factors, SOX18, single-cell RNA-seq, RNA-seq, machine learning, one-class logistic regression models, dual origins

Summary: Androgen-receptor (AR) inhibitors, including enzalutamide, are used for treatment of all metastatic castration-resistant prostate cancers (mCRPCs). However, some patients develop resistance or never respond. We find that the transcription factor CREB5 confers enzalutamide resistance in an open reading frame (ORF) expression screen and in tumor xenografts. CREB5 overexpression is essential for an enzalutamide-resistant patient-derived organoid. In AR-expressing prostate cancer cells, CREB5 interactions enhance AR activity at a subset of promoters and enhancers upon enzalutamide treatment, including MYC and genes involved in the cell cycle. In mCRPC, we found recurrent amplification and overexpression of CREB5. Our observations identify CREB5 as one mechanism that drives resistance to AR antagonists in prostate cancers. : Advanced prostate cancers develop resistance to androgen receptor (AR)-targeting therapies. Hwang et al. show that resistant prostate cancers overexpress or amplify CREB5, mediating resistance to AR inhibition. CREB5 suppression reduces the viability of therapy-resistant patient-derived models, suggesting that CREB5 is a target in prostate cancer patients with limited treatment options. Keywords: metastatic castration-resistant prostate cancer, androgen deprivation therapy, therapy resistance, CREB5, androgen receptor, ORF screen

PurposeGenome-wide-association studies (GWAS) have identified numerous single nucleotide polymorphisms (SNPs) that are associated with an increased risk of breast cancer. Most of these studies were conducted primarily in postmenopausal breast cancer patients. Therefore, we set out to assess whether or not these breast cancer variants are also associated with an elevated risk of breast cancer in young premenopausal patients.MethodsIn 451 women of European ancestry who had prospectively enrolled in a longitudinal cohort study for women diagnosed with breast cancer at or under age 40, we genotyped 44 SNPs that were previously associated with breast cancer risk. A control group was comprised of 1142 postmenopausal healthy women from the Nurses' Health Study (NHS). We assessed if the frequencies of the adequately genotyped SNPs differed significantly (p≤0.05) between the cohort of young breast cancer patients and postmenopausal controls, and then we corrected for multiple testing.ResultsGenotyping of the controls or cases was inadequate for comparisons between the groups for seven of the 44 SNPs. 9 of the remaining 37 were associated with breast cancer risk in young women with a p-value ConclusionAfter correction for multiple testing, none of the alleles for postmenopausal breast cancer were clearly associated with risk of premenopausal breast cancer in this relatively small study.