Cistrome Partitioning Reveals Convergence of Somatic Mutations and Risk Variants on Master Transcription Regulators in Primary Prostate Tumors

Hidden secrets of the cancer genome: unlocking the impact of non-coding mutations in gene regulatory elements

Discoveries in the field of genomics have revealed that non-coding genomic regions are not merely "junk DNA", but rather comprise critical elements involved in gene expression. These gene regulatory elements (GREs) include enhancers, insulators, silencers, and gene promoters. Notably, new evidence shows how mutations within these regions substantially influence gene expression programs, especially in the context of cancer. Advances in high-throughput sequencing technologies have accelerated the identification of somatic and germline single nucleotide mutations in non-coding genomic regions. This review provides an overview of somatic and germline non-coding single nucleotide alterations affecting transcription factor binding sites in GREs, specifically involved in cancer biology. It also summarizes the technologies available for exploring GREs and the challenges associated with studying and characterizing non-coding single nucleotide mutations. Understanding the role of GRE alterations in cancer is essential for improving diagnostic and prognostic capabilities in the precision medicine era, leading to enhanced patient-centered clinical outcomes.

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.

It Takes Two to Tango: The Interplay between Prostate Cancer and Its Microenvironment from an Epigenetic Perspective

Prostate cancer is the second most common cancer in men worldwide and is associated with high morbidity and mortality. Consequently, there is an urgent unmet need for novel treatment avenues. In addition to somatic genetic alterations, deviations in the epigenetic landscape of cancer cells and their tumor microenvironment (TME) are critical drivers of prostate cancer initiation and progression. Unlike genomic mutations, epigenetic modifications are potentially reversible. Therefore, the inhibition of aberrant epigenetic modifications represents an attractive and exciting novel treatment strategy for castration-resistant prostate cancer patients. Moreover, drugs targeting the epigenome also exhibit synergistic interactions with conventional therapeutics by directly enhancing their anti-tumorigenic properties by "priming" the tumor and tumor microenvironment to increase drug sensitivity. This review summarizes the major epigenetic alterations in prostate cancer and its TME, and their involvement in prostate tumorigenesis, and discusses the impact of epigenome-targeted therapies.

Transposable Elements Are Co-opted as Oncogenic Regulatory Elements by Lineage-Specific Transcription Factors in Prostate Cancer

Transposable elements hold regulatory functions that impact cell fate determination by controlling gene expression. However, little is known about the transcriptional machinery engaged at transposable elements in pluripotent and mature versus oncogenic cell states. Through positional analysis over repetitive DNA sequences of H3K27ac chromatin immunoprecipitation sequencing data from 32 normal cell states, we report pluripotent/stem and mature cell state-specific "regulatory transposable elements." Pluripotent/stem elements are binding sites for pluripotency factors (e.g., NANOG, SOX2, OCT4). Mature cell elements are docking sites for lineage-specific transcription factors, including AR and FOXA1 in prostate epithelium. Expanding the analysis to prostate tumors, we identify a subset of regulatory transposable elements shared with pluripotent/stem cells, including Tigger3a. Using chromatin editing technology, we show how such elements promote prostate cancer growth by regulating AR transcriptional activity. Collectively, our results suggest that oncogenesis arises from lineage-specific transcription factors hijacking pluripotent/stem cell regulatory transposable elements.

SIGNIFICANCE

We show that oncogenesis relies on co-opting transposable elements from pluripotent stem cells as regulatory elements altering the recruitment of lineage-specific transcription factors. We further discover how co-option is dependent on active chromatin states with important implications for developing treatment options against drivers of oncogenesis across the repetitive DNA. This article is featured in Selected Articles from This Issue, p. 2293.

Integrative identification of non-coding regulatory regions driving metastatic prostate cancer

Large-scale sequencing efforts of thousands of tumor samples have been undertaken to understand the mutational landscape of the coding genome. However, the vast majority of germline and somatic variants occur within non-coding portions of the genome. These genomic regions do not directly encode for specific proteins, but can play key roles in cancer progression, for example by driving aberrant gene expression control. Here, we designed an integrative computational and experimental framework to identify recurrently mutated non-coding regulatory regions that drive tumor progression. Application of this approach to whole-genome sequencing (WGS) data from a large cohort of metastatic castration-resistant prostate cancer (mCRPC) revealed a large set of recurrently mutated regions. We used (i) in silico prioritization of functional non-coding mutations, (ii) massively parallel reporter assays, and (iii) in vivo CRISPR-interference (CRISPRi) screens in xenografted mice to systematically identify and validate driver regulatory regions that drive mCRPC. We discovered that one of these enhancer regions, GH22I030351, acts on a bidirectional promoter to simultaneously modulate expression of U2-associated splicing factor SF3A1 and chromosomal protein CCDC157. We found that both SF3A1 and CCDC157 are promoters of tumor growth in xenograft models of prostate cancer. We nominated a number of transcription factors, including SOX6, to be responsible for higher expression of SF3A1 and CCDC157. Collectively, we have established and confirmed an integrative computational and experimental approach that enables the systematic detection of non-coding regulatory regions that drive the progression of human cancers.

Epigenetics regulation of prostate cancer: Biomarker and therapeutic potential

