TAD hierarchy restricts poised LTR activation and loss of TAD hierarchy promotes LTR co-option in cancer.

Transposable elements (TEs) are abundant in the human genome, and they provide the sources for genetic and functional diversity. The regulation of TEs expression and their functional consequences in physiological conditions and cancer development remain to be fully elucidated. Previous studies suggested TEs are repressed by DNA methylation and chromatin modifications. The effect of 3D chromatin topology on TE regulation remains elusive. Here, by integrating transcriptome and 3D genome architecture studies, we showed that haploinsufficient loss of NIPBL selectively activates alternative promoters at the long terminal repeats (LTRs) of the TE subclasses. This activation occurs through the reorganization of topologically associating domain (TAD) hierarchical structures and recruitment of proximal enhancers. These observations indicate that TAD hierarchy restricts transcriptional activation of LTRs that already possess open chromatin features. In cancer, perturbation of the hierarchical chromatin topology can lead to co-option of LTRs as functional alternative promoters in a context-dependent manner and drive aberrant transcriptional activation of novel oncogenes and other divergent transcripts. These data uncovered a new layer of regulatory mechanism of TE expression beyond DNA and chromatin modification in human genome. They also posit the TAD hierarchy dysregulation as a novel mechanism for alternative promoter-mediated oncogene activation and transcriptional diversity in cancer, which may be exploited therapeutically.

3D Enhancer-promoter networks provide predictive features for gene expression and coregulation in early embryonic lineages.

Mammalian embryogenesis commences with two pivotal and binary cell fate decisions that give rise to three essential lineages: the trophectoderm, the epiblast and the primitive endoderm. Although key signaling pathways and transcription factors that control these early embryonic decisions have been identified, the non-coding regulatory elements through which transcriptional regulators enact these fates remain understudied. Here, we characterize, at a genome-wide scale, enhancer activity and 3D connectivity in embryo-derived stem cell lines that represent each of the early developmental fates. We observe extensive enhancer remodeling and fine-scale 3D chromatin rewiring among the three lineages, which strongly associate with transcriptional changes, although distinct groups of genes are irresponsive to topological changes. In each lineage, a high degree of connectivity, or 'hubness', positively correlates with levels of gene expression and enriches for cell-type specific and essential genes. Genes within 3D hubs also show a significantly stronger probability of coregulation across lineages compared to genes in linear proximity or within the same contact domains. By incorporating 3D chromatin features, we build a predictive model for transcriptional regulation (3D-HiChAT) that outperforms models using only 1D promoter or proximal variables to predict levels and cell-type specificity of gene expression. Using 3D-HiChAT, we identify, in silico, candidate functional enhancers and hubs in each cell lineage, and with CRISPRi experiments, we validate several enhancers that control gene expression in their respective lineages. Our study identifies 3D regulatory hubs associated with the earliest mammalian lineages and describes their relationship to gene expression and cell identity, providing a framework to comprehensively understand lineage-specific transcriptional behaviors.

Systematic mapping and modeling of 3D enhancer-promoter interactions in early mouse embryonic lineages reveal regulatory principles that determine the levels and cell-type specificity of gene expression.

Mammalian embryogenesis commences with two pivotal and binary cell fate decisions that give rise to three essential lineages, the trophectoderm (TE), the epiblast (EPI) and the primitive endoderm (PrE). Although key signaling pathways and transcription factors that control these early embryonic decisions have been identified, the non-coding regulatory elements via which transcriptional regulators enact these fates remain understudied. To address this gap, we have characterized, at a genome-wide scale, enhancer activity and 3D connectivity in embryo-derived stem cell lines that represent each of the early developmental fates. We observed extensive enhancer remodeling and fine-scale 3D chromatin rewiring among the three lineages, which strongly associate with transcriptional changes, although there are distinct groups of genes that are irresponsive to topological changes. In each lineage, a high degree of connectivity or "hubness" positively correlates with levels of gene expression and enriches for cell-type specific and essential genes. Genes within 3D hubs also show a significantly stronger probability of coregulation across lineages, compared to genes in linear proximity or within the same contact domains. By incorporating 3D chromatin features, we build a novel predictive model for transcriptional regulation (3D-HiChAT), which outperformed models that use only 1D promoter or proximal variables in predicting levels and cell-type specificity of gene expression. Using 3D-HiChAT, we performed genome-wide in silico perturbations to nominate candidate functional enhancers and hubs in each cell lineage, and with CRISPRi experiments we validated several novel enhancers that control expression of one or more genes in their respective lineages. Our study comprehensively identifies 3D regulatory hubs associated with the earliest mammalian lineages and describes their relationship to gene expression and cell identity, providing a framework to understand lineage-specific transcriptional behaviors.

