In Situ Capture of Chromatin Interactions by Biotinylated dCas9.

Cis-regulatory elements (CREs) are commonly recognized by correlative chromatin features, yet the molecular composition of the vast majority of CREs in chromatin remains unknown. Here, we describe a CRISPR affinity purification in situ of regulatory elements (CAPTURE) approach to unbiasedly identify locus-specific chromatin-regulating protein complexes and long-range DNA interactions. Using an in vivo biotinylated nuclease-deficient Cas9 protein and sequence-specific guide RNAs, we show high-resolution and selective isolation of chromatin interactions at a single-copy genomic locus. Purification of human telomeres using CAPTURE identifies known and new telomeric factors. In situ capture of individual constituents of the enhancer cluster controlling human ß-globin genes establishes evidence for composition-based hierarchical organization. Furthermore, unbiased analysis of chromatin interactions at disease-associated cis-elements and developmentally regulated super-enhancers reveals spatial features that causally control gene transcription. Thus, comprehensive and unbiased analysis of locus-specific regulatory composition provides mechanistic insight into genome structure and function in development and disease.

Structural variation cooperates with permissive chromatin to control enhancer hijacking-mediated oncogenic transcription.

Structural variants (SVs) involving enhancer hijacking can rewire chromatin topologies to cause oncogene activation in human cancers, including hematologic malignancies; however, because of the lack of tools to assess their effects on gene regulation and chromatin organization, the molecular determinants for the functional output of enhancer hijacking remain poorly understood. Here, we developed a multimodal approach to integrate genome sequencing, chromosome conformation, chromatin state, and transcriptomic alteration for quantitative analysis of transcriptional effects and structural reorganization imposed by SVs in leukemic genomes. We identified known and new pathogenic SVs, including recurrent t(5;14) translocations that cause the hijacking of BCL11B enhancers for the allele-specific activation of TLX3 in a subtype of pediatric leukemia. Epigenetic perturbation of SV-hijacked BCL11B enhancers impairs TLX3 transcription, which are required for the growth of t(5;14) leukemia cells. By CRISPR engineering of patient-derived t(5;14) in isogenic leukemia cells, we uncovered a new mechanism whereby the transcriptional output of SV-induced BCL11B enhancer hijacking is dependent on the loss of DNA hypermethylation at the TLX3 promoter. Our results highlight the importance of the cooperation between genetic alteration and permissive chromatin as a critical determinant of SV-mediated oncogene activation, with implications for understanding aberrant gene transcription after epigenetic therapies in patients with leukemia. Hence, leveraging the interdependency of genetic alteration on chromatin variation may provide new opportunities to reprogram gene regulation as targeted interventions in human disease.

T/myeloid mixed phenotype acute leukaemia harbouring TLX3::BCL11B with TLX3 activation.

T/myeloid mixed phenotype acute leukaemia (MPAL) is a rare aggressive acute leukaemia with poorly understood pathogenesis. Herein, we report two cases of T/myeloid MPAL harbouring BCL11B-associated structural variants that activate TLX3 (TLX3::BCL11B-TLX3-activation) by genome sequencing and transcriptomic analyses. Both patients were young males with extramedullary involvement. Cooperative gene alterations characteristic of T/myeloid MPAL and T-lymphoblastic leukaemia (T-ALL) were detected. Both patients achieved initial remission following lineage-matched ALL-based therapy with one patient requiring a lineage-switched myeloid-based therapy. Our study is the first to demonstrate the clinicopathological and genomic features of TLX3::BCL11B-TLX3-activated T/myeloid MPAL and provide insights into leukaemogenesis.

Epo-IGF1R cross talk expands stress-specific progenitors in regenerative erythropoiesis and myeloproliferative neoplasm.

