Cellular reprogramming by the conjoint action of ERα, FOXA1, and GATA3 to a ligand-inducible growth state

Selective Estrogen Receptor Modulators' (SERMs) Influence on TET3 Expression in Breast Cancer Cell Lines with Distinct Biological Subtypes

Tamoxifen, a selective estrogen receptor modulator (SERM), exhibits dual agonist or antagonist effects contingent upon its binding to either G-protein-coupled estrogen receptor (GPER) or estrogen nuclear receptor (ESR). Estrogen signaling plays a pivotal role in initiating epigenetic alterations and regulating estrogen-responsive genes in breast cancer. Employing three distinct breast cancer cell lines-MCF-7 (ESR+; GPER+), MDA-MB-231 (ESR-; GPER-), and SkBr3 (ESR-; GPER+)-this study subjected them to treatment with two tamoxifen derivatives: 4-hydroxytamoxifen (4-HT) and endoxifen (Endox). Through 2D high-performance liquid chromatography with tandem mass spectrometry detection (HPLC-MS/MS), varying levels of 5-methylcytosine (5-mC) were found, with MCF-7 displaying the highest levels. Furthermore, TET3 mRNA expression levels varied among the cell lines, with MCF-7 exhibiting the lowest expression. Notably, treatment with 4-HT induced significant changes in TET3 expression across all cell lines, with the most pronounced increase seen in MCF-7 and the least in MDA-MB-231. These findings underscore the influence of tamoxifen derivatives on DNA methylation patterns, particularly through modulating TET3 expression, which appears to be contingent on the presence of estrogen receptors. This study highlights the potential of targeting epigenetic modifications for personalized anti-cancer therapy, offering a novel avenue to improve treatment outcomes.

The androgen receptor interacts with GATA3 to transcriptionally regulate a luminal epithelial cell phenotype in breast cancer

BACKGROUND

The androgen receptor (AR) is a tumor suppressor in estrogen receptor (ER) positive breast cancer, a role sustained in some ER negative breast cancers. Key factors dictating AR genomic activity in a breast context are largely unknown. Herein, we employ an unbiased chromatin immunoprecipitation-based proteomic technique to identify endogenous AR interacting co-regulatory proteins in ER positive and negative models of breast cancer to gain new insight into mechanisms of AR signaling in this disease.

RESULTS

The DNA-binding factor GATA3 is identified and validated as a novel AR interacting protein in breast cancer cells irrespective of ER status. AR activation by the natural ligand 5α-dihydrotestosterone (DHT) increases nuclear AR-GATA3 interactions, resulting in AR-dependent enrichment of GATA3 chromatin binding at a sub-set of genomic loci. Silencing GATA3 reduces but does not prevent AR DNA binding and transactivation of genes associated with AR/GATA3 co-occupied loci, indicating a co-regulatory role for GATA3 in AR signaling. DHT-induced AR/GATA3 binding coincides with upregulation of luminal differentiation genes, including EHF and KDM4B, established master regulators of a breast epithelial cell lineage. These findings are validated in a patient-derived xenograft model of breast cancer. Interaction between AR and GATA3 is also associated with AR-mediated growth inhibition in ER positive and ER negative breast cancer.

CONCLUSIONS

AR and GATA3 interact to transcriptionally regulate luminal epithelial cell differentiation in breast cancer regardless of ER status. This interaction facilitates the tumor suppressor function of AR and mechanistically explains why AR expression is associated with less proliferative, more differentiated breast tumors and better overall survival in breast cancer.

Enhancer-driven 3D chromatin domain folding modulates transcription in human mammary tumor cells

The genome is organized in functional compartments and structural domains at the sub-megabase scale. How within these domains interactions between numerous cis-acting enhancers and promoters regulate transcription remains an open question. Here, we determined chromatin folding and composition over several hundred kb around estrogen-responsive genes in human breast cancer cell lines after hormone stimulation. Modeling of 5C data at 1.8 kb resolution was combined with quantitative 3D analysis of multicolor FISH measurements at 100 nm resolution and integrated with ChIP-seq data on transcription factor binding and histone modifications. We found that rapid estradiol induction of the progesterone gene expression occurs in the context of preexisting, cell type-specific chromosomal architectures encompassing the 90 kb progesterone gene coding region and an enhancer-spiked 5' 300 kb upstream genomic region. In response to estradiol, interactions between estrogen receptor α (ERα) bound regulatory elements are reinforced. Whereas initial enhancer-gene contacts coincide with RNA Pol 2 binding and transcription initiation, sustained hormone stimulation promotes ERα accumulation creating a regulatory hub stimulating transcript synthesis. In addition to implications for estrogen receptor signaling, we uncover that preestablished chromatin architectures efficiently regulate gene expression upon stimulation without the need for de novo extensive rewiring of long-range chromatin interactions.

Meta-analysis of integrated ChIP-seq and transcriptome data revealed genomic regions affected by estrogen receptor alpha in breast cancer

BACKGROUND

The largest group of patients with breast cancer are estrogen receptor-positive (ER+) type. The estrogen receptor acts as a transcription factor and triggers cell proliferation and differentiation. Hence, investigating ER-DNA interaction genomic regions can help identify genes directly regulated by ER and understand the mechanism of ER action in cancer progression.

METHODS

In the present study, we employed a workflow to do a meta-analysis of ChIP-seq data of ER+ cell lines stimulated with 10 nM and 100 nM of E2. All publicly available data sets were re-analyzed with the same platform. Then, the known and unknown batch effects were removed. Finally, the meta-analysis was performed to obtain meta-differentially bound sites in estrogen-treated MCF7 cell lines compared to vehicles (as control). Also, the meta-analysis results were compared with the results of T47D cell lines for more precision. Enrichment analyses were also employed to find the functional importance of common meta-differentially bound sites and associated genes among both cell lines.

RESULTS

Remarkably, POU5F1B, ZNF662, ZNF442, KIN, ZNF410, and SGSM2 transcription factors were recognized in the meta-analysis but not in individual studies. Enrichment of the meta-differentially bound sites resulted in the candidacy of pathways not previously reported in breast cancer. PCGF2, HNF1B, and ZBED6 transcription factors were also predicted through the enrichment analysis of associated genes. In addition, comparing the meta-analysis results of both ChIP-seq and RNA-seq data showed that many transcription factors affected by ER were up-regulated.

CONCLUSION

The meta-analysis of ChIP-seq data of estrogen-treated MCF7 cell line leads to the identification of new binding sites of ER that have not been previously reported. Also, enrichment of the meta-differentially bound sites and their associated genes revealed new terms and pathways involved in the development of breast cancer which should be examined in future in vitro and in vivo studies.

A Tale of Three Proteomes: Visualizing Protein and Transcript Abundance Relationships in the Breast Cancer Proteome Portal

Molecular characterization is transforming research on novel therapeutics in breast cancer. High-throughput methodologies are unbiased to hypotheses; thus, data produced are relevant to address unlimited questions and provide resources for the experimental design process. However, the opportunity is often overlooked because data are not readily accessed or analyzed. Herein, the Breast Cancer Proteome Portal, the only online tool for analyzing protein and transcript abundances across the three breast cancer proteomics studies, is presented. The tool is applied to demonstrate that cofunctioning protein abundances are highly correlated and, conversely, high abundance correlation may be an indicator of cofunction. Furthermore, the cofunction-correlation relationship is less resolved at the transcript level. By applying analysis and visualization tools within the Breast Cancer Proteome Portal, insights are garnered about serine synthesis and the compartmentalization of one-carbon metabolism in breast cancer, and a transcription factor tumorigenic regulatory network of glutamine deamination and oxidation is proposed, illustrating that the Breast Cancer Proteome Portal provides an interface for garnering insights from the information-rich studies of the breast cancer proteome.

