FOXA1 Directs H3K4 Monomethylation at Enhancers via Recruitment of the Methyltransferase MLL3

Methylation of the chromatin modifier KMT2D by SMYD2 contributes to therapeutic response in hormone-dependent breast cancer

Activating mutations in PIK3CA are frequently found in estrogen-receptor-positive (ER+) breast cancer, and the combination of the phosphatidylinositol 3-kinase (PI3K) inhibitor alpelisib with anti-ER inhibitors is approved for therapy. We have previously demonstrated that the PI3K pathway regulates ER activity through phosphorylation of the chromatin modifier KMT2D. Here, we discovered a methylation site on KMT2D, at K1330 directly adjacent to S1331, catalyzed by the lysine methyltransferase SMYD2. SMYD2 loss attenuates alpelisib-induced KMT2D chromatin binding and alpelisib-mediated changes in gene expression, including ER-dependent transcription. Knockdown or pharmacological inhibition of SMYD2 sensitizes breast cancer cells, patient-derived organoids, and tumors to PI3K/AKT inhibition and endocrine therapy in part through KMT2D K1330 methylation. Together, our findings uncover a regulatory crosstalk between post-translational modifications that fine-tunes KMT2D function at the chromatin. This provides a rationale for the use of SMYD2 inhibitors in combination with PI3Kα/AKT inhibitors in the treatment of ER+/PIK3CA mutant breast cancer.

H3K4me1 facilitates promoter-enhancer interactions and gene activation during embryonic stem cell differentiation

Histone H3 lysine 4 mono-methylation (H3K4me1) marks poised or active enhancers. KMT2C (MLL3) and KMT2D (MLL4) catalyze H3K4me1, but their histone methyltransferase activities are largely dispensable for transcription during early embryogenesis in mammals. To better understand the role of H3K4me1 in enhancer function, we analyze dynamic enhancer-promoter (E-P) interactions and gene expression during neural differentiation of the mouse embryonic stem cells. We found that KMT2C/D catalytic activities were only required for H3K4me1 and E-P contacts at a subset of candidate enhancers, induced upon neural differentiation. By contrast, a majority of enhancers retained H3K4me1 in KMT2C/D catalytic mutant cells. Surprisingly, H3K4me1 signals at these KMT2C/D-independent sites were reduced after acute depletion of KMT2B, resulting in aggravated transcriptional defects. Our observations therefore implicate KMT2B in the catalysis of H3K4me1 at enhancers and provide additional support for an active role of H3K4me1 in enhancer-promoter interactions and transcription in mammalian cells.

KMT2 Family of H3K4 Methyltransferases: Enzymatic Activity-dependent and -independent Functions

Histone-lysine N-methyltransferase 2 (KMT2) methyltransferases are critical for gene regulation, cell differentiation, animal development, and human diseases. KMT2 biological roles are often attributed to their methyltransferase activities on lysine 4 of histone H3 (H3K4). However, recent data indicate that KMT2 proteins also possess non-enzymatic functions. In this review, we discuss the current understanding of KMT2 family, with a focus on their enzymatic activity-dependent and -independent functions. Six mammalian KMT2 proteins of three subgroups, KMT2A/B (MLL1/2), KMT2C/D (MLL3/4), and KMT2F/G (SETD1A/B or SET1A/B), have shared and distinct protein domains, catalytic substrates, genomic localizations, and associated complex subunits. Recent studies have revealed the importance of KMT2C/D in enhancer regulation, differentiation, development, tumor suppression and highlighted KMT2C/D enzymatic activity-dependent and -independent roles in mouse embryonic development and cell differentiation. Catalytic dependent and independent functions for KMT2A/B and KMT2F/G in gene regulation, differentiation, and development are less understood. Finally, we provide our perspectives and lay out future research directions that may help advance the investigation on enzymatic activity-dependent and -independent biological roles and working mechanisms of KMT2 methyltransferases.

Transcriptomic dysregulation and autistic-like behaviors in Kmt2c haploinsufficient mice rescued by an LSD1 inhibitor

Recent studies have consistently demonstrated that the regulation of chromatin and gene transcription plays a pivotal role in the pathogenesis of neurodevelopmental disorders. Among many genes involved in these pathways, KMT2C, encoding one of the six known histone H3 lysine 4 (H3K4) methyltransferases in humans and rodents, was identified as a gene whose heterozygous loss-of-function variants are causally associated with autism spectrum disorder (ASD) and the Kleefstra syndrome phenotypic spectrum. However, little is known about how KMT2C haploinsufficiency causes neurodevelopmental deficits and how these conditions can be treated. To address this, we developed and analyzed genetically engineered mice with a heterozygous frameshift mutation of Kmt2c (Kmt2c+/fs mice) as a disease model with high etiological validity. In a series of behavioral analyses, the mutant mice exhibit autistic-like behaviors such as impairments in sociality, flexibility, and working memory, demonstrating their face validity as an ASD model. To investigate the molecular basis of the observed abnormalities, we performed a transcriptomic analysis of their bulk adult brains and found that ASD risk genes were specifically enriched in the upregulated differentially expressed genes (DEGs), whereas KMT2C peaks detected by ChIP-seq were significantly co-localized with the downregulated genes, suggesting an important role of putative indirect effects of Kmt2c haploinsufficiency. We further performed single-cell RNA sequencing of newborn mouse brains to obtain cell type-resolved insights at an earlier stage. By integrating findings from ASD exome sequencing, genome-wide association, and postmortem brain studies to characterize DEGs in each cell cluster, we found strong ASD-associated transcriptomic changes in radial glia and immature neurons with no obvious bias toward upregulated or downregulated DEGs. On the other hand, there was no significant gross change in the cellular composition. Lastly, we explored potential therapeutic agents and demonstrate that vafidemstat, a lysine-specific histone demethylase 1 (LSD1) inhibitor that was effective in other models of neuropsychiatric/neurodevelopmental disorders, ameliorates impairments in sociality but not working memory in adult Kmt2c+/fs mice. Intriguingly, the administration of vafidemstat was shown to alter the vast majority of DEGs in the direction to normalize the transcriptomic abnormalities in the mutant mice (94.3 and 82.5% of the significant upregulated and downregulated DEGs, respectively, P < 2.2 × 10-16, binomial test), which could be the molecular mechanism underlying the behavioral rescuing. In summary, our study expands the repertoire of ASD models with high etiological and face validity, elucidates the cell-type resolved molecular alterations due to Kmt2c haploinsufficiency, and demonstrates the efficacy of an LSD1 inhibitor that might be generalizable to multiple categories of psychiatric disorders along with a better understanding of its presumed mechanisms of action.

KMT2C and KMT2D aberrations in breast cancer

KMT2C and KMT2D are histone lysine methyltransferases responsible for the monomethylation of histone 3 lysine 4 (H3K4) residues at gene enhancer sites. KMT2C/D are the most frequently mutated histone methyltransferases (HMTs) in breast cancer, occurring at frequencies of 10-20% collectively. Frequent damaging and truncating somatic mutations indicate a tumour-suppressive role of KMT2C/D in breast oncogenesis. Recent studies using cell lines and mouse models to replicate KMT2C/D loss show that these genes contribute to oestrogen receptor (ER)-driven transcription in ER+ breast cancers through the priming of gene enhancer regions. This review provides an overview of the functions of KMT2C/D and outlines the recent clinical and experimental evidence of the roles of KMT2C and KMT2D in breast cancer development.

Out with the old, in with the new: Pioneer transcription factors as activators and repressors of lineage specification genes

