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

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

Inhibition of CSF-1R and IL-6R prevents conversion of cDC2s into immune incompetent tumor-induced DC3s boosting DC-driven therapy potential.

The human dendritic cell (DC) family has recently been expanded by CD1c+CD14+CD163+ DCs, introduced as DC3s. DC3s are found in tumors and peripheral blood of cancer patients. Here, we report elevated frequencies of CD14+ cDC2s, which restore to normal frequencies after tumor resection, in non-small cell lung cancer patients. These CD14+ cDC2s phenotypically resemble DC3s and exhibit increased PD-L1, MERTK, IL-10, and IDO expression, consistent with inferior T cell activation ability compared with CD14- cDC2s. In melanoma patients undergoing CD1c+ DC vaccinations, increased CD1c+CD14+ DC frequencies correlate with reduced survival. We demonstrate conversion of CD5+/-CD1c+CD14- cDC2s to CD14+ cDC2s by tumor-associated factors, whereas monocytes failed to express CD1c under similar conditions. Targeted proteomics identified IL-6 and M-CSF as dominant drivers, and we show that IL-6R and CSF1R inhibition prevents tumor-induced CD14+ cDC2s. Together, this indicates cDC2s as direct pre-cursors of DC3-like CD1c+CD14+ DCs and provides insights into the importance and modulation of CD14+ DC3s in anti-tumor immune responses.

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

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

Myeloid-derived suppressor cells and tolerogenic dendritic cells are distinctively induced by PI3K and Wnt signaling pathways.

Imbalanced immune responses are a prominent hallmark of cancer and autoimmunity. Myeloid cells can be overly suppressive, inhibiting protective immune responses or inactive not controlling autoreactive immune cells. Understanding the mechanisms that induce suppressive myeloid cells, such as myeloid-derived suppressor cells (MDSCs) and tolerogenic dendritic cells (TolDCs), can facilitate the development of immune-restoring therapeutic approaches. MDSCs are a major barrier for effective cancer immunotherapy by suppressing antitumor immune responses in cancer patients. TolDCs are administered to patients to promote immune tolerance with the intent to control autoimmune disease. Here, we investigated the development and suppressive/tolerogenic activity of human MDSCs and TolDCs to gain insight into signaling pathways that drive immunosuppression in these different myeloid subsets. Moreover, monocyte-derived MDSCs (M-MDSCs) generated in vitro were compared to M-MDSCs isolated from head-and-neck squamous cell carcinoma patients. PI3K-AKT signaling was identified as being crucial for the induction of human M-MDSCs. PI3K inhibition prevented the downregulation of HLA-DR and the upregulation of reactive oxygen species and MerTK. In addition, we show that the suppressive activity of dexamethasone-induced TolDCs is induced by ß-catenin-dependent Wnt signaling. The identification of PI3K-AKT and Wnt signal transduction pathways as respective inducers of the immunomodulatory capacity of M-MDSCs and TolDCs provides opportunities to overcome suppressive myeloid cells in cancer patients and optimize therapeutic application of TolDCs. Lastly, the observed similarities between generated- and patient-derived M-MDSCs support the use of in vitro-generated M-MDSCs as powerful model to investigate the functionality of human MDSCs.

Grade Group 1 Prostate Cancers Exhibit Tumor-defining Androgen Receptor-driven Programs.

Grade group 1 (GG1) primary prostate cancers with a pathologic Gleason score of 6 are considered indolent and generally not associated with fatal outcomes, so treatment is not indicated for most cases. These low-grade cancers have an overall negligible risk of locoregional progression and metastasis to distant organs, which is why there is an ongoing debate about whether these lesions should be reclassified as "noncancerous". However, the underlying molecular activity of key disease drivers, such as the androgen receptor (AR), have thus far not been thoroughly characterized in low-grade tumors. Therefore, we set out to delineate the AR chromatin-binding landscape in low-grade GG1 prostate cancers to gain insights into whether these AR-driven programs are actually tumor-specific or are normal prostate epithelium-like. These analyses showed that GG1 tumors do not harbor a distinct AR cistrome and, similar to higher-grade cancers, AR preferentially binds to tumor-defining cis-regulatory elements. Furthermore, the enhancer activity of these regions and the expression of their respective target genes were not significantly different in GG1 tumors. From an epigenetic perspective, this finding supports the cancer designation currently given to these low-grade tumors and clearly distinguishes them from noncancerous benign tissue. PATIENT SUMMARY: We characterized the molecular activity of the androgen receptor protein, which drives prostate cancer disease, in low-grade tumors. Our results show that these tumors are true cancers and are clearly separate from benign prostate tissue despite their low clinical aggressiveness.

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.