Prostate cancer (CaP) is the second leading cause of cancer death and displays a broad range of clinical behavior from relatively indolent to aggressive metastatic disease. The etiology of most cases of CaP is not understood completely, which makes it imperative to search for the molecular basis of CaP and markers for early diagnosis. Epigenetic modifications, including changes in DNA methylation patterns, histone modifications, miRNAs, and lncRNAs are key drivers of prostate tumorigenesis. These epigenetic defects might be due to deregulated expression of the epigenetic machinery, affecting the expression of several important genes like GSTP1, RASSF1, CDKN2, RARRES1, IGFBP3, RARB, TMPRSS2-ERG, ITGB4, AOX1, HHEX, WT1, HSPE, PLAU, FOXA1, ASC, GPX3, EZH2, LSD1, etc. In this review, we highlighted the most important epigenetic gene alterations and their variations as a diagnostic marker and target for therapeutic intervention of CaP in the future. Characterization of epigenetic changes involved in CaP is obscure and adequate validation studies are still required to corroborate the present results that would be the impending future of transforming basic research settings into clinical practice.

Associations of CTCF and FOXA1 with androgen and IGF pathways in men with localized prostate cancer

AIMS

To examine associations between the transcription factors CCCTC-binding factor (CTCF) and forkhead box protein A1 (FOXA1) and the androgen receptor (AR) and their association with components of the insulin-like growth factor (IGF)-pathway in a cohort of men with localized prostate cancer.

METHODS

Using prostate tissue samples collected during the Prostate cancer: Evidence of Exercise and Nutrition Trial (PrEvENT) trial (N = 70 to 92, depending on section availability), we assessed the abundance of CTCF, FOXA1, AR, IGFIR, p-mTOR, PTEN and IGFBP-2 proteins using a modified version of the Allred scoring system. Validation studies were performed using large, publicly available datasets (TCGA) (N = 489).

RESULTS

We identified a strong correlation between CTCF and AR staining with benign prostate tissue. CTCF also strongly associated with the IGFIR, with PTEN and with phospho-mTOR. FOXA1 was also correlated with staining for the IGF-IR, with IGFBP-2 and with staining for activated phosphor-mTOR. The staining for the IGF-IR was strongly correlated with the AR.

CONCLUSION

Our findings emphasise the close and complex links between the endocrine controls, well known to play an important role in prostate cancer, and the transcription factors implicated by the recent genetic evidence.

MYC reshapes CTCF-mediated chromatin architecture in prostate cancer

MYC is a well characterized oncogenic transcription factor in prostate cancer, and CTCF is the main architectural protein of three-dimensional genome organization. However, the functional link between the two master regulators has not been reported. In this study, we find that MYC rewires prostate cancer chromatin architecture by interacting with CTCF protein. Through combining the H3K27ac, AR and CTCF HiChIP profiles with CRISPR deletion of a CTCF site upstream of MYC gene, we show that MYC activation leads to profound changes of CTCF-mediated chromatin looping. Mechanistically, MYC colocalizes with CTCF at a subset of genomic sites, and enhances CTCF occupancy at these loci. Consequently, the CTCF-mediated chromatin looping is potentiated by MYC activation, resulting in the disruption of enhancer-promoter looping at neuroendocrine lineage plasticity genes. Collectively, our findings define the function of MYC as a CTCF co-factor in three-dimensional genome organization.

Enhancer profiling identifies epigenetic markers of endocrine resistance and reveals therapeutic options for metastatic castration-resistant prostate cancer patients

Androgen Receptor (AR) signaling inhibitors, including enzalutamide, are treatment options for patients with metastatic castration-resistant prostate cancer (mCRPC), but resistance inevitably develops. Using metastatic samples from a prospective phase II clinical trial, we epigenetically profiled enhancer/promoter activities with H3K27ac chromatin immunoprecipitation followed by sequencing, before and after AR-targeted therapy. We identified a distinct subset of H3K27ac-differentially marked regions that associated with treatment responsiveness. These data were successfully validated in mCRPC patient-derived xenograft models (PDX). In silico analyses revealed HDAC3 as a critical factor that can drive resistance to hormonal interventions, which we validated in vitro . Using cell lines and mCRPC PDX tumors in vitro , we identified drug-drug synergy between enzalutamide and the pan-HDAC inhibitor vorinostat, providing therapeutic proof-of-concept. These findings demonstrate rationale for new therapeutic strategies using a combination of AR and HDAC inhibitors to improve patient outcome in advanced stages of mCRPC.

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

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.

Extensive germline-somatic interplay contributes to prostate cancer progression through HNF1B co-option of TMPRSS2-ERG

Genome-wide association studies have identified 270 loci conferring risk for prostate cancer (PCa), yet the underlying biology and clinical impact remain to be investigated. Here we observe an enrichment of transcription factor genes including HNF1B within PCa risk-associated regions. While focused on the 17q12/HNF1B locus, we find a strong eQTL for HNF1B and multiple potential causal variants involved in the regulation of HNF1B expression in PCa. An unbiased genome-wide co-expression analysis reveals PCa-specific somatic TMPRSS2-ERG fusion as a transcriptional mediator of this locus and the HNF1B eQTL signal is ERG fusion status dependent. We investigate the role of HNF1B and find its involvement in several pathways related to cell cycle progression and PCa severity. Furthermore, HNF1B interacts with TMPRSS2-ERG to co-occupy large proportion of genomic regions with a remarkable enrichment of additional PCa risk alleles. We finally show that HNF1B co-opts ERG fusion to mediate mechanistic and biological effects of the PCa risk-associated locus 17p13.3/VPS53/FAM57A/GEMIN4. Taken together, we report an extensive germline-somatic interaction between TMPRSS2-ERG fusion and genetic variations underpinning PCa risk association and progression.

