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

Deciphering lineage specification during early embryogenesis in mouse gastruloids using multilayered proteomics

Gastrulation is a critical stage in embryonic development during which the germ layers are established. Advances in sequencing technologies led to the identification of gene regulatory programs that control the emergence of the germ layers and their derivatives. However, proteome-based studies of early mammalian development are scarce. To overcome this, we utilized gastruloids and a multilayered mass spectrometry-based proteomics approach to investigate the global dynamics of (phospho) protein expression during gastruloid differentiation. Our findings revealed many proteins with temporal expression and unique expression profiles for each germ layer, which we also validated using single-cell proteomics technology. Additionally, we profiled enhancer interaction landscapes using P300 proximity labeling, which revealed numerous gastruloid-specific transcription factors and chromatin remodelers. Subsequent degron-based perturbations combined with single-cell RNA sequencing (scRNA-seq) identified a critical role for ZEB2 in mouse and human somitogenesis. Overall, this study provides a rich resource for developmental and synthetic biology communities endeavoring to understand mammalian embryogenesis.

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.

Rescuing biologically relevant consensus regions across replicated samples

BACKGROUND

Protein-DNA binding sites of ChIP-seq experiments are identified where the binding affinity is significant based on a given threshold. The choice of the threshold is a trade-off between conservative region identification and discarding weak, but true binding sites.

RESULTS

We rescue weak binding sites using MSPC, which efficiently exploits replicates to lower the threshold required to identify a site while keeping a low false-positive rate, and we compare it to IDR, a widely used post-processing method for identifying highly reproducible peaks across replicates. We observe several master transcription regulators (e.g., SP1 and GATA3) and HDAC2-GATA1 regulatory networks on rescued regions in K562 cell line.

CONCLUSIONS

We argue the biological relevance of weak binding sites and the information they add when rescued by MSPC. An implementation of the proposed extended MSPC methodology and the scripts to reproduce the performed analysis are freely available at https://genometric.github.io/MSPC/ ; MSPC is distributed as a command-line application and an R package available from Bioconductor ( https://doi.org/doi:10.18129/B9.bioc.rmspc ).

Low Gut Microbial Diversity Augments Estrogen-Driven Pulmonary Fibrosis in Female-Predominant Interstitial Lung Disease

Although profibrotic cytokines, such as IL-17A and TGF-β1, have been implicated in the pathogenesis of interstitial lung disease (ILD), the interactions between gut dysbiosis, gonadotrophic hormones and molecular mediators of profibrotic cytokine expression, such as the phosphorylation of STAT3, have not been defined. Here, through chromatin immunoprecipitation sequencing (ChIP-seq) analysis of primary human CD4+ T cells, we show that regions within the STAT3 locus are significantly enriched for binding by the transcription factor estrogen receptor alpha (ERa). Using the murine model of bleomycin-induced pulmonary fibrosis, we found significantly increased regulatory T cells compared to Th17 cells in the female lung. The genetic absence of ESR1 or ovariectomy in mice significantly increased pSTAT3 and IL-17A expression in pulmonary CD4+ T cells, which was reduced after the repletion of female hormones. Remarkably, there was no significant reduction in lung fibrosis under either condition, suggesting that factors outside of ovarian hormones also contribute. An assessment of lung fibrosis among menstruating females in different rearing environments revealed that environments favoring gut dysbiosis augment fibrosis. Furthermore, hormone repletion following ovariectomy further augmented lung fibrosis, suggesting pathologic interactions between gonadal hormones and gut microbiota in relation to lung fibrosis severity. An analysis of female sarcoidosis patients revealed a significant reduction in pSTAT3 and IL-17A levels and a concomitant increase in TGF-β1 levels in CD4+ T cells compared to male sarcoidosis patients. These studies reveal that estrogen is profibrotic in females and that gut dysbiosis in menstruating females augments lung fibrosis severity, supporting a critical interaction between gonadal hormones and gut flora in lung fibrosis pathogenesis.

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.

Low Gut Microbial Diversity Augments Estrogen-driven Pulmonary Fibrosis in Female-Predominant Interstitial Lung Disease

