Loss of Kmt2c or Kmt2d drives brain metastasis via KDM6A-dependent upregulation of MMP3.

KMT2C and KMT2D, encoding histone H3 lysine 4 methyltransferases, are among the most commonly mutated genes in triple-negative breast cancer (TNBC). However, how these mutations may shape epigenomic and transcriptomic landscapes to promote tumorigenesis is largely unknown. Here we describe that deletion of Kmt2c or Kmt2d in non-metastatic murine models of TNBC drives metastasis, especially to the brain. Global chromatin profiling and chromatin immunoprecipitation followed by sequencing revealed altered H3K4me1, H3K27ac and H3K27me3 chromatin marks in knockout cells and demonstrated enhanced binding of the H3K27me3 lysine demethylase KDM6A, which significantly correlated with gene expression. We identified Mmp3 as being commonly upregulated via epigenetic mechanisms in both knockout models. Consistent with these findings, samples from patients with KMT2C-mutant TNBC have higher MMP3 levels. Downregulation or pharmacological inhibition of KDM6A diminished Mmp3 upregulation induced by the loss of histone-lysine N-methyltransferase 2 (KMT2) and prevented brain metastasis similar to direct downregulation of Mmp3. Taken together, we identified the KDM6A-matrix metalloproteinase 3 axis as a key mediator of KMT2C/D loss-driven metastasis in TNBC.

Transcription factor dynamics, oscillation, and functions in human enteroendocrine cell differentiation.

Enteroendocrine cells (EECs) secrete serotonin (enterochromaffin [EC] cells) or specific peptide hormones (non-EC cells) that serve vital metabolic functions. The basis for terminal EEC diversity remains obscure. By forcing activity of the transcription factor (TF) NEUROG3 in 2D cultures of human intestinal stem cells, we replicated physiologic EEC differentiation and examined transcriptional and cis-regulatory dynamics that culminate in discrete cell types. Abundant EEC precursors expressed stage-specific genes and TFs. Before expressing pre-terminal NEUROD1, post-mitotic precursors oscillated between transcriptionally distinct ASCL1+ and HES6hi cell states. Loss of either factor accelerated EEC differentiation substantially and disrupted EEC individuality; ASCL1 or NEUROD1 deficiency had opposing consequences on EC and non-EC cell features. These TFs mainly bind cis-elements that are accessible in undifferentiated stem cells, and they tailor subsequent expression of TF combinations that underlie discrete EEC identities. Thus, early TF oscillations retard EEC maturation to enable accurate diversity within a medically important cell lineage.

Transcription factor dynamics, oscillation, and functions in human enteroendocrine cell differentiation.

Enteroendocrine cells (EECs), which secrete serotonin (enterochromaffin cells, EC) or a dominant peptide hormone, serve vital physiologic functions. As with any adult human lineage, the basis for terminal cell diversity remains obscure. We replicated human EEC differentiation in vitro , mapped transcriptional and chromatin dynamics that culminate in discrete cell types, and studied abundant EEC precursors expressing selected transcription factors (TFs) and gene programs. Before expressing the pre-terminal factor NEUROD1, non-replicating precursors oscillated between epigenetically similar but transcriptionally distinct ASCL1 + and HES6 hi cell states. Loss of either factor substantially accelerated EEC differentiation and disrupted EEC individuality; ASCL1 or NEUROD1 deficiency had opposing consequences on EC and hormone-producing cell features. Expressed late in EEC differentiation, the latter TFs mainly bind cis -elements that are accessible in undifferentiated stem cells and tailor the subsequent expression of TF combinations that specify EEC types. Thus, TF oscillations retard EEC maturation to enable accurate EEC diversification.

scATAnno: Automated Cell Type Annotation for single-cell ATAC Sequencing Data.