In this issue of Molecular Cell, Matsui et al.1 examine lineage determination by pioneer transcription factors, finding that they control cell fate in cooperation with PRDM family members by repressing alternative-lineage and precocious gene expression through establishment of bivalent enhancers.

Pioneer and PRDM transcription factors coordinate bivalent epigenetic states to safeguard cell fate

Pioneer transcription factors (TFs) regulate cell fate by establishing transcriptionally primed and active states. However, cell fate control requires the coordination of both lineage-specific gene activation and repression of alternative-lineage programs, a process that is poorly understood. Here, we demonstrate that the pioneer TF FOXA coordinates with PRDM1 TF to recruit nucleosome remodeling and deacetylation (NuRD) complexes and Polycomb repressive complexes (PRCs), which establish highly occupied, accessible nucleosome conformation with bivalent epigenetic states, thereby preventing precocious and alternative-lineage gene expression during human endoderm differentiation. Similarly, the pioneer TF OCT4 coordinates with PRDM14 to form bivalent enhancers and repress cell differentiation programs in human pluripotent stem cells, suggesting that this may be a common and critical function of pioneer TFs. We propose that pioneer and PRDM TFs coordinate to safeguard cell fate through epigenetic repression mechanisms.

Pioneer factors: roles and their regulation in development

Pioneer factors are a subclass of transcription factors that can bind and initiate opening of silent chromatin regions. Pioneer factors subsequently regulate lineage-specific genes and enhancers and, thus, activate the zygotic genome after fertilization, guide cell fate transitions during development, and promote various forms of human cancers. As such, pioneer factors are useful in directed cell reprogramming. In this review, we define the structural and functional characteristics of pioneer factors, how they bind and initiate opening of closed chromatin regions, and the consequences for chromatin dynamics and gene expression during cell differentiation. We also discuss emerging mechanisms that modulate pioneer factors during development.

Novel germline variants in KMT2C in Chinese patients with Kleefstra syndrome-2

Kleefstra syndrome (KLEFS) refers to a rare inherited neurodevelopmental disorder characterized by intellectual disability (ID), language and motor delays, behavioral abnormalities, abnormal facial appearance, and other variable clinical features. KLEFS is subdivided into two subtypes: Kleefstra syndrome-1 (KLEFS1, OMIM: 610253), caused by a heterozygous microdeletion encompassing the Euchromatic Histone Lysine Methyltransferase 1 (EHMT1) gene on chromosome 9q34.3 or pathogenic variants in the EHMT1 gene, and Kleefstra syndrome-2 (KLEFS2, OMIM: 617768), caused by pathogenic variants in the KMT2C gene. More than 100 cases of KLEFS1 have been reported with pathogenic variants in the EHMT1 gene. However, only 13 patients with KLEFS2 have been reported to date. In the present study, five unrelated Chinese patients were diagnosed with KLEFS2 caused by KMT2C variants through whole-exome sequencing (WES). We identified five different variants of the KMT2C gene in these patients: c.9166C>T (p.Gln3056*), c.9232_9247delCAGCGATCAGAACCGT (p.Gln3078fs*13), c.5068dupA (p.Arg1690fs*10), c.10815_10819delAAGAA (p.Lys3605fs*7), and c.6911_6912insA (p.Met2304fs*8). All five patients had a clinical profile similar to that of patients with KLEFS2. To analyze the correlation between the genotype and phenotype of KLEFS2, we examined 18 variants and their associated phenotypes in 18 patients with KLEFS2. Patients carrying KMT2C variants presented with a wide range of phenotypic defects and an extremely variable phenotype. We concluded that the core phenotypes associated with KMT2C variants were intellectual disability, facial dysmorphisms, language and motor delays, behavioral abnormalities, hypotonia, short stature, and weight loss. Additionally, sex may be one factor influencing the outcome. Our findings expand the phenotypic and genetic spectrum of KLEFS2 and help to clarify the genotype-phenotype correlation.

Pioneer factor Pax7 initiates two-step cell-cycle-dependent chromatin opening

Pioneer transcription factors direct cell differentiation by deploying new enhancer repertoires through their unique ability to target and initiate remodelling of closed chromatin. The initial steps of their action remain undefined, although pioneers have been shown to interact with nucleosomal target DNA and with some chromatin-remodeling complexes. We now define the sequence of events that enables the pioneer Pax7 with its unique abilities. Chromatin condensation exerted by linker histone H1 is the first constraint on Pax7 recruitment, and this establishes the initial speed of chromatin remodeling. The first step of pioneer action involves recruitment of the KDM1A (LSD1) H3K9me2 demethylase for removal of this repressive mark, as well as recruitment of the MLL complex for deposition of the activating H3K4me1 mark. Further progression of pioneer action requires passage through cell division, and this involves dissociation of pioneer targets from perinuclear lamin B. Only then are the SWI-SNF remodeling complex and the coactivator p300 recruited, leading to nucleosome displacement and enhancer activation. Thus, the unique features of pioneer actions are those occurring in the lamin-associated compartment of the nucleus. This model is consistent with previous work that showed a dependence on cell division for establishment of new cell fates.

The chromatin signatures of enhancers and their dynamic regulation

Enhancers are cis-regulatory elements that can stimulate gene expression from distance, and drive precise spatiotemporal gene expression profiles during development. Functional enhancers display specific features including an open chromatin conformation, Histone H3 lysine 27 acetylation, Histone H3 lysine 4 mono-methylation enrichment, and enhancer RNAs production. These features are modified upon developmental cues which impacts their activity. In this review, we describe the current state of knowledge about enhancer functions and the diverse chromatin signatures found on enhancers. We also discuss the dynamic changes of enhancer chromatin signatures, and their impact on lineage specific gene expression profiles, during development or cellular differentiation.

FoxA1/2-dependent epigenomic reprogramming drives lineage switching in lung adenocarcinoma

The ability of cancer cells to alter their identity is essential for tumor survival and progression. Loss of the pulmonary lineage specifier NKX2-1 within KRAS-driven lung adenocarcinoma (LUAD) enhances tumor progression and results in a pulmonary-to-gastric lineage switch that is dependent upon the activity of pioneer factors FoxA1 and FoxA2; however, the underlying mechanism remains largely unknown. Here, we show that FoxA1/2 reprogram the epigenetic landscape of NKX2-1-negative LUAD to facilitate a gastric identity. After Nkx2-1 deletion, FoxA1/2 mediate demethylation of gastric-defining genes through recruitment of TET3, an enzyme that induces DNA demethylation. H3K27ac ChIP-seq and HiChIP show that FoxA1/2 also control the activity of regulatory elements and their 3D interactions at gastric loci. Furthermore, oncogenic KRAS is required for the FoxA1/2-dependent epigenetic reprogramming. This work demonstrates the role of FoxA1/2 in rewiring the methylation and histone landscape and cis-regulatory dynamics of NKX2-1-negative LUAD to drive cancer cell lineage switching.

BRWD3 promotes KDM5 degradation to maintain H3K4 methylation levels