Although profibrotic cytokines such as IL-17A and TGF-β1 have been implicated in interstitial lung disease (ILD) pathogenesis, interactions between gut dysbiosis, gonadotrophic hormones and molecular mediators of profibrotic cytokine expression, such as phosphorylation of STAT3, have not been defined. Here we show by chromatin immunoprecipitation sequencing (ChIP-seq) analysis of primary human CD4+ T cells that regions within the STAT3 locus are significantly enriched for binding by the transcription factor estrogen receptor alpha (ERa). Using the murine model of bleomycin-induced pulmonary fibrosis, we found significantly increased regulatory T cells compared to Th17 cells in the female lung. Genetic absence of ESR1 or ovariectomy in mice significantly increased pSTAT3 and IL-17A expression in pulmonary CD4+ T cells, which was reduced after repletion of female hormones. Remarkably, there was no significant reduction in lung fibrosis under either condition, suggesting that factors outside of ovarian hormones also contribute. Assessment of lung fibrosis among menstruating females in different rearing environments revealed that environments favoring gut dysbiosis augment fibrosis. Furthermore, hormone repletion following ovariectomy further augmented lung fibrosis, suggesting pathologic interactions between gonadal hormones and gut microbiota on lung fibrosis severity. Analysis in female sarcoidosis patients revealed a significant reduction in pSTAT3 and IL-17A levels and a concomitant increase in TGF-β1 levels in CD4+ T cells, compared to male sarcoidosis patients. These studies reveal that estrogen is profibrotic in females and that gut dysbiosis in menstruating females augments lung fibrosis severity, supporting a critical interaction between gonadal hormones and gut flora in lung fibrosis pathogenesis.

Structural variants drive context-dependent oncogene activation in cancer

Higher-order chromatin structure is important for the regulation of genes by distal regulatory sequences^1,2. Structural variants (SVs) that alter three-dimensional (3D) genome organization can lead to enhancer-promoter rewiring and human disease, particularly in the context of cancer³. However, only a small minority of SVs are associated with altered gene expression^4,5, and it remains unclear why certain SVs lead to changes in distal gene expression and others do not. To address these questions, we used a combination of genomic profiling and genome engineering to identify sites of recurrent changes in 3D genome structure in cancer and determine the effects of specific rearrangements on oncogene activation. By analysing Hi-C data from 92 cancer cell lines and patient samples, we identified loci affected by recurrent alterations to 3D genome structure, including oncogenes such as MYC, TERT and CCND1. By using CRISPR-Cas9 genome engineering to generate de novo SVs, we show that oncogene activity can be predicted by using 'activity-by-contact' models that consider partner region chromatin contacts and enhancer activity. However, activity-by-contact models are only predictive of specific subsets of genes in the genome, suggesting that different classes of genes engage in distinct modes of regulation by distal regulatory elements. These results indicate that SVs that alter 3D genome organization are widespread in cancer genomes and begin to illustrate predictive rules for the consequences of SVs on oncogene activation.

MYC promotes immune-suppression in triple-negative breast cancer via inhibition of interferon signaling

The limited efficacy of immune checkpoint inhibitor treatment in triple-negative breast cancer (TNBC) patients is attributed to sparse or unresponsive tumor-infiltrating lymphocytes, but the mechanisms that lead to a therapy resistant tumor immune microenvironment are incompletely known. Here we show a strong correlation between MYC expression and loss of immune signatures in human TNBC. In mouse models of TNBC proficient or deficient of breast cancer type 1 susceptibility gene (BRCA1), MYC overexpression dramatically decreases lymphocyte infiltration in tumors, along with immune signature remodelling. MYC-mediated suppression of inflammatory signalling induced by BRCA1/2 inactivation is confirmed in human TNBC cell lines. Moreover, MYC overexpression prevents the recruitment and activation of lymphocytes in both human and mouse TNBC co-culture models. Chromatin-immunoprecipitation-sequencing reveals that MYC, together with its co-repressor MIZ1, directly binds promoters of multiple interferon-signalling genes, resulting in their downregulation. MYC overexpression thus counters tumor growth inhibition by a Stimulator of Interferon Genes (STING) agonist via suppressing induction of interferon signalling. Together, our data reveal that MYC suppresses innate immunity and facilitates tumor immune escape, explaining the poor immunogenicity of MYC-overexpressing TNBCs.

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.

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.

UALCAN: An update to the integrated cancer data analysis platform