Recent advances in single-cell epigenomic techniques have created a growing demand for scATAC-seq analysis. One key analysis task is to determine cell type identity based on the epigenetic data. We introduce scATAnno, a python package designed to automatically annotate scATAC-seq data using large-scale scATAC-seq reference atlases. This workflow generates the reference atlases from publicly available datasets enabling accurate cell type annotation by integrating query data with reference atlases, without the use of scRNA-seq data. To enhance annotation accuracy, we have incorporated KNN-based and weighted distance-based uncertainty scores to effectively detect cell populations within the query data that are distinct from all cell types in the reference data. We compare and benchmark scATAnno against 7 other published approaches for cell annotation and show superior performance in multiple data sets and metrics. We showcase the utility of scATAnno across multiple datasets, including peripheral blood mononuclear cell (PBMC), Triple Negative Breast Cancer (TNBC), and basal cell carcinoma (BCC), and demonstrate that scATAnno accurately annotates cell types across conditions. Overall, scATAnno is a useful tool for scATAC-seq reference building and cell type annotation in scATAC-seq data and can aid in the interpretation of new scATAC-seq datasets in complex biological systems.

KDM6A epigenetically regulates subtype plasticity in small cell lung cancer.

Small cell lung cancer (SCLC) exists broadly in four molecular subtypes: ASCL1, NEUROD1, POU2F3 and Inflammatory. Initially, SCLC subtypes were thought to be mutually exclusive, but recent evidence shows intra-tumoural subtype heterogeneity and plasticity between subtypes. Here, using a CRISPR-based autochthonous SCLC genetically engineered mouse model to study the consequences of KDM6A/UTX inactivation, we show that KDM6A inactivation induced plasticity from ASCL1 to NEUROD1 resulting in SCLC tumours that express both ASCL1 and NEUROD1. Mechanistically, KDM6A normally maintains an active chromatin state that favours the ASCL1 subtype with its loss decreasing H3K4me1 and increasing H3K27me3 at enhancers of neuroendocrine genes leading to a cell state that is primed for ASCL1-to-NEUROD1 subtype switching. This work identifies KDM6A as an epigenetic regulator that controls ASCL1 to NEUROD1 subtype plasticity and provides an autochthonous SCLC genetically engineered mouse model to model ASCL1 and NEUROD1 subtype heterogeneity and plasticity, which is found in 35-40% of human SCLCs.

The Relationship between the Expression of GATA4 and GATA6 with the Clinical Characteristics and Prognosis of Resectable Pancreatic Adenocarcinoma.

GATA4 and GATA6 are transcription factors involved in the differentiation and development of PDAC. GATA6 expression is related to the classic molecular subtype, while its absence is related to the basal-like molecular subtype. The aim was to determine the clinical utility of IHC determination of GATA4 and GATA6 in a series of patients with resected PDAC. GATA4 and GATA6 expression was studied by IHC in TMA samples of normal tissue, PanIN, tumor tissue and lymph node metastases from a series of 89 patients with resected PDAC. Its relationship with clinicopathologic variables and the outcome was investigated. Seventy-two (81%) tumors were GATA6+ and 37 (42%) were GATA4+. While GATA4 expression was reduced during tumor progression, GATA6 expression remained highly conserved, except in lymph node metastases. All patients with early stages and well-differentiated tumors were GATA6+. The absence of GATA4 expression was related to smoking. Patients with GATA4+ or GATA6+ tumors had significantly lower Ca 19.9 levels. The expression of GATA4 and GATA6 was related to DFS, being more favorable in the GATA4+/GATA6+ group. The determination of the expression of GATA4 and GATA6 by IHC is feasible and provides complementary clinical and prognostic information that can help improve the stratification of patients with PDAC.

Universal NicE-Seq: A Simple and Quick Method for Accessible Chromatin Detection in Fixed Cells.

Genome-wide accessible chromatin sequencing and identification has enabled deciphering the epigenetic information encoded in chromatin, revealing accessible promoters, enhancers, nucleosome positioning, transcription factor occupancy, and other chromosomal protein binding. The starting biological materials are often fixed using formaldehyde crosslinking. Here, we describe accessible chromatin library preparation from low numbers of formaldehyde-crosslinked cells using a modified nick translation method, where a nicking enzyme nicks one strand of DNA and DNA polymerase incorporates biotin-conjugated dATP, dCTP, and methyl-dCTP. Once the DNA is labeled, it can be isolated for NGS library preparation. We termed this method as universal NicE-seq (nicking enzyme-assisted sequencing). We also demonstrate a single tube method that enables direct NGS library preparation from low cell numbers without DNA purification. Furthermore, we demonstrated universal NicE-seq on FFPE tissue section sample.