Histone modifications are critical for regulating chromatin structure and gene expression. Dysregulation of histone modifications likely contributes to disease states and cancer. Depletion of the chromatin-binding protein BRWD3 (Bromodomain and WD repeat-containing protein 3), a known substrate-specificity factor of the Cul4-DDB1 E3 ubiquitin ligase complex, results in increased H3K4me1 (H3 lysine 4 monomethylation) levels. The underlying mechanism linking BRWD3 and H3K4 methylation, however, has yet to be defined. Here, we show that depleting BRWD3 not only causes an increase in H3K4me1 levels but also causes a decrease in H3K4me3 (H3 lysine 4 trimethylation) levels, indicating that BRWD3 influences H3K4 methylation more broadly. Using immunoprecipitation coupled to quantitative mass spectrometry, we identified an interaction between BRWD3 and the H3K4-specific lysine demethylase 5 (KDM5/Lid), an enzyme that removes tri- and dimethyl marks from H3K4. Moreover, analysis of ChIP-seq (chromatin immunoprecipitation sequencing) data revealed that BRWD3 and KDM5 are significantly colocalized throughout the genome and H3K4me3 are highly enriched at BRWD3 binding sites. We show that BRWD3 promotes K48-linked polyubiquitination and degradation of KDM5 and that KDM5 degradation is dependent on both BRWD3 and Cul4. Critically, depleting KDM5 fully restores altered H3K4me3 levels and partially restores H3K4me1 levels upon BRWD3 depletion. Together, our results demonstrate that BRWD3 regulates KDM5 activity to balance H3K4 methylation levels.

A Druggable FOXA1-Glucocorticoid Receptor Transcriptional Axis Drives Tumor Growth in a Subset of Non-Small Cell Lung Cancer

The FOXA1 pioneer factor is an essential mediator of steroid receptor function in multiple hormone-dependent cancers, including breast and prostate cancers, enabling nuclear receptors such as estrogen receptor (ER) and androgen receptor (AR) to activate lineage-specific growth programs. FOXA1 is also highly expressed in non-small cell lung cancer (NSCLC), but whether and how it regulates tumor growth in this context is not known. Analyzing data from loss-of-function screens, we identified a subset of NSCLC tumor lines where proliferation is FOXA1 dependent. Using rapid immunoprecipitation and mass spectrometry of endogenous protein, we identified chromatin-localized interactions between FOXA1 and glucocorticoid receptor (GR) in these tumor cells. Knockdown of GR inhibited proliferation of FOXA1-dependent, but not FOXA1-independent NSCLC cells. In these FOXA1-dependent models, FOXA1 and GR cooperate to regulate gene targets involved in EGF signaling and G1-S cell-cycle progression. To investigate the therapeutic potential for targeting this complex, we examined the effects of highly selective inhibitors of the GR ligand-binding pocket and found that GR antagonism with ORIC-101 suppressed FOXA1/GR target expression, activation of EGF signaling, entry into the S-phase, and attendant proliferation in vitro and in vivo. Taken together, our findings point to a subset of NSCLCs harboring a dependence on the FOXA1/GR growth program and provide rationale for its therapeutic targeting.

Significance

NSCLC is the leading cause of cancer deaths worldwide. There is a need to identify novel druggable dependencies. We identify a subset of NSCLCs dependent on FOXA1-GR and sensitive to GR antagonism.

Refined read-out: The hUHRF1 Tandem-Tudor domain prefers binding to histone H3 tails containing K4me1 in the context of H3K9me2/3

UHRF1 is an essential chromatin protein required for DNA methylation maintenance, mammalian development, and gene regulation. We investigated the Tandem-Tudor domain (TTD) of human UHRF1 that is known to bind H3K9me2/3 histones and is a major driver of UHRF1 localization in cells. We verified binding to H3K9me2/3 but unexpectedly discovered stronger binding to H3 peptides and mononucleosomes containing K9me2/3 with additional K4me1. We investigated the combined binding of TTD to H3K4me1-K9me2/3 versus H3K9me2/3 alone, engineered mutants with specific and differential changes of binding, and discovered a novel read-out mechanism for H3K4me1 in an H3K9me2/3 context that is based on the interaction of R207 with the H3K4me1 methyl group and on counting the H-bond capacity of H3K4. Individual TTD mutants showed up to a 10,000-fold preference for the double-modified peptides, suggesting that after a conformational change, WT TTD could exhibit similar effects. The frequent appearance of H3K4me1-K9me2 regions in human chromatin demonstrated in our TTD chromatin pull-down and ChIP-western blot data suggests that it has specific biological roles. Chromatin pull-down of TTD from HepG2 cells and full-length murine UHRF1 ChIP-seq data correlate with H3K4me1 profiles indicating that the H3K4me1-K9me2/3 interaction of TTD influences chromatin binding of full-length UHRF1. We demonstrate the H3K4me1-K9me2/3 specific binding of UHRF1-TTD to enhancers and promoters of cell-type-specific genes at the flanks of cell-type-specific transcription factor binding sites, and provided evidence supporting an H3K4me1-K9me2/3 dependent and TTD mediated downregulation of these genes by UHRF1. All these findings illustrate the important physiological function of UHRF1-TTD binding to H3K4me1-K9me2/3 double marks in a cellular context.

KMT2C knockout generates ASD-like behaviors in mice

Neurodevelopmental disorders have been associated with genetic mutations that affect cellular function, including chromatin regulation and epigenetic modifications. Recent studies in humans have identified mutations in KMT2C, an enzyme responsible for modifying histone tails and depositing H3K4me1 and H3K4me3, as being associated with Kleefstra syndrome 2 and autism spectrum disorder (ASD). However, the precise role of KMT2C mutations in brain disorders remains poorly understood. Here we employed CRISPR/Cas9 gene editing to analyze the effects of KMT2C brain specific knockout on animal behavior. Knocking out KMT2C expression in cortical neurons and the mouse brain resulted in decreased KMT2C levels. Importantly, KMT2C brain specific knockout animals exhibited repetitive behaviors, social deficits, and intellectual disability resembling ASD. Our findings shed light on the involvement of KMT2C in neurodevelopmental processes and establish a valuable model for elucidating the cellular and molecular mechanisms underlying KMT2C mutations and their relationship to Kleefstra syndrome 2 and ASD.

FOXA1 Reprogramming Dictates Retinoid X Receptor Response in ESR1-Mutant Breast Cancer

Estrogen receptor alpha (ER/ESR1) mutations occur in 30% to 40% of endocrine resistant ER-positive (ER+) breast cancer. Forkhead box A1 (FOXA1) is a key pioneer factor mediating ER-chromatin interactions and endocrine response in ER+ breast cancer, but its role in ESR1-mutant breast cancer remains unclear. Our previous FOXA1 chromatin immunoprecipitation sequencing (ChIP-seq) identified a large portion of redistributed binding sites in T47D genome-edited Y537S and D538G ESR1-mutant cells. Here, we further integrated FOXA1 genomic binding profile with the isogenic ER cistrome, accessible genome, and transcriptome data of T47D cell model. FOXA1 redistribution was significantly associated with transcriptomic alterations caused by ESR1 mutations. Furthermore, in ESR1-mutant cells, FOXA1-binding sites less frequently overlapped with ER, and differential gene expression was less associated with the canonical FOXA1-ER axis. Motif analysis revealed a unique enrichment of retinoid X receptor (RXR) motifs in FOXA1-binding sites of ESR1-mutant cells. Consistently, ESR1-mutant cells were more sensitive to growth stimulation with the RXR agonist LG268. The mutant-specific response was dependent on two RXR isoforms, RXR-α and RXR-β, with a stronger dependency on the latter. In addition, T3, the agonist of thyroid receptor (TR) also showed a similar growth-promoting effect in ESR1-mutant cells. Importantly, RXR antagonist HX531 blocked growth of ESR1-mutant cells and a patient-derived xenograft (PDX)-derived organoid with an ESR1 D538G mutation. Collectively, our data support the evidence for a stronger RXR response associated with FOXA1 reprograming in ESR1-mutant cells, suggesting development of therapeutic strategies targeting RXR pathways in breast tumors with ESR1 mutation.