We found that in regenerative erythropoiesis, the erythroid progenitor landscape is reshaped, and a previously undescribed progenitor population with colony-forming unit-erythroid (CFU-E) activity (stress CFU-E [sCFU-E]) is expanded markedly to restore the erythron. sCFU-E cells are targets of erythropoietin (Epo), and sCFU-E expansion requires signaling from the Epo receptor (EpoR) cytoplasmic tyrosines. Molecularly, Epo promotes sCFU-E expansion via JAK2- and STAT5-dependent expression of IRS2, thus engaging the progrowth signaling from the IGF1 receptor (IGF1R). Inhibition of IGF1R and IRS2 signaling impairs sCFU-E cell growth, whereas exogenous IRS2 expression rescues cell growth in sCFU-E expressing truncated EpoR-lacking cytoplasmic tyrosines. This sCFU-E pathway is the major pathway involved in erythrocytosis driven by the oncogenic JAK2 mutant JAK2(V617F) in myeloproliferative neoplasm. Inability to expand sCFU-E cells by truncated EpoR protects against JAK2(V617F)-driven erythrocytosis. In samples from patients with myeloproliferative neoplasm, the number of sCFU-E-like cells increases, and inhibition of IGR1R and IRS2 signaling blocks Epo-hypersensitive erythroid cell colony formation. In summary, we identified a new stress-specific erythroid progenitor cell population that links regenerative erythropoiesis to pathogenic erythrocytosis.

Regulation of embryonic haematopoietic multipotency by EZH1.

All haematopoietic cell lineages that circulate in the blood of adult mammals derive from multipotent haematopoietic stem cells (HSCs). By contrast, in the blood of mammalian embryos, lineage-restricted progenitors arise first, independently of HSCs, which only emerge later in gestation. As best defined in the mouse, 'primitive' progenitors first appear in the yolk sac at 7.5 days post-coitum. Subsequently, erythroid-myeloid progenitors that express fetal haemoglobin, as well as fetal lymphoid progenitors, develop in the yolk sac and the embryo proper, but these cells lack HSC potential. Ultimately, 'definitive' HSCs with long-term, multilineage potential and the ability to engraft irradiated adults emerge at 10.5 days post-coitum from arterial endothelium in the aorta-gonad-mesonephros and other haemogenic vasculature. The molecular mechanisms of this reverse progression of haematopoietic ontogeny remain unexplained. We hypothesized that the definitive haematopoietic program might be actively repressed in early embryogenesis through epigenetic silencing, and that alleviating this repression would elicit multipotency in otherwise lineage-restricted haematopoietic progenitors. Here we show that reduced expression of the Polycomb group protein EZH1 enhances multi-lymphoid output from human pluripotent stem cells. In addition, Ezh1 deficiency in mouse embryos results in precocious emergence of functional definitive HSCs in vivo. Thus, we identify EZH1 as a repressor of haematopoietic multipotency in the early mammalian embryo.

Temporal dynamics of cardiac hypertrophic growth in response to pressure overload.

Hypertension is one of the most important risk factors of heart failure. In response to high blood pressure, the left ventricle manifests hypertrophic growth to ameliorate wall stress, which may progress into decompensation and trigger pathological cardiac remodeling. Despite the clinical importance, the temporal dynamics of pathological cardiac growth remain elusive. Here, we took advantage of the puromycin labeling approach to measure the relative rates of protein synthesis as a way to delineate the temporal regulation of cardiac hypertrophic growth. We first identified the optimal treatment conditions for puromycin in neonatal rat ventricular myocyte culture. We went on to demonstrate that myocyte growth reached its peak rate after 8-10 h of growth stimulation. At the in vivo level, with the use of an acute surgical model of pressure-overload stress, we observed the maximal growth rate to occur at day 7 after surgery. Moreover, RNA sequencing analysis supports that the most profound transcriptomic changes occur during the early phase of hypertrophic growth. Our results therefore suggest that cardiac myocytes mount an immediate growth response in reply to pressure overload followed by a gradual return to basal levels of protein synthesis, highlighting the temporal dynamics of pathological cardiac hypertrophic growth.NEW & NOTEWORTHY We determined the optimal conditions of puromycin incorporation in cardiac myocyte culture. We took advantage of this approach to identify the growth dynamics of cardiac myocytes in vitro. We went further to discover the protein synthesis rate in vivo, which provides novel insights about cardiac temporal growth dynamics in response to pressure overload.

Electron transport chain inhibition increases cellular dependence on purine transport and salvage.