A Distinct Chromatin State Drives Therapeutic Resistance in Invasive Lobular Breast Cancer.

Most invasive lobular breast cancers (ILC) are of the luminal A subtype and are strongly hormone receptor-positive. Yet, ILC is relatively resistant to tamoxifen and associated with inferior long-term outcomes compared with invasive ductal cancers (IDC). In this study, we sought to gain mechanistic insights into these clinical findings that are not explained by the genetic landscape of ILC and to identify strategies to improve patient outcomes. A comprehensive analysis of the epigenome of ILC in preclinical models and clinical samples showed that, compared with IDC, ILC harbored a distinct chromatin state linked to gained recruitment of FOXA1, a lineage-defining pioneer transcription factor. This resulted in an ILC-unique FOXA1-estrogen receptor (ER) axis that promoted the transcription of genes associated with tumor progression and poor outcomes. The ILC-unique FOXA1-ER axis led to retained ER chromatin binding after tamoxifen treatment, which facilitated tamoxifen resistance while remaining strongly dependent on ER signaling. Mechanistically, gained FOXA1 binding was associated with the autoinduction of FOXA1 in ILC through an ILC-unique FOXA1 binding site. Targeted silencing of this regulatory site resulted in the disruption of the feed-forward loop and growth inhibition in ILC. In summary, ILC is characterized by a unique chromatin state and FOXA1-ER axis that is associated with tumor progression, offering a novel mechanism of tamoxifen resistance. These results underscore the importance of conducting clinical trials dedicated to patients with ILC in order to optimize treatments in this breast cancer subtype. SIGNIFICANCE: A unique FOXA1-ER axis in invasive lobular breast cancer promotes disease progression and tamoxifen resistance, highlighting a potential therapeutic avenue for clinical investigations dedicated to this disease. See related commentary by Blawski and Toska, p. 3668.

FOXR2 Is an Epigenetically Regulated Pan-Cancer Oncogene That Activates ETS Transcriptional Circuits.

Forkhead box R2 (FOXR2) is a forkhead transcription factor located on the X chromosome whose expression is normally restricted to the testis. In this study, we performed a pan-cancer analysis of FOXR2 activation across more than 10,000 adult and pediatric cancer samples and found FOXR2 to be aberrantly upregulated in 70% of all cancer types and 8% of all individual tumors. The majority of tumors (78%) aberrantly expressed FOXR2 through a previously undescribed epigenetic mechanism that involves hypomethylation of a novel promoter, which was functionally validated as necessary for FOXR2 expression and proliferation in FOXR2-expressing cancer cells. FOXR2 promoted tumor growth across multiple cancer lineages and co-opted ETS family transcription circuits across cancers. Taken together, this study identifies FOXR2 as a potent and ubiquitous oncogene that is epigenetically activated across the majority of human cancers. The identification of hijacking of ETS transcription circuits by FOXR2 extends the mechanisms known to active ETS transcription factors and highlights how transcription factor families cooperate to enhance tumorigenesis. SIGNIFICANCE: This work identifies a novel promoter that drives aberrant FOXR2 expression and delineates FOXR2 as a pan-cancer oncogene that specifically activates ETS transcriptional circuits across human cancers. See related commentary by Liu and Northcott, p. 2977.

MYC drives aggressive prostate cancer by disrupting transcriptional pause release at androgen receptor targets.

c-MYC (MYC) is a major driver of prostate cancer tumorigenesis and progression. Although MYC is overexpressed in both early and metastatic disease and associated with poor survival, its impact on prostate transcriptional reprogramming remains elusive. We demonstrate that MYC overexpression significantly diminishes the androgen receptor (AR) transcriptional program (the set of genes directly targeted by the AR protein) in luminal prostate cells without altering AR expression. Analyses of clinical specimens reveal that concurrent low AR and high MYC transcriptional programs accelerate prostate cancer progression toward a metastatic, castration-resistant disease. Data integration of single-cell transcriptomics together with ChIP-seq uncover an increase in RNA polymerase II (Pol II) promoter-proximal pausing at AR-dependent genes following MYC overexpression without an accompanying deactivation of AR-bound enhancers. Altogether, our findings suggest that MYC overexpression antagonizes the canonical AR transcriptional program and contributes to prostate tumor initiation and progression by disrupting transcriptional pause release at AR-regulated genes.