BAP18 facilitates CTCF-mediated chromatin accessible to regulate enhancer activity in breast cancer

The estrogen receptor alpha (ERα) signaling pathway is a crucial target for ERα-positive breast cancer therapeutic strategies. Co-regulators and other transcription factors cooperate for effective ERα-related enhancer activation. Recent studies demonstrate that the transcription factor CTCF is essential to participate in ERα/E2-induced enhancer transactivation. However, the mechanism of how CTCF is achieved remains unknown. Here, we provided evidence that BAP18 is required for CTCF recruitment on ERα-enriched enhancers, facilitating CTCF-mediated chromatin accessibility to promote enhancer RNAs transcription. Consistently, GRO-seq demonstrates that the enhancer activity is positively correlated with BAP18 enrichment. Furthermore, BAP18 interacts with SMARCA1/BPTF to accelerate the recruitment of CTCF to ERα-related enhancers. Interestingly, BAP18 is involved in chromatin accessibility within enhancer regions, thereby increasing enhancer transactivation and enhancer-promoter looping. BAP18 depletion increases the sensitivity of anti-estrogen and anti-enhancer treatment in MCF7 cells. Collectively, our study indicates that BAP18 coordinates with CTCF to enlarge the transactivation of ERα-related enhancers, providing a better understanding of BAP18/CTCF coupling chromatin remodeling and E-P looping in the regulation of enhancer transcription.

CRISPR/Cas9 screen uncovers functional translation of cryptic lncRNA-encoded open reading frames in human cancer

Emerging evidence suggests that cryptic translation within long noncoding RNAs (lncRNAs) may produce novel proteins with important developmental/physiological functions. However, the role of this cryptic translation in complex diseases (e.g., cancer) remains elusive. Here, we applied an integrative strategy combining ribosome profiling and CRISPR/Cas9 screening with large-scale analysis of molecular/clinical data for breast cancer (BC) and identified estrogen receptor α-positive (ER+) BC dependency on the cryptic ORFs encoded by lncRNA genes that were upregulated in luminal tumors. We confirmed the in vivo tumor-promoting function of an unannotated protein, GATA3-interacting cryptic protein (GT3-INCP) encoded by LINC00992, the expression of which was associated with poor prognosis in luminal tumors. GTE-INCP was upregulated by estrogen/ER and regulated estrogen-dependent cell growth. Mechanistically, GT3-INCP interacted with GATA3, a master transcription factor key to mammary gland development/BC cell proliferation, and coregulated a gene expression program that involved many BC susceptibility/risk genes and impacted estrogen response/cell proliferation. GT3-INCP/GATA3 bound to common cis regulatory elements and upregulated the expression of the tumor-promoting and estrogen-regulated BC susceptibility/risk genes MYB and PDZK1. Our study indicates that cryptic lncRNA-encoded proteins can be an important integrated component of the master transcriptional regulatory network driving aberrant transcription in cancer, and suggests that the "hidden" lncRNA-encoded proteome might be a new space for therapeutic target discovery.

Estrogen receptor α K303R mutation reorganizes its binding to forkhead box protein A1 regions and induces chromatin opening

BACKGROUND

Estrogen receptor alpha (ERα) is a frequently mutated gene in breast cancer (BC). While many studies have investigated molecular dysregulation by hotspot mutations at Y537 and D538, which exhibit an estrogen-independent constitutively active phenotype, the functional abnormalities of other mutations remain obscure. The K303R mutation in primary invasive BC has been implicated with endocrine resistance, tumor size, and lymph node positivity. However, the impact of the K303R mutation on the cell epigenome is yet unknown.

METHODS AND RESULTS

We introduced the K303R ERα mutant in ERα-negative MDA-MB-453 cells to monitor ERα-dependent transactivation and to perform epigenomic analyses. ATAC-seq and ChIP-Seq analyses indicated that both wild-type (WT) and the K303R mutant associated with Forkhead box (Fox) protein family motif regions at similar rates, even without an ERα-binding sequence, but only the K303R mutant induced chromatin opening at those regions. Biochemical analyses demonstrated that the WT and the K303R mutant can be tethered on DNA by FoxA1 indirectly, but only the K303R/FoxA1/DNA complex can induce associations with the nuclear receptor cofactor 2 (NCOA2).

CONCLUSIONS

These findings suggest that the K303R mutant induces chromatin opening at the Fox binding region through the FoxA1-dependent associations of the K303R mutant to NCOA2 and then probably disrupts the regulation of Fox-target genes, resulting in K303R-related BC events.

Association between the CASC16 rs4784227 polymorphism and breast cancer risk and prognosis in a northeast Chinese Han population

Background

Breast cancer (BC) poses a serious threat to women worldwide. This research was designed to explore the association between the rs4784227 polymorphism of cancer susceptibility candidate gene 16 (CASC16) and BC susceptibility and prognosis, aiming to provide further information for the early detection of BC and to accelerate comprehensive cancer management.

Methods

A total of 1,733 subjects were recruited for this case-control study, of which 828 are BC patients and 905 are healthy individuals. The relevance between SNP rs4784227 and BC risk in diverse genetic models was analyzed by using the SNPStats analysis program and was assessed by odds ratios (ORs) and 95% confidence intervals (CIs) using the binary logistic regression model. Pearson's χ ² test was used to determine the correlation between the polymorphism and clinical characteristics of BC patients. Additionally, univariate survival analysis was performed by the Kaplan-Meier method and log-rank test, and multivariate survival analysis was performed by Cox regression.

Results

SNP rs4784227 was significantly associated with susceptibility to BC in the dominant model (CT/TT versus CC, OR = 1.237, 95% CI = 1.012-1.513, P = 0.038). The minor allele of SNP rs4784227 was significantly linked to an increased risk of BC (OR = 1.197, 95% CI = 1.022-1.401, P = 0.026). In addition, the rs4784227 polymorphism of CASC16 was associated with perineural invasion (P = 0.030), menstrual status (P = 0.016) and histological grade (P = 0.001, P = 0.003, P = 0.025; respectively) of BC patients. There was no significant association between the genotypes of rs4784227 and disease-free survival (DFS) or overall survival (OS) of breast cancer patients (P > 0.05).

Conclusions

The rs4784227 polymorphism of CASC16 may affect susceptibility to breast cancer and is associated with perineural invasion, menstrual status and histological grade in BC patients. Additionally, our results could not confirm that this polymorphism was related to breast cancer prognosis.

14-3-3τ drives estrogen receptor loss via ERα36 induction and GATA3 inhibition in breast cancer

About one-fourth of recurrent estrogen receptor-positive (ER+) breast cancers lose ER expression, leading to endocrine therapy failure. However, the mechanisms underlying ER loss remain to be fully explored. We now show that 14-3-3τ, up-regulated in ∼60% of breast cancer, drives the conversion of ER+ to ER- and epithelial-to-mesenchymal transition (EMT). We identify ERα36, an isoform of ERα66, as a downstream effector of 14-3-3τ. Overexpression of 14-3-3τ induces ERα36 in xenografts and tumor spheroids. The regulation is further supported by a positive correlation between ERα36 and 14-3-3τ expression in human breast cancers. ERα36 can antagonize ERα66 and inhibit ERα66 expression. Isoform-specific depletion of ERα36 blocks the ER conversion and EMT induced by 14-3-3τ overexpression in tumor spheroids, thus establishing ERα36 as a key mediator in 14-3-3τ-driven ER loss and EMT. ERα36 promoter is repressed by GATA3, which can be phosphorylated by AKT at consensus binding sites for 14-3-3. Upon AKT activation, 14-3-3τ binds phosphorylated GATA3 and facilitates the degradation of GATA3 causing GATA3 to lose transcriptional control over its target genes ERα66 and ERα36. We also demonstrate a role for the collaboration between 14-3-3τ and AKT in ERα36 induction and endocrine therapy resistance by three-dimensional spheroid and tamoxifen treatment models in MCF7 and T47D ER+ breast cancer cells. Thus, the 14-3-3τ-ERα36 regulation provides a previously unrecognized mechanism for ER loss and endocrine therapy failure.