Homing in on a Moving Target: Androgen Receptor Cistromic Plasticity in Prostate Cancer

The androgen receptor (AR) is the critical driver in prostate cancer and exerts its function mainly through transcriptional control. Recent advances in clinical studies and cell line models have illustrated that AR chromatin binding features are not static; rather they are highly variable yet reproducibly altered between clinical stages. Extensive genomic analyses of AR chromatin binding features in different disease stages have revealed a high degree of plasticity of AR chromatin interactions in clinical samples. Mechanistically, AR chromatin binding patterns are associated with specific somatic mutations on AR and other permutations, including mutations of AR-interacting proteins. Here we summarize the most recent studies on how the AR cistrome is dynamically altered in prostate cancer models and patient samples, and what implications this has for the identification of therapeutic targets to avoid the emergence of treatment resistance.

Chromatin Variants Reveal the Genetic Determinants of Oncogenesis in Breast Cancer

Breast cancer presents as multiple distinct disease entities. Each tumor harbors diverse cell populations defining a phenotypic heterogeneity that impinges on our ability to treat patients. To date, efforts mainly focused on genetic variants to find drivers of inter- and intratumor phenotypic heterogeneity. However, these efforts have failed to fully capture the genetic basis of breast cancer. Through recent technological and analytical approaches, the genetic basis of phenotypes can now be decoded by characterizing chromatin variants. These variants correspond to polymorphisms in chromatin states at DNA sequences that serve a distinct role across cell populations. Here, we review the function and causes of chromatin variants as they relate to breast cancer inter- and intratumor heterogeneity and how they can guide the development of treatment alternatives to fulfill the goal of precision cancer medicine.

Recurrent noncoding somatic and germline WT1 variants converge to disrupt MYB binding in acute promyelocytic leukemia

Genetic alternations can occur at noncoding regions, but how they contribute to cancer pathogenesis is poorly understood. Here, we established a mutational landscape of cis-regulatory regions (CREs) in acute promyelocytic leukemia (APL) based on whole-genome sequencing analysis of paired tumor and germline samples from 24 patients and epigenetic profiling of 16 patients. Mutations occurring in CREs occur preferentially in active enhancers bound by the complex of master transcription factors in APL. Among significantly enriched mutated CREs, we found a recurrently mutated region located within the third intron of WT1, an essential regulator of normal and malignant hematopoiesis. Focusing on noncoding mutations within this WT1 intron, an analysis on 169 APL patients revealed that somatic mutations were clustered into a focal hotspot region, including one site identified as a germline polymorphism contributing to APL risk. Significantly decreased WT1 expression was observed in APL patients bearing somatic and/or germline noncoding WT1 variants. Furthermore, biallelic WT1 inactivation was recurrently found in APL patients with noncoding WT1 variants, which resulted in the complete loss of WT1. The high incidence of biallelic inactivation suggested the tumor suppressor activity of WT1 in APL. Mechanistically, noncoding WT1 variants disrupted MYB binding on chromatin and suppressed the enhancer activity and WT1 expression through destroying the chromatin looping formation. Our study highlights the important role of noncoding variants in the leukemogenesis of APL.

Drug-induced epigenomic plasticity reprograms circadian rhythm regulation to drive prostate cancer towards androgen-independence

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 multi-omics analyses on tissues isolated before and after 3 months of AR-targeting enzalutamide monotherapy from high-risk prostate cancer patients enrolled in a neoadjuvant clinical trial. Transcriptomic analyses demonstrated that AR inhibition drove tumors towards a neuroendocrine-like disease state. Additionally, epigenomic profiling revealed massive enzalutamide-induced reprogramming of pioneer factor FOXA1 - from inactive chromatin sites towards active cis-regulatory elements that dictate pro-survival signals. Notably, treatment-induced FOXA1 sites were enriched for circadian clock component ARNTL. Post-treatment ARNTL levels associated with poor outcome, 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 circadian regulator ARNTL, a novel candidate therapeutic target.

From Omics to Multi-Omics Approaches for In-Depth Analysis of the Molecular Mechanisms of Prostate Cancer

Cancer arises following alterations at different cellular levels, including genetic and epigenetic modifications, transcription and translation dysregulation, as well as metabolic variations. High-throughput omics technologies that allow one to identify and quantify processes involved in these changes are now available and have been instrumental in generating a wealth of steadily increasing data from patient tumors, liquid biopsies, and from tumor models. Extensive investigation and integration of these data have led to new biological insights into the origin and development of multiple cancer types and helped to unravel the molecular networks underlying this complex pathology. The comprehensive and quantitative analysis of a molecule class in a biological sample is named omics and large-scale omics studies addressing different prostate cancer stages have been performed in recent years. Prostate tumors represent the second leading cancer type and a prevalent cause of cancer death in men worldwide. It is a very heterogenous disease so that evaluating inter- and intra-tumor differences will be essential for a precise insight into disease development and plasticity, but also for the development of personalized therapies. There is ample evidence for the key role of the androgen receptor, a steroid hormone-activated transcription factor, in driving early and late stages of the disease, and this led to the development and approval of drugs addressing diverse targets along this pathway. Early genomic and transcriptomic studies have allowed one to determine the genes involved in prostate cancer and regulated by androgen signaling or other tumor-relevant signaling pathways. More recently, they have been supplemented by epigenomic, cistromic, proteomic and metabolomic analyses, thus, increasing our knowledge on the intricate mechanisms involved, the various levels of regulation and their interplay. The comprehensive investigation of these omics approaches and their integration into multi-omics analyses have led to a much deeper understanding of the molecular pathways involved in prostate cancer progression, and in response and resistance to therapies. This brings the hope that novel vulnerabilities will be identified, that existing therapies will be more beneficial by targeting the patient population likely to respond best, and that bespoke treatments with increased efficacy will be available soon.