The nuclear receptor THRB facilitates differentiation of human PSCs into more mature hepatocytes.

To understand the mechanisms regulating the in vitro maturation of hPSC-derived hepatocytes, we developed a 3D differentiation system and compared gene regulatory elements in human primary hepatocytes with those in hPSC-hepatocytes that were differentiated in 2D or 3D conditions by RNA-seq, ATAC-seq, and H3K27Ac ChIP-seq. Regulome comparisons showed a reduced enrichment of thyroid receptor THRB motifs in accessible chromatin and active enhancers without a reduced transcription of THRB. The addition of thyroid hormone T3 increased the binding of THRB to the CYP3A4 proximal enhancer, restored the super-enhancer status and gene expression of NFIC, and reduced the expression of AFP. The resultant hPSC-hepatocytes showed gene expression, epigenetic status, and super-enhancer landscape closer to primary hepatocytes and activated regulatory regions including non-coding SNPs associated with liver-related diseases. Transplanting the hPSC-hepatocytes resulted in the engraftment of human hepatocytes into the mouse liver without disrupting normal liver histology. This work implicates the environmental factor-nuclear receptor axis in regulating the maturation of hPSC-hepatocytes.

Non-muscle-invasive micropapillary bladder cancer has a distinct lncRNA profile associated with unfavorable prognosis.

BACKGROUND: Molecular subtyping of bladder cancer has revealed luminal tumors generally have a more favourable prognosis. However, some aggressive forms of variant histology, including micropapillary, are often classified luminal. In previous work, we found long non-coding RNA (lncRNA) expression profiles could identify a subgroup of luminal bladder tumors with less aggressive biology and better outcomes. OBJECTIVE: In the present study, we aimed to investigate whether lncRNA expression profiles could identify high-grade T1 micropapillary bladder cancer with differential outcome. DESIGN, SETTING, AND PARTICIPANTS: LncRNAs were quantified from RNA-seq data from a HGT1 bladder cancer cohort that was enriched for primary micropapillary cases (15/84). Unsupervised consensus clustering of variant lncRNAs identified a three-cluster solution, which was further characterised using a panel of micropapillary-associated biomarkers, molecular subtypes, gene signatures, and survival analysis. A single-sample genomic signature was trained using lasso-penalized logistic regression to classify micropapillary-like gene-expression, as characterised by lncRNA clustering. The genomic classifier (GC) was tested on luminal tumors derived from the TCGA cohort (N = 202). OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Patient and tumor characteristics were compared between subgroups by using X2 tests and two-sided Wilcoxon rank-sum tests. Primary endpoints were overall, progression-free and high-grade recurrence-free survival, calculated as the date of high-grade T1 disease at TURBT till date of death from any cause, progression, or recurrence, respectively. Survival rates were estimated using weighted Kaplan-Meier (KM) curves. RESULTS AND LIMITATIONS: Primary micropapillary HGT1 showed decreased FGFR3, SHH, and p53 pathway activity relative to tumors with conventional urothelial carcinoma. Many bladder cancer-associated lncRNAs were downregulated in micropapillary tumors, including UCA1, LINC00152, and MALAT1. Unsupervised consensus clustering resulted in a lncRNA cluster 1 (LC1) with worse prognosis that was enriched for primary micropapillary histology and the Luminal Unstable (LumU) molecular subtype. Interestingly, LC1 appeared to better identify aggressive HGT1 disease, compared to stratifying outcomes using primary histologic characteristics. A signature trained to identify LC1 cases showed good performance in the testing cohort, identifying seven cases with significantly worse survival (p < 0.001). Limitations include the retrospective nature of the study and the lack of a validation cohort. CONCLUSIONS: Using the lncRNA transcriptome we identified a subgroup of aggressive HGT1 bladder cancer that was enriched with micropapillary histology. These data suggest that lncRNAs can facilitate the identification of aggressive micropapillary-like tumors, potentially improving patient management.