Single-cell analysis of bidirectional reprogramming between early embryonic states reveals mechanisms of differential lineage plasticities.

Two distinct fates, pluripotent epiblast (EPI) and primitive (extra-embryonic) endoderm (PrE), arise from common progenitor cells, the inner cell mass (ICM), in mammalian embryos. To study how these sister identities are forged, we leveraged embryonic (ES) and eXtraembryonic ENdoderm (XEN) stem cells - in vitro counterparts of the EPI and PrE. Bidirectional reprogramming between ES and XEN coupled with single-cell RNA and ATAC-seq analyses uncovered distinct rates, efficiencies and trajectories of state conversions, identifying drivers and roadblocks of reciprocal conversions. While GATA4-mediated ES-to-iXEN conversion was rapid and nearly deterministic, OCT4, KLF4 and SOX2-induced XEN-to-iPS reprogramming progressed with diminished efficiency and kinetics. The dominant PrE transcriptional program, safeguarded by Gata4, and globally elevated chromatin accessibility of EPI underscored the differential plasticities of the two states. Mapping in vitro trajectories to embryos revealed reprogramming in either direction tracked along, and toggled between, EPI and PrE in vivo states without transitioning through the ICM.

A Knock-In Mouse Model of Thymoma With the GTF2I L424H Mutation.

INTRODUCTION: The pathogenesis of thymic epithelial tumors remains largely unknown. We previously identified GTF2I L424H as the most frequently recurrent mutation in thymic epithelial tumors. Nevertheless, the precise role of this mutation in tumorigenesis of thymic epithelial cells is unclear. METHODS: To investigate the role of GTF2I L424H mutation in thymic epithelial cells in vivo, we generated and characterized a mouse model in which the Gtf2i L424H mutation was conditionally knocked-in in the Foxn1+ thymic epithelial cells. Digital spatial profiling was performed on thymomas and normal thymic tissues with GeoMx-mouse whole transcriptome atlas. Immunohistochemistry staining was performed using both mouse tissues and human thymic epithelial tumors. RESULTS: We observed that the Gtf2i mutation impairs development of the thymic medulla and maturation of medullary thymic epithelial cells in young mice and causes tumor formation in the thymus of aged mice. Cell cycle-related pathways, such as E2F targets and MYC targets, are enriched in the tumor epithelial cells. Results of gene set variation assay analysis revealed that gene signatures of cortical thymic epithelial cells and thymic epithelial progenitor cells are also enriched in the thymomas of the knock-in mice, which mirrors the human counterparts in The Cancer Genome Atlas database. Immunohistochemistry results revealed similar expression pattern of epithelial cell markers between mouse and human thymomas. CONCLUSIONS: We have developed and characterized a novel thymoma mouse model. This study improves knowledge of the molecular drivers in thymic epithelial cells and provides a tool for further study of the biology of thymic epithelial tumors and for development of novel therapies.

Assessing Specific Networks of Chromatin Interactions with HiChIP.

The introduction of chromosome conformation capture (3C)-based technologies coupled with next-generation sequencing have significantly advanced our understanding of how the genetic material is organized within the eukaryotic nucleus. Three-dimensional (3D) genomic organization occurs at hierarchical levels, ranging from chromosome territories and subnuclear compartments to smaller self-associated domains and fine-scale chromatin interactions. The latter can be further categorized into different subtypes, such as structural or regulatory, based either on their presumed functionality and/or the factors that mediate their formation. Various enrichment strategies coupled with 3C-based technologies have been developed to prospectively isolate and quantify chromatin interactions around regions occupied by specific proteins or marks of interest. These approaches not only enable high-resolution characterization of the selected chromatin contacts at a cost-effective manner, but also offer important biological insights into their organizational principles and regulatory function. In this chapter, we will focus on the recently developed HiChIP technology with an emphasis on the discovery of putative active enhancers and promoter interactions in cell types of interest. We will describe the specific steps for designing, performing and analyzing successful HiChIP experiments as well as important limitations and considerations.