Drug-Induced Epigenomic Plasticity Reprograms Circadian Rhythm Regulation to Drive Prostate Cancer toward 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 multiomics analyses on tissues isolated before and after 3 months of AR-targeting enzalutamide monotherapy from patients with high-risk prostate cancer enrolled in a neoadjuvant clinical trial. Transcriptomic analyses demonstrated that AR inhibition drove tumors toward a neuroendocrine-like disease state. Additionally, epigenomic profiling revealed massive enzalutamide-induced reprogramming of pioneer factor FOXA1 from inactive chromatin sites toward active cis-regulatory elements that dictate prosurvival signals. Notably, treatment-induced FOXA1 sites were enriched for the circadian clock component ARNTL. Posttreatment ARNTL levels were associated with patients' clinical outcomes, and ARNTL knockout strongly decreased prostate cancer cell growth. Our data highlight a remarkable cistromic plasticity of FOXA1 following AR-targeted therapy and revealed an acquired dependency on the circadian regulator ARNTL, a novel candidate therapeutic target. SIGNIFICANCE: Understanding how prostate cancers adapt to AR-targeted interventions is critical for identifying novel drug targets to improve the clinical management of treatment-resistant disease. Our study revealed an enzalutamide-induced epigenomic plasticity toward prosurvival signaling and uncovered the circadian regulator ARNTL as an acquired vulnerability after AR inhibition, presenting a novel lead for therapeutic development. See related commentary by Zhang et al., p. 2017. This article is highlighted in the In This Issue feature, p. 2007.

An androgen receptor switch underlies lineage infidelity in treatment-resistant prostate cancer.

Cancers adapt to increasingly potent targeted therapies by reprogramming their phenotype. Here we investigated such a phenomenon in prostate cancer, in which tumours can escape epithelial lineage confinement and transition to a high-plasticity state as an adaptive response to potent androgen receptor (AR) antagonism. We found that AR activity can be maintained as tumours adopt alternative lineage identities, with changes in chromatin architecture guiding AR transcriptional rerouting. The epigenetic regulator enhancer of zeste homologue 2 (EZH2) co-occupies the reprogrammed AR cistrome to transcriptionally modulate stem cell and neuronal gene networks-granting privileges associated with both fates. This function of EZH2 was associated with T350 phosphorylation and establishment of a non-canonical polycomb subcomplex. Our study provides mechanistic insights into the plasticity of the lineage-infidelity state governed by AR reprogramming that enabled us to redirect cell fate by modulating EZH2 and AR, highlighting the clinical potential of reversing resistance phenotypes.

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.

The androgen receptor is a tumor suppressor in estrogen receptor-positive breast cancer.

The role of the androgen receptor (AR) in estrogen receptor (ER)-α-positive breast cancer is controversial, constraining implementation of AR-directed therapies. Using a diverse, clinically relevant panel of cell-line and patient-derived models, we demonstrate that AR activation, not suppression, exerts potent antitumor activity in multiple disease contexts, including resistance to standard-of-care ER and CDK4/6 inhibitors. Notably, AR agonists combined with standard-of-care agents enhanced therapeutic responses. Mechanistically, agonist activation of AR altered the genomic distribution of ER and essential co-activators (p300, SRC-3), resulting in repression of ER-regulated cell cycle genes and upregulation of AR target genes, including known tumor suppressors. A gene signature of AR activity positively predicted disease survival in multiple clinical ER-positive breast cancer cohorts. These findings provide unambiguous evidence that AR has a tumor suppressor role in ER-positive breast cancer and support AR agonism as the optimal AR-directed treatment strategy, revealing a rational therapeutic opportunity.

Androgen modulation of XBP1 is functionally driving part of the AR transcriptional program.

Prostate cancer development and progression is largely dependent on androgen receptor (AR) signaling. AR is a hormone-dependent transcription factor, which binds to thousands of sites throughout the human genome to regulate expression of directly responsive genes, including pro-survival genes that enable tumor cells to cope with increased cellular stress. ERN1 and XBP1 - two key players of the unfolded protein response (UPR) - are among such stress-associated genes. Here, we show that XBP1 levels in primary prostate cancer are associated with biochemical recurrence in five independent cohorts. Patients who received AR-targeted therapies had significantly lower XBP1 expression, whereas expression of the active form of XBP1 (XBP1s) was elevated. In vitro results show that AR-induced ERN1 expression led to increased XBP1s mRNA and protein levels. Furthermore, ChIP-seq analysis revealed that XBP1s binds enhancers upon stress stimuli regulating genes involved in UPR processes, eIF2 signaling and protein ubiquitination. We further demonstrate genomic overlap of AR- and XBP1s-binding sites, suggesting genomic conversion of the two signaling cascades. Transcriptomic effects of XBP1 were further studied by knockdown experiments, which lead to decreased expression of androgen-responsive genes and UPR genes. These results suggest a two-step mechanism of gene regulation, which involves androgen-induced expression of ERN1, thereby enhancing XBP1 splicing and transcriptional activity. This signaling cascade may prepare the cells for the increased protein folding, mRNA decay and translation that accompanies AR-regulated tumor cell proliferation.