Nutritional Epigenetics in Cancer.

Alterations in the epigenome are well known to affect cancer development and progression. Epigenetics is highly influenced by the environment, including diet, which is a source of metabolic substrates that influence the synthesis of cofactors or substrates for chromatin and RNA modifying enzymes. In addition, plants are a common source of bioactives that can directly modify the activity of these enzymes. Here, we review and discuss the impact of diet on epigenetic mechanisms, including chromatin and RNA regulation, and its potential implications for cancer prevention and treatment.

High-Resolution ATAC-Seq Analysis of Frozen Clinical Tissues.

The ATAC-seq method enables the genome-wide analysis of accessible chromatin revealing transcriptionally active and poised regulatory elements. The ATAC-seq analysis of clinical specimens at a single-cell resolution reveals the cellular composition of the tissue contributing to the understanding of intra-tissue heterogeneity. Here we describe our method for nuclei isolation from frozen specimens with wide applicability across tissue types, producing nuclei suitable for a number of molecular profiling methods including ATAC-seq in bulk and at a single-cell resolution.

BCOR and BCORL1 Mutations Drive Epigenetic Reprogramming and Oncogenic Signaling by Unlinking PRC1.1 from Target Genes.

Polycomb repressive epigenetic complexes are recurrently dysregulated in cancer. Unlike polycomb repressive complex 2 (PRC2), the role of PRC1 in oncogenesis and therapy resistance is not well-defined. Here, we demonstrate that highly recurrent mutations of the PRC1 subunits BCOR and BCORL1 in leukemia disrupt assembly of a noncanonical PRC1.1 complex, thereby selectively unlinking the RING-PCGF enzymatic core from the chromatin-targeting auxiliary subcomplex. As a result, BCOR-mutated PRC1.1 is localized to chromatin but lacks repressive activity, leading to epigenetic reprogramming and transcriptional activation at target loci. We define a set of functional targets that drive aberrant oncogenic signaling programs in PRC1.1-mutated cells and primary patient samples. Activation of these PRC1.1 targets in BCOR-mutated cells confers acquired resistance to treatment while sensitizing to targeted kinase inhibition. Our study thus reveals a novel epigenetic mechanism that explains PRC1.1 tumor-suppressive activity and identifies a therapeutic strategy in PRC1.1-mutated cancer. SIGNIFICANCE: We demonstrate that BCOR and BCORL1 mutations in leukemia unlink PRC1.1 repressive function from target genes, resulting in epigenetic reprogramming and activation of aberrant cell signaling programs that mediate treatment resistance. Our study provides mechanistic insights into the pathogenesis of PRC1.1-mutated leukemia that inform novel therapeutic approaches. This article is highlighted in the In This Issue feature, p. 85.

An Enhancer-Driven Stem Cell-Like Program Mediated by SOX9 Blocks Intestinal Differentiation in Colorectal Cancer.

BACKGROUND AND AIMS: Genomic alterations that encourage stem cell activity and hinder proper maturation are central to the development of colorectal cancer (CRC). Key molecular mediators that promote these malignant properties require further elucidation to galvanize translational advances. We therefore aimed to characterize a key factor that blocks intestinal differentiation, define its transcriptional and epigenetic program, and provide preclinical evidence for therapeutic targeting in CRC. METHODS: Intestinal tissue from transgenic mice and patients were analyzed by means of histopathology and immunostaining. Human CRC cells and neoplastic murine organoids were genetically manipulated for functional studies. Gene expression profiling was obtained through RNA sequencing. Histone modifications and transcription factor binding were determined with the use of chromatin immunoprecipitation sequencing. RESULTS: We demonstrate that SRY-box transcription factor 9 (SOX9) promotes CRC by activating a stem cell-like program that hinders intestinal differentiation. Intestinal adenomas and colorectal adenocarcinomas from mouse models and patients demonstrate ectopic and elevated expression of SOX9. Functional experiments indicate a requirement for SOX9 in human CRC cell lines and engineered neoplastic organoids. Disrupting SOX9 activity impairs primary CRC tumor growth by inducing intestinal differentiation. By binding to genome wide enhancers, SOX9 directly activates genes associated with Paneth and stem cell activity, including prominin 1 (PROM1). SOX9 up-regulates PROM1 via a Wnt-responsive intronic enhancer. A pentaspan transmembrane protein, PROM1 uses its first intracellular domain to support stem cell signaling, at least in part through SOX9, reinforcing a PROM1-SOX9 positive feedback loop. CONCLUSIONS: These studies establish SOX9 as a central regulator of an enhancer-driven stem cell-like program and carry important implications for developing therapeutics directed at overcoming differentiation defects in CRC.