Cancer genomic, transcriptomic, and proteomic profiling has generated extensive data that necessitate the development of tools for its analysis and dissemination. We developed UALCAN to provide a portal for easy exploring, analyzing, and visualizing these data, allowing users to integrate the data to better understand the gene, proteins, and pathways perturbed in cancer and make discoveries. UALCAN web portal enables analyzing and delivering cancer transcriptome, proteomics, and patient survival data to the cancer research community. With data obtained from The Cancer Genome Atlas (TCGA) project, UALCAN has enabled users to evaluate protein-coding gene expression and its impact on patient survival across 33 types of cancers. The web portal has been used extensively since its release and received immense popularity, underlined by its usage from cancer researchers in more than 100 countries. The present manuscript highlights the task we have undertaken and updates that we have made to UALCAN since its release in 2017. Extensive user feedback motivated us to expand the resource by including data on a) microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and promoter DNA methylation from TCGA and b) mass spectrometry-based proteomics from the Clinical Proteomic Tumor Analysis Consortium (CPTAC). UALCAN provides easy access to pre-computed, tumor subgroup-based gene/protein expression, promoter DNA methylation status, and Kaplan-Meier survival analyses. It also provides new visualization features to comprehend and integrate observations and aids in generating hypotheses for testing. UALCAN is accessible at http://ualcan.path.uab.edu

Chromatin Immunoprecipitation Sequencing (ChIP-seq) Protocol for Small Amounts of Frozen Biobanked Cardiac Tissue

Chromatin immunoprecipitation and sequencing (ChIP-seq) is a well-established method to study the epigenetic profile at the genome-wide scale, including histone modifications and DNA-protein interactions. It provides valuable insights to better understand disease mechanisms. Here we present an optimized ChIP-seq protocol suitable for human cardiac tissues, especially the frozen biobanked small biopsy samples.

Opposing transcriptional programs of KLF5 and AR emerge during therapy for advanced prostate cancer

Endocrine therapies for prostate cancer inhibit the androgen receptor (AR) transcription factor. In most cases, AR activity resumes during therapy and drives progression to castration-resistant prostate cancer (CRPC). However, therapy can also promote lineage plasticity and select for AR-independent phenotypes that are uniformly lethal. Here, we demonstrate the stem cell transcription factor Krüppel-like factor 5 (KLF5) is low or absent in prostate cancers prior to endocrine therapy, but induced in a subset of CRPC, including CRPC displaying lineage plasticity. KLF5 and AR physically interact on chromatin and drive opposing transcriptional programs, with KLF5 promoting cellular migration, anchorage-independent growth, and basal epithelial cell phenotypes. We identify ERBB2 as a point of transcriptional convergence displaying activation by KLF5 and repression by AR. ERBB2 inhibitors preferentially block KLF5-driven oncogenic phenotypes. These findings implicate KLF5 as an oncogene that can be upregulated in CRPC to oppose AR activities and promote lineage plasticity.

In vivo CRISPR screens identify the E3 ligase Cop1 as a modulator of macrophage infiltration and cancer immunotherapy target

Despite remarkable clinical efficacy of immune checkpoint blockade (ICB) in cancer treatment, ICB benefits for triple-negative breast cancer (TNBC) remain limited. Through pooled in vivo CRISPR knockout (KO) screens in syngeneic TNBC mouse models, we found that deletion of the E3 ubiquitin ligase Cop1 in cancer cells decreases secretion of macrophage-associated chemokines, reduces tumor macrophage infiltration, enhances anti-tumor immunity, and strengthens ICB response. Transcriptomics, epigenomics, and proteomics analyses revealed that Cop1 functions through proteasomal degradation of the C/ebpδ protein. The Cop1 substrate Trib2 functions as a scaffold linking Cop1 and C/ebpδ, which leads to polyubiquitination of C/ebpδ. In addition, deletion of the E3 ubiquitin ligase Cop1 in cancer cells stabilizes C/ebpδ to suppress expression of macrophage chemoattractant genes. Our integrated approach implicates Cop1 as a target for improving cancer immunotherapy efficacy in TNBC by regulating chemokine secretion and macrophage infiltration in the tumor microenvironment.

Reshaping of the androgen-driven chromatin landscape in normal prostate cells by early cancer drivers and effect on therapeutic sensitivity

The normal androgen receptor (AR) cistrome and transcriptional program are fundamentally altered in prostate cancer (PCa). Here, we profile the chromatin landscape and AR-directed transcriptional program in normal prostate cells and show the impact of SPOP mutations, an early event in prostate tumorigenesis. In genetically normal mouse prostate organoids, SPOP mutation results in accessibility and AR binding patterns similar to that of human PCa. Consistent with dependence on AR signaling, castration of SPOP mutant mouse models results in the loss of neoplastic phenotypes, and human SPOP mutant PCa shows a favorable response to AR-targeted therapies. Together, these data validate mouse prostate organoids as a robust model for studying epigenomic and transcriptional alterations in normal prostate, provide valuable datasets for further studies, and show that a single genomic alteration may be sufficient to reprogram the chromatin of normal prostate cells toward oncogenic phenotypes, with potential therapeutic implications for AR-targeting therapies.

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.