HOXB13 suppresses de novo lipogenesis through HDAC3-mediated epigenetic reprogramming in prostate cancer

HOXB13, a homeodomain transcription factor, critically regulates androgen receptor (AR) activities and androgen-dependent prostate cancer (PCa) growth. However, its functions in AR-independent contexts remain elusive. Here we report HOXB13 interaction with histone deacetylase HDAC3, which is disrupted by the HOXB13 G84E mutation that has been associated with early-onset PCa. Independently of AR, HOXB13 recruits HDAC3 to lipogenic enhancers to catalyze histone deacetylation and suppress lipogenic regulators such as fatty acid synthase. Analysis of human tissues reveals that the HOXB13 gene is hypermethylated and downregulated in approximately 30% of metastatic castration-resistant PCa. HOXB13 loss or G84E mutation leads to lipid accumulation in PCa cells, thereby promoting cell motility and xenograft tumor metastasis, which is mitigated by pharmaceutical inhibition of fatty acid synthase. In summary, we present evidence that HOXB13 recruits HDAC3 to suppress de novo lipogenesis and inhibit tumor metastasis and that lipogenic pathway inhibitors may be useful to treat HOXB13-low PCa.

Androgen Receptor-Mediated Transcription in Prostate Cancer

Androgen receptor (AR)-mediated transcription is critical in almost all stages of prostate cancer (PCa) growth and differentiation. This process involves a complex interplay of coregulatory proteins, chromatin remodeling complexes, and other transcription factors that work with AR at cis-regulatory enhancer regions to induce the spatiotemporal transcription of target genes. This enhancer-driven mechanism is remarkably dynamic and undergoes significant alterations during PCa progression. In this review, we discuss the AR mechanism of action in PCa with a focus on how cis-regulatory elements modulate gene expression. We explore emerging evidence of genetic variants that can impact AR regulatory regions and alter gene transcription in PCa. Finally, we highlight several outstanding questions and discuss potential mechanisms of this critical transcription factor.

Mutations in Noncoding Cis-Regulatory Elements Reveal Cancer Driver Cistromes in Luminal Breast Cancer

Whole-genome sequencing of primary breast tumors enabled the identification of cancer driver genes and noncoding cancer driver plexuses from somatic mutations. However, differentiating driver from passenger events among noncoding genetic variants remains a challenge. Herein, we reveal cancer-driver cis-regulatory elements linked to transcription factors previously shown to be involved in development of luminal breast cancers by defining a tumor-enriched catalogue of approximately 100,000 unique cis-regulatory elements from 26 primary luminal estrogen receptor (ER)⁺ progesterone receptor (PR)⁺ breast tumors. Integrating this catalog with somatic mutations from 350 publicly available breast tumor whole genomes, we uncovered cancer driver cistromes, defined as the sum of binding sites for a transcription factor, for ten transcription factors in luminal breast cancer such as FOXA1 and ER, nine of which are essential for growth in breast cancer with four exclusive to the luminal subtype. Collectively, we present a strategy to find cancer driver cistromes relying on quantifying the enrichment of noncoding mutations over cis-regulatory elements concatenated into a functional unit. IMPLICATIONS: Mapping the accessible chromatin of luminal breast cancer led to discovery of an accumulation of mutations within cistromes of transcription factors essential to luminal breast cancer. This demonstrates coopting of regulatory networks to drive cancer and provides a framework to derive insight into the noncoding space of cancer.

Androgen receptor and MYC equilibration centralizes on developmental super-enhancer

Androgen receptor (AR) in prostate cancer (PCa) can drive transcriptional repression of multiple genes including MYC, and supraphysiological androgen is effective in some patients. Here, we show that this repression is independent of AR chromatin binding and driven by coactivator redistribution, and through chromatin conformation capture methods show disruption of the interaction between the MYC super-enhancer within the PCAT1 gene and the MYC promoter. Conversely, androgen deprivation in vitro and in vivo increases MYC expression. In parallel, global AR activity is suppressed by MYC overexpression, consistent with coactivator redistribution. These suppressive effects of AR and MYC are mitigated at shared AR/MYC binding sites, which also have markedly higher levels of H3K27 acetylation, indicating enrichment for functional enhancers. These findings demonstrate an intricate balance between AR and MYC, and indicate that increased MYC in response to androgen deprivation contributes to castration-resistant PCa, while decreased MYC may contribute to responses to supraphysiological androgen therapy.

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.

Mechanistic insights into genetic susceptibility to prostate cancer