Mitochondria house many metabolic pathways required for homeostasis and growth. To explore how human cells respond to mitochondrial dysfunction, we performed metabolomics in fibroblasts from patients with various mitochondrial disorders and cancer cells with electron transport chain (ETC) blockade. These analyses revealed extensive perturbations in purine metabolism, and stable isotope tracing demonstrated that ETC defects suppress de novo purine synthesis while enhancing purine salvage. In human lung cancer, tumors with markers of low oxidative mitochondrial metabolism exhibit enhanced expression of the salvage enzyme hypoxanthine phosphoribosyl transferase 1 (HPRT1) and high levels of the HPRT1 product inosine monophosphate. Mechanistically, ETC blockade activates the pentose phosphate pathway, providing phosphoribosyl diphosphate to drive purine salvage supplied by uptake of extracellular bases. Blocking HPRT1 sensitizes cancer cells to ETC inhibition. These findings demonstrate how cells remodel purine metabolism upon ETC blockade and uncover a new metabolic vulnerability in tumors with low respiration.

Light-activated macromolecular phase separation modulates transcription by reconfiguring chromatin interactions.

Biomolecular condensates participate in the regulation of gene transcription, yet the relationship between nuclear condensation and transcriptional activation remains elusive. Here, we devised a biotinylated CRISPR-dCas9-based optogenetic method, light-activated macromolecular phase separation (LAMPS), to enable inducible formation, affinity purification, and multiomic dissection of nuclear condensates at the targeted genomic loci. LAMPS-induced condensation at enhancers and promoters activates endogenous gene transcription by chromatin reconfiguration, causing increased chromatin accessibility and de novo formation of long-range chromosomal loops. Proteomic profiling of light-induced condensates by dCas9-mediated affinity purification uncovers multivalent interaction-dependent remodeling of macromolecular composition, resulting in the selective enrichment of transcriptional coactivators and chromatin structure proteins. Our findings support a model whereby the formation of nuclear condensates at native genomic loci reconfigures chromatin architecture and multiprotein assemblies to modulate gene transcription. Hence, LAMPS facilitates mechanistic interrogation of the relationship between nuclear condensation, genome structure, and gene transcription in living cells.

Electron transport chain inhibition increases cellular dependence on purine transport and salvage.

Cancer cells reprogram their metabolism to support cell growth and proliferation in harsh environments. While many studies have documented the importance of mitochondrial oxidative phosphorylation (OXPHOS) in tumor growth, some cancer cells experience conditions of reduced OXPHOS in vivo and induce alternative metabolic pathways to compensate. To assess how human cells respond to mitochondrial dysfunction, we performed metabolomics in fibroblasts and plasma from patients with inborn errors of mitochondrial metabolism, and in cancer cells subjected to inhibition of the electron transport chain (ETC). All these analyses revealed extensive perturbations in purine-related metabolites; in non-small cell lung cancer (NSCLC) cells, ETC blockade led to purine metabolite accumulation arising from a reduced cytosolic NAD + /NADH ratio (NADH reductive stress). Stable isotope tracing demonstrated that ETC deficiency suppressed de novo purine nucleotide synthesis while enhancing purine salvage. Analysis of NSCLC patients infused with [U- 13 C]glucose revealed that tumors with markers of low oxidative mitochondrial metabolism exhibited high expression of the purine salvage enzyme HPRT1 and abundant levels of the HPRT1 product inosine monophosphate (IMP). ETC blockade also induced production of ribose-5' phosphate (R5P) by the pentose phosphate pathway (PPP) and import of purine nucleobases. Blocking either HPRT1 or nucleoside transporters sensitized cancer cells to ETC inhibition, and overexpressing nucleoside transporters was sufficient to drive growth of NSCLC xenografts. Collectively, this study mechanistically delineates how cells compensate for suppressed purine metabolism in response to ETC blockade, and uncovers a new metabolic vulnerability in tumors experiencing NADH excess.

Disabling Uncompetitive Inhibition of Oncogenic IDH Mutations Drives Acquired Resistance.