CREBBP/EP300 acetyltransferase inhibition disrupts FOXA1-bound enhancers to inhibit the proliferation of ER+ breast cancer cells

Therapeutic targeting of the estrogen receptor (ER) is a clinically validated approach for estrogen receptor positive breast cancer (ER+ BC), but sustained response is limited by acquired resistance. Targeting the transcriptional coactivators required for estrogen receptor activity represents an alternative approach that is not subject to the same limitations as targeting estrogen receptor itself. In this report we demonstrate that the acetyltransferase activity of coactivator paralogs CREBBP/EP300 represents a promising therapeutic target in ER+ BC. Using the potent and selective inhibitor CPI-1612, we show that CREBBP/EP300 acetyltransferase inhibition potently suppresses in vitro and in vivo growth of breast cancer cell line models and acts in a manner orthogonal to directly targeting ER. CREBBP/EP300 acetyltransferase inhibition suppresses ER-dependent transcription by targeting lineage-specific enhancers defined by the pioneer transcription factor FOXA1. These results validate CREBBP/EP300 acetyltransferase activity as a viable target for clinical development in ER+ breast cancer.

Epigenetic Factors as Etiological Agents, Diagnostic Markers, and Therapeutic Targets for Luminal Breast Cancer

Luminal breast cancer, an etiologically heterogeneous disease, is characterized by high steroid hormone receptor activity and aberrant gene expression profiles. Endocrine therapy and chemotherapy are promising therapeutic approaches to mitigate breast cancer proliferation and recurrence. However, the treatment of therapy-resistant breast cancer is a major challenge. Recent studies on breast cancer etiology have revealed the critical roles of epigenetic factors in luminal breast cancer tumorigenesis and drug resistance. Tumorigenic epigenetic factor-induced aberrant chromatin dynamics dysregulate the onset of gene expression and consequently promote tumorigenesis and metastasis. Epigenetic dysregulation, a type of somatic mutation, is a high-risk factor for breast cancer progression and therapy resistance. Therefore, epigenetic modulators alone or in combination with other therapies are potential therapeutic agents for breast cancer. Several clinical trials have analyzed the therapeutic efficacy of potential epi-drugs for breast cancer and reported beneficial clinical outcomes, including inhibition of tumor cell adhesion and invasiveness and mitigation of endocrine therapy resistance. This review focuses on recent findings on the mechanisms of epigenetic factors in the progression of luminal breast cancer. Additionally, recent findings on the potential of epigenetic factors as diagnostic biomarkers and therapeutic targets for breast cancer are discussed.

GATA3 Truncation Mutants Alter EMT Related Gene Expression via Partial Motif Recognition in Luminal Breast Cancer Cells

GATA3 is known to be one of the most frequently mutated genes in breast cancer. More than 10% of breast tumors carry mutations in this gene. However, the functional consequence of GATA3 mutations is still largely unknown. Clinical data suggest that different types of GATA3 mutations may have distinct roles in breast cancer characterization. In this study, we have established three luminal breast cancer cell lines that stably express different truncation mutants (X308 splice site deletion, C321 frameshift, and A333 frameshift mutants) found in breast cancer patients. Transcriptome analysis identified common and distinct gene expression patterns in these GATA3 mutant cell lines. In particular, the impacts on epithelial-to-mesenchymal transition (EMT) related genes are similar across these mutant cell lines. Chromatin localization of the mutants is highly overlapped and exhibits non-canonical motif enrichment. Interestingly, the A333 frameshift mutant expressed cells displayed the most significant impact on the GATA3 binding compared to X308 splice site deletion and C321fs mutants expressed cells. Our results suggest the common and different roles of GATA3 truncation mutations during luminal breast cancer development.

ESR1 fusions and therapeutic resistance in metastatic breast cancer

Breast cancer is the most frequent female malignant tumor, and the leading cause of cancer death in women worldwide. The most common subtype of breast cancer is hormone receptor positive that expresses the estrogen receptor (ER). Targeting ER with endocrine therapy (ET) is the current standard of care for ER positive (ER+) breast cancer, reducing mortality by up to 40% in early- stage disease. However, resistance to ET represents a major clinical challenge for ER+ breast cancer patients leading to disease recurrence or progression of metastatic disease. Salient drivers of ET resistance are missense mutations in the ER gene (ESR1) leading to constitutive transcriptional activity and reduced ET sensitivity. These mutations are particularly prominent and deleterious in metastatic breast cancer (MBC). In addition to activating ESR1 point mutations, emerging evidence imposes that chromosomal translocation involving the ESR1 gene can also drive ET resistance through the formation of chimeric transcription factors with constitutive transcriptional activity. Although these ESR1 gene fusions are relatively rare, they are enriched in ET resistant metastatic disease. This review discusses the characteristics of ER fusion proteins and their association with clinical outcomes in more aggressive and metastatic breast cancer. The structure and classification of ER fusion proteins based on function and clinical significance are also addressed. Finally, this review summarizes the metastatic phenotypes exhibited by the ER fusion proteins and their role in intrinsic ET resistance.

HiC1Dmetrics: framework to extract various one-dimensional features from chromosome structure data

Eukaryotic genomes are organized in a three-dimensional spatial structure. In this regard, the development of chromosome conformation capture methods has enabled studies of chromosome organization on a genomic scale. Hi-C, the high-throughput chromosome conformation capture method, can reveal a population-averaged, hierarchical chromatin structure. The typical Hi-C analysis uses a two-dimensional (2D) contact matrix that indicates contact frequencies between all possible genomic position pairs. Oftentimes, however, such a 2D matrix is not amenable to handling quantitative comparisons, visualizations and integrations across multiple datasets. Although several one-dimensional (1D) metrics have been proposed to depict structural information in Hi-C data, their effectiveness is still underappreciated. Here, we first review the currently available 1D metrics for individual Hi-C samples or two-sample comparisons and then discuss their validity and suitable analysis scenarios. We also propose several new 1D metrics to identify additional unique features of chromosome structures. We highlight that the 1D metrics are reproducible and robust for comparing and visualizing multiple Hi-C samples. Moreover, we show that 1D metrics can be easily combined with epigenome tracks to annotate chromatin states in greater details. We develop a new framework, called HiC1Dmetrics, to summarize all 1D metrics discussed in this study. HiC1Dmetrics is open-source (github.com/wangjk321/HiC1Dmetrics) and can be accessed from both command-line and web-based interfaces. Our tool constitutes a useful resource for the community of chromosome-organization researchers.

E239K mutation abolishes the suppressive effects of lysine-specific demethylase 1 on migration and invasion of MCF7 cells

Lysine‐specific demethylase 1 (LSD1) is an important histone demethylase that mediates epithelial to mesenchymal transition (EMT). The E239K mutation of LSD1 was identified in a luminal breast cancer patient from the COSMIC Breast Cancer dataset. To investigate the functional effects of the E239K mutation of LSD1, a stable LSD1 knockdown MCF7 cell line was generated. Rescue with WT LSD1, but not E239K mutated LSD1, suppressed the invasion and migration of the LSD1 knockdown cells, indicating that the E239K mutation abolished the suppressive effects of LSD1 on the invasion and migration of MCF7 cells. Further analysis showed that the E239K mutation abolished LSD1‐mediated invasion and migration of MCF7 cells through downregulation of estrogen receptor α (ERα). Most importantly, the E239K mutation disrupted the interaction between LSD1 and GATA3, which reduced the enrichment of LSD1 at the promoter region of the ERα gene; the reduced enrichment of LSD1 at the promoter region of the ERα gene caused enhanced histone H3K9 methylation, which subsequently suppressed the transcription of the ERα gene. In summary, the E239K mutation abolishes the suppressive function of LSD1 on migration and invasion of breast cancer cells by disrupting the interaction between LSD1 and GATA3.