Prostate cancer (PCa) is the second most common cancer in men and is a highly heritable disease that affects millions of individuals worldwide. Genome-wide association studies have to date discovered nearly 270 genetic loci harboring hundreds of single nucleotide polymorphisms (SNPs) that are associated with PCa susceptibility. In contrast, the functional characterization of the mechanisms underlying PCa risk association is still growing. Given that PCa risk-associated SNPs are highly enriched in noncoding cis-regulatory genomic regions, accumulating evidence suggests a widespread modulation of transcription factor chromatin binding and allelic enhancer activity by these noncoding SNPs, thereby dysregulating gene expression. Emerging studies have shown that a proportion of noncoding variants can modulate the formation of transcription factor complexes at enhancers and CTCF-mediated 3D genome architecture. Interestingly, DNA methylation-regulated CTCF binding could orchestrate a long-range chromatin interaction between PCa risk enhancer and causative genes. Additionally, one-causal-variant-two-risk genes or multiple-risk-variant-multiple-genes are prevalent in some PCa risk-associated loci. In this review, we will discuss the current understanding of the general principles of SNP-mediated gene regulation, experimental advances, and functional evidence supporting the mechanistic roles of several PCa genetic loci with potential clinical impact on disease prevention and treatment.

Understanding and targeting prostate cancer cell heterogeneity and plasticity

Prostate cancer (PCa) is a prevalent malignancy that occurs primarily in old males. Prostate tumors in different patients manifest significant inter-patient heterogeneity with respect to histo-morphological presentations and molecular architecture. An individual patient tumor also harbors genetically distinct clones in which PCa cells display intra-tumor heterogeneity in molecular features and phenotypic marker expression. This inherent PCa cell heterogeneity, e.g., in the expression of androgen receptor (AR), constitutes a barrier to the long-term therapeutic efficacy of AR-targeting therapies. Furthermore, tumor progression as well as therapeutic treatments induce PCa cell plasticity such that AR-positive PCa cells may turn into AR-negative cells and prostate tumors may switch lineage identity from adenocarcinomas to neuroendocrine-like tumors. This induced PCa cell plasticity similarly confers resistance to AR-targeting and other therapies. In this review, I first discuss PCa from the perspective of an abnormal organ development and deregulated cellular differentiation, and discuss the luminal progenitor cells as the likely cells of origin for PCa. I then focus on intrinsic PCa cell heterogeneity in treatment-naïve tumors with the presence of prostate cancer stem cells (PCSCs). I further elaborate on PCa cell plasticity induced by genetic alterations and therapeutic interventions, and present potential strategies to therapeutically tackle PCa cell heterogeneity and plasticity. My discussions will make it clear that, to achieve enduring clinical efficacy, both intrinsic PCa cell heterogeneity and induced PCa cell plasticity need to be targeted with novel combinatorial approaches.

Reorganization of the 3D genome pinpoints non-coding drivers of primary prostate tumors

Prostate cancer is a heterogeneous disease whose progression is linked to genome instability. However, the impact of this instability on the non-coding genome and its three-dimensional organization to aid progression is unclear. Using primary benign and tumor tissue, we find a high concordance in higher order three-dimensional genome organization. This concordance argues for constraints to the topology of prostate tumor genomes. Nonetheless, we identified changes in focal chromatin interactions, typical of loops bridging non-coding cis-regulatory elements, and showed how structural variants can induce these changes to guide cis-regulatory element hijacking. Such events resulted in opposing differential expression of genes found at antipodes of rearrangements. Collectively, these results argue that changes to focal chromatin interactions, as opposed to higher order genome organization, allow for aberrant gene regulation and are repeatedly mediated by structural variants in primary prostate cancer.

Somatic Alterations Impact AR Transcriptional Activity and Efficacy of AR-Targeting Therapies in Prostate Cancer

Simple Summary

For patients whose prostate cancer spreads beyond the confines of the prostate, treatment options continue to increase. However, we are missing the information that is needed to choose for each patient the best treatment at each step of his cancer progression so we can ensure that maximal remissions and prolonged survival are achieved. In this review, we examine whether a better understanding of how the activity of the target for the default first treatment, the androgen receptor, is regulated in prostate cancer tissues can improve prostate cancer treatment plans. We consider the evidence for variability of androgen receptor activity among patients and examine the molecular basis for this variable action. We summarize clinical evidence supporting that information on a prostate cancer’s genomic composition may inform on its level of androgen receptor action, which may facilitate choice for the most effective first-line therapy and ultimately improve prostate cancer treatment plans overall.

Abstract

Inhibiting the activity of the ligand-activated transcription factor androgen receptor (AR) is the default first-line treatment for metastatic prostate cancer (CaP). Androgen deprivation therapy (ADT) induces remissions, however, their duration varies widely among patients. The reason for this heterogeneity is not known. A better understanding of its molecular basis may improve treatment plans and patient survival. AR’s transcriptional activity is regulated in a context-dependent manner and relies on an interplay between its associated transcriptional regulators, DNA recognition motifs, and ligands. Alterations in one or more of these factors induce shifts in the AR cistrome and transcriptional output. Significant variability in AR activity is seen in both castration-sensitive (CS) and castration-resistant CaP (CRPC). Several AR transcriptional regulators undergo somatic alterations that impact their function in clinical CaPs. Some alterations occur in a significant fraction of cases, resulting in CaP subtypes, while others affect only a few percent of CaPs. Evidence is emerging that these alterations may impact the response to CaP treatments such as ADT, radiation therapy, and chemotherapy. Here, we review the contribution of recurring somatic alterations on AR cistrome and transcriptional output and the efficacy of CaP treatments and explore strategies to use these insights to improve treatment plans and outcomes for CaP patients.