TLE3 loss confers AR inhibitor resistance by facilitating GR-mediated human prostate cancer cell growth.

Androgen receptor (AR) inhibitors represent the mainstay of prostate cancer treatment. In a genome-wide CRISPR-Cas9 screen using LNCaP prostate cancer cells, loss of co-repressor TLE3 conferred resistance to AR antagonists apalutamide and enzalutamide. Genes differentially expressed upon TLE3 loss share AR as the top transcriptional regulator, and TLE3 loss rescued the expression of a subset of androgen-responsive genes upon enzalutamide treatment. GR expression was strongly upregulated upon AR inhibition in a TLE3-negative background. This was consistent with binding of TLE3 and AR at the GR locus. Furthermore, GR binding was observed proximal to TLE3/AR-shared genes. GR inhibition resensitized TLE3KO cells to enzalutamide. Analyses of patient samples revealed an association between TLE3 and GR levels that reflected our findings in LNCaP cells, of which the clinical relevance is yet to be determined. Together, our findings reveal a mechanistic link between TLE3 and GR-mediated resistance to AR inhibitors in human prostate cancer.

CHD1 Loss Alters AR Binding at Lineage-Specific Enhancers and Modulates Distinct Transcriptional Programs to Drive Prostate Tumorigenesis.

Deletion of the gene encoding the chromatin remodeler CHD1 is among the most common alterations in prostate cancer (PCa); however, the tumor-suppressive functions of CHD1 and reasons for its tissue-specific loss remain undefined. We demonstrated that CHD1 occupied prostate-specific enhancers enriched for the androgen receptor (AR) and lineage-specific cofactors. Upon CHD1 loss, the AR cistrome was redistributed in patterns consistent with the oncogenic AR cistrome in PCa samples and drove tumor formation in the murine prostate. Notably, this cistrome shift was associated with a unique AR transcriptional signature enriched for pro-oncogenic pathways unique to this tumor subclass. Collectively, these data credential CHD1 as a tumor suppressor in the prostate that constrains AR binding/function to limit tumor progression.

Androgen receptor enhancer usage and the chromatin regulatory landscape in human prostate cancers.

The androgen receptor (AR) is commonly known as a key transcription factor in prostate cancer development, progression and therapy resistance. Genome-wide chromatin association studies revealed that transcriptional regulation by AR mainly depends on binding to distal regulatory enhancer elements that control gene expression through chromatin looping to gene promoters. Changes in the chromatin epigenetic landscape and DNA sequence can locally alter AR-DNA-binding capacity and consequently impact transcriptional output and disease outcome. The vast majority of reports describing AR chromatin interactions have been limited to cell lines, identifying numerous other factors and interacting transcription factors that impact AR chromatin interactions. Do these factors also impact AR cistromics - the genome-wide chromatin-binding landscape of AR - in vivo? Recent technological advances now enable researchers to identify AR chromatin-binding sites and their target genes in human specimens. In this review, we provide an overview of the different factors that influence AR chromatin binding in prostate cancer specimens, which is complemented with knowledge from cell line studies. Finally, we discuss novel perspectives on studying AR cistromics in clinical samples.

Optimized ChIP-seq method facilitates transcription factor profiling in human tumors.

Chromatin immunoprecipitation (ChIP)-seq analyses of transcription factors in clinical specimens are challenging due to the technical limitations and low quantities of starting material, often resulting in low enrichments and poor signal-to-noise ratio. Here, we present an optimized protocol for transcription factor ChIP-seq analyses in human tissue, yielding an ∼100% success rate for all transcription factors analyzed. As proof of concept and to illustrate general applicability of the approach, human tissue from the breast, prostate, and endometrial cancers were analyzed. In addition to standard formaldehyde fixation, disuccinimidyl glutarate was included in the procedure, greatly increasing data quality. To illustrate the sensitivity of the optimized protocol, we provide high-quality ChIP-seq data for three independent factors (AR, FOXA1, and H3K27ac) from a single core needle prostate cancer biopsy specimen. In summary, double-cross-linking strongly improved transcription factor ChIP-seq quality on human tumor samples, further facilitating and enhancing translational research on limited amounts of tissue.