FGFR-inhibitor-mediated dismissal of SWI/SNF complexes from YAP-dependent enhancers induces adaptive therapeutic resistance.

How cancer cells adapt to evade the therapeutic effects of drugs targeting oncogenic drivers is poorly understood. Here we report an epigenetic mechanism leading to the adaptive resistance of triple-negative breast cancer (TNBC) to fibroblast growth factor receptor (FGFR) inhibitors. Prolonged FGFR inhibition suppresses the function of BRG1-dependent chromatin remodelling, leading to an epigenetic state that derepresses YAP-associated enhancers. These chromatin changes induce the expression of several amino acid transporters, resulting in increased intracellular levels of specific amino acids that reactivate mTORC1. Consistent with this mechanism, addition of mTORC1 or YAP inhibitors to FGFR blockade synergistically attenuated the growth of TNBC patient-derived xenograft models. Collectively, these findings reveal a feedback loop involving an epigenetic state transition and metabolic reprogramming that leads to adaptive therapeutic resistance and provides potential therapeutic strategies to overcome this mechanism of resistance.

Predicting master transcription factors from pan-cancer expression data.

Critical developmental "master transcription factors" (MTFs) can be subverted during tumorigenesis to control oncogenic transcriptional programs. Current approaches to identifying MTFs rely on ChIP-seq data, which is unavailable for many cancers. We developed the CaCTS (Cancer Core Transcription factor Specificity) algorithm to prioritize candidate MTFs using pan-cancer RNA sequencing data. CaCTS identified candidate MTFs across 34 tumor types and 140 subtypes including predictions for cancer types/subtypes for which MTFs are unknown, including e.g. PAX8, SOX17, and MECOM as candidates in ovarian cancer (OvCa). In OvCa cells, consistent with known MTF properties, these factors are required for viability, lie proximal to superenhancers, co-occupy regulatory elements globally, co-bind loci encoding OvCa biomarkers, and are sensitive to pharmacologic inhibition of transcription. Our predictions of MTFs, especially for tumor types with limited understanding of transcriptional drivers, pave the way to therapeutic targeting of MTFs in a broad spectrum of cancers.

Subtype heterogeneity and epigenetic convergence in neuroendocrine prostate cancer.

Neuroendocrine carcinomas (NEC) are tumors expressing markers of neuronal differentiation that can arise at different anatomic sites but have strong histological and clinical similarities. Here we report the chromatin landscapes of a range of human NECs and show convergence to the activation of a common epigenetic program. With a particular focus on treatment emergent neuroendocrine prostate cancer (NEPC), we analyze cell lines, patient-derived xenograft (PDX) models and human clinical samples to show the existence of two distinct NEPC subtypes based on the expression of the neuronal transcription factors ASCL1 and NEUROD1. While in cell lines and PDX models these subtypes are mutually exclusive, single-cell analysis of human clinical samples exhibits a more complex tumor structure with subtypes coexisting as separate sub-populations within the same tumor. These tumor sub-populations differ genetically and epigenetically contributing to intra- and inter-tumoral heterogeneity in human metastases. Overall, our results provide a deeper understanding of the shared clinicopathological characteristics shown by NECs. Furthermore, the intratumoral heterogeneity of human NEPCs suggests the requirement of simultaneous targeting of coexisting tumor populations as a therapeutic strategy.