1. Introduction
Chromatin structure and epigenetics are intertwined but, nonetheless, distinct entities that have been implicated in prostate cancer (PC) disease initiation and progression. Here, we consider epigenetics to be the collection of DNA modifications such as DNA methylation. This definition of epigenetics can also include histone modifications and the binding of transcription factors (TFs) to DNA. Histone modifications are often referred to as epigenetic profiles as they determine chromatin states and nucleosome positioning, which, in turn, allows for DNA accessibility (Figure 1). Therefore, we additionally define chromatin structure as histone modifications and nucleosome positioning, as well as the three-dimensional (3D) organization of the chromatin within the nucleus (Figure 1).
Figure 1. Epigenetic- and chromatin-related mechanisms with potential for dysregulation in prostate cancer cells. Epigenetic dysregulation can occur at multiple levels, including changes in chromatin accessibility, histone and DNA modification through processes such as methylation, chromatin remodeling, modification of transcription factors and changes in their availability, cis-regulatory elements, chromatin loops, and topologically associated domains. These chromatin and epigenetic features can be analyzed via the integration of multiple high-throughput sequencing data types. These data include chromatin conformation capture (Hi-C) to understand the 3D chromatin structure and topologically associated domains, chromatin immunoprecipitation (ChIP-seq) to study histone markers, assay for transposase-accessible chromatin (ATAC-seq) to show chromatin accessibility patterns, and DNA methylation sequencing. Figure created with
BioRender.com.
PC continues to be a major cause of cancer-related death in men worldwide
[1]. Although primary intervention with radiotherapy or surgery and androgen deprivation therapy (ADT) have a curative intent in hormone-sensitive PC (HSPC), metastatic disease remains incurable, despite the introduction of combination approaches
[2]. Targeted systemic therapies with androgen receptor (AR) signaling inhibitors (ARSIs) such as abiraterone or enzalutamide are primarily used to treat relapsed castration-resistant prostate cancer (CRPC)
[3]. Employment of these agents has been shown to be effective, especially in high risk primary metastatic HSPCs patients, and their use is becoming more common
[4,5][4][5]. Although these combination approaches have demonstrated survival benefits, they have also been shown to contribute to the emergence of more aggressive castration-resistant tumors
[6]. The majority of castration- and some ARSI-resistant PCs are characterized by increased AR signaling
[7,8][7][8]. However, alternative forms of castration-resistance have also been identified, including forms that are AR-negative with neuroendocrine-like features
[9], forms indifferent to AR
[10], and forms that are more dependent on alternative signaling pathways such as fibroblast growth factor (FGF)
[11]. Additionally, forms of castration-resistance dependent on other TFs such as the glucocorticoid receptor (GR) or the pluripotent stem cell TF SOX2
[12,13,14,15][12][13][14][15] have been described.
Along with others, we previously reported that the emergence of castration-resistance and AR overexpression are associated with chromatin reprogramming
[16,17,18,19][16][17][18][19]. As with other malignancies, PC is thought to arise from and be driven by oncogenic genetic alterations. However, as with many cancers, PC cannot be explained solely on the basis of genetic alterations
[20]. PC in particular has a relatively low mutational load at presentation
[21] and prostate carcinogenesis is not clearly driven by any particular genetic alterations
[22]. Nevertheless, HSPCs and CRPCs are characterized by the typical cancer hallmarks
[23,24][23][24] that mediate carcinogenesis, disease progression under treatment pressure, and cancer growth beyond the tumor microenvironment (TME).
Flavahan and colleagues first proposed the concept of epigenetic plasticity, by which alterations in the structure of the chromatin, or chromatin states, and epigenetic alterations would be able to confer the full range of cancer hallmarks by altering transcriptional regulation
[25,26][25][26]. Indeed, many of the epigenetic and chromatin regulators that drive normal tissue development and differentiation are well-known oncogenes and tumor suppressors recurrently mutated or aberrantly expressed in different malignancies.
Similarly, it is well known that cellular lineage identity is defined by tightly regulated chromatin-related processes and the epigenetic landscape
[27]. The concept of lineage plasticity is of clinical relevance for PC as it is a common mechanism of resistance following the increasing usage of more potent ARSI in primary disease
[11,28][11][28].
Due to the dynamic sets of chromatin and epigenetic alterations and their reversibility, such alterations can play a significant role in driving both carcinogenesis and progression to treatment-resistant disease. These alterations can also be stochastically different in individual cells
[25]. Therefore, epigenetic plasticity can be the basis of the heterogeneity observed among PC patients, but also between different tumor foci of the same prostate
[29,30,31][29][30][31]. It is apt that epigenetic plasticity can consist of genetic changes that allow normal cells to transform through alteration of lineage commitment and further allow malignant cells to adapt under diverse treatment pressure by mediating lineage plasticity.
In this review, we highlight the contribution of chromatin- and epigenetics-related processes to prostate carcinogenesis and progression to treatment resistance. We explore the role of epigenetic regulation and TFs in lineage commitment and plasticity in the normal prostate and PC cells, the mutations and the altered expression of key genes coding for chromatin-associated proteins, alterations in DNA methylation patterns, and changes in the structure and 3D organization of the chromatin (Figure 2).
Figure 2. Epigenetic plasticity in prostate cancer. Epigenetic dysregulation (light blue boxes) is in the forefront of lineage plasticity as well as in carcinogenesis and therapy resistance. Normal prostate epithelium is renewing at a steady state as terminally differentiated luminal cells are slowly replaced by progenitor cells. As genetic alterations accumulate due to cell division and the normal aging process, driver alterations (red boxes) such as ETS gene fusions or SPOP mutations emerge. The mutational processes lead to less ordered chromatin structure, as characterized by chromatin relaxation at distal regulatory regions, alterations in DNA methylation and histone modifications, and dysregulation of higher order chromatin structures, which alters the binding of key TFs such as AR. Some cells gain stem-like properties, leading to increased proliferation capacity and reduced apoptotic rates, which leads to tumor formation over time. The plasticity of the cellular identity is also in a key role during the emergence of treatment resistance as the cancer cells can repurpose differentiation-promoting transcription factors such as AR into regulatory regions supporting developmental gene expression (seen in castration-resistant prostate adenocarcinoma, CRPC), or transdifferentiate into non-luminal, small cell prostate carcinoma or neuroendocrine prostate cancer (NEPC).
2. Prostate Lineage Commitment and Prostate Cancer Plasticity
Here, we collect evidence of the involvement of lineage commitment and cancer cell plasticity in PC initiation and resistance to treatment, respectively.
The normal epithelium of the prostate secretory acini is composed of a layer of luminal cells oriented towards the acinar lumen and a basal cell layer that also includes rare neuroendocrine cells. The determinants of prostatic epithelial differentiation include master regulators such as AR, FOXA1, and NKX3-1. Talos and colleagues showed that these TFs were sufficient and required for differentiation of iPS cells of mouse fibroblast origin in prostatic cells engrafted in the renal capsule
[32]. Xie et al. showed that AR may also be required in the basal-luminal intermediate cells to produce fully differentiated luminal progeny in adult mice
[33]. The same study found that AR expression is not needed for survival of the luminal cells but is essential for normal luminal differentiation and morphology
[33]. This may be, in part, due to an indirect effect by
NKX3-1, as it is an AR target gene and was not expressed after deletion of AR
[33]. Dutta et al. had previously shown that NKX3-1 is a prostate-specific master transcriptional regulator that can transdifferentiate the seminal vesicle epithelium into prostate epithelium
[34]. Of note, the seminal vesicle epithelium expresses AR and FOXA1 following the introduction of NKX3-1
[34]. Differentiation by NKX3-1 was found to be mediated by histone demethylase UTY and histone methylase G9a
[34]. The stromal component seems to contribute to the re-differentiation process as the urogenital sinus mesenchyme is always co-transplanted with the epithelial cells in these experiments
[32,34][32][34].
AR and FOXA1 have been found to also be important for the regulation of the Homeobox (HOX) A genes
[35]. The HOX genes regulate prostate development in mice
[36] and the paralogous HOX13 genes (
HOXA13,
HOXB13, and
HOXD13) are still expressed in the luminal epithelium of the human adult prostate
[36]. In recent single cell analyses, Guo et al. showed that terminally differentiated luminal cells express
NKX3-1 and
HOXB13 together with
AR [37].
Experiments in mice have shown that the prostate epithelium displays a regenerative capacity following repeated cycles of androgen deprivation
[38,39][38][39]. This led to the hypothesis that the prostate epithelium may also harbor stem cells responsible for tissue renewal
[38]. Recent advances in single cell sequencing have facilitated a more detailed characterization of the complexity of prostate cell types without having to rely on restricted sets of cell surface markers for their classification
[40]. These single cell approaches have led to the identification of club cells (KRT5
−, KRT8
−, and SCGB1A1
+) and hillock cells (KRT5
+, KRT14
−, KRT13
+) with stem cell potential
[40]. Similar to lung club and hillock cells, progenitor-like cells that are able to differentiate into goblet cells and ciliated cells, these newly described prostate cells are efficient at reconstituting the prostate epithelium in in vivo studies following androgen deprivation, as shown by Karthaus et al.
[41], although their transdifferentiation to basal cells is rather limited
[40,41,42][40][41][42].
Maitland and Collins have reviewed how the prostate epithelium constantly renews and the terminal differentiation of luminal cells is associated with a higher rate of apoptosis
[43]. The existence of intermediate cells, which are cells with both basal and luminal features in the prostate epithelium, has been acknowledged for a long time
[44]. In lineage tracing studies with mice, both basal and luminal cells have been shown to contribute to the renewal of prostate epithelium during androgen deprivation and add back cycles, indicating that there are several degrees of stemness within the adult prostate epithelium
[37,41,45][37][41][45]. In previous studies, both luminal and basal stem or progenitor cells have been shown to act as the PC-initiating cells, but PC is generally characterized by the absence of basal cells
[45,46,47,48,49,50][45][46][47][48][49][50].
Lineage plasticity refers to the reversal of the process of lineage commitment, either by dedifferentiation of the more differentiated cells or, in extreme cases, transdifferentiation directly (or via an intermediate) to another epithelial cell lineage. This is a process that is unlikely to occur in the normal prostate under physiological conditions. However, pathological processes or stressors such as inflammation may alter this scheme
[51]. Alterations in key epigenetic- and chromatin-associated or environmental stress factors can disrupt the epigenetic homeostasis of normal cells, which may lead to lineage plasticity, differentiation arrest, and accumulation of undifferentiated cells in transition. The role of the TME in determining such alterations is not well known, but factors such as hypoxia
[52] and other metabolic stressors
[53] have been associated with more aggressive PC phenotypes. A well-established example of metabolic stress leading to epigenetic changes is the dysregulation of one-carbon metabolism and its effects on both DNA and histone methylation
[53,54][53][54].
Lineage plasticity during prostate carcinogenesis is poorly understood. However, for example, the upregulation of c-MYC is a common early event
[55] and has been implicated in the gain of stem cell properties and repression of differentiation
[56]. Normal luminal cells repress MYC expression via AR/β-catenin/TCF-4 signaling in the presence of androgens, leading to growth arrest, but overexpressing c-MYC rescues cell proliferation
[57,58][57][58]. This highlights the opposing roles of AR as the main differentiation factor and growth suppressor in the non-transformed secretory luminal prostate cells and as a prominent driver of PC cells’ proliferation, as discussed below.
Cells with regenerative potential within the normal prostate tissue have been hypothesized to be the cells of origin for PC
[59]. Indeed, Song et al. recently described club-like PC cells in primary PC specimens
[60]. These cells are transcriptionally similar but have higher AR expression and an enhanced androgen signaling signature when compared to club cells from normal prostates
[60] (
Table S1). The overexpression of AR in these cells is consistent with the modulation and expansion of the AR cistrome, which is a well-documented feature of PC initiation
[61]. Under this scenario, AR overexpression induces changes in the AR transcriptional program, leading to cell survival and proliferation, possibly via the alteration of the pioneer activity of cooperative TFs of the AR such as HOXB13
[61]. Interestingly, recent analyses with mice expressing the F133V mutant form of Speckle Type BTB/POZ Protein (
SPOP) showed that this mutation is sufficient to modify chromatin accessibility and binding of AR and FOXA1 at PC specific genomic sites
[61,62][61][62]. Further expansion of the AR cistrome has also been reported during progression to CRPC
[63,64][63][64].
The role of some luminal progenitor cells as possible cells of origin for PC is also supported by studies in Tmprss2-CreER;Pten
flox/flox mice
[37]. The TACSTD2/Trop2 (encoded by
Tacstd2)-expressing luminal progenitor cells characterized by Guo et al. have similar transcriptional features to the club cells characterized by Henry et al., Karthaus et al., and Song et al., based on the presence of markers such as
PIGR,
PSCA, and
KRT4 [37,40,41,60][37][40][41][60] and high expression of TACSTD2
[41,60][41][60] (
Table S1). Similarly, Kwon et al. found that TACSTD2+ luminal prostate cells were more efficiently transformed in vitro in an organoid-forming assay than TACSTD2−luminal cells
[65]. The authors also showed that TACSTD2+ luminal cells express SOX2 and display remarkable plasticity by transdifferentiation to de novo neuroendocrine PC (NEPC), even in the absence of selection pressure from treatment with antiandrogens
[65].
Although the above-described studies suggest that multiple subpopulations of cells can give rise to PC
[37,45[37][45][46][47][48][50],
46,47,48,50], their transcriptional features seem to converge toward features of luminal progenitor or club-like cells
[37,41,48,49,60][37][41][48][49][60]. Undoubtedly, the cell type of origin of PC can have clinically relevant consequences in terms of tumor trajectory and prognosis
[50].
Treatment of metastatic HSPCs with ARSI in combination with ADT imposes a strong negative selection pressure on the cancer cell population. Multiple treatment resistance mechanisms arise from genetic and epigenetic alterations, leading to a rewiring of alternative bypass pathways driving tumor cell growth
[66], but the majority of CRPCs remain AR-dependent through maintenance of AR signaling
[67]. This imposes chromatin reconfiguration and epigenetic plasticity in ARSI-resistant PC cells, leading to heterogeneous CRPC phenotypes more or less dependent on AR signaling, with both luminal and basal features
[68,69,70][68][69][70]. A recent single cell analysis by He and colleagues showed that different AR transcript variants are ubiquitously expressed prior to treatment with ARSI
[71]. They also found that ARSI drives resistance pathways such as epithelial to mesenchymal transition (EMT), a form of lineage plasticity generally associated with metastatic disease
[71].
Chromatin reprogramming is a key feature of lineage plasticity during treatment resistance, for example by allowing reactivation of pathways that are normally active only during prostate tissue development
[19,64,72][19][64][72]. Zhang et al. found that loss of the chromatin-modifying helicase CHD1 confers enzalutamide resistance by inducing marked changes in chromatin accessibility, along with transcriptional rewiring and upregulation of TFs
NR3C1 (GR),
POU3F2 (BRN2),
TBX2, and
NR2F1, leading to gene expression changes including a reduction in luminal and an increase in EMT markers
[68].
Cellular plasticity has also been associated with alterations in tumor suppressor genes and chromatin reconfiguration, exemplified by the deletion of
TP53 and
RB1 that can lead to treatment resistance by allowing for diverse transcriptional programs
[6,69,70][6][69][70]. P53 and pRb (encoded by
TP53 and
RB1, respectively) cooperate to suppress expression of
SOX2, a well-known factor of pluripotency, which, in part, explains the association of these factors with lineage plasticity
[15,73][15][73]. Using an isogenic model of HSPC (LNCaP) and CRPC (C4-2), Mandingo et al. showed that
RB1 loss leads to a reconfiguration of E2F activity for increasing the production of antioxidants, protecting cells against doxorubicin in the CRPC cells that is not observed in the HSPC cells
[74]. These findings highlight how transcriptional reprogramming by the same TF can be modulated in different disease stages
[74].
In the context of ARSI resistance, so-called treatment-induced neuroendocrine prostate cancers (t-NEPCs) (reviewed by Rubin et al. and Kaarijärvi et al.
[28,75][28][75]) represent a tumor phenotype mediated by chromatin and transcriptional plasticity, which, in turn, drives cellular plasticity. During therapies with ARSIs, AR-dependent cells are depleted, while transcriptionally and epigenetically heterogeneous AR-indifferent cells increase their growth and survival potential
[11,68,69][11][68][69]. Increased mutation frequency in tumor suppressor genes such as
TP53 and
RB1 suggests that t-NEPCs require these alterations for tumor selective reprogramming
[6,15,73,76][6][15][73][76]. The t-NEPC tumors acquire neuroendocrine features and suppress luminal transcriptional programs via the upregulation of SOX2 and EZH2 in mouse models
[15,73][15][73] and the upregulation of
LHX2 and
ISL1 in human cells, which goes hand in hand with increased chromatin accessibility and the transcriptional output of neuroendocrine lineage-specific genes
[76]. FOXA1, a well-known AR pioneering factor and a mediator of the transcriptional output in prostate adenocarcinoma, has also been shown to mediate NE-specific transcription in NEPC
[77]. NE-specific genes are repressed in normal differentiated prostate epithelial cells, but during transdifferentiation following the activation of NE TFs such as ASCL1 and NKX2-1 in these cells, the chromatin structure in their proximal regulatory regions is rewired and allows for their expression
[77]. The loss of the transcriptional repressor REST, which normally represses
ASCL1 and
NKX2-1, has been implicated in this process
[78]. Active AR signaling keeps REST stable via inhibition of E-ubiquitin ligase β-TrcP, so ARSIs may contribute to the downregulation of REST
[79].
Ultimately, PC progression is a continuous process resulting from the selection pressures of treatment upon the tumor, the TME, and the host immune system
[80].
3. Mutations and Expression Dysregulation of Genes Coding for Chromatin-Associated Factors
Many studies have illustrated how chromatin-associated factors affect both lineage commitment and lineage plasticity, which is likely to drive tumor progression
[64,68,77][64][68][77]. PC cell lineage plasticity is driven by genetic alterations, gene expression changes, and the altered activity of chromatin-associated and epigenetic regulators
[6,25,81,82,83,84][6][25][81][82][83][84]. These chromatin-associated factors can be broadly categorized by their protein function into TFs, transcriptional co-regulators, chromatin modifiers, and genes involved in mRNA transcript synthesis or processing. In this section, we performed an analysis of their mutational and expression patterns in prostate carcinogenesis and the development of treatment resistance. To gain a comprehensive view on chromatin-associated factors, we utilized a list of 2754 genes previously annotated to the above-mentioned functional groupings
[85] and queried their alteration status in publicly available PC patient datasets.
Somatic mutation data from the International Cancer Genome Consortium (ICGC) prostate adenocarcinoma datasets
[86] showed that only 34 were genes mutated in more than 1% of patients with potential protein function-altering effects, reflecting the overall low mutation frequency of these tumors (
Table S2). As the ICGC cohort largely consists of early stage primary prostate tumors, we repeated the analysis using two metastatic CRPC patient cohorts from Robinson et al. and Grasso et al. to include 150 pre- or post-ARSI mCRPC biopsies
[87] and 50 heavily pre-treated lethal CRPCs exposed only to first generation ADT
[88], respectively. We identified an additional 78 genes coding for chromatin-associated proteins in more than 2% of these patients with protein-altering mutations (
Table S2). Altogether, 18 genes, including
TP53,
FOXA1,
SPOP, and
CDK12, were recurrently mutated in both the ICGC early stage and advanced CRPC tumors. The proportions of mutated genes from different functional categories did not differ significantly in the early versus treatment-resistant disease stages.
We further assessed the expression of the 94 recurrently mutated genes (
Table S2) during prostate carcinogenesis and the development of treatment resistance in a previously published RNA-sequencing dataset of benign prostatic hyperplasia (BPH,
n = 10), untreated PC (
n = 16), and CRPC (
n = 11) samples (
Figure 3)
[89]. We observed expression changes in several known PC driver genes including the loss of expression of
TP53 and
SPOP during the transition from untreated PC to CRPC, as well as the increased expression of
AR,
BRCA2,
KDM6A, and
CDK12 in CRPC. Interestingly, TFs that were mutated in PC and CRPC or only in CRPC showed higher expression in those categories, with the exception of decreased
TP53 expression in CRPC, which is potentially a consequence of its wide-ranging tumor suppressive functions
[70]. TFs mutated only in CRPC typically had lower expression in CRPC, with the exception of
AR and
EHMT1, a histone methyltransferase with potential tumor suppressive function in the prostate
[90]. Nearly all recurrently mutated chromatin modifiers had the highest expression in CRPC irrespective of the mutation frequency (
Figure 3). Among the genes involved in epigenetic, chromatin, or gene regulation that were not recurrently mutated, 68 (2.6%) genes were significantly upregulated during prostate carcinogenesis (PC vs. BPH), while 80 (3.0%) genes were upregulated during the development of treatment resistance (CRPC vs. PC). In contrast, 34 (1.3%) non-recurrently mutated genes were downregulated during prostate carcinogenesis and 77 (2.9%) were downregulated during treatment resistance. TFs were the most common group of genes to be aberrantly expressed, particularly amongst the genes upregulated during carcinogenesis (46 of 68 genes, 68%) or downregulated during carcinogenesis (18 of 34 genes, 53%).
Figure 3. Expression changes of 94 recurrently mutated genes coding for chromatin-associated proteins in prostate carcinogenesis and development of treatment resistance. Row-scaled log2 mean expression values for each gene (the rows) are shown in a heatmap for benign prostatic hyperplasia (BPH), untreated prostate cancer (PC), and castration-resistant prostate cancer (CRPC) patient samples. The names of frequently studied genes are shown in bold. Each row is annotated with the disease stage, in which the gene is found to be recurrently mutated (either early stage (PC), late stage (CRPC), or both early and late stage). The rows are also annotated with the functional category of each gene. On the right, two columns show the mutation frequency of the gene in early- and late-stage disease in four categories (0–1%, 1–5%, 5–10%, and >10%).
Our analysis further highlighted alterations in groups of genes with functions related to PC lineage plasticity. These included the master transcriptional regulators
FOXA1 and
AR [91,92,93,94][91][92][93][94], the
KMT2A-D transcriptional co-activators involved in development
[95], and a number of chromodomain (CHD) genes, including
CHD1, involved in chromatin remodeling and transcription activation
[96]. Loss of
CHD1 has been reported in 15% of HSPCs and 17% of CRPCs
[97] and
CHD1 loss has been implicated in PC cell chromatin rewiring with tumor-suppressing functions
[68] and increased sensitivity of PC tumors to DNA damage
[98]. The simultaneous dysregulation of tumor suppressors such as
TP53 and
RB1 in our recurrently mutated gene list has also been shown to promote PC plasticity
[99].
In addition to genetic variation in coding regions as explored here, recent studies have also shown that alterations within noncoding regions affect chromatin conformation and transcriptional regulation in PC. By studying somatic single nucleotide variants and germline single nucleotide polymorphisms in the
cis-regulatory elements of prostate tumors, Mazrooei et al. found that these variants are specifically enriched in the cistromes of master transcriptional regulators AR, FOXA1, HOXB13, and SOX9
[100]. This implicates noncoding variation in the dysregulation of the chromatin binding activity of these factors and therefore, in prostate carcinogenesis.
AR expression has further been found to be modulated through a somatically gained upstream enhancer in CRPC
[101] and its binding sites have been shown to be lost or gained due to the activity of FOXA1 and HOXB13 during carcinogenesis
[61]. The chromatin regions linked to altered gene regulation in PC have also been shown to be enriched in PC predisposing genetic variants. Pomerantz et al. analyzed prostate lineage-specific enhancers and promoters marked by a combination of histone modifications and PC-specific TFs AR, FOXA1, and HOXB13, and showed that these were enriched for genetic variants linked to increased PC risk, depicting an active epigenetic state
[64]. In addition, they showed that prostate lineage-specific regulatory regions exhibit active epigenetic states that are predicted to increase somatic mutational burden
[64]. Collectively, these studies highlight the importance of the noncoding genome in PC lineage commitment and plasticity.
4. Dysregulation of Chromatin States through Histone Modification in Prostate Cancer
The organization and structure of the chromatin is the result of the activity of various proteins interacting with DNA. The negatively charged 2 nm DNA double helix wraps around positively charged histone octamers, creating the 10 nm fiber of adjacent nucleosomes. By continuous folding and looping, this primary structure is further compacted to higher order structures to form the 3D organization of the chromatin. Both DNA replication and RNA transcription require accessibility to DNA by reversible decondensation of these chromatin structures. Chromatin compaction and decondensation at certain loci are tightly regulated and cell type-specific. This is primarily achieved through the ability of chromatin-associated proteins to write (e.g., with acetyl transferase activity), read (e.g., with bromodomains), and erase (e.g., with demethylases) the histone code.
The histone code is the set of histone modifications (methylation, acetylation, etc.) defining the state of the chromatin. For instance, a chromatin region marked by methylation in lysine 4 of the histone 3 tail (typically marked as H3K4me1) and H3K36me1 could be defined as a promoter
[131,132,133][102][103][104], but additional histone markers such as H3K27 Acetylation (H2K27Ac), H3K4me3, H3K27me3, or the presence of RNAPol2 could determine whether the promoter is actually active or poised
[131,132,133,134][102][103][104][105]. ChromHMM is a tool that converts the histone code to an annotated chromatin state
[133,134][104][105]. The presence of histone marks H3K4me1, H3K4me3, H3K27me3, H3K27Ac, H3K36me3, and H3K9me3 are used to annotate the chromatin state based on the Roadmap Epigenomics 18-state expanded model
[135][106]. For example, in this model, active transcription is marked with H3K4me3, whereas repressed chromatin is marked with H3K27me3 and H3K9me3 heterochromatin marks
[135][106]
PC cells exhibit significant dysregulation of chromatin states (
Figure 4). Firstly, androgen stimulation can trigger several changes in the chromatin states of PC cells. H3 acetylation has been demonstrated to accumulate at the
KLK3 (encoding PSA) promoter and enhancer regions upon androgen treatment in LNCaP cells
[136,137,138][107][108][109]. In addition, the presence of androgens triggers the recruitment of histone acetylases (HATs) CBP and p300 to the regulatory regions of
KLK3, thus changing the chromatin state to a more active one via their acetyl transferase activity
[138,139][109][110]. p300 has also been shown to protect AR from degradation in PTEN-depleted PC mouse models, in which p300 was essential for the expression of AR target genes
[140][111]. Moreover, combined inhibition of p300 and CBP significantly decreases the expression of programmed death-ligand 1 (PD-L1), which increases the efficacy of immunotherapy that aims to reactivate T-cells by PD-L1 blockade
[141][112].
Figure 4. Forms of epigenetic and chromatin dysregulation in PC. DNA methylation is increased at specific regulatory regions, such as the
GSTP1 promoter, but is generally reduced genome-wide. Aberrant histone modifications especially at distal regulatory regions, often harboring binding motifs for key transcription factors (TFs) such as the androgen receptor (AR), shift to a more active state. This also occurs due to the action of transcriptional coactivators and chromatin readers such as BRD4, leading to recruitment of these TFs to previously repressed regions. Loss of
CHD1 induces chromatin rewiring, and pioneer factors such as FOXA1 are able to bind to repressed regions and recruit other TFs and histone remodeling complexes. Dysregulation of chromatin accessibility is partly due to increased activity of the SWI/SNF remodeling complex, which shifts the chromatin towards more permissive states, a process termed chromatin relaxation. Figure created with
BioRender.com.
In addition to HATs, histone demethylases (JHDM2A, JMJD2C, and LSD1) have been reported to exhibit androgen stimulation-dependent recruitment to AR target genes such as
KLK3,
TMPRSS2, and
NKX3-1, resulting in less methylated H3K9 and an active chromatin state
[142,143,144][113][114][115]. Low levels of H3K9me2 have been shown to predict poor prognosis in PC patients
[145][116], whereas higher IHC staining of H3K4me1 is associated with a higher risk of biochemical recurrence
[146][117]. High levels of H3K4me2 are also associated with shorter time to biochemical recurrence after radical prostatectomy
[147][118]. Despite active histone marks being associated with prognosis, suppressive chromatin states marked by H3K27me3 seem to be more abundant than active chromatin states marked by H3K4me3 in PC cells
[148,149][119][120]. Ke et al. performed a genome-wide profiling of H3K27me3 and H3K4me3 in PC and normal epithelial cells and showed that while the number of each modification is relatively similar in both cells, the genome-wide profiles of these marks vary so that the genes marked by H3K27me3 or by both H3K27me3 and H3K4me3 are approximately 70% different
[148][119]. Genes marked by H3K27me3 in PC cells are enriched with developmental functions, which was not the case for normal epithelial cells. This suggests that epigenetic reprogramming of developmental genes occurs during PC carcinogenesis. Examining genes that exhibited differential chromatin state in PC compared to normal epithelial cells, Ke et al. observed that frequent switches between H3K4me3 and H3K27me3 caused the strongest expression changes in genes, including most of the HOXA genes cluster
[148][119].
Changes in chromatin states have also been associated with PC progression. The expression of the aforementioned HATs p300 and CBP was shown to be significantly higher in metastatic CRPC lesions than in local tumors
[150][121]. p300 functions to stabilize a histone demethylase JMJD1A, which is upregulated in CRPC in comparison to PC, increasing JMJD1A recruitment to AR target genes and thus, facilitating the conversion of the chromatin state into an active one
[151][122]. The levels of H3K4me1, H3K4me2, and H3K4me3 based on IHC increase during PC progression
[146][117]. This is interesting as the binding of FOXA1 to its targets strongly depends on the distribution of H3K4me1 and H3K4me2
[93,152][93][123]. Selective accumulation of these histone marks is found at the AR-bound enhancers of M-phase cell cycle genes, such as
CDK1 and
UBEC2, in CRPC tissues and cell models when compared to androgen-dependent cell lines and tumor tissue samples
[153][124]. Consistently analyzing global levels of H3K27me3 by IHC, Pellakuru et al. showed that the levels decrease as markers of aggressive disease increase
[154][125]. This decline in H3K27me3 levels is also observed when more differentiated normal prostatic luminal epithelial cells are compared to more stem-like basal cells
[154][125].
These findings suggest that overall active chromatin state increases during PC progression and that the epigenome is reprogrammed to resemble stem-like cells. Supporting this, genes repressed and marked by H3K27me3 levels in normal prostate tissue and derepressed during PC progression are developmental regulators and homeobox proteins
[155,156][126][127]. Yu et al. revealed that such epigenetic reprogramming might partly originate from the aberrant function of EZH2
[155,156][126][127]. Many studies have now shown that EZH2 is overexpressed in PC and even more so in CRPC
[157,158,159,160][128][129][130][131]. However, it has been shown to also have a PRC2-independent role in CRPC
[161][132]. Mechanistically, EZH2 functions as a co-activator for important transcription factors such as AR upon phosphorylation
[161][132], and as a part of the PRC2 complex to induce silenced chromatin states.
The combination of aberrant functions of HATs, demethylases, and EZH2 all seem to contribute to increased active chromatin states in PC progression and treatment resistance. In line with this, Pomerantz et al. showed that the profiles of active enhancers defined by H3K27ac are similar in CRPC and prostate-specific fetal tissue, implicating reactivation of the developmental transcriptional programs in late-stage disease
[64].
5. Dysregulation of Chromatin Accessibility
Chromatin accessibility is facilitated by the chromatin remodeling complexes by either ejecting histones or full nucleosomes or by sliding nucleosomes to different positions in the DNA
[162,163][133][134]. Several publications have linked the chromatin remodeling complex SWI/SNF to PC carcinogenesis and progression. BRG1 (also known as SMARCA4), a core subunit of the SWI/SNF complex, is required for chromatin accessibility
[163][134]. Previous studies have reported the increased expression of BRG1 in primary PC
[164,165,166,167][135][136][137][138]. Higher BRG1 expression has been associated with larger tumor mass and increased invasion of PC-3 cells
[166][137]. Inhibition of BRG1 in PC-3 cell lines has also been demonstrated to decrease cell proliferation and induce apoptosis
[168][139]. Inhibition of BRG1 in mouse xenografts increased the survival of the mice and decreased tumor growth
[168][139]. Moreover, its increased expression has been associated with more malignant features
[164,166,167,168][135][137][138][139]. Expression of BRG1 has a positive correlation with Gleason score, as it increases as the tumor progresses and the highest expression is observed in NEPC
[164,166,167,168][135][137][138][139].
PTEN loss is associated with PC progression to CRPC
[169][140] and has been demonstrated to stabilize, and thus, increase, the amount of BRG1 protein
[167][138]. IHC staining of radical prostatectomy tissues demonstrated that in tumors exhibiting low PTEN expression, high BRG1 expression correlated with poor outcome
[167][138]. PTEN and BRG1 were indeed shown to be synthetic lethal in a panel of PC cell lines and in a mouse model
[167][138]. Open chromatin regions were reduced by approximately 60% when BRG1 was depleted in cells lacking PTEN, highlighting the importance of this SWI/SNF subunit to chromatin accessibility
[167][138].
In addition to BRG1, increased expression of SWI/SNF complex subunits BAF57 and BAF155 has been described in PC
[170,171,172][141][142][143]. Inhibition of BAF57 results in the inhibition of AR-dependent genes in PC cell lines
[172][143]. IHC studies of primary tumors and metastatic specimens have revealed that BAF57 expression positively correlates with Gleason score
[170,172][141][143]. Moreover, cell line studies with LNCaP and C4-2 (a CRPC cell line) showed that BAF57 mediates ligand-independent AR activity, thus inducing the expression of migration- and metastasis-related genes in the absence of androgens
[170][141]. TMA IHC analysis of BAF155 revealed the highest staining in recurrent and metastatic prostate cancer specimens and a positive correlation with Gleason score
[171][142].
Interestingly, the differential expression of certain SWI/SNF complex subunits has been linked to the progression of PC into NEPC
[164][135]. For example, the RNA and protein expression of BRG1 and neuron-specific SWI/SNF subunit proteins, such as BAF53B and BAF45B, exhibited significantly higher expression in NEPC than in CRPC
[164][135]. This suggests that differential chromatin remodeling is required for the progression of prostate cancer to CRPC and NEPC.
Abnormal activation of chromatin-related proteins such as chromatin readers and the SWI/SNF complex and its subunits suggests an overall increase in accessible chromatin in gene regulatory regions in CRPC. Indeed, multiple studies have shown that relaxed chromatin correlates with PC progression
[19,72,145,154,173,174][19][72][116][125][144][145] and that different markers for open chromatin can function as predictors for tumor recurrence
[175,176][146][147]. Moreover, Bromodomain and Extraterminal domain (BET) proteins BRD2 and BRD4, chromatin readers that recognize acetylated histones and function to promote transcription, have been shown to maintain accessible chromatin in PC
[19,177,178,179][19][148][149][150] and to be overexpressed in CRPC
[19]. These proteins have been demonstrated to mediate tumor growth in both CRPC cell lines and in vivo experiments with xenograft models
[180,181,182,183][151][152][153][154]. In addition, BET proteins are associated with the increased expression of migration- and invasion-related genes in CRPC
[184,185][155][156]. Although BRD4 can co-regulate invasion-related genes together with ERG
[184][155], it seems that the main mechanism of function for these proteins is to enhance AR-mediated transcriptional regulation
[180[151][152][153][154][157],
181,182,183,186], highlighted by the fact that AR-negative cells are not sensitive to BET inhibitors
[180][151]. AR overexpression has been associated with the chromatin accessibility maintenance activities of BET proteins BRD4 and BRD2 together with AR-regulated bromodomain-containing protein ATAD2
[19], indicating that these proteins remodel chromatin to promote AR-mediated gene regulation.
Similarly, the SWI/SNF complex has been shown to be an important transcriptional cofactor for AR activity
[172,187,188,189,190][143][158][159][160][161]. Multiple studies have shown that SWI/SNF activity is critical for the AR-mediated activation of
KLK3 and
TMPRSS2 transcription
[172,188,190][143][159][161]. In addition, SWI/SNF activity positively regulates the expression of other androgen regulated genes, such as
FKBP5 and
KLK2 [188][159]. Inhibition of different subunits of the SWI/SNF complex inhibits the proliferation of AR-dependent PC cells
[188[159][160],
189], but different subunits of the SWI/SNF complex seem to promote AR-mediated activation on different target genes
[190][161]. This highlights the dynamic nature of SWI/SNF function in mediating transcriptional plasticity in response to potential PC treatments.
The importance of bromodomain activity in modifying the structure of the chromatin in PC has been reviewed previously
[191][162]. Here, we add that Welti et al. recently reported that targeting p300/CBP, histone acetyltransferases with known AR transcriptional coactivator function, could be a new therapeutic strategy for CRPC, and describe a novel small-molecule inhibitor (CCS1477) of the p300/CBP conserved bromodomain
[150][121].
These findings demonstrate the abnormal function of chromatin accessibility regulators in PC and in the progression to treatment resistance. These changes precede the reprogramming of TF chromatin binding, especially the binding of the AR.