Purification of DNA repair protein complexes from mammalian cells

Cells possess multiple DNA repair pathways to tackle a variety of DNA lesions. Often, DNA repair proteins function as large protein complexes. Here, we describe a protocol to purify DNA repair protein complexes from nuclei of mammalian cells. The method permits purification of protein complexes containing stable as well as transiently associated proteins, which subsequently can be identified by mass-spectrometry analysis. This protocol can be applied to uncover the functions and mechanism of DNA repair pathways. For complete information on the use and execution of this protocol, please refer to Socha et al. (2020).

Systematic characterization of mutations altering protein degradation in human cancers

The ubiquitin-proteasome system (UPS) is the primary route for selective protein degradation in human cells. The UPS is an attractive target for novel cancer therapies, but the precise UPS genes and substrates important for cancer growth are incompletely understood. Leveraging multi-omics data across more than 9,000 human tumors and 33 cancer types, we found that over 19% of all cancer driver genes affect UPS function. We implicate transcription factors as important substrates and show that c-Myc stability is modulated by CUL3. Moreover, we developed a deep learning model (deepDegron) to identify mutations that result in degron loss and experimentally validated the prediction that gain-of-function truncating mutations in GATA3 and PPM1D result in increased protein stability. Last, we identified UPS driver genes associated with prognosis and the tumor microenvironment. This study demonstrates the important role of UPS dysregulation in human cancer and underscores the potential therapeutic utility of targeting the UPS.

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.

A comprehensive enhancer screen identifies TRAM2 as a key and novel mediator of YAP oncogenesis

BACKGROUND

Frequent activation of the co-transcriptional factor YAP is observed in a large number of solid tumors. Activated YAP associates with enhancer loci via TEAD4-DNA-binding protein and stimulates cancer aggressiveness. Although thousands of YAP/TEAD4 binding-sites are annotated, their functional importance is unknown. Here, we aim at further identification of enhancer elements that are required for YAP functions.

RESULTS

We first apply genome-wide ChIP profiling of YAP to systematically identify enhancers that are bound by YAP/TEAD4. Next, we implement a genetic approach to uncover functions of YAP/TEAD4-associated enhancers, demonstrate its robustness, and use it to reveal a network of enhancers required for YAP-mediated proliferation. We focus on EnhancerTRAM2, as its target gene TRAM2 shows the strongest expression-correlation with YAP activity in nearly all tumor types. Interestingly, TRAM2 phenocopies the YAP-induced cell proliferation, migration, and invasion phenotypes and correlates with poor patient survival. Mechanistically, we identify FSTL-1 as a major direct client of TRAM2 that is involved in these phenotypes. Thus, TRAM2 is a key novel mediator of YAP-induced oncogenic proliferation and cellular invasiveness.

CONCLUSIONS

YAP is a transcription co-factor that binds to thousands of enhancer loci and stimulates tumor aggressiveness. Using unbiased functional approaches, we dissect YAP enhancer network and characterize TRAM2 as a novel mediator of cellular proliferation, migration, and invasion. Our findings elucidate how YAP induces cancer aggressiveness and may assist diagnosis of cancer metastasis.

MethylationToActivity: a deep-learning framework that reveals promoter activity landscapes from DNA methylomes in individual tumors

Although genome-wide DNA methylomes have demonstrated their clinical value as reliable biomarkers for tumor detection, subtyping, and classification, their direct biological impacts at the individual gene level remain elusive. Here we present MethylationToActivity (M2A), a machine learning framework that uses convolutional neural networks to infer promoter activities based on H3K4me3 and H3K27ac enrichment, from DNA methylation patterns for individual genes. Using publicly available datasets in real-world test scenarios, we demonstrate that M2A is highly accurate and robust in revealing promoter activity landscapes in various pediatric and adult cancers, including both solid and hematologic malignant neoplasms.

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.

Resolving the 3D Landscape of Transcription-Linked Mammalian Chromatin Folding