IMPLICATIONS

It provides comprehensive characterization of the role of FOXA1 in ESR1-mutant breast cancer and potential therapeutic strategy through blocking RXR activation.

BRWD3 promotes KDM5 degradation to maintain H3K4 methylation levels

Histone modifications are critical for regulating chromatin structure and gene expression. Dysregulation of histone modifications likely contributes to disease states and cancer. Depletion of the chromatin-binding protein BRWD3, a known substrate-specificity factor of the Cul4-DDB1 E3 ubiquitin ligase complex, results in increased in H3K4me1 levels. The underlying mechanism linking BRWD3 and H3K4 methylation, however, has yet to be defined. Here, we show that depleting BRWD3 not only causes an increase in H3K4me1 levels, but also causes a decrease in H3K4me3 levels, indicating that BRWD3 influences H3K4 methylation more broadly. Using immunoprecipitation coupled to quantitative mass spectrometry, we identified an interaction between BRWD3 and the H3K4-specific demethylase 5 (KDM5/Lid), an enzyme that removes tri- and di- methyl marks from H3K4. Moreover, analysis of ChIP-seq data revealed that BRWD3 and KDM5 are significantly co- localized throughout the genome and that sites of H3K4me3 are highly enriched at BRWD3 binding sites. We show that BRWD3 promotes K48-linked polyubiquitination and degradation of KDM5 and that KDM5 degradation is dependent on both BRWD3 and Cul4. Critically, depleting KDM5 fully restores altered H3K4me3 levels and partially restores H3K4me1 levels upon BRWD3 depletion. Together, our results demonstrate that BRWD3 regulates KDM5 activity to balance H3K4 methylation levels.

Loss of Kmt2c in vivo leads to EMT, mitochondrial dysfunction and improved response to lapatinib in breast cancer

Deep sequencing of human tumours has uncovered a previously unappreciated role for epigenetic regulators in tumorigenesis. H3K4 methyltransferase KMT2C/MLL3 is mutated in several solid malignancies, including more than 10% of breast tumours. To study the tumour suppressor role of KMT2C in breast cancer, we generated mouse models of Erbb2/Neu, Myc or PIK3CA-driven tumorigenesis, in which the Kmt2c locus is knocked out specifically in the luminal lineage of mouse mammary glands using the Cre recombinase. Kmt2c knock out mice develop tumours earlier, irrespective of the oncogene, assigning a bona fide tumour suppressor role for KMT2C in mammary tumorigenesis. Loss of Kmt2c induces extensive epigenetic and transcriptional changes, which lead to increased ERK1/2 activity, extracellular matrix re-organization, epithelial-to-mesenchymal transition and mitochondrial dysfunction, the latter associated with increased reactive oxygen species production. Loss of Kmt2c renders the Erbb2/Neu-driven tumours more responsive to lapatinib. Publicly available clinical datasets revealed an association of low Kmt2c gene expression and better long-term outcome. Collectively, our findings solidify the role of KMT2C as a tumour suppressor in breast cancer and identify dependencies that could be therapeutically amenable.

KDM6A Loss Triggers an Epigenetic Switch That Disrupts Urothelial Differentiation and Drives Cell Proliferation in Bladder Cancer

Disruption of KDM6A, a histone lysine demethylase, is one of the most common somatic alternations in bladder cancer. Insights into how KDM6A mutations affect the epigenetic landscape to promote carcinogenesis could help reveal potential new treatment approaches. Here, we demonstrated that KDM6A loss triggers an epigenetic switch that disrupts urothelial differentiation and induces a neoplastic state characterized by increased cell proliferation. In bladder cancer cells with intact KDM6A, FOXA1 interacted with KDM6A to activate genes instructing urothelial differentiation. KDM6A-deficient cells displayed simultaneous loss of FOXA1 target binding and genome-wide redistribution of the bZIP transcription factor ATF3, which in turn repressed FOXA1-target genes and activated cell-cycle progression genes. Importantly, ATF3 depletion reversed the cell proliferation phenotype induced by KDM6A deficiency. These data establish that KDM6A loss engenders an epigenetic state that drives tumor growth in an ATF3-dependent manner, creating a potentially targetable molecular vulnerability.

SIGNIFICANCE

A gain-of-function epigenetic switch that disrupts differentiation is triggered by inactivating KDM6A mutations in bladder cancer and can serve as a potential target for novel therapies.

JMJ28 guides sequence-specific targeting of ATX1/2-containing COMPASS-like complex in Arabidopsis

Despite extensive investigations in mammals and yeasts, the importance and specificity of COMPASS-like complex, which catalyzes histone 3 lysine 4 methylation (H3K4me), are not fully understood in plants. Here, we report that JMJ28, a Jumonji C domain-containing protein in Arabidopsis, recognizes specific DNA motifs through a plant-specific WRC domain and acts as an interacting factor to guide the chromatin targeting of ATX1/2-containing COMPASS-like complex. JMJ28 associates with COMPASS-like complex in vivo via direct interaction with RBL. The DNA-binding activity of JMJ28 is essential for both the targeting specificity of ATX1/2-COMPASS and the deposition of H3K4me at specific loci but exhibit functional redundancy with alternative COMPASS-like complexes at other loci. Finally, we demonstrate that JMJ28 is a negative regulator of plant immunity. In summary, our findings reveal a plant-specific recruitment mechanism of COMPASS-like complex. These findings help to gain deeper insights into the regulatory mechanism of COMPASS-like complex in plants.

Transcription networks rewire gene repertoire to coordinate cellular reprograming in prostate cancer

Transcription factors (TFs) represent the most commonly deregulated DNA-binding class of proteins associated with multiple human cancers. They can act as transcriptional activators or repressors that rewire the cistrome, resulting in cellular reprogramming during cancer progression. Deregulation of TFs is associated with the onset and maintenance of various cancer types including prostate cancer. An emerging subset of TFs has been implicated in the regulation of multiple cancer hallmarks during tumorigenesis. Here, we discuss the role of key TFs which modulate transcriptional cicuitries involved in the development and progression of prostate cancer. We further highlight the role of TFs associated with key cancer hallmarks, including, chromatin remodeling, genome instability, DNA repair, invasion, and metastasis. We also discuss the pluripotent function of TFs in conferring lineage plasticity, that aids in disease progression to neuroendocrine prostate cancer. At the end, we summarize the current understanding and approaches employed for the therapeutic targeting of TFs and their cofactors in the clinical setups to prevent disease progression.

GRHL2-controlled gene expression networks in luminal breast cancer