A bipartite element with allele-specific functions safeguards DNA methylation imprints at the Dlk1-Dio3 locus.

Loss of imprinting (LOI) results in severe developmental defects, but the mechanisms preventing LOI remain incompletely understood. Here, we dissect the functional components of the imprinting control region of the essential Dlk1-Dio3 locus (called IG-DMR) in pluripotent stem cells. We demonstrate that the IG-DMR consists of two antagonistic elements: a paternally methylated CpG island that prevents recruitment of TET dioxygenases and a maternally unmethylated non-canonical enhancer that ensures expression of the Gtl2 lncRNA by counteracting de novo DNA methyltransferases. Genetic or epigenetic editing of these elements leads to distinct LOI phenotypes with characteristic alternations of allele-specific gene expression, DNA methylation, and 3D chromatin topology. Although repression of the Gtl2 promoter results in dysregulated imprinting, the stability of LOI phenotypes depends on the IG-DMR, suggesting a functional hierarchy. These findings establish the IG-DMR as a bipartite control element that maintains imprinting by allele-specific restriction of the DNA (de)methylation machinery.

H3K27ac bookmarking promotes rapid post-mitotic activation of the pluripotent stem cell program without impacting 3D chromatin reorganization.

During self-renewal, cell-type-defining features are drastically perturbed in mitosis and must be faithfully reestablished upon G1 entry, a process that remains largely elusive. Here, we characterized at a genome-wide scale the dynamic transcriptional and architectural resetting of mouse pluripotent stem cells (PSCs) upon mitotic exit. We captured distinct waves of transcriptional reactivation with rapid induction of stem cell genes and transient activation of lineage-specific genes. Topological reorganization at different hierarchical levels also occurred in an asynchronous manner and showed partial coordination with transcriptional resetting. Globally, rapid transcriptional and architectural resetting associated with mitotic retention of H3K27 acetylation, supporting a bookmarking function. Indeed, mitotic depletion of H3K27ac impaired the early reactivation of bookmarked, stem-cell-associated genes. However, 3D chromatin reorganization remained largely unaffected, suggesting that these processes are driven by distinct forces upon mitotic exit. This study uncovers principles and mediators of PSC molecular resetting during self-renewal.

Identification of a dynamic gene regulatory network required for pluripotency factor-induced reprogramming of mouse fibroblasts and hepatocytes.

The generation of induced pluripotent stem cells (iPSCs) from somatic cells provides an excellent model to study mechanisms of transcription factor-induced global alterations of the epigenome and genome function. Here, we have investigated the early transcriptional events of cellular reprogramming triggered by the co-expression of Oct4, Sox2, Klf4, and c-Myc (OSKM) in mouse embryonic fibroblasts (MEFs) and mouse hepatocytes (mHeps). In this analysis, we identified a gene regulatory network composed of nine transcriptional regulators (9TR; Cbfa2t3, Gli2, Irf6, Nanog, Ovol1, Rcan1, Taf1c, Tead4, and Tfap4), which are directly targeted by OSKM, in vivo. Functional studies using single and double shRNA knockdowns of any of these factors caused disruption of the network and dramatic reductions in reprogramming efficiency, indicating that this network is essential for the induction and establishment of pluripotency. We demonstrate that the stochastic co-expression of 9TR network components occurs in a remarkably small number of cells, approximating the percentage of terminally reprogrammed cells as a result of dynamic molecular events. Thus, the early DNA-binding patterns of OSKM and the subsequent probabilistic co-expression of essential 9TR components in subpopulations of cells undergoing reprogramming steer the reconstruction of a gene regulatory network marking the transition to pluripotency.

Context-Dependent Requirement of Euchromatic Histone Methyltransferase Activity during Reprogramming to Pluripotency.