Whereas folding of genomes at the large scale of epigenomic compartments and topologically associating domains (TADs) is now relatively well understood, how chromatin is folded at finer scales remains largely unexplored in mammals. Here, we overcome some limitations of conventional 3C-based methods by using high-resolution Micro-C to probe links between 3D genome organization and transcriptional regulation in mouse stem cells. Combinatorial binding of transcription factors, cofactors, and chromatin modifiers spatially segregates TAD regions into various finer-scale structures with distinct regulatory features including stripes, dots, and domains linking promoters-to-promoters (P-P) or enhancers-to-promoters (E-P) and bundle contacts between Polycomb regions. E-P stripes extending from the edge of domains predominantly link co-expressed loci, often in the absence of CTCF and cohesin occupancy. Acute inhibition of transcription disrupts these gene-related folding features without altering higher-order chromatin structures. Our study uncovers previously obscured finer-scale genome organization, establishing functional links between chromatin folding and gene regulation.

TRPS1 acts as a context-dependent regulator of mammary epithelial cell growth/differentiation and breast cancer development

In this study, Cornelissen et al. set out to elucidate the role of the GATA-type zinc finger transcription factor TRPS1 in breast cancer. Using in vitro and in vivo loss-of-function approaches, the authors demonstrate that TRPS1 can function as a context-dependent tumor suppressor in breast cancer, while being essential for growth and differentiation of normal mammary epithelial cells.

The GATA-type zinc finger transcription factor TRPS1 has been implicated in breast cancer. However, its precise role remains unclear, as both amplifications and inactivating mutations in TRPS1 have been reported. Here, we used in vitro and in vivo loss-of-function approaches to dissect the role of TRPS1 in mammary gland development and invasive lobular breast carcinoma, which is hallmarked by functional loss of E-cadherin. We show that TRPS1 is essential in mammary epithelial cells, since TRPS1-mediated suppression of interferon signaling promotes in vitro proliferation and lactogenic differentiation. Similarly, TRPS1 expression is indispensable for proliferation of mammary organoids and in vivo survival of luminal epithelial cells during mammary gland development. However, the consequences of TRPS1 loss are dependent on E-cadherin status, as combined inactivation of E-cadherin and TRPS1 causes persistent proliferation of mammary organoids and accelerated mammary tumor formation in mice. Together, our results demonstrate that TRPS1 can function as a context-dependent tumor suppressor in breast cancer, while being essential for growth and differentiation of normal mammary epithelial cells.

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 TLE3^KO 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.

Eomes and Brachyury control pluripotency exit and germ-layer segregation by changing the chromatin state

The first lineage specification of pluripotent mouse epiblast segregates neuroectoderm (NE) from mesoderm and definitive endoderm (ME) by mechanisms that are not well understood. Here we demonstrate that the induction of ME gene programs critically relies on the T-box transcription factors Eomesodermin (also known as Eomes) and Brachyury, which concomitantly repress pluripotency and NE gene programs. Cells deficient in these T-box transcription factors retain pluripotency and differentiate to NE lineages despite the presence of ME-inducing signals transforming growth factor β (TGF-β)/Nodal and Wnt. Pluripotency and NE gene networks are additionally repressed by ME factors downstream of T-box factor induction, demonstrating a redundancy in program regulation to safeguard mutually exclusive lineage specification. Analyses of chromatin revealed that accessibility of ME enhancers depends on T-box factor binding, whereas NE enhancers are accessible and already activation primed at pluripotency. This asymmetry of the chromatin landscape thus explains the default differentiation of pluripotent cells to NE in the absence of ME induction that depends on activating and repressive functions of Eomes and Brachyury.

Exogenous ERα Expression in the Mammary Epithelium Decreases Over Time and Does Not Contribute to p53-Deficient Mammary Tumor Formation in Mice

Approximately 75% of all breast cancers express the nuclear hormone receptor estrogen receptor α (ERα). However, the majority of mammary tumors from genetically engineered mouse models (GEMMs) are ERα-negative. To model ERα-positive breast cancer in mice, we exogenously introduced expression of mouse and human ERα in an existing GEMM of p53-deficient breast cancer. After initial ERα expression during mammary gland development, expression was reduced or lost in adult glands and p53-deficient mammary tumors. Chromatin immunoprecipitation (ChIP)-sequencing analysis of primary mouse mammary epithelial cells (MMECs) derived from these models, in which expression of the ERα constructs was induced in vitro, confirmed interaction of ERα with the DNA. In human breast and endometrial cancer, and also in healthy breast tissue, DNA binding of ERα is facilitated by the pioneer factor FOXA1. Surprisingly, the ERα binding sites identified in primary MMECs, but also in mouse mammary gland and uterus, showed an high enrichment of ERE motifs, but were devoid of Forkhead motifs. Furthermore, exogenous introduction of FOXA1 and GATA3 in ERα-expressing MMECs was not sufficient to promote ERα-responsiveness of these cells. Together, this suggests that species-specific differences in pioneer factor usage between mouse and human are dictated by the DNA sequence, resulting in ERα-dependencies in mice that are not FOXA1 driven. These species-specific differences in ERα-biology may limit the utility of mice for in vivo modeling of ERα-positive breast cancer.