Mutations in IDH genes occur frequently in acute myeloid leukemia (AML) and other human cancers to generate the oncometabolite R-2HG. Allosteric inhibition of mutant IDH suppresses R-2HG production in a subset of patients with AML; however, acquired resistance emerges as a new challenge, and the underlying mechanisms remain incompletely understood. Here we establish isogenic leukemia cells containing common IDH oncogenic mutations by CRISPR base editing. By mutational scanning of IDH single amino acid variants in base-edited cells, we describe a repertoire of IDH second-site mutations responsible for therapy resistance through disabling uncompetitive enzyme inhibition. Recurrent mutations at NADPH binding sites within IDH heterodimers act in cis or trans to prevent the formation of stable enzyme-inhibitor complexes, restore R-2HG production in the presence of inhibitors, and drive therapy resistance in IDH-mutant AML cells and patients. We therefore uncover a new class of pathogenic mutations and mechanisms for acquired resistance to targeted cancer therapies. SIGNIFICANCE: Comprehensive scanning of IDH single amino acid variants in base-edited leukemia cells uncovers recurrent mutations conferring resistance to IDH inhibition through disabling NADPH-dependent uncompetitive inhibition. Together with targeted sequencing, structural, and functional studies, we identify a new class of pathogenic mutations and mechanisms for acquired resistance to IDH-targeting cancer therapies. This article is highlighted in the In This Issue feature, p. 1.

Metabolon formation regulates branched-chain amino acid oxidation and homeostasis.

The branched-chain aminotransferase isozymes BCAT1 and BCAT2, segregated into distinct subcellular compartments and tissues, initiate the catabolism of branched-chain amino acids (BCAAs). However, whether and how BCAT isozymes cooperate with downstream enzymes to control BCAA homeostasis in an intact organism remains largely unknown. Here, we analyse system-wide metabolomic changes in BCAT1- and BCAT2-deficient mouse models. Loss of BCAT2 but not BCAT1 leads to accumulation of BCAAs and branched-chain α-keto acids (BCKAs), causing morbidity and mortality that can be ameliorated by dietary BCAA restriction. Through proximity labelling, isotope tracing and enzymatic assays, we provide evidence for the formation of a mitochondrial BCAA metabolon involving BCAT2 and branched-chain α-keto acid dehydrogenase. Disabling the metabolon contributes to BCAT2 deficiency-induced phenotypes, which can be reversed by BCAT1-mediated BCKA reamination. These findings establish a role for metabolon formation in BCAA metabolism in vivo and suggest a new strategy to modulate this pathway in diseases involving dysfunctional BCAA metabolism.

Convergence of oncogenic cooperation at single-cell and single-gene levels drives leukemic transformation.

Cancers develop from the accumulation of somatic mutations, yet it remains unclear how oncogenic lesions cooperate to drive cancer progression. Using a mouse model harboring NRasG12D and EZH2 mutations that recapitulates leukemic progression, we employ single-cell transcriptomic profiling to map cellular composition and gene expression alterations in healthy or diseased bone marrows during leukemogenesis. At cellular level, NRasG12D induces myeloid lineage-biased differentiation and EZH2-deficiency impairs myeloid cell maturation, whereas they cooperate to promote myeloid neoplasms with dysregulated transcriptional programs. At gene level, NRasG12D and EZH2-deficiency independently and synergistically deregulate gene expression. We integrate results from histopathology, leukemia repopulation, and leukemia-initiating cell assays to validate transcriptome-based cellular profiles. We use this resource to relate developmental hierarchies to leukemia phenotypes, evaluate oncogenic cooperation at single-cell and single-gene levels, and identify GEM as a regulator of leukemia-initiating cells. Our studies establish an integrative approach to deconvolute cancer evolution at single-cell resolution in vivo.

A unified model of human hemoglobin switching through single-cell genome editing.

Key mechanisms of fetal hemoglobin (HbF) regulation and switching have been elucidated through studies of human genetic variation, including mutations in the HBG1/2 promoters, deletions in the ß-globin locus, and variation impacting BCL11A. While this has led to substantial insights, there has not been a unified understanding of how these distinct genetically-nominated elements, as well as other key transcription factors such as ZBTB7A, collectively interact to regulate HbF. A key limitation has been the inability to model specific genetic changes in primary isogenic human hematopoietic cells to uncover how each of these act individually and in aggregate. Here, we describe a single-cell genome editing functional assay that enables specific mutations to be recapitulated individually and in combination, providing insights into how multiple mutation-harboring functional elements collectively contribute to HbF expression. In conjunction with quantitative modeling and chromatin capture analyses, we illustrate how these genetic findings enable a comprehensive understanding of how distinct regulatory mechanisms can synergistically modulate HbF expression.