Grainyhead like 2 (GRHL2) is an essential transcription factor for development and function of epithelial tissues. It has dual roles in cancer by supporting tumor growth while suppressing epithelial to mesenchymal transitions (EMT). GRHL2 cooperates with androgen and estrogen receptors (ER) to regulate gene expression. We explore genome wide GRHL2 binding sites conserved in three ER⍺/GRHL2 positive luminal breast cancer cell lines by ChIP-Seq. Interaction with the ER⍺/FOXA1/GATA3 complex is observed, however, only for a minor fraction of conserved GRHL2 peaks. We determine genome wide transcriptional dynamics in response to loss of GRHL2 by nascent RNA Bru-seq using an MCF7 conditional knockout model. Integration of ChIP- and Bru-seq pinpoints candidate direct GRHL2 target genes in luminal breast cancer. Multiple connections between GRHL2 and proliferation are uncovered, including transcriptional activation of ETS and E2F transcription factors. Among EMT-related genes, direct regulation of CLDN4 is corroborated but several targets identified in other cells (including CDH1 and ZEB1) are ruled out by both ChIP- and Bru-seq as being directly controlled by GRHL2 in luminal breast cancer cells. Gene clusters correlating positively (including known GRHL2 targets such as ErbB3, CLDN4/7) or negatively (including TGFB1 and TGFBR2) with GRHL2 in the MCF7 knockout model, display similar correlation with GRHL2 in ER positive as well as ER negative breast cancer patients. Altogether, this study uncovers gene sets regulated directly or indirectly by GRHL2 in luminal breast cancer, identifies novel GRHL2-regulated genes, and points to distinct GRHL2 regulation of EMT in luminal breast cancer cells. Video Abstract.

MLL3 loss drives metastasis by promoting a hybrid epithelial-mesenchymal transition state

Phenotypic plasticity associated with the hybrid epithelial-mesenchymal transition (EMT) is crucial to metastatic seeding and outgrowth. However, the mechanisms governing the hybrid EMT state remain poorly defined. Here we showed that deletion of the epigenetic regulator MLL3, a tumour suppressor frequently altered in human cancer, promoted the acquisition of hybrid EMT in breast cancer cells. Distinct from other EMT regulators that mediate only unidirectional changes, MLL3 loss enhanced responses to stimuli inducing EMT and mesenchymal-epithelial transition in epithelial and mesenchymal cells, respectively. Consequently, MLL3 loss greatly increased metastasis by enhancing metastatic colonization. Mechanistically, MLL3 loss led to increased IFNγ signalling, which contributed to the induction of hybrid EMT cells and enhanced metastatic capacity. Furthermore, BET inhibition effectively suppressed the growth of MLL3-mutant primary tumours and metastases. These results uncovered MLL3 mutation as a key driver of hybrid EMT and metastasis in breast cancer that could be targeted therapeutically.

FOXQ1 recruits the MLL complex to activate transcription of EMT and promote breast cancer metastasis

Aberrant expression of the Forkhead box transcription factor, FOXQ1, is a prevalent mechanism of epithelial-mesenchymal transition (EMT) and metastasis in multiple carcinoma types. However, it remains unknown how FOXQ1 regulates gene expression. Here, we report that FOXQ1 initiates EMT by recruiting the MLL/KMT2 histone methyltransferase complex as a transcriptional coactivator. We first establish that FOXQ1 promoter recognition precedes MLL complex assembly and histone-3 lysine-4 trimethylation within the promoter regions of critical genes in the EMT program. Mechanistically, we identify that the Forkhead box in FOXQ1 functions as a transactivation domain directly binding the MLL core complex subunit RbBP5 without interrupting FOXQ1 DNA binding activity. Moreover, genetic disruption of the FOXQ1-RbBP5 interaction or pharmacologic targeting of KMT2/MLL recruitment inhibits FOXQ1-dependent gene expression, EMT, and in vivo tumor progression. Our study suggests that targeting the FOXQ1-MLL epigenetic axis could be a promising strategy to combat triple-negative breast cancer metastatic progression.

GRHL2 Enhances Phosphorylated Estrogen Receptor (ER) Chromatin Binding and Regulates ER-Mediated Transcriptional Activation and Repression

Phosphorylation of estrogen receptor α (ER) at serine 118 (pS118-ER) is induced by estrogen and is the most abundant posttranslational mark associated with a transcriptionally active receptor. Cistromic analysis of pS118-ER from our group revealed enrichment of the GRHL2 motif near pS118-ER binding sites. In this study, we used cistromic and transcriptomic analyses to interrogate the relationship between GRHL2 and pS118-ER. We found that GRHL2 is bound to chromatin at pS118-ER/GRHL2 co-occupancy sites prior to ligand treatment, and GRHL2 binding is required for maximal pS118-ER recruitment. pS118-ER/GRHL2 co-occupancy sites were enriched at active enhancers marked by H3K27ac and H3K4me1, along with FOXA1 and p300, compared to sites where each factor binds independently. Transcriptomic analysis yielded four subsets of ER/GRHL2-coregulated genes revealing that GRHL2 can both enhance and antagonize E2-mediated ER transcriptional activity. Gene ontology analysis indicated that coregulated genes are involved in cell migration. Accordingly, knockdown of GRHL2, combined with estrogen treatment, resulted in increased cell migration but no change in proliferation. These results support a model in which GRHL2 binds to selected enhancers and facilitates pS118-ER recruitment to chromatin, which then results in differential activation and repression of genes that control estrogen-regulated ER-positive breast cancer cell migration.

Interpretable deep learning for chromatin-informed inference of transcriptional programs driven by somatic alterations across cancers

Cancer is a disease of gene dysregulation, where cells acquire somatic and epigenetic alterations that drive aberrant cellular signaling. These alterations adversely impact transcriptional programs and cause profound changes in gene expression. Interpreting somatic alterations within context-specific transcriptional programs will facilitate personalized therapeutic decisions but is a monumental task. Toward this goal, we develop a partially interpretable neural network model called Chromatin-informed Inference of Transcriptional Regulators Using Self-attention mechanism (CITRUS). CITRUS models the impact of somatic alterations on transcription factors and downstream transcriptional programs. Our approach employs a self-attention mechanism to model the contextual impact of somatic alterations. Furthermore, CITRUS uses a layer of hidden nodes to explicitly represent the state of transcription factors (TFs) to learn the relationships between TFs and their target genes based on TF binding motifs in the open chromatin regions of tumor samples. We apply CITRUS to genomic, transcriptomic, and epigenomic data from 17 cancer types profiled by The Cancer Genome Atlas. CITRUS predicts patient-specific TF activities and reveals transcriptional program variations between and within tumor types. We show that CITRUS yields biological insights into delineating TFs associated with somatic alterations in individual tumors. Thus, CITRUS is a promising tool for precision oncology.

Chromatin Variants Reveal the Genetic Determinants of Oncogenesis in Breast Cancer

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

Histone modification and histone modification-targeted anti-cancer drugs in breast cancer: Fundamentals and beyond

Dysregulated epigenetic enzymes and resultant abnormal epigenetic modifications (EMs) have been suggested to be closely related to tumor occurrence and progression. Histone modifications (HMs) can assist in maintaining genome stability, DNA repair, transcription, and chromatin modulation within breast cancer (BC) cells. In addition, HMs are reversible, dynamic processes involving the associations of different enzymes with molecular compounds. Abnormal HMs (e.g. histone methylation and histone acetylation) have been identified to be tightly related to BC occurrence and development, even though their underlying mechanisms remain largely unclear. EMs are reversible, and as a result, epigenetic enzymes have aroused wide attention as anti-tumor therapeutic targets. At present, treatments to restore aberrant EMs within BC cells have entered preclinical or clinical trials. In addition, no existing studies have comprehensively analyzed aberrant HMs within BC cells; in addition, HM-targeting BC treatments remain to be further investigated. Histone and non-histone protein methylation is becoming an attractive anti-tumor epigenetic therapeutic target; such methylation-related enzyme inhibitors are under development at present. Consequently, the present work focuses on summarizing relevant studies on HMs related to BC and the possible mechanisms associated with abnormal HMs. Additionally, we also aim to analyze existing therapeutic agents together with those drugs approved and tested through pre-clinical and clinical trials, to assess their roles in HMs. Moreover, epi-drugs that target HMT inhibitors and HDAC inhibitors should be tested in preclinical and clinical studies for the treatment of BC. Epi-drugs that target histone methylation (HMT inhibitors) and histone acetylation (HDAC inhibitors) have now entered clinical trials or are approved by the US Food and Drug Administration (FDA). Therefore, the review covers the difficulties in applying HM-targeting treatments in clinics and proposes feasible approaches for overcoming such difficulties and promoting their use in treating BC cases.