Chromatin and Epigenetic Dysregulation of Prostate Cancer Development, Progression, and Therapeutic Response

The dysregulation of chromatin and epigenetics has been defined as the overarching cancer hallmark. By disrupting transcriptional regulation in normal cells and mediating tumor progression by promoting cancer cell plasticity, this process has the ability to mediate all defined hallmarks of cancer. In this review, we collect and assess evidence on the contribution of chromatin and epigenetic dysregulation in prostate cancer. We highlight important mechanisms leading to prostate carcinogenesis, the emergence of castration-resistance upon treatment with androgen deprivation therapy, and resistance to antiandrogens. We examine in particular the contribution of chromatin structure and epigenetics to cell lineage commitment, which is dysregulated during tumorigenesis, and cell plasticity, which is altered during tumor progression.

Evolution of Castration-Resistant Prostate Cancer in ctDNA during Sequential Androgen Receptor Pathway Inhibition

PURPOSE

Cross-resistance renders multiple lines of androgen receptor (AR) signaling inhibitors increasingly futile in metastatic castration-resistant prostate cancer (mCRPC). We sought to determine acquired genomic contributors to cross-resistance.

EXPERIMENTAL DESIGN

We collected 458 serial plasma cell-free DNA samples at baseline and progression timepoints from 202 patients with mCRPC receiving sequential AR signaling inhibitors (abiraterone and enzalutamide) in a randomized phase II clinical trial (NCT02125357). We utilized deep targeted and whole-exome sequencing to compare baseline and posttreatment somatic genomic profiles in circulating tumor DNA (ctDNA).

RESULTS

Patient ctDNA abundance was correlated across plasma collections and independently prognostic for sequential therapy response and overall survival. Most driver alterations in established prostate cancer genes were consistently detected in ctDNA over time. However, shifts in somatic populations after treatment were identified in 53% of patients, particularly after strong treatment responses. Treatment-associated changes converged upon the AR gene, with an average 50% increase in AR copy number, changes in AR mutation frequencies, and a 2.5-fold increase in the proportion of patients carrying AR ligand binding domain truncating rearrangements.

CONCLUSIONS

Our data show that the dominant AR genotype continues to evolve during sequential lines of AR inhibition and drives acquired resistance in patients with mCRPC.

Functional mapping of androgen receptor enhancer activity

Background

Androgen receptor (AR) is critical to the initiation, growth, and progression of prostate cancer. Once activated, the AR binds to cis-regulatory enhancer elements on DNA that drive gene expression. Yet, there are 10–100× more binding sites than differentially expressed genes. It is unclear how or if these excess binding sites impact gene transcription.

Results

To characterize the regulatory logic of AR-mediated transcription, we generated a locus-specific map of enhancer activity by functionally testing all common clinical AR binding sites with Self-Transcribing Active Regulatory Regions sequencing (STARRseq). Only 7% of AR binding sites displayed androgen-dependent enhancer activity. Instead, the vast majority of AR binding sites were either inactive or constitutively active enhancers. These annotations strongly correlated with enhancer-associated features of both in vitro cell lines and clinical prostate cancer samples. Evaluating the effect of each enhancer class on transcription, we found that AR-regulated enhancers frequently interact with promoters and form central chromosomal loops that are required for transcription. Somatic mutations of these critical AR-regulated enhancers often impact enhancer activity.

Conclusions

Using a functional map of AR enhancer activity, we demonstrated that AR-regulated enhancers act as a regulatory hub that increases interactions with other AR binding sites and gene promoters.

Functional and genetic determinants of mutation rate variability in regulatory elements of cancer genomes

Background

Cancer genomes are shaped by mutational processes with complex spatial variation at multiple scales. Entire classes of regulatory elements are affected by local variations in mutation frequency. However, the underlying mechanisms with functional and genetic determinants remain poorly understood.

Results

We characterise the mutational landscape of 1.3 million gene-regulatory and chromatin architectural elements in 2419 whole cancer genomes with transcriptional and pathway activity, functional conservation and recurrent driver events. We develop RM2, a statistical model that quantifies mutational enrichment or depletion in classes of genomic elements through genetic, trinucleotide and megabase-scale effects. We report a map of localised mutational processes affecting CTCF binding sites, transcription start sites (TSS) and tissue-specific open-chromatin regions. Increased mutation frequency in TSSs associates with mRNA abundance in most cancer types, while open-chromatin regions are generally enriched in mutations. We identify ~ 10,000 CTCF binding sites with core DNA motifs and constitutive binding in 66 cell types that represent focal points of mutagenesis. We detect site-specific mutational signature enrichments, such as SBS40 in open-chromatin regions in prostate cancer and SBS17b in CTCF binding sites in gastrointestinal cancers. Candidate drivers of localised mutagenesis are also apparent: BRAF mutations associate with mutational enrichments at CTCF binding sites in melanoma, and ARID1A mutations with TSS-specific mutagenesis in pancreatic cancer.

Conclusions

Our method and catalogue of localised mutational processes provide novel perspectives to cancer genome evolution, mutagenesis, DNA repair and driver gene discovery. The functional and genetic correlates of mutational processes suggest mechanistic hypotheses for future studies.

CRISPRi screens reveal a DNA methylation-mediated 3D genome dependent causal mechanism in prostate cancer