Silencing of LINE-1 retrotransposons is a selective dependency of myeloid leukemia.

Transposable elements or transposons are major players in genetic variability and genome evolution. Aberrant activation of long interspersed element-1 (LINE-1 or L1) retrotransposons is common in human cancers, yet their tumor-type-specific functions are poorly characterized. We identified MPHOSPH8/MPP8, a component of the human silencing hub (HUSH) complex, as an acute myeloid leukemia (AML)-selective dependency by epigenetic regulator-focused CRISPR screening. Although MPP8 is dispensable for steady-state hematopoiesis, MPP8 loss inhibits AML development by reactivating L1s to induce the DNA damage response and cell cycle exit. Activation of endogenous or ectopic L1s mimics the phenotype of MPP8 loss, whereas blocking retrotransposition abrogates MPP8-deficiency-induced phenotypes. Expression of AML oncogenic mutations promotes L1 suppression, and enhanced L1 silencing is associated with poor prognosis in human AML. Hence, while retrotransposons are commonly recognized for their cancer-promoting functions, we describe a tumor-suppressive role for L1 retrotransposons in myeloid leukemia.

Noncoding Variants Connect Enhancer Dysregulation with Nuclear Receptor Signaling in Hematopoietic Malignancies.

Mutations in protein-coding genes are well established as the basis for human cancer, yet how alterations within noncoding genome, a substantial fraction of which contain cis-regulatory elements (CRE), contribute to cancer pathophysiology remains elusive. Here, we developed an integrative approach to systematically identify and characterize noncoding regulatory variants with functional consequences in human hematopoietic malignancies. Combining targeted resequencing of hematopoietic lineage-associated CREs and mutation discovery, we uncovered 1,836 recurrently mutated CREs containing leukemia-associated noncoding variants. By enhanced CRISPR/dCas9-based CRE perturbation screening and functional analyses, we identified 218 variant-associated oncogenic or tumor-suppressive CREs in human leukemia. Noncoding variants at KRAS and PER2 enhancers reside in proximity to nuclear receptor (NR) binding regions and modulate transcriptional activities in response to NR signaling in leukemia cells. NR binding sites frequently colocalize with noncoding variants across cancer types. Hence, recurrent noncoding variants connect enhancer dysregulation with nuclear receptor signaling in hematopoietic malignancies. SIGNIFICANCE: We describe an integrative approach to identify noncoding variants in human leukemia, and reveal cohorts of variant-associated oncogenic and tumor-suppressive cis-regulatory elements including KRAS and PER2 enhancers. Our findings support a model in which noncoding regulatory variants connect enhancer dysregulation with nuclear receptor signaling to modulate gene programs in hematopoietic malignancies.See related commentary by van Galen, p. 646.This article is highlighted in the In This Issue feature, p. 627.

Multiplexed capture of spatial configuration and temporal dynamics of locus-specific 3D chromatin by biotinylated dCas9.

The spatiotemporal control of 3D genome is fundamental for gene regulation, yet it remains challenging to profile high-resolution chromatin structure at cis-regulatory elements (CREs). Using C-terminally biotinylated dCas9, endogenous biotin ligases, and pooled sgRNAs, we describe the dCas9-based CAPTURE method for multiplexed analysis of locus-specific chromatin interactions. The redesigned system allows for quantitative analysis of the spatial configuration of a few to hundreds of enhancers or promoters in a single experiment, enabling comparisons across CREs within and between gene clusters. Multiplexed analyses of the spatiotemporal configuration of erythroid super-enhancers and promoter-centric interactions reveal organizational principles of genome structure and function.

Discovering How Heme Controls Genome Function Through Heme-omics.