Insights into the multi-faceted role of Pioneer transcription factors in glioma formation and progression with targeting options

Pioneer transcription factors (TFs) present an important subtype of transcription factors which are vital for cell programming during embryonic development and cellular memory during mitotic growth, as well as cell fate reprogramming. Pioneer TFs can engage specific target binding sites on nucleosomal DNA to attract chromatin remodeling complexes, cofactors, and other transcription factors, ultimately controlling gene expression by shaping locally the epigenome. The priority of binding that they exhibit in contrast to other transcription factors and their involvement in crucial events regarding cell fate, has implicated their aberrant function in the pathogenesis of several disorders including carcinogenesis. Emerging experimental data indicate that certain Pioneer TFs are highly implicated in gliomas development, in neoplastic cell proliferation, angiogenic processes, resistance to therapy, and patient survival. Herein, we describe the main structural characteristics and functional mechanisms of pioneer TFs, focusing on their central role in the pathogenesis and progression of gliomas. We further highlight the current treatment options ranging from natural agents (oleanolic acid) to a variety of chemical compounds (APR-246, COTI-2) and discuss potential delivery systems, including nanoparticles, viral vectors, and intracellular protein delivery techniques.

Structural insight into ASH1L PHD finger recognizing methylated histone H3K4 and promoting cell growth in prostate cancer

ASH1L is a member of the Trithorax-group protein and acts as a histone methyltransferase for gene transcription activation. It is known that ASH1L modulates H3K4me3 and H3K36me2/3 at its gene targets, but its specific mechanism of histone recognition is insufficiently understood. In this study, we found that the ASH1L plant homeodomain (PHD) finger interacts with mono-, di-, and trimethylated states of H3K4 peptides with comparable affinities, indicating that ASH1L PHD non-selectively binds to all three methylation states of H3K4. We solved nuclear magnetic resonance structures picturing the ASH1L PHD finger binding to the dimethylated H3K4 peptide and found that a narrow binding groove and residue composition in the methylated-lysine binding pocket restricts the necessary interaction with the dimethyl-ammonium moiety of K4. In addition, we found that the ASH1L protein is overexpressed in castrate-resistant prostate cancer (PCa) PC3 and DU145 cells in comparison to PCa LNCaP cells. The knockdown of ASH1L modulated gene expression and cellular pathways involved in apoptosis and cell cycle regulation and consequently induced cell cycle arrest, cell apoptosis, and reduced colony-forming abilities in PC3 and DU145 cells. The overexpression of the C-terminal core of ASH1L but not the PHD deletion mutant increased the overall H3K36me2 level but had no effect on the H3K4me2/3 level. Overall, our study identifies the ASH1L PHD finger as the first native reader that non-selectively recognizes the three methylation states of H3K4. Additionally, ASH1L is required for the deregulation of cell cycle and survival in PCas.

Circadian clock function does not require the histone methyltransferase MLL3

The circadian clock controls the physiological function of tissues through the regulation of thousands of genes in a cell-type-specific manner. The core cellular circadian clock is a transcription-translation negative feedback loop, which can recruit epigenetic regulators to facilitate temporal control of gene expression. Histone methyltransferase, mixed lineage leukemia gene 3 (MLL3) was reported to be required for the maintenance of circadian oscillations in cultured cells. Here, we test the role of MLL3 in circadian organization in whole animals. Using mice expressing catalytically inactive MLL3, we show that MLL3 methyltransferase activity is in fact not required for circadian oscillations in vitro in a range of tissues, nor for the maintenance of circadian behavioral rhythms in vivo. In contrast to a previous report, loss of MLL3-dependent methylation did not affect the global levels of H3K4 methylation in liver, indicating substantial compensation from other methyltransferases. Furthermore, we found little evidence of genomic repositioning of H3K4me3 marks. We did, however, observe repositioning of H3K4me1 from intronic regions to intergenic regions and gene promoters; however, there were no changes in H3K4me1 mark abundance around core circadian clock genes. Output functions of the circadian clock, such as control of inflammation, were largely intact in MLL3-methyltransferase-deficient mice, although some gene-specific changes were observed, with sexually dimorphic loss of circadian regulation of specific cytokines. Taken together, these observations indicate that MLL3-directed histone methylation is not essential for core circadian clock function; however, it may influence the inflammatory response.

Genomic Markers of CDK 4/6 Inhibitor Resistance in Hormone Receptor Positive Metastatic Breast Cancer

Cyclin-dependent kinase 4/6 inhibitors are the standard of care for hormone receptor-positive metastatic breast cancer. This retrospective study reports on genomic biomarkers of CDK 4/6i resistance utilizing genomic data acquired through routine clinical practice. Patients with HR+ MBC treated with palbociclib, ribociclib, or abemaciclib and antiestrogen therapy were identified. Patients were grouped into early (<6 months); intermediate (6−24 months for 0−1 lines; 6−9 months for ≥2 lines); or late progressors (>24 months for 0−1 lines; >9 months PFS for ≥2 lines). NGS and RNA sequencing data were analyzed in association with PFS, and survival analysis was stratified by prior lines of chemotherapy. A total of 795 patients with HR+ MBC treated with CDK 4/6i were identified. Of these, 144 (18%) patients had genomic data and 29 (3.6%) had RNA data. Among the 109 patients who received CDK4/6i as 1st- or 2nd-line therapy, 17 genes showed associations with PFS (p-value ≤ 0.15 and HR ≥ 1.5 or HR < 0.5). Whole transcriptome RNAseq was analyzed for 24/109 (22%) patients with 0−1 prior lines of therapy and 56 genes associated with PFS (HR ≥ 4 or HR ≤ 0.25 and FDR ≤ 0.15). In this retrospective analysis, genomic biomarkers including FGFR1 amplification, PTEN loss, and DNA repair pathway gene mutations showed significant associations with shorter PFS for patients receiving CDK4/6 inhibitor therapy.

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.

Single-Cell Chromatin and Gene-Regulatory Dynamics of Mouse Nephron Progenitors

BACKGROUND

We reasoned that unraveling the dynamic changes in accessibility of genomic regulatory elements and gene expression at single-cell resolution will inform the basic mechanisms of nephrogenesis.

METHODS

We performed single-cell ATAC-seq and RNA-seq both individually (Singleomes; Six2GFPcells) and jointly in the same cells (Multiomes; kidneys) to generate integrated chromatin and transcriptional maps in mouse embryonic and neonatal nephron progenitor cells (NPCs).

RESULTS

WWe demonstrate that singleomes and multiomes are comparable in assigning most cell states, identification of new cell type markers, and defining the transcription factors driving cell identity. However, multiomes are more precise in defining the progenitor population. Multiomes identified a "pioneer" bHLH/Fox motif signature in NPCs. Moreover, we identified a subset of Fox factors exhibiting high chromatin activity in podocytes. One of these Fox factors, Foxp1, is important for nephrogenesis. Key nephrogenic factors are distinguished by strong correlation between linked generegulatory elements and gene expression.

CONCLUSION