Methylation of histone 3 at lysine 9 (H3K9) constitutes a roadblock for cellular reprogramming. Interference with methyltransferases or activation of demethylases by the cofactor ascorbic acid (AA) facilitates the derivation of induced pluripotent stem cells (iPSCs), but possible interactions between specific methyltransferases and AA treatment remain insufficiently explored. We show that chemical inhibition of the methyltransferases EHMT1 and EHMT2 counteracts iPSC formation in an enhanced reprogramming system in the presence of AA, an effect that is dependent on EHMT1. EHMT inhibition during enhanced reprogramming is associated with rapid loss of H3K9 dimethylation, inefficient downregulation of somatic genes, and failed mesenchymal-to-epithelial transition. Furthermore, transient EHMT inhibition during reprogramming yields iPSCs that fail to efficiently give rise to viable mice upon blastocyst injection. Our observations establish novel functions of H3K9 methyltransferases and suggest that a functional balance between AA-stimulated enzymes and EHMTs supports efficient and less error-prone iPSC reprogramming to pluripotency.

Transcription factors: building hubs in the 3D space.

The hierarchical three-dimensional folding of the mammalian genome constitutes an important regulatory layer of gene expression and cell fate control during processes such as development and tumorigenesis. Accumulating evidence supports the existence of complex topological assemblies in which multiple genes and regulatory elements are frequently interacting with each other in the 3D nucleus. Here, we will discuss the nature, organizational principles, and potential function of such assemblies, including the recently reported enhancer "hubs," "cliques," and FIREs (frequently interacting regions) as well as multi-contact hubs. We will also review recent studies that investigate the role of transcription factors (TFs) in driving the topological genome reorganization and hub formation in the context of cell fate transitions and cancer. Finally, we will highlight technological advances that enabled these studies, current limitations, and future directions necessary to advance our understating in the field.

A Deep Learning Framework for Predicting Response to Therapy in Cancer.

A major challenge in cancer treatment is predicting clinical response to anti-cancer drugs on a personalized basis. Using a pharmacogenomics database of 1,001 cancer cell lines, we trained deep neural networks for prediction of drug response and assessed their performance on multiple clinical cohorts. We demonstrate that deep neural networks outperform the current state in machine learning frameworks. We provide a proof of concept for the use of deep neural network-based frameworks to aid precision oncology strategies.

CHUK/IKK-α loss in lung epithelial cells enhances NSCLC growth associated with HIF up-regulation.

Through the progressive accumulation of genetic and epigenetic alterations in cellular physiology, non-small-cell lung cancer (NSCLC) evolves in distinct steps involving mutually exclusive oncogenic mutations in K-Ras or EGFR along with inactivating mutations in the p53 tumor suppressor. Herein, we show two independent in vivo lung cancer models in which CHUK/IKK-α acts as a major NSCLC tumor suppressor. In a novel transgenic mouse strain, wherein IKKα ablation is induced by tamoxifen (Tmx) solely in alveolar type II (AT-II) lung epithelial cells, IKKα loss increases the number and size of lung adenomas in response to the chemical carcinogen urethane, whereas IKK-ß instead acts as a tumor promoter in this same context. IKKα knockdown in three independent human NSCLC lines (independent of K-Ras or p53 status) enhances their growth as tumor xenografts in immune-compromised mice. Bioinformatics analysis of whole transcriptome profiling followed by quantitative protein and targeted gene expression validation experiments reveals that IKKα loss can result in the up-regulation of activated HIF-1-α protein to enhance NSCLC tumor growth under hypoxic conditions in vivo.

KLF4 is involved in the organization and regulation of pluripotency-associated three-dimensional enhancer networks.

Cell fate transitions are accompanied by global transcriptional, epigenetic and topological changes driven by transcription factors, as is exemplified by reprogramming somatic cells to pluripotent stem cells through the expression of OCT4, KLF4, SOX2 and cMYC. How transcription factors orchestrate the complex molecular changes around their target gene loci remains incompletely understood. Here, using KLF4 as a paradigm, we provide a transcription-factor-centric view of chromatin reorganization and its association with three-dimensional enhancer rewiring and transcriptional changes during the reprogramming of mouse embryonic fibroblasts to pluripotent stem cells. Inducible depletion of KLF factors in PSCs caused a genome-wide decrease in enhancer connectivity, whereas disruption of individual KLF4 binding sites within pluripotent-stem-cell-specific enhancers was sufficient to impair enhancer-promoter contacts and reduce the expression of associated genes. Our study provides an integrative view of the complex activities of a lineage-specifying transcription factor and offers novel insights into the nature of the molecular events that follow transcription factor binding.

Machine learning and data mining frameworks for predicting drug response in cancer: An overview and a novel in silico screening process based on association rule mining.