ARSD, a novel ERα downstream target gene, inhibits proliferation and migration of breast cancer cells via activating Hippo/YAP pathway

Advanced breast cancer (BC), especially basal like triple-negative BC (TNBC), is a highly malignant tumor without viable treatment option, highlighting the urgent need to seek novel therapeutic targets. Arylsulfatase D (ARSD), localized at Xp22.3, is a female-biased gene due to its escaping from X chromosome inactivation (XCI). Unfortunately, no systematic investigation of ARSD on BC has been reported. In this study, we observed that ARSD expression was positively related to ERα status either in BC cells or tissue specimens, which were associated with good prognosis. Furthermore, we found a set of hormone-responsive lineage-specific transcription factors, FOXA1, GATA3, ERα, directly drove high expression of ARSD through chromatin looping in luminal subtype BC cells. Opposingly, ARSD still subjected to XCI in TNBC cells mediated by Xist, CpG islands methylation, and inhibitory histone modification. Unexpectedly, we also found that ectopic ARSD overexpression could inhibit proliferation and migration of TNBC cells by activating Hippo/YAP pathway, indicating that ARSD may be a molecule brake on ERα signaling pathway, which restricted ERα to be an uncontrolled active status. Combined with other peoples' researches that Hippo signaling maintained ER expression and ER + BC growth, we believed that there should exist a regulative feedback loop formation among ERα, ARSD, and Hippo/YAP pathway. Collectively, our findings will help filling the knowledge gap about the influence of ARSD on BC and providing evidence that ARSD may serve as a potential marker to predict prognosis and as a therapeutic target.

Positive Regulation of Estrogen Receptor Alpha in Breast Tumorigenesis

Estrogen receptor alpha (ERα, NR3A1) contributes through its expression in different tissues to a spectrum of physiological processes, including reproductive system development and physiology, bone mass maintenance, as well as cardiovascular and central nervous system functions. It is also one of the main drivers of tumorigenesis in breast and uterine cancer and can be targeted by several types of hormonal therapies. ERα is expressed in a subset of luminal cells corresponding to less than 10% of normal mammary epithelial cells and in over 70% of breast tumors (ER+ tumors), but the basis for its selective expression in normal or cancer tissues remains incompletely understood. The mapping of alternative promoters and regulatory elements has delineated the complex genomic structure of the ESR1 gene and shed light on the mechanistic basis for the tissue-specific regulation of ESR1 expression. However, much remains to be uncovered to better understand how ESR1 expression is regulated in breast cancer. This review recapitulates the current body of knowledge on the structure of the ESR1 gene and the complex mechanisms controlling its expression in breast tumors. In particular, we discuss the impact of genetic alterations, chromatin modifications, and enhanced expression of other luminal transcription regulators on ESR1 expression in tumor cells.

The estrogen receptor/GATA3/FOXA1 transcriptional network: lessons learned from breast cancer

Cellular identity and physiologic function in mammary epithelial cells and in many breast cancers flow from the action of a network of master transcriptional regulators including estrogen receptor alpha, GATA3, and FOXA1. The last decade has seen the completion of multiple large sequencing projects focusing on breast cancer. These massive compendia of sequence data have provided a wealth of new information linking mutation in these transcription factors to alterations in tumor biology and transcriptional program. The emerging details on mutation in cancer, and direct experimental exploration of hypotheses based on it, are now providing a wealth of new information on the roles played by estrogen receptor alpha, GATA3, and FOXA1 in regulating gene transcription and how their combined action contributes to a network shaping cell function in both physiologic and disease states.

Estrogen receptor-α signaling in post-natal mammary development and breast cancers

17β-estradiol controls post-natal mammary gland development and exerts its effects through Estrogen Receptor ERα, a member of the nuclear receptor family. ERα is also critical for breast cancer progression and remains a central therapeutic target for hormone-dependent breast cancers. In this review, we summarize the current understanding of the complex ERα signaling pathways that involve either classical nuclear “genomic” or membrane “non-genomic” actions and regulate in concert with other hormones the different stages of mammary development. We describe the cellular and molecular features of the luminal cell lineage expressing ERα and provide an overview of the transgenic mouse models impacting ERα signaling, highlighting the pivotal role of ERα in mammary gland morphogenesis and function and its implication in the tumorigenic processes. Finally, we describe the main features of the ERα-positive luminal breast cancers and their modeling in mice.

Enhancer RNAs Mediate Estrogen-Induced Decommissioning of Selective Enhancers by Recruiting ERα and Its Cofactor

The function of enhancer RNAs (eRNAs) in transcriptional regulation remains obscure. By analyzing the genome-wide nascent transcript profiles in breast cancer cells, we identify a special group of eRNAs that are essential for estrogen-induced transcriptional repression. Using eRNAs of TM4SF1 and EFEMP1 as the paradigms, we find that these RNA molecules not only stabilize promoter-enhancer interactions but also recruit liganded estrogen receptor α (ERα) to particular enhancer regions, facilitate the formation of a functional transcriptional complex, and cause gene silencing. Interestingly, ERα is shown to directly bind with eRNAs by its DNA-binding domain. These eRNAs help with the formation of a specific ERα-centered transcriptional complex and promote the association of the histone demethylase KDM2A, which dismisses RNA polymerase II from designated enhancers and suppresses the transcription of target genes. Our work demonstrates a complete mechanism underlying the action of eRNAs in modulating and refining the locus-specific transcriptional program.

Yang et al. identified a group of eRNAs that are essential for estrogen-induced transcriptional repression by assisting with the chromatin recruitment of ERα through binding to its DNA-binding domain and facilitating the interaction of ERα with its cofactors, which leads to the dismissal of RNA polymerase II.

Subtype-associated epigenomic landscape and 3D genome structure in bladder cancer

Muscle-invasive bladder cancers are characterized by their distinct expression of luminal and basal genes, which could be used to predict key clinical features such as disease progression and overall survival. Transcriptionally, FOXA1, GATA3, and PPARG are shown to be essential for luminal subtype-specific gene regulation and subtype switching, while TP63, STAT3, and TFAP2 family members are critical for regulation of basal subtype-specific genes. Despite these advances, the underlying epigenetic mechanisms and 3D chromatin architecture responsible for subtype-specific regulation in bladder cancer remain unknown. RESULT: We determine the genome-wide transcriptome, enhancer landscape, and transcription factor binding profiles of FOXA1 and GATA3 in luminal and basal subtypes of bladder cancer. Furthermore, we report the first-ever mapping of genome-wide chromatin interactions by Hi-C in both bladder cancer cell lines and primary patient tumors. We show that subtype-specific transcription is accompanied by specific open chromatin and epigenomic marks, at least partially driven by distinct transcription factor binding at distal enhancers of luminal and basal bladder cancers. Finally, we identify a novel clinically relevant transcription factor, Neuronal PAS Domain Protein 2 (NPAS2), in luminal bladder cancers that regulates other subtype-specific genes and influences cancer cell proliferation and migration. CONCLUSION: In summary, our work identifies unique epigenomic signatures and 3D genome structures in luminal and basal urinary bladder cancers and suggests a novel link between the circadian transcription factor NPAS2 and a clinical bladder cancer subtype.

TET2 is a component of the estrogen receptor complex and controls 5mC to 5hmC conversion at estrogen receptor cis-regulatory regions