Mapping the regulatory landscape at single-cell resolution informs the regulatory hierarchy of nephrogenesis. Paired single-cell epigenomes and transcriptomes of nephron progenitors should provide a foundation to understand prenatal programming, regeneration following injury, and ex vivo nephrogenesis.

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.

Pioneer factors as master regulators of the epigenome and cell fate

Pioneer factors are transcription factors with the unique ability to initiate opening of closed chromatin. The stability of cell identity relies on robust mechanisms that maintain the epigenome and chromatin accessibility to transcription factors. Pioneer factors counter these mechanisms to implement new cell fates through binding of DNA target sites in closed chromatin and introduction of active-chromatin histone modifications, primarily at enhancers. As master regulators of enhancer activation, pioneers are thus crucial for the implementation of correct cell fate decisions in development, and as such, they hold tremendous potential for therapy through cellular reprogramming. The power of pioneer factors to reshape the epigenome also presents an Achilles heel, as their misexpression has major pathological consequences, such as in cancer. In this Review, we discuss the emerging mechanisms of pioneer factor functions and their roles in cell fate specification, cellular reprogramming and cancer.

Androgen Receptor-Mediated Transcription in Prostate Cancer

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

HNF4A defines tissue-specific circadian rhythms by beaconing BMAL1::CLOCK chromatin binding and shaping the rhythmic chromatin landscape

Transcription modulated by the circadian clock is diverse across cell types, underlying circadian control of peripheral metabolism and its observed perturbation in human diseases. We report that knockout of the lineage-specifying Hnf4a gene in mouse liver causes associated reductions in the genome-wide distribution of core clock component BMAL1 and accessible chromatin marks (H3K4me1 and H3K27ac). Ectopically expressing HNF4A remodels chromatin landscape and nucleates distinct tissue-specific BMAL1 chromatin binding events, predominantly in enhancer regions. Circadian rhythms are disturbed in Hnf4a knockout liver and HNF4A-MODY diabetic model cells. Additionally, the epigenetic state and accessibility of the liver genome dynamically change throughout the day, synchronized with chromatin occupancy of HNF4A and clustered expression of circadian outputs. Lastly, Bmal1 knockout attenuates HNF4A genome-wide binding in the liver, likely due to downregulated Hnf4a transcription. Our results may provide a general mechanism for establishing circadian rhythm heterogeneity during development and disease progression, governed by chromatin structure.

FOXA1: A Pioneer of Nuclear Receptor Action in Breast Cancer

The pioneering function of FOXA1 establishes estrogen-responsive transcriptomes in luminal breast cancer. Dysregulated FOXA1 chromatin occupancy through focal amplification, mutation, or cofactor recruitment modulates estrogen receptor (ER) transcriptional programs and drives endocrine-resistant disease. However, ER is not the sole nuclear receptor (NR) expressed in breast cancers, nor is it the only NR for which FOXA1 serves as a licensing factor. Receptors for androgens, glucocorticoids, and progesterone are also found in the majority of breast cancers, and their functions are also impacted by FOXA1. These NRs interface with ER transcriptional programs and, depending on their activation level, can reprogram FOXA1-ER cistromes. Thus, NR interplay contributes to endocrine therapy response and resistance and may provide a vulnerability for future therapeutic benefit in patients. Herein, we review what is known regarding FOXA1 regulation of NR function in breast cancer in the context of cell identity, endocrine resistance, and NR crosstalk in breast cancer progression and treatment.

MLL3 is a de novo cause of endocrine therapy resistance

Background

Cancer resequencing studies have revealed epigenetic enzymes as common targets for recurrent mutations. The monomethyltransferase MLL3 is among the most recurrently mutated enzymes in ER+ breast cancer. The H3K4me1 marks created by MLL3 can define enhancers. In ER+ breast cancer, ERα genome‐binding sites are primarily distal enhancers. Thus, we hypothesize that mutation of MLL3 will alter the genomic binding and transcriptional regulatory activity of ERα.

Methods

We investigated the genomic consequences of knocking down MLL3 in an MLL3/PIK3CA WT ER+ breast cancer cell line.

Results

Loss of MLL3 led to a large loss of H3K4me1 across the genome, and a shift in genomic location of ERα‐binding sites, which was accompanied by a re‐organization of the breast cancer transcriptome. Gene set enrichment analyses of ERα‐binding sites in MLL3 KD identified endocrine therapy resistance terms, and we showed that MLL3 KD cells are resistant to tamoxifen and fulvestrant. Many differentially expressed genes are controlled by the small collection of new locations of H3K4me1 deposition and ERα binding, suggesting that loss of functional MLL3 leads to new transcriptional regulation of essential genes. Motif analysis of RNA‐seq and ChIP‐seq data highlighted SP1 as a critical transcription factor in the MLL3 KD cells. Differentially expressed genes that display a loss of ERα binding upon MLL3 KD also harbor increased SP1 binding.

Conclusions

Our data show that a decrease in functional MLL3 leads to endocrine therapy resistance. This highlights the importance of genotyping patient tumor samples for MLL3 mutation upon initial resection, prior to deciding upon treatment plans.

Epigenetic remodelling of enhancers in response to estrogen deprivation and re-stimulation

Estrogen hormones are implicated in a majority of breast cancers and estrogen receptor alpha (ER), the main nuclear factor mediating estrogen signaling, orchestrates a complex molecular circuitry that is not yet fully elucidated. Here, we investigated genome-wide DNA methylation, histone acetylation and transcription after estradiol (E2) deprivation and re-stimulation to better characterize the ability of ER to coordinate gene regulation. We found that E2 deprivation mostly resulted in DNA hypermethylation and histone deacetylation in enhancers. Transcriptome analysis revealed that E2 deprivation leads to a global down-regulation in gene expression, and more specifically of TET2 demethylase that may be involved in the DNA hypermethylation following short-term E2 deprivation. Further enrichment analysis of transcription factor (TF) binding and motif occurrence highlights the importance of ER connection mainly with two partner TF families, AP-1 and FOX. These interactions take place in the proximity of E2 deprivation-mediated differentially methylated and histone acetylated enhancers. Finally, while most deprivation-dependent epigenetic changes were reversed following E2 re-stimulation, DNA hypermethylation and H3K27 deacetylation at certain enhancers were partially retained. Overall, these results show that inactivation of ER mediates rapid and mostly reversible epigenetic changes at enhancers, and bring new insight into early events, which may ultimately lead to endocrine resistance.

Single allele loss-of-function mutations select and sculpt conditional cooperative networks in breast cancer

The most common events in breast cancer (BC) involve chromosome arm losses and gains. Here we describe identification of 1089 gene-centric common insertion sites (gCIS) from transposon-based screens in 8 mouse models of BC. Some gCIS are driver-specific, others driver non-specific, and still others associated with tumor histology. Processes affected by driver-specific and histology-specific mutations include well-known cancer pathways. Driver non-specific gCIS target the Mediator complex, Ca++ signaling, Cyclin D turnover, RNA-metabolism among other processes. Most gCIS show single allele disruption and many map to genomic regions showing high-frequency hemizygous loss in human BC. Two gCIS, Nf1 and Trps1, show synthetic haploinsufficient tumor suppressor activity. Many gCIS act on the same pathway responsible for tumor initiation, thereby selecting and sculpting just enough and just right signaling. These data highlight ~1000 genes with predicted conditional haploinsufficient tumor suppressor function and the potential to promote chromosome arm loss in BC.

Integrated Chromatin Accessibility and Transcriptome Landscapes of Doxorubicin-Resistant Breast Cancer Cells

BACKGROUND