A major challenge in cancer treatment is predicting the clinical response to anti-cancer drugs on a personalized basis. The success of such a task largely depends on the ability to develop computational resources that integrate big "omic" data into effective drug-response models. Machine learning is both an expanding and an evolving computational field that holds promise to cover such needs. Here we provide a focused overview of: 1) the various supervised and unsupervised algorithms used specifically in drug response prediction applications, 2) the strategies employed to develop these algorithms into applicable models, 3) data resources that are fed into these frameworks and 4) pitfalls and challenges to maximize model performance. In this context we also describe a novel in silico screening process, based on Association Rule Mining, for identifying genes as candidate drivers of drug response and compare it with relevant data mining frameworks, for which we generated a web application freely available at: https://compbio.nyumc.org/drugs/. This pipeline explores with high efficiency large sample-spaces, while is able to detect low frequency events and evaluate statistical significance even in the multidimensional space, presenting the results in the form of easily interpretable rules. We conclude with future prospects and challenges of applying machine learning based drug response prediction in precision medicine.

The lysine-specific methyltransferase KMT2C/MLL3 regulates DNA repair components in cancer.

Genome-wide studies in tumor cells have indicated that chromatin-modifying proteins are commonly mutated in human cancers. The lysine-specific methyltransferase 2C (KMT2C/MLL3) is a putative tumor suppressor in several epithelia and in myeloid cells. Here, we show that downregulation of KMT2C in bladder cancer cells leads to extensive changes in the epigenetic status and the expression of DNA damage response and DNA repair genes. More specifically, cells with low KMT2C activity are deficient in homologous recombination-mediated double-strand break DNA repair. Consequently, these cells suffer from substantially higher endogenous DNA damage and genomic instability. Finally, these cells seem to rely heavily on PARP1/2 for DNA repair, and treatment with the PARP1/2 inhibitor olaparib leads to synthetic lethality, suggesting that cancer cells with low KMT2C expression are attractive targets for therapies with PARP1/2 inhibitors.

IFNα Impairs Autophagic Degradation of mtDNA Promoting Autoreactivity of SLE Monocytes in a STING-Dependent Fashion.

Interferon α (IFNα) is a prompt and efficient orchestrator of host defense against nucleic acids but upon chronicity becomes a potent mediator of autoimmunity. Sustained IFNα signaling is linked to pathogenesis of systemic lupus erythematosus (SLE), an incurable autoimmune disease characterized by aberrant self-DNA sensing that culminates in anti-DNA autoantibody-mediated pathology. IFNα instructs monocytes differentiation into autoinflammatory dendritic cells (DCs) than potentiates the survival and expansion of autoreactive lymphocytes, but the molecular mechanism bridging sterile IFNα-danger alarm with adaptive response against self-DNA remains elusive. Herein, we demonstrate IFNα-mediated deregulation of mitochondrial metabolism and impairment of autophagic degradation, leading to cytosolic accumulation of mtDNA that is sensed via stimulator of interferon genes (STING) to promote induction of autoinflammatory DCs. Identification of mtDNA as a cell-autonomous enhancer of IFNα signaling underlines the significance of efficient mitochondrial recycling in the maintenance of peripheral tolerance. Antioxidant treatment and metabolic rescue of autolysosomal degradation emerge as drug targets in SLE and other IFNα-related pathologies.

A metabolic interplay coordinated by HLX regulates myeloid differentiation and AML through partly overlapping pathways.

The H2.0-like homeobox transcription factor (HLX) regulates hematopoietic differentiation and is overexpressed in Acute Myeloid Leukemia (AML), but the mechanisms underlying these functions remain unclear. We demonstrate here that HLX overexpression leads to a myeloid differentiation block both in zebrafish and human hematopoietic stem and progenitor cells (HSPCs). We show that HLX overexpression leads to downregulation of genes encoding electron transport chain (ETC) components and upregulation of PPARδ gene expression in zebrafish and human HSPCs. HLX overexpression also results in AMPK activation. Pharmacological modulation of PPARδ signaling relieves the HLX-induced myeloid differentiation block and rescues HSPC loss upon HLX knockdown but it has no effect on AML cell lines. In contrast, AMPK inhibition results in reduced viability of AML cell lines, but minimally affects myeloid progenitors. This newly described role of HLX in regulating the metabolic state of hematopoietic cells may have important therapeutic implications.