Estrogen receptor-α (ER) drives tumor development in ER-positive (ER+) breast cancer. The transcription factor GATA3 has been closely linked to ER function, but its precise role in this setting remains unclear. Quantitative proteomics was used to assess changes to the ER complex in response to GATA3 depletion. Unexpectedly, few proteins were lost from the ER complex in the absence of GATA3, with the only major change being depletion of the dioxygenase TET2. TET2 binding constituted a near-total subset of ER binding in multiple breast cancer models, with loss of TET2 associated with reduced activation of proliferative pathways. TET2 knockdown did not appear to change global methylated cytosine (5mC) levels; however, oxidation of 5mC to 5-hydroxymethylcytosine (5hmC) was significantly reduced, and these events occurred at ER enhancers. These findings implicate TET2 in the maintenance of 5hmC at ER sites, providing a potential mechanism for TET2-mediated regulation of ER target genes.

GATA transcription factors, SOX17 and TFAP2C, drive the human germ-cell specification program

This work shows that GATA transcription factors transduce the BMP signaling and, with SOX17 and TFAP2C, induce the human germ-cell fate, delineating the mechanism for human germ-cell specification.

The in vitro reconstitution of human germ-cell development provides a robust framework for clarifying key underlying mechanisms. Here, we explored transcription factors (TFs) that engender the germ-cell fate in their pluripotent precursors. Unexpectedly, SOX17, TFAP2C, and BLIMP1, which act under the BMP signaling and are indispensable for human primordial germ-cell-like cell (hPGCLC) specification, failed to induce hPGCLCs. In contrast, GATA3 or GATA2, immediate BMP effectors, combined with SOX17 and TFAP2C, generated hPGCLCs. GATA3/GATA2 knockouts dose-dependently impaired BMP-induced hPGCLC specification, whereas GATA3/GATA2 expression remained unaffected in SOX17, TFAP2C, or BLIMP1 knockouts. In cynomolgus monkeys, a key model for human development, GATA3, SOX17, and TFAP2C were co-expressed exclusively in early PGCs. Crucially, the TF-induced hPGCLCs acquired a hallmark of bona fide hPGCs to undergo epigenetic reprogramming and mature into oogonia/gonocytes in xenogeneic reconstituted ovaries. By uncovering a TF circuitry driving the germ line program, our study provides a paradigm for TF-based human gametogenesis.

Emerging role of pioneer transcription factors in targeted ERα positive breast cancer

Transcription factors (TFs) are modular protein groups that preferably bind to DNA sequences and guide genomic expression through transcription. Among these key regulators, “pioneer factors” are an emerging class of TFs that specifically interact with nucleosomal DNA and facilitate accessible genomic binding sites for the additional TFs. There is growing evidence of these specialized modulators in particular malignancies, as highlighted by agents’ clinical efficacy, specifically targeting nuclear hormone receptors. They have been implicated in multiple cancers more recently, with a high proportion inculpating on hormone influential cancers. Moreover, extended crosstalk and cooperation between ERα pioneering factors in estrogen-dependent breast cancer (BC) remain elucidated. This review discusses on the recent advances in our understanding of pioneer TFs in cancer, especially highlighting its potentiality to modulate chromatin condensation to permit ERα recruitment in BC cells. Through the study it was concluded that the highly prospected pioneer TFs in BC, including FOXA1, TLE1, PBX1, and GATA3, possess the potential therapeutic significance and further innovations in the field could yield targeted therapy in cancer treatment.

Integrated genomic and transcriptomic analysis reveals unique characteristics of hepatic metastases and pro-metastatic role of complement C1q in pancreatic ductal adenocarcinoma

BACKGROUND

Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal cancers due to its high metastasis rate in the liver. However, little is known about the molecular features of hepatic metastases due to difficulty in obtaining fresh tissues and low tumor cellularity.

RESULTS

We conduct exome sequencing and RNA sequencing for synchronous surgically resected primary tumors and the paired hepatic metastases from 17 hepatic oligometastatic pancreatic ductal adenocarcinoma and validate our findings in specimens from 35 of such cases. The comprehensive analysis of somatic mutations, copy number alterations, and gene expressions show high similarity between primary tumors and hepatic metastases. However, hepatic metastases also show unique characteristics, such as a higher degree of 3p21.1 loss, stronger abilities of proliferation, downregulation of epithelial to mesenchymal transition activity, and metabolic rewiring. More interesting, altered tumor microenvironments are observed in hepatic metastases, especially a higher proportion of tumor infiltrating M2 macrophage and upregulation of complement cascade. Further experiments demonstrate that expression of C1q increases in primary tumors and hepatic metastases, C1q is mainly produced by M2 macrophage, and C1q promotes migration and invasion of PDAC cells.

CONCLUSION

Taken together, we find potential factors that contribute to different stages of PDAC metastasis. Our study broadens the understanding of molecular mechanisms driving PDAC metastasis.

Chromatin Modifiers in Transcriptional Regulation: New Findings and Prospects

Histone-modifying and remodeling complexes are considered the main coregulators that affect transcription by changing the chromatin structure. Coordinated action by these complexes is essential for the transcriptional activation of any eukaryotic gene. In this review, we discuss current trends in the study of histone modifiers and chromatin remodelers, including the functional impact of transcriptional proteins/ complexes i.e., "pioneers"; remodeling and modification of non-histone proteins by transcriptional complexes; the supplementary functions of the non-catalytic subunits of remodelers, and the participation of histone modifiers in the "pause" of RNA polymerase II. The review also includes a scheme illustrating the mechanisms of recruitment of the main classes of remodelers and chromatin modifiers to various sites in the genome and their functional activities.

Role of GATA3 exon 6 germline mutations in breast cancer progression in Egyptian female patients

IMPACT STATEMENT

GATA3 mutations are known to play an important role in breast cancer progression. The exact role and mechanisms of these mutations remain controversial as some studies suggest a relation to breast tumor growth, while others suggest a relation to longer survival. GATA3 germline mutations are not well studied in breast cancer. In this study, it was hypothesized that different types of GATA3 mutations could contribute to the breast cancer progression in different ways. GATA3 exon 6, which is important for GATA3 protein functions, was reported to have hotspots, and hence it was selected for study. Intronic GATA3 germline mutations were found to be related to favorable prognosis, while protein coding mutations were found to be related to unfavorable prognosis. Bioinformatics study of large publically available datasets showed that GATA3 mutations lead to dysregulation of pathways related to T-cells activation, inflammation, and breast cancer development.

TET2 directs mammary luminal cell differentiation and endocrine response

Epigenetic regulation plays an important role in governing stem cell fate and tumorigenesis. Lost expression of a key DNA demethylation enzyme TET2 is associated with human cancers and has been linked to stem cell traits in vitro; however, whether and how TET2 regulates mammary stem cell fate and mammary tumorigenesis in vivo remains to be determined. Here, using our recently established mammary specific Tet2 deletion mouse model, the data reveals that TET2 plays a pivotal role in mammary gland development and luminal lineage commitment. We show that TET2 and FOXP1 form a chromatin complex that mediates demethylation of ESR1, GATA3, and FOXA1, three key genes that are known to coordinately orchestrate mammary luminal lineage specification and endocrine response, and also are often silenced by DNA methylation in aggressive breast cancers. Furthermore, Tet2 deletion-PyMT breast cancer mouse model exhibits enhanced mammary tumor development with deficient ERα expression that confers tamoxifen resistance in vivo. As a result, this study elucidates a role for TET2 in governing luminal cell differentiation and endocrine response that underlies breast cancer resistance to anti-estrogen treatments.

TET2 loss is associated with human cancers but its role in the mammary gland development and tumorigenesis is unclear. Here, the authors show that TET2–FOXP1 complex mediates demethylation of genes involved in luminal lineage commitment and endocrine response, underlying a role of TET2 loss in endocrine resistant breast cancer.