Prostate cancer (PCa) risk-associated SNPs are enriched in noncoding cis-regulatory elements (rCREs), yet their modi operandi and clinical impact remain elusive. Here, we perform CRISPRi screens of 260 rCREs in PCa cell lines. We find that rCREs harboring high risk SNPs are more essential for cell proliferation and H3K27ac occupancy is a strong indicator of essentiality. We also show that cell-line-specific essential rCREs are enriched in the 8q24.21 region, with the rs11986220-containing rCRE regulating MYC and PVT1 expression, cell proliferation and tumorigenesis in a cell-line-specific manner, depending on DNA methylation-orchestrated occupancy of a CTCF binding site in between this rCRE and the MYC promoter. We demonstrate that CTCF deposition at this site as measured by DNA methylation level is highly variable in prostate specimens, and observe the MYC eQTL in the 8q24.21 locus in individuals with low CTCF binding. Together our findings highlight a causal mechanism synergistically driven by a risk SNP and DNA methylation-mediated 3D genome architecture, advocating for the integration of genetics and epigenetics in assessing risks conferred by genetic predispositions.

Epigenetic and transcriptional analysis reveals a core transcriptional program conserved in clonal prostate cancer metastases

The epigenomic regulation of transcriptional programs in metastatic prostate cancer is poorly understood. We studied the epigenomic landscape of prostate cancer drivers using transcriptional profiling and ChIP‐seq in four clonal metastatic tumors derived from a single prostate cancer patient. Our epigenomic analyses focused on androgen receptor (AR), which is a key oncogenic driver in prostate cancer, the AR pioneer factor FOXA1, chromatin insulator CCCTC‐Binding Factor, as well as for modified histones H3K27ac and H3K27me3. The vast majority of AR binding sites were shared among healthy prostate, primary prostate cancer, and metastatic tumor samples, signifying core AR‐driven transcriptional regulation within the prostate cell lineage. Genes associated with core AR‐binding events were significantly enriched for essential genes in prostate cancer cell proliferation. Remarkably, the metastasis‐specific active AR binding sites showed no differential transcriptional output, indicating a robust transcriptional program across metastatic samples. Combined, our data reveal a core transcriptional program in clonal metastatic prostate cancer, despite epigenomic differences in the AR cistrome.

The circadian cryptochrome, CRY1, is a pro-tumorigenic factor that rhythmically modulates DNA repair

Mechanisms regulating DNA repair processes remain incompletely defined. Here, the circadian factor CRY1, an evolutionally conserved transcriptional coregulator, is identified as a tumor specific regulator of DNA repair. Key findings demonstrate that CRY1 expression is androgen-responsive and associates with poor outcome in prostate cancer. Functional studies and first-in-field mapping of the CRY1 cistrome and transcriptome reveal that CRY1 regulates DNA repair and the G2/M transition. DNA damage stabilizes CRY1 in cancer (in vitro, in vivo, and human tumors ex vivo), which proves critical for efficient DNA repair. Further mechanistic investigation shows that stabilized CRY1 temporally regulates expression of genes required for homologous recombination. Collectively, these findings reveal that CRY1 is hormone-induced in tumors, is further stabilized by genomic insult, and promotes DNA repair and cell survival through temporal transcriptional regulation. These studies identify the circadian factor CRY1 as pro-tumorigenic and nominate CRY1 as a new therapeutic target.

Cryptochrome 1 (CRY1) is a transcriptional coregulator associated with the circadian clock. Here the authors reveal that CRY1 is hormone-regulated, stabilized by genomic insult, and promotes DNA repair and cell survival through temporal transcriptional regulation.

The Nexus of Endocrine Signaling and Cancer: How Steroid Hormones Influence Genomic Stability

Endocrine-driven malignancies, including breast and prostate cancer, are among the most common human cancers. The relationship between sex steroid hormones (eg, androgen, estrogen, and progesterone), their cognate receptors, and genomic stability lie at the center of endocrine-driven cancer development, progression, and therapeutic resistance. A variety of direct and indirect mechanisms have been described that link steroid hormone signaling to the loss of genomic integrity that drives early carcinogenesis. These effects are often enriched within endocrine receptor cistromes, accounting for the high proportion of mutations and rearrangements in the region of hormone response elements. In other cases, the effects are generalized and rely on a complex array of genetic, epigenetic, and metabolic interactions. Both androgen and estrogen receptors directly modulate the DNA damage response by trans-activating DNA damage response genes and redirecting the cellular repair machinery in the wake of genotoxic stress. Here we review the key mechanistic underpinnings of the relationship between sex steroid hormone receptors and genomic stability. In addition, we summarize emerging research in this area and discuss important implications for cancer prevention and treatment.

Pioneer of prostate cancer: past, present and the future of FOXA1

Prostate cancer is the most commonly diagnosed non-cutaneous cancers in North American men. While androgen deprivation has remained as the cornerstone of prostate cancer treatment, resistance ensues leading to lethal disease. Forkhead box A1 (FOXA1) encodes a pioneer factor that induces open chromatin conformation to allow the binding of other transcription factors. Through direct interactions with the Androgen Receptor (AR), FOXA1 helps to shape AR signaling that drives the growth and survival of normal prostate and prostate cancer cells. FOXA1 also possesses an AR-independent role of regulating epithelial-to-mesenchymal transition (EMT). In prostate cancer, mutations converge onto the coding sequence and cis-regulatory elements (CREs) of FOXA1, leading to functional alterations. In addition, FOXA1 activity in prostate cancer can be modulated post-translationally through various mechanisms such as LSD1-mediated protein demethylation. In this review, we describe the latest discoveries related to the function and regulation of FOXA1 in prostate cancer, pointing to their relevance to guide future clinical interventions.