Protein ensembles control genome function by establishing, maintaining, and deconstructing cell-type-specific chromosomal landscapes. A plethora of small molecules orchestrate cellular functions and therefore may link physiological processes with genome biology. The metabolic enzyme and hemoglobin cofactor heme induces proteolysis of a transcriptional repressor, Bach1, and regulates gene expression post-transcriptionally. However, whether heme controls genome function broadly or through prescriptive actions is unclear. Using assay for transposase-accessible chromatin sequencing (ATAC-seq), we establish a heme-dependent chromatin atlas in wild-type and mutant erythroblasts lacking enhancers that confer normal heme synthesis. Amalgamating chromatin landscapes and transcriptomes in cells with sub-physiological heme and post-heme rescue reveals parallel Bach1-dependent and Bach1-independent mechanisms that target heme-sensing chromosomal hotspots. The hotspots harbor a DNA motif demarcating heme-regulated chromatin and genes encoding proteins not known to be heme regulated, including metabolic enzymes. The heme-omics analysis establishes how an essential biochemical cofactor controls genome function and cellular physiology.

Interrogation of enhancer function by enhancer-targeting CRISPR epigenetic editing.

Tissue-specific gene expression requires coordinated control of gene-proximal and -distal cis-regulatory elements (CREs), yet functional analysis of gene-distal CREs such as enhancers remains challenging. Here we describe CRISPR/dCas9-based enhancer-targeting epigenetic editing systems, enCRISPRa and enCRISPRi, for efficient analysis of enhancer function in situ and in vivo. Using dual effectors capable of re-writing enhancer-associated chromatin modifications, we show that enCRISPRa and enCRISPRi modulate gene transcription by remodeling local epigenetic landscapes at sgRNA-targeted enhancers and associated genes. Comparing with existing methods, the improved systems display more robust perturbations of enhancer activity and gene transcription with minimal off-targets. Allele-specific targeting of enCRISPRa to oncogenic TAL1 super-enhancer modulates TAL1 expression and cancer progression in xenotransplants. Single or multi-loci perturbations of lineage-specific enhancers using an enCRISPRi knock-in mouse establish in vivo evidence for lineage-restricted essentiality of developmental enhancers during hematopoiesis. Hence, enhancer-targeting CRISPR epigenetic editing provides opportunities for interrogating enhancer function in native biological contexts.

Lactate Dehydrogenase A Governs Cardiac Hypertrophic Growth in Response to Hemodynamic Stress.

The heart manifests hypertrophic growth in response to high blood pressure, which may decompensate and progress to heart failure under persistent stress. Metabolic remodeling is an early event in this process. However, its role remains to be fully characterized. Here, we show that lactate dehydrogenase A (LDHA), a critical glycolytic enzyme, is elevated in the heart in response to hemodynamic stress. Cardiomyocyte-restricted deletion of LDHA leads to defective cardiac hypertrophic growth and heart failure by pressure overload. Silencing of LDHA in cultured cardiomyocytes suppresses cell growth from pro-hypertrophic stimulation in vitro, while overexpression of LDHA is sufficient to drive cardiomyocyte growth. Furthermore, we find that lactate is capable of rescuing the growth defect from LDHA knockdown. Mechanistically, lactate stabilizes NDRG3 (N-myc downregulated gene family 3) and stimulates ERK (extracellular signal-regulated kinase). Our results together suggest that the LDHA/NDRG3 axis may play a critical role in adaptive cardiomyocyte growth in response to hemodynamic stress.

MAP: model-based analysis of proteomic data to detect proteins with significant abundance changes.

Isotope-labeling-based mass spectrometry (MS) is widely used in quantitative proteomic studies. With this technique, the relative abundance of thousands of proteins can be efficiently profiled in parallel, greatly facilitating the detection of proteins differentially expressed across samples. However, this task remains computationally challenging. Here we present a new approach, termed Model-based Analysis of Proteomic data (MAP), for this task. Unlike many existing methods, MAP does not require technical replicates to model technical and systematic errors, and instead utilizes a novel step-by-step regression analysis to directly assess the significance of observed protein abundance changes. We applied MAP to compare the proteomic profiles of undifferentiated and differentiated mouse embryonic stem cells (mESCs), and found it has superior performance compared with existing tools in detecting proteins differentially expressed during mESC differentiation. A web-based application of MAP is provided for online data processing at http://bioinfo.sibs.ac.cn/shaolab/MAP.