Cancer-specific mutation of GATA3 disrupts the transcriptional regulatory network governed by Estrogen Receptor alpha, FOXA1 and GATA3

Estrogen receptors (ER) are activated by the steroid hormone 17β-estradiol. Estrogen receptor alpha (ER-α) forms a regulatory network in mammary epithelial cells and in breast cancer with the transcription factors FOXA1 and GATA3. GATA3 is one of the most frequently mutated genes in breast cancer and is capable of specifying chromatin localization of FOXA1 and ER-α. How GATA3 mutations found in breast cancer impact genomic localization of ER-α and the transcriptional network downstream of ER-α and FOXA1 remains unclear. Here, we investigate the function of a recurrent patient-derived GATA3 mutation (R330fs) on this regulatory network. Genomic analysis indicates that the R330fs mutant can disrupt localization of ER-α and FOXA1. Loci co-bound by all three factors are enriched for genes integral to mammary gland development as well as epithelial cell biology. This gene set is differentially regulated in GATA3 mutant cells in culture and in tumors bearing similar mutations in vivo. The altered distribution of ER-α and FOXA1 in GATA3-mutant cells is associated with altered chromatin architecture, which leads to differential gene expression. These results suggest an active role for GATA3 zinc finger 2 mutants in ER-α positive breast tumors.

STAT3 and GR Cooperate to Drive Gene Expression and Growth of Basal-Like Triple-Negative Breast Cancer

Breast cancers are divided into subtypes with different prognoses and treatment responses based on global differences in gene expression. Luminal breast cancer gene expression and proliferation are driven by estrogen receptor alpha, and targeting this transcription factor is the most effective therapy for this subtype. By contrast, it remains unclear which transcription factors drive the gene expression signature that defines basal-like triple-negative breast cancer, and there are no targeted therapies approved to treat this aggressive subtype. In this study, we utilized integrated genomic analysis of DNA methylation, chromatin accessibility, transcription factor binding, and gene expression in large collections of breast cancer cell lines and patient tumors to identify transcription factors responsible for the basal-like gene expression program. Glucocorticoid receptor (GR) and STAT3 bind to the same genomic regulatory regions, which were specifically open and unmethylated in basal-like breast cancer. These transcription factors cooperated to regulate expression of hundreds of genes in the basal-like gene expression signature, which were associated with poor prognosis. Combination treatment with small-molecule inhibitors of both transcription factors resulted in synergistic decreases in cell growth in cell lines and patient-derived organoid models. This study demonstrates that GR and STAT3 cooperate to regulate the basal-like breast cancer gene expression program and provides the basis for improved therapy for basal-like triple-negative breast cancer through rational combination of STAT3 and GR inhibitors. SIGNIFICANCE: This study demonstrates that GR and STAT3 cooperate to activate the canonical gene expression signature of basal-like triple-negative breast cancer and that combination treatment with STAT3 and GR inhibitors could provide synergistic therapeutic efficacy.

The GATA3 X308_Splice breast cancer mutation is a hormone context-dependent oncogenic driver

As the catalog of oncogenic driver mutations is expanding, it becomes clear that alterations in a given gene might have different functions and should not be lumped into one class. The transcription factor GATA3 is a paradigm of this. We investigated the functions of the most common GATA3 mutation (X308_Splice) and five additional mutations, which converge into a neoprotein that we called "neoGATA3," associated with excellent prognosis in patients. Analysis of available molecular data from >3000 breast cancer patients revealed a dysregulation of the ER-dependent transcriptional response in tumors carrying neoGATA3-generating mutations. Mechanistic studies in vitro showed that neoGATA3 interferes with the transcriptional programs controlled by estrogen and progesterone receptors, without fully abrogating them. ChIP-Seq analysis indicated that ER binding is reduced in neoGATA3-expressing cells, especially at distal regions, suggesting that neoGATA3 interferes with the fine tuning of ER-dependent gene expression. This has opposite outputs in distinct hormonal context, having pro- or anti-proliferative effects, depending on the estrogen/progesterone ratio. Our data call for functional analyses of putative cancer drivers to guide clinical application.

ASH2L is involved in promotion of endometrial cancer progression via upregulation of PAX2 transcription

Absent, small or homeotic 2‐like protein (ASH2L) is a core component of a multimeric histone methyltransferase complex that is involved in the maintenance of active transcription, participating in several cancers, however the biological function and molecular mechanism of ASH2L in endometrial cancer (ECa) are largely unknown. Endometrial cancer is a common malignant tumor in women and the incidence of this cancer is on the rise. Estrogen‐ERα signaling, as an oncogenic pathway, plays a crucial role in endometrial carcinogenesis. Therefore, further exploration of the molecular mechanisms around ERα‐mediated gene transcription in ECa would be helpful to the understanding of tumor development and to finding a new therapeutic target for ECa. Here, our study demonstrated that ASH2L was highly expressed in ECa samples, and higher expression of ASH2L was positively correlated with a poor prognosis. Moreover, we identified that ASH2L associated with ERα and that knockdown of ASH2L resulted in decreased expression of a subset of the estrogen‐induced target genes, including paired box 2 (PAX2), an oncogenic gene in ECa. ASH2L was recruited to cis‐regulatory elements in PAX2, thereby altering histone H3K4me3 and H3K27me3 levels, to enhance ERα‐mediated transactivation. Finally, depletion of ASH2L suppressed endometrial cancer cell proliferation and migration. Our findings suggest that ASH2L participates in the promotion of ECa progression, if not totally at least partially, via upregulation of PAX2 transcription.

NR2F2 Orphan Nuclear Receptor is Involved in Estrogen Receptor Alpha-Mediated Transcriptional Regulation in Luminal A Breast Cancer Cells

Nuclear Receptor Subfamily 2 Group F Member 2 (NR2F2) is a member of the steroid/thyroid hormone receptor superfamily with a crucial role in organogenesis, angiogenesis, cardiovascular development and tumorigenesis. However, there is limited knowledge about the cistrome and transcriptome of NR2F2 in breast cancer. In this study, we mapped the regulatory mechanism by NR2F2 using functional genomic methods. To investigate the clinical significance of NR2F2 in breast cancer, The Cancer Genome Atlas (TCGA) data were used. These results show that a high NR2F2 is associated with better survival of a specific subset of patients, namely those with luminal A breast cancer. Therefore, genome-wide NR2F2 and estrogen receptor alpha (ERα) binding sites were mapped in luminal A breast cancer cells using chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-Seq), revealing that most NR2F2 overlap with ERα that are co-occupied by forkhead box A1 (FOXA1) and GATA binding protein 3 (GATA3) in active enhancer regions. NR2F2 overlaps with highly frequent ERα chromatin interactions, which are essential for the formation of ERα-bound super-enhancers. In the process of the transcriptome profiling of NR2F2-depleted breast cancer cells such differentially expressed genes have been identified that are involved in endocrine therapy resistance and are also ERα target genes. Overall, these findings demonstrate that the NR2F2 nuclear receptor has a key role in ERα-mediated transcription and it can offer a potential therapeutic target in patients with luminal A breast cancer.

Liver Tropism in Cancer: The Hepatic Metastatic Niche

The liver is the largest organ in the human body and is prone for cancer metastasis. Although the metastatic pattern can differ depending on the cancer type, the liver is the organ to which cancer cells most frequently metastasize for the majority of prevalent malignancies. The liver is unique in several aspects: the vascular structure is highly permeable and has unparalleled dual blood connectivity, and the hepatic tissue microenvironment presents a natural soil for the seeding of disseminated tumor cells. Although 70% of the liver is composed of the parenchymal hepatocytes, the remaining 30% is composed of nonparenchymal cells including Kupffer cells, liver sinusoidal endothelial cells, and hepatic stellate cells. Recent discoveries show that both the parenchymal and the nonparenchymal cells can modulate each step of the hepatic metastatic cascade, including the initial seeding and colonization as well as the decision to undergo dormancy versus outgrowth. Thus, a better understanding of the molecular mechanisms orchestrating the formation of a hospitable hepatic metastatic niche and the identification of the drivers supporting this process is critical for the development of better therapies to stop or at least decrease liver metastasis. The focus of this perspective is on the bidirectional interactions between the disseminated cancer cells and the unique hepatic metastatic niche.