Prostate cancer reactivates developmental epigenomic programs during metastatic progression

Epigenetic processes govern prostate cancer (PCa) biology, as evidenced by the dependency of PCa cells on the androgen receptor (AR), a prostate master transcription factor. We generated 268 epigenomic datasets spanning two state transitions-from normal prostate epithelium to localized PCa to metastases-in specimens derived from human tissue. We discovered that reprogrammed AR sites in metastatic PCa are not created de novo; rather, they are prepopulated by the transcription factors FOXA1 and HOXB13 in normal prostate epithelium. Reprogrammed regulatory elements commissioned in metastatic disease hijack latent developmental programs, accessing sites that are implicated in prostate organogenesis. Analysis of reactivated regulatory elements enabled the identification and functional validation of previously unknown metastasis-specific enhancers at HOXB13, FOXA1 and NKX3-1. Finally, we observed that prostate lineage-specific regulatory elements were strongly associated with PCa risk heritability and somatic mutation density. Examining prostate biology through an epigenomic lens is fundamental for understanding the mechanisms underlying tumor progression.

Functional and genetic determinants of mutation rate variability in regulatory elements of cancer genomes.. [DOI: 10.1101/2020.07.29.226373] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Background

Cancer genomes are shaped by mutational processes with complex spatial variation at multiple scales. Entire classes of regulatory elements are affected by local variations in mutation frequency. However, the underlying mutational mechanisms with functional and genetic determinants remain poorly understood.

Results

We characterised the mutational landscape of 1.3 million gene regulatory and chromatin architectural elements in 2,419 whole cancer genomes with transcriptional and pathway activity, functional conservation and recurrent driver events. We developed RM2, a statistical model that quantifies mutational enrichment or depletion in classes of genomic elements through genetic, trinucleotide and megabase-scale effects. We report a map of localised mutational processes affecting CTCF binding sites, transcription start sites (TSS) and tissue-specific open-chromatin regions. We show that increased mutational frequency in TSSs correlates with mRNA abundance in most cancer types, while open-chromatin regions are generally enriched in mutations. We identified ∼10,000 CTCF binding sites with core DNA motifs and constitutive binding in 66 cell types that represent focal points of local mutagenesis. We detected site-specific mutational signatures, such as SBS40 in open-chromatin regions in prostate cancer and SBS17b in CTCF binding sites in gastrointestinal cancers. We also proposed candidate drivers of localised mutagenesis: BRAF mutations associate with mutational enrichments at CTCF binding sites in melanoma, and ARID1A mutations with TSS-specific mutations in pancreatic cancer.

Conclusions

Our method and catalogue of localised mutational processes provide novel perspectives to cancer genome evolution, mutagenesis, DNA repair and driver discovery. Functional and genetic correlates of localised mutagenesis provide mechanistic hypotheses for future studies.

Molecular and Clinical Relevance of ZBTB38 Expression Levels in Prostate Cancer

Prostate cancer is one of the most commonly diagnosed cancers in men. A number of genomic and clinical studies have led to a better understanding of prostate cancer biology. Still, the care of patients as well as the prediction of disease aggressiveness, recurrence and outcome remain challenging. Here, we showed that expression of the gene ZBTB38 is associated with poor prognosis in localised prostate cancer and could help discriminate aggressive localised prostate tumours from those who can benefit only from observation. Analysis of different prostate cancer cohorts indicates that low expression levels of ZBTB38 associate with increased levels of chromosomal abnormalities and more aggressive pathological features, including higher rate of biochemical recurrence of the disease. Importantly, gene expression profiling of these tumours, complemented with cellular assays on prostate cancer cell lines, unveiled that tumours with low levels of ZBTB38 expression might be targeted by doxorubicin, a compound generating reactive oxygen species. Our study shows that ZBTB38 is involved in prostate cancer pathogenesis and may represent a useful marker to identify high risk and highly rearranged localised prostate cancer susceptible to doxorubicin.

Noncoding mutations target cis-regulatory elements of the FOXA1 plexus in prostate cancer

Prostate cancer is the second most commonly diagnosed malignancy among men worldwide. Recurrently mutated in primary and metastatic prostate tumors, FOXA1 encodes a pioneer transcription factor involved in disease onset and progression through both androgen receptor-dependent and androgen receptor-independent mechanisms. Despite its oncogenic properties however, the regulation of FOXA1 expression remains unknown. Here, we identify a set of six cis-regulatory elements in the FOXA1 regulatory plexus harboring somatic single-nucleotide variants in primary prostate tumors. We find that deletion and repression of these cis-regulatory elements significantly decreases FOXA1 expression and prostate cancer cell growth. Six of the ten single-nucleotide variants mapping to FOXA1 regulatory plexus significantly alter the transactivation potential of cis-regulatory elements by modulating the binding of transcription factors. Collectively, our results identify cis-regulatory elements within the FOXA1 plexus mutated in primary prostate tumors as potential targets for therapeutic intervention.

FOXA1 pioneer transcription factor is recurrently mutated in primary and metastatic prostate tumors. Here, authors identify a set of six cis-regulatory elements in the FOXA1 regulatory plexus harboring somatic SNVs in primary prostate tumors and characterize their role in regulating FOXA1 expression and prostate cancer cell growth.