Integrative epigenetic taxonomy of primary prostate cancer.

The Androgen Receptor (AR) is the key-driving transcription factor in prostate cancer, tightly controlled by epigenetic regulation. To date, most epigenetic profiling has been performed in cell lines or limited tissue samples. Here, to comprehensively study the epigenetic landscape, we perform RNA-seq with ChIP-seq for AR and histone modification marks (H3K27ac, H3K4me3, H3K27me3) in 100 primary prostate carcinomas. Integrative molecular subtyping of the five data streams revealed three major subtypes of which two were clearly TMPRSS2-ERG dictated. Importantly, we identify a third subtype with low chromatin binding and activity of AR, but with high activity of FGF and WNT signaling. While positive for neuroendocrine-hallmark genes, these tumors were copy number-neutral with low mutational burden, significantly depleted for genes characteristic of poor-outcome associated luminal B-subtype. We present a unique resource on transcriptional and epigenetic control in prostate cancer, revealing tight control of gene regulation differentially dictated by AR over three subtypes.

Loss of androgen receptor signaling in prostate cancer-associated fibroblasts (CAFs) promotes CCL2- and CXCL8-mediated cancer cell migration.

Fibroblasts are abundantly present in the prostate tumor microenvironment (TME), including cancer-associated fibroblasts (CAFs) which play a key role in cancer development. Androgen receptor (AR) signaling is the main driver of prostate cancer (PCa) progression, and stromal cells in the TME also express AR. High-grade tumor and poor clinical outcome are associated with low AR expression in the TME, which suggests a protective role of AR signaling in the stroma against PCa development. However, the mechanism of this relation is not clear. In this study, we isolated AR-expressing CAF-like cells. Testosterone (R1881) exposure did not affect CAF-like cell morphology, proliferation, or motility. PCa cell growth was not affected by culturing in medium from R1881-exposed CAF-like cells; however, migration of PCa cells was inhibited. AR chromatin immune precipitation sequencing (ChIP-seq) was performed and motif search suggested that AR in CAF-like cells bound the chromatin through AP-1-elements upon R1881 exposure, inducing enhancer-mediated AR chromatin interactions. The vast majority of chromatin binding sites in CAF-like cells were unique and not shared with AR sites observed in PCa cell lines or tumors. AR signaling in CAF-like cells decreased expression of multiple cytokines; most notably CCL2 and CXCL8 and both cytokines increased migration of PCa cells. These results suggest direct paracrine regulation of PCa cell migration by CAFs through AR signaling.

mTOR pathway activation is a favorable prognostic factor in human prostate adenocarcinoma.

Prostate cancer patients with localized disease are treated with curative intent. However, the disease will recur in approximately 30% of patients with a high incidence of morbidity and mortality. Prognostic biomarkers are needed to identify patients with high risk of relapse. mTOR pathway activation is reported in prostate cancer, but clinical trials testing efficacy of mTOR inhibitors were unsuccessful. To explain this clinical observation, we studied the expression and prognostic impact of mTOR-S2448 phosphorylation in localized prostate carcinomas. mTOR-S2448 phosphorylation is indicative for an activated mTOR pathway in prostate cancer. Surprisingly, the mTOR signaling pathway is activated specifically in prostate cancer patients with a favorable outcome. Since tumors from poor-outcome patients have low levels of mTOR-S2448 phosphorylation, this may explain why mTOR inhibitors proved unsuccessful in prostate cancer trials.

Androgen receptor DNA binding and chromatin accessibility profiling in prostate cancer.

Prostate cancer (PCa) is the second most common cancer in men. The Androgen Receptor (AR) is the major driver of PCa and the main target of therapy in the advanced setting. AR is a nuclear receptor that binds the chromatin and regulates transcription of genes involved in cancer cell proliferation and survival. In a study by Stelloo et al. (1) we explored prostate cancer on the level of transcriptional regulation by means of Formaldehyde-Assisted Isolation of Regulatory Elements and Chromatin Immunoprecipitation coupled with massive parallel sequencing (FAIRE-seq and ChIP-seq, respectively). We employed these data for the assessment of differences in transcriptional regulation at distinct stages of PCa progression and to construct a prognostic gene expression classifier. Genomics data includes FAIRE-seq data from normal prostate tissue as well as primary, hormone therapy resistant and metastatic PCa. Furthermore, ChIP-seq data from primary and resistant PCa were generated, along with multiple input controls. The data are publicly available through NCBI GEO database with accession number GSE65478. Here we describe the genomics and clinical data in detail and provide comparative analysis of FAIRE-seq and ChIP-seq data.