Amphiregulin deletion strongly attenuates the development of estrogen receptor-positive tumors in p53 mutant mice

PURPOSE

The epidermal growth factor receptor ligand, Amphiregulin, is a transcriptional target of estrogen receptor alpha and is required for pubertal mammary gland development. Previous studies using immortalized human breast cancer cell line xenografts have suggested that Amphiregulin may be an important effector of estrogen receptor alpha during breast cancer development, at least in immune-compromised animals. Here, we evaluate the requirement for Amphiregulin in an immune-competent mouse model which is prone to developing estrogen receptor-positive tumors.

METHODS

We have intercrossed mice with mammary-specific mutation of p53 with mice deficient in Amphiregulin in order to assess the requirement for Amphiregulin in the initiation and progression of both estrogen receptor-positive and estrogen receptor-negative mammary tumors.

RESULTS

Deletion of Amphiregulin significantly delayed the onset of palpable mammary tumors and also strongly reduced the proportion of estrogen receptor alpha-positive tumors formed. Upon necropsy, no substantial differences in the prevalence of non-palpable lesions were observed between cohorts, suggesting that the importance of Amphiregulin in mammary tumorigenesis is limited to the post-initiation phase.

CONCLUSIONS

This study underlines the importance of the EGFR ligand, Amphiregulin, as a key mediator of estrogen receptor action in breast cancer.

FOXA1 upregulation promotes enhancer and transcriptional reprogramming in endocrine-resistant breast cancer

Forkhead box A1 (FOXA1) is a pioneer factor that facilitates chromatin binding and function of lineage-specific and oncogenic transcription factors. Hyperactive FOXA1 signaling due to gene amplification or overexpression has been reported in estrogen receptor-positive (ER⁺) endocrine-resistant metastatic breast cancer. However, the molecular mechanisms by which FOXA1 up-regulation promotes these processes and the key downstream targets of the FOXA1 oncogenic network remain elusive. Here, we demonstrate that FOXA1 overexpression in ER⁺ breast cancer cells drives genome-wide enhancer reprogramming to activate prometastatic transcriptional programs. Up-regulated FOXA1 employs superenhancers (SEs) to synchronize transcriptional reprogramming in endocrine-resistant breast cancer cells, reflecting an early embryonic development process. We identify the hypoxia-inducible transcription factor hypoxia-inducible factor-2α (HIF-2α) as the top high FOXA1-induced SE target, mediating the impact of high FOXA1 in activating prometastatic gene sets and pathways associated with poor clinical outcome. Using clinical ER⁺/HER2^- metastatic breast cancer datasets, we show that the aberrant FOXA1/HIF-2α transcriptional axis is largely nonconcurrent with the ESR1 mutations, suggesting different mechanisms of endocrine resistance and treatment strategies. We further demonstrate the selective efficacy of an HIF-2α antagonist, currently in clinical trials for advanced kidney cancer and recurrent glioblastoma, in reducing the clonogenicity, migration, and invasion of endocrine-resistant breast cancer cells expressing high FOXA1. Our study has uncovered high FOXA1-induced enhancer reprogramming and HIF-2α-dependent transcriptional programs as vulnerable targets for treating endocrine-resistant and metastatic breast cancer.

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.

GATA3 recruits UTX for gene transcriptional activation to suppress metastasis of breast cancer

GATA3 has emerged as a prominent transcription factor required for maintaining mammary-gland homeostasis. GATA3 loss is associated with aggressive breast cancer development, but the mechanism by which breast cancer is affected by the loss of GATA3 function remains unclear. Here, we report that GATA3 expression is positively correlated with the expression of UTX, a histone H3K27 demethylase contained in the MLL4 methyltransferase complex, and that GATA3 recruits the chromatin-remodeling MLL4 complex and interacts directly with UTX, ASH2L, and RBBP5. Using RNA sequencing and chromatin immunoprecipitation and sequencing, we demonstrate that the GATA3/UTX complex synergistically regulates a cohort of genes including Dicer and UTX, which are critically involved in the epithelial-to-mesenchymal transition (EMT). Our results further show that the GATA3-UTX-Dicer axis inhibits EMT, invasion, and metastasis of breast cancer cells in vitro and the dissemination of breast cancer in vivo. Our study implicates the GATA3-UTX-Dicer axis in breast cancer metastasis and provides new mechanistic insights into the pathophysiological function of GATA3.

The Lineage Determining Factor GRHL2 Collaborates with FOXA1 to Establish a Targetable Pathway in Endocrine Therapy-Resistant Breast Cancer

Notwithstanding the positive clinical impact of endocrine therapies in estrogen receptor-alpha (ERα)-positive breast cancer, de novo and acquired resistance limits the therapeutic lifespan of existing drugs. Taking the position that resistance is nearly inevitable, we undertook a study to identify and exploit targetable vulnerabilities that were manifest in endocrine therapy-resistant disease. Using cellular and mouse models of endocrine therapy-sensitive and endocrine therapy-resistant breast cancer, together with contemporary discovery platforms, we identified a targetable pathway that is composed of the transcription factors FOXA1 and GRHL2, a coregulated target gene, the membrane receptor LYPD3, and the LYPD3 ligand, AGR2. Inhibition of the activity of this pathway using blocking antibodies directed against LYPD3 or AGR2 inhibits the growth of endocrine therapy-resistant tumors in mice, providing the rationale for near-term clinical development of humanized antibodies directed against these proteins.

REGGAE: a novel approach for the identification of key transcriptional regulators

Motivation

Transcriptional regulators play a major role in most biological processes. Alterations in their activities are associated with a variety of diseases and in particular with tumor development and progression. Hence, it is important to assess the effects of deregulated regulators on pathological processes.

Results

Here, we present REGulator-Gene Association Enrichment (REGGAE), a novel method for the identification of key transcriptional regulators that have a significant effect on the expression of a given set of genes, e.g. genes that are differentially expressed between two sample groups. REGGAE uses a Kolmogorov-Smirnov-like test statistic that implicitly combines associations between regulators and their target genes with an enrichment approach to prioritize the influence of transcriptional regulators. We evaluated our method in two different application scenarios, which demonstrate that REGGAE is well suited for uncovering the influence of transcriptional regulators and is a valuable tool for the elucidation of complex regulatory mechanisms.

Availability and implementation

REGGAE is freely available at https://regulatortrail.bioinf.uni-sb.de.

Supplementary information

Supplementary data are available at Bioinformatics online.

GATA factor transcriptional activity: Insights from genome-wide binding profiles

The members of the GATA family of transcription factors have homologous zinc fingers and bind to similar sequence motifs. Recent advances in genome-wide technologies and the integration of bioinformatics data have led to a better understanding of how GATA factors regulate gene expression; GATA-factor-induced transcriptional and epigenetic changes have now been analyzed at unprecedented levels of detail. Here, we review the results of genome-wide studies of GATA factor occupancy in human and murine cell lines and primary cells (as determined by chromatin immunoprecipitation sequencing), and then discuss the molecular mechanisms underlying the mediation of transcriptional and epigenetic regulation by GATA factors.

Meta-analysis of Chromatin Programming by Steroid Receptors