Principles and methods of integrative chromatin analysis in primary tissues and tumors

Recent methodological advances have enabled the genome-wide interrogation of chromatin from primary tumor tissues. Integrative analysis of histone post-translational modifications, transcription factor (TF) binding and open chromatin sites in tumors across cancer stages can elucidate the aberrant epigenetic states accompanying tumor progression. Cancer-associated chromatin alterations can activate or inactivate enhancers at genes involved in cancer while still respecting cell-of-origin constrictions. Accordingly, enhancer analysis in cancer could have uses for biomarker discovery to further refine patient diagnosis and potentially sub-classify patients for tailored therapy. Methodologies used for chromatin analyses of primary tissues need to address issues distinct from cell line studies including the specific sources of variability coming from the heterogeneous cellular composition of tissues and from inter-individual (epi)genetic differences. This leads to requirements for careful histological analysis to select the specific samples and cells of interest. In analyzing tumors somatic changes should be taken into account to distinguish the genuine epigenetic changes across tumor specimens from any genetic alterations such as copy number variations (CNV). In this contribution we review a selection of current results from chromatin profiling, examine experimental methodologies and discuss specific analysis approaches. We also review specific considerations regarding tissue preparation for epigenetic analysis and conclude with our perspectives on emerging approaches that will impact studies of chromatin landscapes of clinical samples in the future.

The genomic landscape of metastatic castration-resistant prostate cancers reveals multiple distinct genotypes with potential clinical impact

Metastatic castration-resistant prostate cancer (mCRPC) has a highly complex genomic landscape. With the recent development of novel treatments, accurate stratification strategies are needed. Here we present the whole-genome sequencing (WGS) analysis of fresh-frozen metastatic biopsies from 197 mCRPC patients. Using unsupervised clustering based on genomic features, we define eight distinct genomic clusters. We observe potentially clinically relevant genotypes, including microsatellite instability (MSI), homologous recombination deficiency (HRD) enriched with genomic deletions and BRCA2 aberrations, a tandem duplication genotype associated with CDK12-/- and a chromothripsis-enriched subgroup. Our data suggests that stratification on WGS characteristics may improve identification of MSI, CDK12-/- and HRD patients. From WGS and ChIP-seq data, we show the potential relevance of recurrent alterations in non-coding regions identified with WGS and highlight the central role of AR signaling in tumor progression. These data underline the potential value of using WGS to accurately stratify mCRPC patients into clinically actionable subgroups.

Identification of mineralocorticoid receptor target genes in the mouse hippocampus

Brain mineralocorticoid receptors (MRs) and glucocorticoid receptors (GRs) respond to the same glucocorticoid hormones but can have differential effects on cellular function. Several lines of evidence suggest that MR-specific target genes must exist and might underlie the distinct effects of the receptors. The present study aimed to identify MR-specific target genes in the hippocampus, a brain region where MR and GR are co-localised and play a role in the stress response. Using genome-wide binding of both receptor types, we previously identified MR-specific, MR-GR overlapping and GR-specific putative target genes. We now report altered gene expression levels of such genes in the hippocampus of forebrain MR knockout (fbMRKO) mice, killed at the time of their endogenous corticosterone peak. Of those genes associated with MR-specific binding, the most robust effect was a 50% reduction in Jun dimerization protein 2 (Jdp2) mRNA levels in fbMRKO mice. Down-regulation was also observed for the MR-specific Nitric oxide synthase 1 adaptor protein (Nos1ap) and Suv3 like RNA helicase (Supv3 l1). Interestingly, the classical glucocorticoid target gene FK506 binding protein 5 (Fkbp5), which is associated with MR and GR chromatin binding, was expressed at substantially lower levels in fbMRKO mice. Subsequently, hippocampal Jdp2 was confirmed to be up-regulated in a restraint stress model, posing Jdp2 as a bona fide MR target that is also responsive in an acute stress condition. Thus, we show that MR-selective DNA binding can reveal functional regulation of genes and further identify distinct MR-specific effector pathways.

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.