Doxorubicin is one of the most effective chemotherapeutic drugs for breast cancer while its common drug resistance leads to poor patient prognosis and survival. Growing evidence indicate dynamically reorganized chromatin allows rapid access of the gene regulatory machinery to open genomic regions facilitating subsequent gene expression through direct transcription factor (TF) activation and regulatory element binding.

METHODS

To better understand the regulatory network underlying doxorubicin resistance in breast cancer cells, we explored the systematic alterations of chromatin accessibility and gene expression by the assay for transposase-accessible chromatin using sequencing (ATAC-seq) in combination with RNA sequencing, followed by integrative analysis to identify potential regulators and their targets associated with differentially accessible regions (DARs) in doxorubicin-resistant MCF7 (MCF7-DR) cells.

RESULTS

A total of 3,963 differentially expressed genes (DEGs) related to doxorubicin resistance were identified, including dramatically up-regulated MT1E, GSTP1, LDHB, significantly down-regulated TFF1, UBB, DSCAM-AS1, and histone-modifying enzyme coding genes HDAC2, EZH2, PRMT5, etc. By integrating with transcriptomic datasets, we identified 18,228 DARs in MCF7-DR cells compared to control, which were positively correlated with their nearest DEGs (r = 0.6). There were 11,686 increased chromatin-accessible regions, which were enriched in up-regulated genes related to diverse KEGG pathways, such as the cell cycle, regulation of actin cytoskeleton, signaling pathways of MAPK, PI3K/Akt and Hippo, which play essential roles in regulating cell apoptosis, proliferation, metabolism, and inflammatory responses. The 6,542 decreased chromatin-accessible regions were identified for the declined doxorubicin-associated biological processes, for instance, endocrine and insulin resistance, central carbon metabolism, signaling pathways of TGF-beta and P53. Combining data from TCGA, analyses of the DAR sequences associated with the DNA-binding motifs of significantly enriched TF families including AP-1, TEAD and FOX, indicated that the loss-function of FOXA1 might play a critical role in doxorubicin-resistant breast cancer cells (DOX-R BCCs).

CONCLUSION

These data exhibit the non-genetic landscape of chromatin accessibility and transcript levels in the DOX-R BCCs, and provide clear insights and resources for the detection of critical TFs and potential cis-regulatory elements-based putative therapeutic targets.

Loss of KMT2C reprograms the epigenomic landscape in hPSCs resulting in NODAL overexpression and a failure of hemogenic endothelium specification

Germline or somatic variation in the family of KMT2 lysine methyltransferases have been associated with a variety of congenital disorders and cancers. Notably, KMT2A-fusions are prevalent in 70% of infant leukaemias but fail to phenocopy short latency leukaemogenesis in mammalian models, suggesting additional factors are necessary for transformation. Given the lack of additional somatic mutation, the role of epigenetic regulation in cell specification, and our prior results of germline KMT2C variation in infant leukaemia patients, we hypothesized that germline dysfunction of KMT2C altered haematopoietic specification. In isogenic KMT2C KO hPSCs, we found genome-wide differences in histone modifications at active and poised enhancers, leading to gene expression profiles akin to mesendoderm rather than mesoderm highlighted by a significant increase in NODAL expression and WNT inhibition, ultimately resulting in a lack of in vitro hemogenic endothelium specification. These unbiased multi-omic results provide new evidence for germline mechanisms increasing risk of early leukaemogenesis.

KMT2C is a potential biomarker of prognosis and chemotherapy sensitivity in breast cancer

PURPOSE

Epigenetic regulation plays critical roles in cancer progression, and high-frequency mutations or expression variations in epigenetic regulators have been frequently observed in tumorigenesis, serving as biomarkers and targets for cancer therapy. Here, we aimed to explore the function of epigenetic regulators in breast cancer.

METHODS

The mutational landscape of epigenetic regulators in breast cancer samples was investigated based on datasets from the Cancer Genome Atlas. The Kaplan-Meier method was used for survival analysis. RNA sequencing (RNA-seq) in MCF-7 cells transfected with control siRNA or KMT2C siRNA was performed. Quantitative reverse transcription-PCR and chromatin immunoprecipitation were used to validate the RNA-seq results.

RESULTS

Among the 450 epigenetic regulators, KMT2C was frequently mutated in breast cancer samples. The tumor mutational burden (TMB) was elevated in breast cancer samples with KMT2C mutations or low KMT2C mRNA levels compared to their counterparts with wild-type KMT2C or high KMT2C mRNA levels. Somatic mutation and low expression of KMT2C were independently correlated with the poor overall survival (OS) and disease-free survival (DFS) of the breast cancer samples, respectively. RNA-seq analysis combined with chromatin immunoprecipitation and qRT-PCR assays revealed that the depletion of KMT2C remarkably affected the expression of DNA damage repair-related genes. More importantly, the low expression of KMT2C was related to breast cancer cell sensitivity to chemotherapy and longer OS of breast cancer patients who underwent chemotherapy.

CONCLUSION

We conclude that KMT2C could serve as a potential biomarker of prognosis and chemotherapy sensitivity by affecting the DNA damage repair-related genes of breast cancer.

Multiomics metabolic and epigenetics regulatory network in cancer: A systems biology perspective

Genetic, epigenetic, and metabolic alterations are all hallmarks of cancer. However, the epigenome and metabolome are both highly complex and dynamic biological networks in vivo. The interplay between the epigenome and metabolome contributes to a biological system that is responsive to the tumor microenvironment and possesses a wealth of unknown biomarkers and targets of cancer therapy. From this perspective, we first review the state of high-throughput biological data acquisition (i.e. multiomics data) and analysis (i.e. computational tools) and then propose a conceptual in silico metabolic and epigenetic regulatory network (MER-Net) that is based on these current high-throughput methods. The conceptual MER-Net is aimed at linking metabolomic and epigenomic networks through observation of biological processes, omics data acquisition, analysis of network information, and integration with validated database knowledge. Thus, MER-Net could be used to reveal new potential biomarkers and therapeutic targets using deep learning models to integrate and analyze large multiomics networks. We propose that MER-Net can serve as a tool to guide integrated metabolomics and epigenomics research or can be modified to answer other complex biological and clinical questions using multiomics data.

The MLL3/4 H3K4 methyltransferase complex in establishing an active enhancer landscape

Enhancers are cis-regulatory elements that play essential roles in tissue-specific gene expression during development. Enhancer function in the expression of developmental genes requires precise regulation, while deregulation of enhancer function could be the main cause of tissue-specific cancer development. MLL3/KMT2C and MLL4/KMT2D are two paralogous histone modifiers that belong to the SET1/MLL (also named COMPASS) family of lysine methyltransferases and play critical roles in enhancer-regulated gene activation. Importantly, large-scale DNA sequencing studies have revealed that they are amongst the most frequently mutated genes associated with human cancers. MLL3 and MLL4 form identical multi-protein complexes for modifying mono-methylation of histone H3 lysine 4 (H3K4) at enhancers, which together with the p300/CBP-mediated H3K27 acetylation can generate an active enhancer landscape for long-range target gene activation. Recent studies have provided a better understanding of the possible mechanisms underlying the roles of MLL3/MLL4 complexes in enhancer regulation. Moreover, accumulating studies offer new insights into our knowledge of the potential role of MLL3/MLL4 in cancer development. In this review, we summarize recent evidence on the molecular mechanisms of MLL3/MLL4 in the regulation of active enhancer landscape and long-range gene expression, and discuss their clinical implications in human cancers.