Transcription factors (TFs) must access chromatin to bind to their response elements and regulate gene expression. A widely accepted model proposes that only a special subset of TFs, pioneer factors, can associate with condensed chromatin and initiate chromatin opening. We previously reported that steroid receptors (SRs), not considered pioneer factors, can assist the binding of an archetypal pioneer, the forkhead box protein 1 (FOXA1), at a subset of receptor-activated enhancers. These findings have been challenged recently, with the suggestion that newly acquired data fully support the prevailing pioneer model. Here, we reexamine our results and confirm the original conclusions. We also analyze and discuss a number of available datasets relevant to chromatin penetration by SRs and find a general consensus supporting our original observations. Hence, we propose that chromatin opening at some sites can be initiated by SRs, with a parallel recruitment of factors often treated as having a unique pioneer function. This Matters Arising paper is in response to Glont et al. (2019), published in Cell Reports.

Epigenetic re-wiring of breast cancer by pharmacological targeting of C-terminal binding protein

The C-terminal binding protein (CtBP) is an NADH-dependent dimeric family of nuclear proteins that scaffold interactions between transcriptional regulators and chromatin-modifying complexes. Its association with poor survival in several cancers implicates CtBP as a promising target for pharmacological intervention. We employed computer-assisted drug design to search for CtBP inhibitors, using quantitative structure-activity relationship (QSAR) modeling and docking. Functional screening of these drugs identified 4 compounds with low toxicity and high water solubility. Micro molar concentrations of these CtBP inhibitors produces significant de-repression of epigenetically silenced pro-epithelial genes, preferentially in the triple-negative breast cancer cell line MDA-MB-231. This epigenetic reprogramming occurs through eviction of CtBP from gene promoters; disrupted recruitment of chromatin-modifying protein complexes containing LSD1, and HDAC1; and re-wiring of activating histone marks at targeted genes. In functional assays, CtBP inhibition disrupts CtBP dimerization, decreases cell migration, abolishes cellular invasion, and improves DNA repair. Combinatorial use of CtBP inhibitors with the LSD1 inhibitor pargyline has synergistic influence. Finally, integrated correlation of gene expression in breast cancer patients with nuclear levels of CtBP1 and LSD1, reveals new potential therapeutic vulnerabilities. These findings implicate a broad role for this class of compounds in strategies for epigenetically targeted therapeutic intervention.

GATA3 Truncating Mutations Promote Cistromic Re-Programming In Vitro, but Not Mammary Tumor Formation in Mice

Heterozygous mutations in the transcription factor GATA3 are identified in 10-15% of all breast cancer cases. Most of these are protein-truncating mutations, concentrated within or downstream of the second GATA-type zinc-finger domain. Here, we investigated the functional consequences of expression of two truncated GATA3 mutants, in vitro in breast cancer cell lines and in vivo in the mouse mammary gland. We found that the truncated GATA3 mutants display altered DNA binding activity caused by preferred tethering through FOXA1. In addition, expression of the truncated GATA3 mutants reduces E-cadherin expression and promotes anchorage-independent growth in vitro. However, we could not identify any effects of truncated GATA3 expression on mammary gland development or mammary tumor formation in mice. Together, our results demonstrate that both truncated GATA3 mutants promote cistromic re-programming of GATA3 in vitro, but these mutants are not sufficient to induce tumor formation in mice.

Cooperativity of co-factor NR2F2 with Pioneer Factors GATA3, FOXA1 in promoting ERα function

Estrogen receptor α (ERα) drives growth in the majority of human breast cancers by binding to regulatory elements and inducing transcriptional events that promote tumor growth. ERα binding activity largely depends on access to binding sites on chromatin, which is facilitated in part by Pioneer Factors (PFs). Transcription factors operate in complexes through thousands of genomic binding sites in a combinatorial fashion to control the expression of genes. However, the extent of crosstalk and cooperation between ERα pioneer factors and more collaborative transcription factors in breast cancer still remains to be elucidated systematically. Methods: Here, we determined the genomic binding information of 40 transcription-related factors and histone modifications with ChIP-seq in ENCODE and integrated it with other genomic information (RNA-seq, ATAC-seq, Gene microarray, 450k methylation chip, GRO-seq), forming a multi-dimension network to illuminate ERα associated transcription. Results: We show that transcription factor, NR2F2 binds to most sites independently of estrogen. Perturbation of NR2F2 expression decreases ERα DNA binding, chromatin openning, and estrogen-dependent cell growth. In the genome-wide analysis, we show that most binding events of NR2F2 and known pioneer factors FOXA1, GATA3 occur together, covering 85% of the ERα binding sites. Regions bound by all the three TFs appeared to be the most active, to have the strongest ERα binding and to be enriched for the super enhancers. Conclusions: The ERα binds to pre-accessible sites containing ERE elements bound by the three transcription factors (NR2F2, FOXA1 and GATA3).The three genes were also identified to correlate with decreased metastatic potential in patient cohorts and co-regulate each other. Together, our results suggest that NR2F2 is a cofactor with FOXA1 and GATA3 in ERα-mediated transcription.

The Histone Demethylase Enzymes KDM3A and KDM4B Co-Operatively Regulate Chromatin Transactions of the Estrogen Receptor in Breast Cancer

Many estrogen receptor (ER)-positive breast cancers develop resistance to endocrine therapy but retain canonical receptor signalling in the presence of selective ER antagonists. Numerous co-regulatory proteins, including enzymes that modulate the chromatin environment, control the transcriptional activity of the ER. Targeting ER co-regulators has therefore been proposed as a novel therapeutic approach. By assessing DNA-binding dynamics in ER-positive breast cancer cells, we have identified that the histone H3 lysine 9 demethylase enzymes, KDM3A and KDM4B, co-operate to regulate ER activity via an auto-regulatory loop that facilitates the recruitment of each co-activating enzyme to chromatin. We also provide evidence that suggests that KDM3A primes chromatin for deposition of the ER pioneer factor FOXA1 and recruitment of the ER-transcriptional complex, all prior to ER recruitment, therefore establishing an important mechanism of chromatin regulation involving histone demethylases and pioneer factors, which controls ER functionality. Importantly, we show via global gene-expression analysis that a KDM3A/KDM4B/FOXA1 co-regulated gene signature is enriched for pro-proliferative and ER-target gene sets, suggesting that abrogation of this network could be an efficacious therapeutic strategy. Finally, we show that depletion of both KDM3A and KDM4B has a greater inhibitory effect on ER activity and cell growth than knockdown of each individual enzyme, suggesting that targeting both enzymes represents a potentially efficacious therapeutic option for ER-driven breast cancer.

Single-cell RNA-seq highlights intra-tumoral heterogeneity and malignant progression in pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is the most common type of pancreatic cancer featured with high intra-tumoral heterogeneity and poor prognosis. To comprehensively delineate the PDAC intra-tumoral heterogeneity and the underlying mechanism for PDAC progression, we employed single-cell RNA-seq (scRNA-seq) to acquire the transcriptomic atlas of 57,530 individual pancreatic cells from primary PDAC tumors and control pancreases, and identified diverse malignant and stromal cell types, including two ductal subtypes with abnormal and malignant gene expression profiles respectively, in PDAC. We found that the heterogenous malignant subtype was composed of several subpopulations with differential proliferative and migratory potentials. Cell trajectory analysis revealed that components of multiple tumor-related pathways and transcription factors (TFs) were differentially expressed along PDAC progression. Furthermore, we found a subset of ductal cells with unique proliferative features were associated with an inactivation state in tumor-infiltrating T cells, providing novel markers for the prediction of antitumor immune response. Together, our findings provide a valuable resource for deciphering the intra-tumoral heterogeneity in PDAC and uncover a connection between tumor intrinsic transcriptional state and T cell activation, suggesting potential biomarkers for anticancer treatment such as targeted therapy and immunotherapy.