Genome-wide identification of TAL1's functional targets: insights into its mechanisms of action in primary erythroid cells

Transcription factors regulating vasculogenesis and angiogenesis

Transcription factors (TFs) play a crucial role in regulating the dynamic and precise patterns of gene expression required for the initial specification of endothelial cells (ECs), and during endothelial growth and differentiation. While sharing many core features, ECs can be highly heterogeneous. Differential gene expression between ECs is essential to pattern the hierarchical vascular network into arteries, veins and capillaries, to drive angiogenic growth of new vessels, and to direct specialization in response to local signals. Unlike many other cell types, ECs have no single master regulator, instead relying on differing combinations of a necessarily limited repertoire of TFs to achieve tight spatial and temporal activation and repression of gene expression. Here, we will discuss the cohort of TFs known to be involved in directing gene expression during different stages of mammalian vasculogenesis and angiogenesis, with a primary focus on development.

Efficient genome editing in erythroid cells unveils novel MYB target genes and regulatory functions

Targeted genome editing holds great promise in biology. However, efficient genome modification, including gene knock-in (KI), remains an unattained goal in multiple cell types and loci due to poor transfection efficiencies and low target genes expression, impeding the positive selection of recombined cells. Here, we describe a genome editing approach to achieve efficient gene targeting using hard to transfect erythroid cell lines. We demonstrate robust fluorescent protein KI efficiency in low expressed transcription factor (TF) genes (e.g., Myb or Zeb1). We further show the ability to target two independent loci in individual cells, exemplified by MYB-GFP and NuMA-Cherry double KI, allowing multicolor labeling of regulatory factors at physiological endogenous levels. Our KI tagging approach allowed us to perform genome-wide TF analysis at increased signal-to-noise ratios, and highlighted previously unidentified MYB target genes and pathways. Overall, we establish a versatile CRISPR-Cas9-based platform, offering attractive opportunities for the dissection of the erythroid differentiation process.

A genome-wide association study of blood cell morphology identifies cellular proteins implicated in disease aetiology

Blood cells contain functionally important intracellular structures, such as granules, critical to immunity and thrombosis. Quantitative variation in these structures has not been subjected previously to large-scale genetic analysis. We perform genome-wide association studies of 63 flow-cytometry derived cellular phenotypes-including cell-type specific measures of granularity, nucleic acid content and reactivity-in 41,515 participants in the INTERVAL study. We identify 2172 distinct variant-trait associations, including associations near genes coding for proteins in organelles implicated in inflammatory and thrombotic diseases. By integrating with epigenetic data we show that many intracellular structures are likely to be determined in immature precursor cells. By integrating with proteomic data we identify the transcription factor FOG2 as an early regulator of platelet formation and α-granularity. Finally, we show that colocalisation of our associations with disease risk signals can suggest aetiological cell-types-variants in IL2RA and ITGA4 respectively mirror the known effects of daclizumab in multiple sclerosis and vedolizumab in inflammatory bowel disease.

Isoforms of the TAL1 transcription factor have different roles in hematopoiesis and cell growth

T-cell acute lymphoblastic leukemia (T-ALL) protein 1 (TAL1) is a central transcription factor in hematopoiesis. The timing and level of TAL1 expression orchestrate the differentiation to specialized blood cells and its overexpression is a common cause of T-ALL. Here, we studied the 2 protein isoforms of TAL1, short and long, which are generated by the use of alternative promoters as well as by alternative splicing. We analyzed the expression of each isoform by deleting an enhancer or insulator, or by opening chromatin at the enhancer location. Our results show that each enhancer promotes expression from a specific TAL1 promoter. Expression from a specific promoter gives rise to a unique 5' UTR with differential regulation of translation. Moreover, our study suggests that the enhancers regulate TAL1 exon 3 alternative splicing by inducing changes in the chromatin at the splice site, which we demonstrate is mediated by KMT2B. Furthermore, our results indicate that TAL1-short binds more strongly to TAL1 E-protein partners and functions as a stronger transcription factor than TAL1-long. Specifically TAL1-short has a unique transcription signature promoting apoptosis. Finally, when we expressed both isoforms in mice bone marrow, we found that while overexpression of both isoforms prevents lymphoid differentiation, expression of TAL-short alone leads to hematopoietic stem cell exhaustion. Furthermore, we found that TAL1-short promoted erythropoiesis and reduced cell survival in the CML cell line K562. While TAL1 and its partners are considered promising therapeutic targets in the treatment of T-ALL, our results show that TAL1-short could act as a tumor suppressor and suggest that altering TAL1 isoform's ratio could be a preferred therapeutic approach.

Pathway-based integration of multi-omics data reveals lipidomics alterations validated in an Alzheimer's disease mouse model and risk loci carriers

Alzheimer's disease (AD) is a highly prevalent neurodegenerative disorder. Despite increasing evidence of the importance of metabolic dysregulation in AD, the underlying metabolic changes that may impact amyloid plaque formation are not understood, particularly for late-onset AD. This study analyzed genome-wide association studies (GWAS), transcriptomics, and proteomics data obtained from several data repositories to obtain differentially expressed (DE) multi-omics elements in mouse models of AD. We characterized the metabolic modulation in these data sets using gene ontology, transcription factor, pathway, and cell-type enrichment analyses. A predicted lipid signature was extracted from genome-scale metabolic networks (GSMN) and subsequently validated in a lipidomic data set derived from cortical tissue of ABCA-7 null mice, a mouse model of one of the genes associated with late-onset AD. Moreover, a metabolome-wide association study (MWAS) was performed to further characterize the association between dysregulated lipid metabolism in human blood serum and genes associated with AD risk. We found 203 DE transcripts, 164 DE proteins, and 58 DE GWAS-derived mouse orthologs associated with significantly enriched metabolic biological processes. Lipid and bioenergetic metabolic pathways were significantly over-represented across the AD multi-omics data sets. Microglia and astrocytes were significantly enriched in the lipid-predominant AD-metabolic transcriptome. We also extracted a predicted lipid signature that was validated and robustly modeled class separation in the ABCA7 mice cortical lipidome, with 11 of these lipid species exhibiting statistically significant modulations. MWAS revealed 298 AD single nucleotide polymorphisms-metabolite associations, of which 70% corresponded to lipid classes. These results support the importance of lipid metabolism dysregulation in AD and highlight the suitability of mapping AD multi-omics data into GSMNs to identify metabolic alterations.

Pathway-based integration of multi-omics data reveals lipidomics alterations validated in an Alzheimer's disease mouse model and risk loci carriers

Alzheimer's disease (AD) is a highly prevalent neurodegenerative disorder. Despite increasing evidence of the importance of metabolic dysregulation in AD, the underlying metabolic changes that may impact amyloid plaque formation are not understood, particularly for late-onset AD. This study analyzed genome-wide association studies (GWAS), transcriptomics, and proteomics data obtained from several data repositories to obtain differentially expressed (DE) multi-omics elements in mouse models of AD. We characterized the metabolic modulation in these data sets using gene ontology, transcription factor, pathway, and cell-type enrichment analyses. A predicted lipid signature was extracted from genome-scale metabolic networks (GSMN) and subsequently validated in a lipidomic data set derived from cortical tissue of ABCA-7 null mice, a mouse model of one of the genes associated with late-onset AD. Moreover, a metabolome-wide association study (MWAS) was performed to further characterize the association between dysregulated lipid metabolism in human blood serum and genes associated with AD risk. We found 203 DE transcripts, 164 DE proteins, and 58 DE GWAS-derived mouse orthologs associated with significantly enriched metabolic biological processes. Lipid and bioenergetic metabolic pathways were significantly over-represented across the AD multi-omics data sets. Microglia and astrocytes were significantly enriched in the lipid-predominant AD-metabolic transcriptome. We also extracted a predicted lipid signature that was validated and robustly modeled class separation in the ABCA7 mice cortical lipidome, with 11 of these lipid species exhibiting statistically significant modulations. MWAS revealed 298 AD single nucleotide polymorphisms-metabolite associations, of which 70% corresponded to lipid classes. These results support the importance of lipid metabolism dysregulation in AD and highlight the suitability of mapping AD multi-omics data into GSMNs to identify metabolic alterations.

Cell landscape of larval and adult Xenopus laevis at single-cell resolution

The rapid development of high-throughput single-cell RNA sequencing technology offers a good opportunity to dissect cell heterogeneity of animals. A large number of organism-wide single-cell atlases have been constructed for vertebrates such as Homo sapiens, Macaca fascicularis, Mus musculus and Danio rerio. However, an intermediate taxon that links mammals to vertebrates of more ancient origin is still lacking. Here, we construct the first Xenopus cell landscape to date, including larval and adult organs. Common cell lineage-specific transcription factors have been identified in vertebrates, including fish, amphibians and mammals. The comparison of larval and adult erythrocytes identifies stage-specific hemoglobin subtypes, as well as a common type of cluster containing both larval and adult hemoglobin, mainly at NF59. In addition, cell lineages originating from all three layers exhibits both antigen processing and presentation during metamorphosis, indicating a common regulatory mechanism during metamorphosis. Overall, our study provides a large-scale resource for research on Xenopus metamorphosis and adult organs.

Differential Etv2 threshold requirement for endothelial and erythropoietic development

Endothelial and erythropoietic lineages arise from a common developmental progenitor. Etv2 is a master transcriptional regulator required for the development of both lineages. However, the mechanisms through which Etv2 initiates the gene-regulatory networks (GRNs) for endothelial and erythropoietic specification and how the two GRNs diverge downstream of Etv2 remain incompletely understood. Here, by analyzing a hypomorphic Etv2 mutant, we demonstrate different threshold requirements for initiation of the downstream GRNs for endothelial and erythropoietic development. We show that Etv2 functions directly in a coherent feedforward transcriptional network for vascular endothelial development, and a low level of Etv2 expression is sufficient to induce and sustain the endothelial GRN. In contrast, Etv2 induces the erythropoietic GRN indirectly via activation of Tal1, which requires a significantly higher threshold of Etv2 to initiate and sustain erythropoietic development. These results provide important mechanistic insight into the divergence of the endothelial and erythropoietic lineages.

Epigenetic and Transcriptional Control of Erythropoiesis

Erythropoiesis is a process of enormous magnitude, with the average person generating two to three million red cells every second. Erythroid progenitors start as large cells with large nuclei, and over the course of three to four cell divisions they undergo a dramatic decrease in cell size accompanied by profound nuclear condensation, which culminates in enucleation. As maturing erythroblasts are undergoing these dramatic phenotypic changes, they accumulate hemoglobin and express high levels of other erythroid-specific genes, while silencing much of the non-erythroid transcriptome. These phenotypic and gene expression changes are associated with distinct changes in the chromatin landscape, and require close coordination between transcription factors and epigenetic regulators, as well as precise regulation of RNA polymerase II activity. Disruption of these processes are associated with inherited anemias and myelodysplastic syndromes. Here, we review the epigenetic mechanisms that govern terminal erythroid maturation, and their role in human disease.

Erythropoietin Non-hematopoietic Tissue Response and Regulation of Metabolism During Diet Induced Obesity

Erythropoietin (EPO) receptor (EPOR) determines EPO response. High level EPOR on erythroid progenitor cells gives rise to EPO regulated production of red blood cells. Animal models provide evidence for EPO activity in non-hematopoietic tissue mediated by EPOR expression. Beyond erythropoiesis, EPO activity includes neuroprotection in brain ischemia and trauma, endothelial nitric oxide production and cardioprotection, skeletal muscle wound healing, and context dependent bone remodeling affecting bone repair or bone loss. This review highlights examples of EPO protective activity in select non-hematopoietic tissue with emphasis on metabolic response mediated by EPOR expression in fat and brain and sex-specific regulation of fat mass and inflammation associated with diet induced obesity. Endogenous EPO maintains glucose and insulin tolerance and protects against fat mass accumulation and inflammation. Accompanying the increase in erythropoiesis with EPO treatment is improved glucose tolerance and insulin response. During high fat diet feeding, EPO also decreases fat mass accumulation in male mice. The increased white adipose tissue inflammation and macrophage infiltration associated with diet induced obesity are also reduced with EPO treatment with a shift toward an anti-inflammatory state and decreased inflammatory cytokine production. In female mice the protective effect of estrogen against obesity supersedes EPO regulation of fat mass and inflammation, and requires estrogen receptor alpha activity. In brain, EPOR expression in the hypothalamus localizes to proopiomelanocortin neurons in the arcuate nucleus that promotes a lean phenotype. EPO stimulation of proopiomelanocortin neurons increases STAT3 signaling and production of proopiomelanocortin. Cerebral EPO contributes to metabolic response, and elevated brain EPO reduces fat mass and hypothalamus inflammation during diet induced obesity in male mice without affecting EPO stimulated erythropoiesis. Ovariectomy abrogates the sex-specific metabolic response of brain EPO. The sex-dimorphic EPO metabolic response associated with fat mass accumulation and inflammation during diet induced obesity provide evidence for crosstalk between estrogen and EPO in their anti-obesity potential in female mice mediated in part via tissue specific response in brain and white adipose tissue. Endogenous and exogenous EPO response in non-hematopoietic tissue demonstrated in animal models suggests additional activity by which EPO treatment may affect human health beyond increased erythropoiesis.

Mechanisms of Binding Specificity among bHLH Transcription Factors

The transcriptome of every cell is orchestrated by the complex network of interaction between transcription factors (TFs) and their binding sites on DNA. Disruption of this network can result in many forms of organism malfunction but also can be the substrate of positive natural selection. However, understanding the specific determinants of each of these individual TF-DNA interactions is a challenging task as it requires integrating the multiple possible mechanisms by which a given TF ends up interacting with a specific genomic region. These mechanisms include DNA motif preferences, which can be determined by nucleotide sequence but also by DNA’s shape; post-translational modifications of the TF, such as phosphorylation; and dimerization partners and co-factors, which can mediate multiple forms of direct or indirect cooperative binding. Binding can also be affected by epigenetic modifications of putative target regions, including DNA methylation and nucleosome occupancy. In this review, we describe how all these mechanisms have a role and crosstalk in one specific family of TFs, the basic helix-loop-helix (bHLH), with a very conserved DNA binding domain and a similar DNA preferred motif, the E-box. Here, we compile and discuss a rich catalog of strategies used by bHLH to acquire TF-specific genome-wide landscapes of binding sites.

A self-attention model for inferring cooperativity between regulatory features

Deep learning has demonstrated its predictive power in modeling complex biological phenomena such as gene expression. The value of these models hinges not only on their accuracy, but also on the ability to extract biologically relevant information from the trained models. While there has been much recent work on developing feature attribution methods that discover the most important features for a given sequence, inferring cooperativity between regulatory elements, which is the hallmark of phenomena such as gene expression, remains an open problem. We present SATORI, a Self-ATtentiOn based model to detect Regulatory element Interactions. Our approach combines convolutional layers with a self-attention mechanism that helps us capture a global view of the landscape of interactions between regulatory elements in a sequence. A comprehensive evaluation demonstrates the ability of SATORI to identify numerous statistically significant TF-TF interactions, many of which have been previously reported. Our method is able to detect higher numbers of experimentally verified TF-TF interactions than existing methods, and has the advantage of not requiring a computationally expensive post-processing step. Finally, SATORI can be used for detection of any type of feature interaction in models that use a similar attention mechanism, and is not limited to the detection of TF-TF interactions.

A Positive Regulatory Feedback Loop between EKLF/KLF1 and TAL1/SCL Sustaining the Erythropoiesis

The erythroid Krüppel-like factor EKLF/KLF1 is a hematopoietic transcription factor binding to the CACCC DNA motif and participating in the regulation of erythroid differentiation. With combined use of microarray-based gene expression profiling and the promoter-based ChIP-chip assay of E14.5 fetal liver cells from wild type (WT) and EKLF-knockout (Eklf^−/−) mouse embryos, we identified the pathways and direct target genes activated or repressed by EKLF. This genome-wide study together with the molecular/cellular analysis of the mouse erythroleukemic cells (MEL) indicate that among the downstream direct target genes of EKLF is Tal1/Scl. Tal1/Scl encodes another DNA-binding hematopoietic transcription factor TAL1/SCL, known to be an Eklf activator and essential for definitive erythroid differentiation. Further identification of the authentic Tal gene promoter in combination with the in vivo genomic footprinting approach and DNA reporter assay demonstrate that EKLF activates the Tal gene through binding to a specific CACCC motif located in its promoter. These data establish the existence of a previously unknow positive regulatory feedback loop between two DNA-binding hematopoietic transcription factors, which sustains mammalian erythropoiesis.

A gain-of-function single nucleotide variant creates a new promoter which acts as an orientation-dependent enhancer-blocker

Many single nucleotide variants (SNVs) associated with human traits and genetic diseases are thought to alter the activity of existing regulatory elements. Some SNVs may also create entirely new regulatory elements which change gene expression, but the mechanism by which they do so is largely unknown. Here we show that a single base change in an otherwise unremarkable region of the human α-globin cluster creates an entirely new promoter and an associated unidirectional transcript. This SNV downregulates α-globin expression causing α-thalassaemia. Of note, the new promoter lying between the α-globin genes and their associated super-enhancer disrupts their interaction in an orientation-dependent manner. Together these observations show how both the order and orientation of the fundamental elements of the genome determine patterns of gene expression and support the concept that active genes may act to disrupt enhancer-promoter interactions in mammals as in Drosophila. Finally, these findings should prompt others to fully evaluate SNVs lying outside of known regulatory elements as causing changes in gene expression by creating new regulatory elements.

The role of promoters as potential insulator elements has been largely unexplored in mammals. Here the authors show that a single nucleotide variant in the α-globin locus forms a new promoter and acts as an orientation-dependent enhancer-blocking insulator element.

Integrative Single-Cell RNA-Seq and ATAC-Seq Analysis of Human Developmental Hematopoiesis

Regulation of hematopoiesis during human development remains poorly defined. Here we applied single-cell RNA sequencing (scRNA-seq) and single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq) to over 8,000 human immunophenotypic blood cells from fetal liver and bone marrow. We inferred their differentiation trajectory and identified three highly proliferative oligopotent progenitor populations downstream of hematopoietic stem cells (HSCs)/multipotent progenitors (MPPs). Along this trajectory, we observed opposing patterns of chromatin accessibility and differentiation that coincided with dynamic changes in the activity of distinct lineage-specific transcription factors. Integrative analysis of chromatin accessibility and gene expression revealed extensive epigenetic but not transcriptional priming of HSCs/MPPs prior to their lineage commitment. Finally, we refined and functionally validated the sorting strategy for the HSCs/MPPs and achieved around 90% enrichment. Our study provides a useful framework for future investigation of human developmental hematopoiesis in the context of blood pathologies and regenerative medicine.

The transcriptional regulation of normal and malignant blood cell development

Development of multicellular organisms requires the differential usage of our genetic information to change one cell fate into another. This process drives the appearance of different cell types that come together to form specialized tissues sustaining a healthy organism. In the last decade, by moving away from studying single genes toward a global view of gene expression control, a revolution has taken place in our understanding of how genes work together and how cells communicate to translate the information encoded in the genome into a body plan. The development of hematopoietic cells has long served as a paradigm of development in general. In this review, we highlight how transcription factors and chromatin components work together to shape the gene regulatory networks controlling gene expression in the hematopoietic system and to drive blood cell differentiation. In addition, we outline how this process goes astray in blood cancers. We also touch upon emerging concepts that place these processes firmly into their associated subnuclear structures adding another layer of the control of differential gene expression.

Gene Expression Analysis of Mevalonate Kinase Deficiency Affected Children Identifies Molecular Signatures Related to Hematopoiesis

Mevalonate kinase deficiency (MKD) is a rare autoinflammatory genetic disorder characterized by recurrent fever attacks and systemic inflammation with potentially severe complications. Although it is recognized that the lack of protein prenylation consequent to mevalonate pathway blockade drives IL1β hypersecretion, and hence autoinflammation, MKD pathogenesis and the molecular mechanisms underlaying most of its clinical manifestations are still largely unknown. In this study, we performed a comprehensive bioinformatic analysis of a microarray dataset of MKD patients, using gene ontology and Ingenuity Pathway Analysis (IPA) tools, in order to identify the most significant differentially expressed genes and infer their predicted relationships into biological processes, pathways, and networks. We found that hematopoiesis linked biological functions and pathways are predominant in the gene ontology of differentially expressed genes in MKD, in line with the observed clinical feature of anemia. We also provided novel information about the molecular mechanisms at the basis of the hematological abnormalities observed, that are linked to the chronic inflammation and to defective prenylation. Considering the broad and unspecific spectrum of MKD clinical manifestations and the difficulty in its diagnosis, a better understanding of MKD molecular bases could be translated to the clinical level to facilitate diagnosis, and improve management and therapy.

The transcription factor TAL1 and miR-17-92 create a regulatory loop in hematopoiesis

A network of gene regulatory factors such as transcription factors and microRNAs establish and maintain gene expression patterns during hematopoiesis. In this network, transcription factors regulate each other and are involved in regulatory loops with microRNAs. The microRNA cluster miR-17-92 is located within the MIR17HG gene and encodes six mature microRNAs. It is important for hematopoietic differentiation and plays a central role in malignant disease. However, the transcription factors downstream of miR-17-92 are largely elusive and the transcriptional regulation of miR-17-92 is not fully understood. Here we show that miR-17-92 forms a regulatory loop with the transcription factor TAL1. The miR-17-92 cluster inhibits expression of TAL1 and indirectly leads to decreased stability of the TAL1 transcriptional complex. We found that TAL1 and its heterodimerization partner E47 regulate miR-17-92 transcriptionally. Furthermore, miR-17-92 negatively influences erythroid differentiation, a process that depends on gene activation by the TAL1 complex. Our data give example of how transcription factor activity is fine-tuned during normal hematopoiesis. We postulate that disturbance of the regulatory loop between TAL1 and the miR-17-92 cluster could be an important step in cancer development and progression.

Erythropoietin regulates metabolic response in mice via receptor expression in adipose tissue, brain, and bone

Erythropoietin (EPO) acts by binding to erythroid progenitor cells to regulate red blood cell production. While EPO receptor (Epor) expression is highest on erythroid tissue, animal models exhibit EPO activity in nonhematopoietic tissues, mediated, in part, by tissue-specific Epor expression. This review describes the metabolic response in mice to endogenous EPO and EPO treatment associated with glucose metabolism, fat mass accumulation, and inflammation in white adipose tissue and brain during diet-induced obesity and with bone marrow fat and bone remodeling. During high-fat diet-induced obesity, EPO treatment improves glucose tolerance, decreases fat mass accumulation, and shifts white adipose tissue from a pro-inflammatory to an anti-inflammatory state. Fat mass regulation by EPO is sex dimorphic, apparent in males and abrogated by estrogen in females. Cerebral EPO also regulates fat mass and hypothalamus inflammation associated with diet-induced obesity in males and ovariectomized female mice. In bone, EPO contributes to the balance between adipogenesis and osteogenesis in both male and female mice. EPO treatment promotes bone loss mediated via Epor in osteoblasts and reduces bone marrow adipocytes before and independent of change in white adipose tissue fat mass. EPO regulation of bone loss and fat mass is independent of EPO-stimulated erythropoiesis. EPO nonhematopoietic tissue response may relate to the long-term consequences of EPO treatment of anemia in chronic kidney disease and to the alternative treatment of oral hypoxia-inducible factor prolyl hydroxylase inhibitors that increase endogenous EPO production.

A Unique Epigenomic Landscape Defines Human Erythropoiesis

Mammalian erythropoiesis yields a highly specialized cell type, the mature erythrocyte, evolved to meet the organismal needs of increased oxygen-carrying capacity. To better understand the regulation of erythropoiesis, we performed genome-wide studies of chromatin accessibility, DNA methylation, and transcriptomics using a recently developed strategy to obtain highly purified populations of primary human erythroid cells. The integration of gene expression, DNA methylation, and chromatin state dynamics reveals that stage-specific gene regulation during erythropoiesis is a stepwise and hierarchical process involving many cis-regulatory elements. Erythroid-specific, nonpromoter sites of chromatin accessibility are linked to erythroid cell phenotypic variation and inherited disease. Comparative analyses of stage-specific chromatin accessibility indicate that there is limited early chromatin priming of erythroid genes during hematopoiesis. The epigenome of terminally differentiating erythroid cells defines a distinct subset of highly specialized cells that are vastly dissimilar from other hematopoietic and nonhematopoietic cell types. These epigenomic and transcriptome data are powerful tools to study human erythropoiesis.

Schulz et al. use genome-wide studies of chromatin accessibility, DNA methylation, and transcriptomes in primary human erythroid cells to reveal important characteristics of erythropoiesis. Chromatin accessibility of terminal erythroid differentiation is markedly dissimilar from other hematopoietic cell types. Epigenomic changes are linked to erythroid cell traits and disease genes.

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.

The Many Facets of Erythropoietin Physiologic and Metabolic Response

In mammals, erythropoietin (EPO), produced in the kidney, is essential for bone marrow erythropoiesis, and hypoxia induction of EPO production provides for the important erythropoietic response to ischemic stress, such as during blood loss and at high altitude. Erythropoietin acts by binding to its cell surface receptor which is expressed at the highest level on erythroid progenitor cells to promote cell survival, proliferation, and differentiation in production of mature red blood cells. In addition to bone marrow erythropoiesis, EPO causes multi-tissue responses associated with erythropoietin receptor (EPOR) expression in non-erythroid cells such neural cells, endothelial cells, and skeletal muscle myoblasts. Animal and cell models of ischemic stress have been useful in elucidating the potential benefit of EPO affecting maintenance and repair of several non-hematopoietic organs including brain, heart and skeletal muscle. Metabolic and glucose homeostasis are affected by endogenous EPO and erythropoietin administration affect, in part via EPOR expression in white adipose tissue. In diet-induced obese mice, EPO is protective for white adipose tissue inflammation and gives rise to a gender specific response in weight control associated with white fat mass accumulation. Erythropoietin regulation of fat mass is masked in female mice due to estrogen production. EPOR is also expressed in bone marrow stromal cells (BMSC) and EPO administration in mice results in reduced bone independent of the increase in hematocrit. Concomitant reduction in bone marrow adipocytes and bone morphogenic protein suggests that high EPO inhibits adipogenesis and osteogenesis. These multi-tissue responses underscore the pleiotropic potential of the EPO response and may contribute to various physiological manifestations accompanying anemia or ischemic response and pharmacological uses of EPO.

Regulation of GATA1 levels in erythropoiesis

GATA1 is considered as the "master" transcription factor in erythropoiesis. It regulates at the transcriptional level all aspects of erythroid maturation and function, as revealed by gene knockout studies in mice and by genome-wide occupancies in erythroid cells. The GATA1 protein contains two zinc finger domains and an N-terminal transactivation domain. GATA1 translation results in the production of the full-length protein and of a shorter variant (GATA1s) lacking the N-terminal transactivation domain, which is functionally deficient in supporting erythropoiesis. GATA1 protein abundance is highly regulated in erythroid cells at different levels, including transcription, mRNA translation, posttranslational modifications, and protein degradation, in a differentiation-stage-specific manner. Maintaining high GATA1 protein levels is essential in the early stages of erythroid maturation, whereas downregulating GATA1 protein levels is a necessary step in terminal erythroid differentiation. The importance of maintaining proper GATA1 protein homeostasis in erythropoiesis is demonstrated by the fact that both GATA1 loss and its overexpression result in lethal anemia. Importantly, alterations in any of those GATA1 regulatory checkpoints have been recognized as an important cause of hematological disorders such as dyserythropoiesis (with or without thrombocytopenia), β-thalassemia, Diamond-Blackfan anemia, myelodysplasia, or leukemia. In this review, we provide an overview of the multilevel regulation of GATA1 protein homeostasis in erythropoiesis and of its deregulation in hematological disease.

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

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

Discovering epistatic feature interactions from neural network models of regulatory DNA sequences

Motivation

Transcription factors bind regulatory DNA sequences in a combinatorial manner to modulate gene expression. Deep neural networks (DNNs) can learn the cis-regulatory grammars encoded in regulatory DNA sequences associated with transcription factor binding and chromatin accessibility. Several feature attribution methods have been developed for estimating the predictive importance of individual features (nucleotides or motifs) in any input DNA sequence to its associated output prediction from a DNN model. However, these methods do not reveal higher-order feature interactions encoded by the models.

Results

We present a new method called Deep Feature Interaction Maps (DFIM) to efficiently estimate interactions between all pairs of features in any input DNA sequence. DFIM accurately identifies ground truth motif interactions embedded in simulated regulatory DNA sequences. DFIM identifies synergistic interactions between GATA1 and TAL1 motifs from in vivo TF binding models. DFIM reveals epistatic interactions involving nucleotides flanking the core motif of the Cbf1 TF in yeast from in vitro TF binding models. We also apply DFIM to regulatory sequence models of in vivo chromatin accessibility to reveal interactions between regulatory genetic variants and proximal motifs of target TFs as validated by TF binding quantitative trait loci. Our approach makes significant strides in improving the interpretability of deep learning models for genomics.

Availability and implementation

Code is available at: https://github.com/kundajelab/dfim.

Supplementary information

Supplementary data are available at Bioinformatics online.

CUT&RUNTools: a flexible pipeline for CUT&RUN processing and footprint analysis

We introduce CUT&RUNTools as a flexible, general pipeline for facilitating the identification of chromatin-associated protein binding and genomic footprinting analysis from antibody-targeted CUT&RUN primary cleavage data. CUT&RUNTools extracts endonuclease cut site information from sequences of short-read fragments and produces single-locus binding estimates, aggregate motif footprints, and informative visualizations to support the high-resolution mapping capability of CUT&RUN. CUT&RUNTools is available at https://bitbucket.org/qzhudfci/cutruntools/ .

Translational regulation and deregulation in erythropoiesis

Translational regulation plays a critical role in erythropoiesis, as it reflects the translational needs of enucleated mature erythroid cells in the absence of transcription and the large translational demands of balanced globin chain synthesis during erythroid maturation. In addition, red blood cells need to respond quickly to changes in their environment and the demands of the organism. Translational regulation occurs at several levels in erythroid cells, including the differential utilization of upstream open reading frames during differentiation and in response to signaling and the employment of RNA-binding proteins in an erythroid cell-specific fashion. Translation initiation is a critical juncture for translational regulation in response to environmental signals such as heme and iron availability, whereas regulatory mechanisms for ribosome recycling are consistent with recent observations highlighting the importance of maintaining adequate ribosome levels in differentiating erythroid cells. Translational deregulation in erythroid cells leads to disease associated with ineffective erythropoiesis, further highlighting the pivotal role translational regulation in erythropoiesis plays in human physiology and homeostasis. Overall, erythropoiesis has served as a unique model that has provided invaluable insight into translational regulation.

SCL/TAL1 cooperates with Polycomb RYBP-PRC1 to suppress alternative lineages in blood-fated cells

During development, it is unclear if lineage-fated cells derive from multilineage-primed progenitors and whether active mechanisms operate to restrict cell fate. Here we investigate how mesoderm specifies into blood-fated cells. We document temporally restricted co-expression of blood (Scl/Tal1), cardiac (Mesp1) and paraxial (Tbx6) lineage-affiliated transcription factors in single cells, at the onset of blood specification, supporting the existence of common progenitors. At the same time-restricted stage, absence of SCL results in expansion of cardiac/paraxial cell populations and increased cardiac/paraxial gene expression, suggesting active suppression of alternative fates. Indeed, SCL normally activates expression of co-repressor ETO2 and Polycomb-PRC1 subunits (RYBP, PCGF5) and maintains levels of Polycomb-associated histone marks (H2AK119ub/H3K27me3). Genome-wide analyses reveal ETO2 and RYBP co-occupy most SCL target genes, including cardiac/paraxial loci. Reduction of Eto2 or Rybp expression mimics Scl-null cardiac phenotype. Therefore, SCL-mediated transcriptional repression prevents mis-specification of blood-fated cells, establishing active repression as central to fate determination processes.

Mechanisms that operate during embryonic development to restrict cell fate are currently under investigation. Here the authors characterise the role of SCL/TAL1 at the onset of blood specification in embryonic development using mouse EB differentiation culture as a model system.

A novel role for Lyl1 in primitive erythropoiesis

Stem cell leukemia (Scl or Tal1) and lymphoblastic leukemia 1 (Lyl1) encode highly related members of the basic helix-loop-helix family of transcription factors that are co-expressed in the erythroid lineage. Previous studies have suggested that Scl is essential for primitive erythropoiesis. However, analysis of single-cell RNA-seq data of early embryos showed that primitive erythroid cells express both Scl and Lyl1 Therefore, to determine whether Lyl1 can function in primitive erythropoiesis, we crossed conditional Scl knockout mice with mice expressing a Cre recombinase under the control of the Epo receptor, active in erythroid progenitors. Embryos with 20% expression of Scl from E9.5 survived to adulthood. However, mice with reduced expression of Scl and absence of Lyl1 (double knockout; DKO) died at E10.5 because of progressive loss of erythropoiesis. Gene expression profiling of DKO yolk sacs revealed loss of Gata1 and many of the known target genes of the SCL-GATA1 complex. ChIP-seq analyses in a human erythroleukemia cell line showed that LYL1 exclusively bound a small subset of SCL targets including GATA1. Together, these data show for the first time that Lyl1 can maintain primitive erythropoiesis.

TCF4 promotes erythroid development

Transcription factor 4 (TCF4) is implicated in lymphoid cell differentiation and its expression predicts outcome in acute myeloid leukemia. Here, we investigated the role of TCF4 in myelopoiesis. Overexpression of TCF4 (TCF4^OE) in umbilical cord blood (UCB) cells resulted in a twofold increase in erythroid colony forming units (CFU-Es), whereas knock-down (KD) of TCF4 (TCF4^KD) caused a dramatic decrease in the number of erythroid colonies. In megakaryocyte CFUs (CFU-MKs), both TCF4^KD and TCF4^OE inhibited MK colony formation. TCF4 did not have an impact on granulocyte, macrophage, or granulocyte-macrophage colonies or on the proportion of MK-erythrocyte progenitors (MEPs) in culture. Because TCF4 affects erythroid/MK development and these lineages are affected in myelodysplastic syndrome (MDS), we studied the impact of TCF4 expression in this disease. MDS patients with high (≥median) TCF4 mRNA expression had higher hemoglobin (Hb) levels than MDS patients with low TCF4 expression (mean 9.0 vs. 8.55 g/dL, p = 0.02). Overall, TCF4 mRNA expression was lower in hematopoietic stem cells, common myeloid progenitors, and MEPs from MDS patients, but not in granulocyte-macrophage progenitors, compared with healthy controls. Therefore, in cell fractions with erythroid lineage potential, TCF4 is expressed less in MDS patients than in healthy controls. This correlates with the low overall Hb levels seen in MDS patients compared with healthy individuals and is consistent with the positive impact of TCF4 on erythroid development while not having impact on white colonies. These results indicate a role for TCF4 as a novel factor in erythroid-megakaryocytic differentiation.

Oncogenic transcriptional program driven by TAL1 in T-cell acute lymphoblastic leukemia

TAL1/SCL is a prime example of an oncogenic transcription factor that is abnormally expressed in acute leukemia due to the replacement of regulator elements. This gene has also been recognized as an essential regulator of hematopoiesis. TAL1 expression is strictly regulated in a lineage- and stage-specific manner. Such precise control is crucial for the switching of the transcriptional program. The misexpression of TAL1 in immature thymocytes leads to a widespread series of orchestrated downstream events that affect several different cellular machineries, resulting in a lethal consequence, namely T-cell acute lymphoblastic leukemia (T-ALL). In this article, we will discuss the transcriptional regulatory network and downstream target genes, including protein-coding genes and non-coding RNAs, controlled by TAL1 in normal hematopoiesis and T-cell leukemogenesis.

In utero gene expression in the Slc39a8(neo/neo) knockdown mouse

Slc39a8 encodes ZIP8, a divalent cation/bicarbonate symporter expressed in pluripotent mouse embryonic stem cells, and therefore ubiquitous in adult tissues; ZIP8 influxes Zn²⁺, Mn²⁺ and Fe²⁺. Slc39a8(neo/neo) knockdown mice exhibit 10-15% of wild-type ZIP8 mRNA and protein levels, and show pleiotropic phenotype of stunted growth, neonatal lethality, multi-organ dysmorphogenesis, and dysregulated hematopoiesis manifested as severe anemia. Herein we performed RNA-seq analysis of gestational day (GD)13.5 yolk sac and placenta, and GD16.5 liver, kidney, lung, heart and cerebellum, comparing Slc39a8(neo/neo) with Slc39a8(+/+) wild-type. Meta-data analysis of differentially-expressed genes revealed 29 unique genes from all tissues - having enriched GO categories associated with hematopoiesis and hypoxia and KEGG categories of complement, response to infection, and coagulation cascade - consistent with dysregulated hematopoietic stem cell fate. Based on transcription factor (TF) profiles in the JASPAR database, and searching for TF-binding sites enriched by Pscan, we identified numerous genes encoding zinc-finger and other TFs associated with hematopoietic stem cell functions. We conclude that, in this mouse model, deficient ZIP8-mediated divalent cation transport affects zinc-finger (e.g. GATA proteins) and other TFs interacting with GATA proteins (e.g. TAL1), predominantly in yolk sac. These data strongly support the phenotype of dysmorphogenesis and anemia seen in Slc39a8(neo/neo) mice in utero.

Evolution of hemoglobin loci and their regulatory elements

Across the expanse of vertebrate evolution, each species produces multiple forms of hemoglobin in erythroid cells at appropriate times and in the proper amounts. The multiple hemoglobins are encoded in two globin gene clusters in almost all species. One globin gene cluster, linked to the gene NPRL3, is preserved in all vertebrates, including a gene cluster encoding the highly divergent globins from jawless vertebrates. This preservation of synteny may reflect the presence of a powerful enhancer of globin gene expression in the NPRL3 gene. Despite substantial divergence in noncoding DNA sequences among mammals, several epigenetic features of the globin gene regulatory regions are preserved across vertebrates. The preserved features include multiple DNase hypersensitive sites, at least one of which is an enhancer, and binding by key lineage-restricted transcription factors such as GATA1 and TAL1, which in turn recruit coactivators such as P300 that catalyze acetylation of histones. The maps of epigenetic features are strongly correlated with activity in gene regulation, and resources for accessing and visualizing such maps are readily available to the community of researchers and students.

The LDB1 Complex Co-opts CTCF for Erythroid Lineage-Specific Long-Range Enhancer Interactions

Lineage-specific transcription factors are critical for long-range enhancer interactions, but direct or indirect contributions of architectural proteins such as CCCTC-binding factor (CTCF) to enhancer function remain less clear. The LDB1 complex mediates enhancer-gene interactions at the β-globin locus through LDB1 self-interaction. We find that an LDB1-bound enhancer upstream of carbonic anhydrase 2 (Car2) activates its expression by interacting directly with CTCF at the gene promoter. Both LDB1 and CTCF are required for enhancer-Car2 looping, and the domain of LDB1 contacted by CTCF is necessary to rescue Car2 transcription in LDB1-deficient cells. Genome-wide studies and CRISPR/Cas9 genome editing indicate that LDB1-CTCF enhancer looping underlies activation of a substantial fraction of erythroid genes. Our results provide a mechanism by which long-range interactions of architectural protein CTCF can be tailored to achieve a tissue-restricted pattern of chromatin loops and gene expression.

Epigenetic regulation of megakaryocytic and erythroid differentiation by PHF2 histone demethylase

Plant homeodomain finger 2 (PHF2) is a JmjC family histone demethylase that demethylates H3K9me2, a repressive gene marker. PHF2 was found to play a role in the differentiation of several tissue types such as osteoblast and adipocyte differentiation. We report here that PHF2 plays a role in the epigenetic regulation of megakaryocytic (MK) and erythroid differentiation. We investigated PHF2 expression during MK and erythroid differentiation in K562 and human CD34⁺ progenitor (hCD34⁺ ) cells. Our data demonstrate that PHF2 expression is down-regulated during megakaryopoiesis and erythropoiesis. PHF2 has a negative role in MK and erythroid differentiation of K562 cells; knockdown of PHF2 promotes MK and erythroid differentiation of hCD34⁺ cells. Similarly, we found that p53 expression is also down-regulated during MK and erythroid differentiation, which parallels PHF2 expression. PHF2 binds to the p53 promoter and regulates the expression of p53 by demethylating H3K9me2 in the promoter region of p53. Taken together, our data show that PHF2 is a negative epigenetic regulator of MK and erythroid differentiation, and that one of the pathways through which PHF2 affects MK and erythroid differentiation is via regulation of p53 expression.

HEDD: Human Enhancer Disease Database

Enhancers, as specialized genomic cis-regulatory elements, activate transcription of their target genes and play an important role in pathogenesis of many human complex diseases. Despite recent systematic identification of them in the human genome, currently there is an urgent need for comprehensive annotation databases of human enhancers with a focus on their disease connections. In response, we built the Human Enhancer Disease Database (HEDD) to facilitate studies of enhancers and their potential roles in human complex diseases. HEDD currently provides comprehensive genomic information for ∼2.8 million human enhancers identified by ENCODE, FANTOM5 and RoadMap with disease association scores based on enhancer-gene and gene-disease connections. It also provides Web-based analytical tools to visualize enhancer networks and score enhancers given a set of selected genes in a specific gene network. HEDD is freely accessible at http://zdzlab.einstein.yu.edu/1/hedd.php.

The Role of TAL1 in Hematopoiesis and Leukemogenesis

TAL1 (SCL/TAL1, T-cell acute leukemia protein 1) is a transcription factor that is involved in the process of hematopoiesis and leukemogenesis. It participates in blood cell formation, forms mesoderm in early embryogenesis, and regulates hematopoiesis in adult organisms. TAL1 is essential in maintaining the multipotency of hematopoietic stem cells (HSC) and keeping them in quiescence (stage G0). TAL1 forms complexes with various transcription factors, regulating hematopoiesis (E2A/HEB, GATA1-3, LMO1-2, Ldb1, ETO₂, RUNX1, ERG, FLI1). In these complexes, TAL1 regulates normal myeloid differentiation, controls the proliferation of erythroid progenitors, and determines the choice of the direction of HSC differentiation. The transcription factors TAL1, E2A, GATA1 (or GATA2), LMO2, and Ldb1 are the major components of the SCL complex. In addition to normal hematopoiesis, this complex may also be involved in the process of blood cell malignant transformation. Upregulation of C-KIT expression is one of the main roles played by the SCL complex. Today, TAL1 and its partners are considered promising therapeutic targets in the treatment of T-cell acute lymphoblastic leukemia.

Affected pathways and transcriptional regulators in gene expression response to an ultra-marathon trail: Global and independent activity approaches

Gene expression (GE) analyses on blood samples from marathon and half-marathon runners have reported significant impacts on the immune and inflammatory systems. An ultra-marathon trail (UMT) represents a greater effort due to its more testing conditions. For the first time, we report the genome-wide GE profiling in a group of 16 runners participating in an 82 km UMT competition. We quantified their differential GE profile before and after the race using HuGene2.0st microarrays (Affymetrix Inc., California, US). The results obtained were decomposed by means of an independent component analysis (ICA) targeting independent expression modes. We observed significant differences in the expression levels of 5,084 protein coding genes resulting in an overrepresentation of 14% of the human biological pathways from the Kyoto Encyclopedia of Genes and Genomes database. These were mainly clustered on terms related with protein synthesis repression, altered immune system and infectious diseases related mechanisms. In a second analysis, 27 out of the 196 transcriptional regulators (TRs) included in the Open Regulatory Annotation database were overrepresented. Among these TRs, we identified transcription factors from the hypoxia-inducible factors (HIF) family EPAS1 (p< 0.01) and HIF1A (p<0.001), and others jointly described in the gluconeogenesis program such as HNF4 (p< 0.001), EGR1 (p<0.001), CEBPA (p< 0.001) and a highly specific TR, YY1 (p<0.01). The five independent components, obtained from ICA, further revealed a down-regulation of 10 genes distributed in the complex I, III and V from the electron transport chain. This mitochondrial activity reduction is compatible with HIF-1 system activation. The vascular endothelial growth factor (VEGF) pathway, known to be regulated by HIF, also emerged (p<0.05). Additionally, and related to the brain rewarding circuit, the endocannabinoid signalling pathway was overrepresented (p<0.05).

Identification of Plagl1/Zac1 binding sites and target genes establishes its role in the regulation of extracellular matrix genes and the imprinted gene network

PLAGL1/ZAC1 undergoes parental genomic imprinting, is paternally expressed, and is a member of the imprinted gene network (IGN). It encodes a zinc finger transcription factor with anti-proliferative activity and is a candidate tumor suppressor gene on 6q24 whose expression is frequently lost in various neoplasms. Conversely, gain of PLAGL1 function is responsible for transient neonatal diabetes mellitus, a rare genetic disease that results from defective pancreas development. In the present work, we showed that Plagl1 up-regulation was not associated with DNA damage-induced cell cycle arrest. It was rather associated with physiological cell cycle exit that occurred with contact inhibition, growth factor withdrawal, or cell differentiation. To gain insights into Plagl1 mechanism of action, we identified Plagl1 target genes by combining chromatin immunoprecipitation and genome-wide transcriptomics in transfected cell lines. Plagl1-elicited gene regulation correlated with multiple binding to the proximal promoter region through a GC-rich motif. Plagl1 target genes included numerous genes involved in signaling, cell adhesion, and extracellular matrix composition, including collagens. Plagl1 targets also included 22% of the 409 genes that make up the IGN. Altogether, this work identified Plagl1 as a transcription factor that coordinated the regulation of a subset of IGN genes and controlled extracellular matrix composition.

Sasquatch: predicting the impact of regulatory SNPs on transcription factor binding from cell- and tissue-specific DNase footprints

In the era of genome-wide association studies (GWAS) and personalized medicine, predicting the impact of single nucleotide polymorphisms (SNPs) in regulatory elements is an important goal. Current approaches to determine the potential of regulatory SNPs depend on inadequate knowledge of cell-specific DNA binding motifs. Here, we present Sasquatch, a new computational approach that uses DNase footprint data to estimate and visualize the effects of noncoding variants on transcription factor binding. Sasquatch performs a comprehensive k-mer-based analysis of DNase footprints to determine any k-mer's potential for protein binding in a specific cell type and how this may be changed by sequence variants. Therefore, Sasquatch uses an unbiased approach, independent of known transcription factor binding sites and motifs. Sasquatch only requires a single DNase-seq data set per cell type, from any genotype, and produces consistent predictions from data generated by different experimental procedures and at different sequence depths. Here we demonstrate the effectiveness of Sasquatch using previously validated functional SNPs and benchmark its performance against existing approaches. Sasquatch is available as a versatile webtool incorporating publicly available data, including the human ENCODE collection. Thus, Sasquatch provides a powerful tool and repository for prioritizing likely regulatory SNPs in the noncoding genome.

TAL1 as a master oncogenic transcription factor in T-cell acute lymphoblastic leukemia

In hematopoietic cell development, the transcriptional program is strictly regulated in a lineage- and stage-specific manner that requires a number of transcription factors to work in a cascade or in a loop, in addition to interactions with nonhematopoietic cells in the microenvironment. Disruption of the transcriptional program alters the cellular state and may predispose cells to the acquisition of genetic abnormalities. Early studies have shown that proteins that promote cell differentiation often serve as tumor suppressors, whereas inhibitors of those proteins act as oncogenes in the context of acute leukemia. A prime example is T-cell acute lymphoblastic leukemia (T-ALL), a malignant disorder characterized by clonal proliferation of immature stage thymocytes. Although a relatively small number of genetic abnormalities are observed in T-ALL, these abnormalities are crucial for leukemogenesis. Many oncogenes and tumor suppressors in T-ALL are transcription factors that are required for normal hematopoiesis. The transformation process in T-ALL is efficient and orchestrated; the oncogene disrupts the transcriptional program directing T-cell differentiation and also uses its native ability as a master transcription factor in hematopoiesis. This imbalance in the transcriptional program is a primary determinant underlying the molecular pathogenesis of T-ALL. In this review, we focus on the oncogenic transcription factor TAL1 and the tumor-suppressor E-proteins and discuss the malignant cell state, the transcriptional circuit, and the consequence of molecular abnormalities in T-ALL.

Hematopoietic transcription factors and differential cofactor binding regulate PRKACB isoform expression

Hematopoietic differentiation is controlled by key transcription factors, which regulate stem cell functions and differentiation. TAL1 is a central transcription factor for hematopoietic stem cell development in the embryo and for gene regulation during erythroid/megakaryocytic differentiation. Knowledge of the target genes controlled by a given transcription factor is important to understand its contribution to normal development and disease. To uncover direct target genes of TAL1 we used high affinity streptavidin/biotin-based chromatin precipitation (Strep-CP) followed by Strep-CP on ChIP analysis using ChIP promoter arrays. We identified 451 TAL1 target genes in K562 cells. Furthermore, we analysed the regulation of one of these genes, the catalytic subunit beta of protein kinase A (PRKACB), during megakaryopoiesis of K562 and primary human CD34+ stem cell/progenitor cells. We found that TAL1 together with hematopoietic transcription factors RUNX1 and GATA1 binds to the promoter of the isoform 3 of PRKACB (Cβ3). During megakaryocytic differentiation a coactivator complex on the Cβ3 promoter, which includes WDR5 and p300, is replaced with a corepressor complex. In this manner, activating chromatin modifications are removed and expression of the PRKACB-Cβ3 isoform during megakaryocytic differentiation is reduced. Our data uncover a role of the TAL1 complex in controlling differential isoform expression of PRKACB. These results reveal a novel function of TAL1, RUNX1 and GATA1 in the transcriptional control of protein kinase A activity, with implications for cellular signalling control during differentiation and disease.

SCL/TAL1: a multifaceted regulator from blood development to disease

SCL/TAL1 (stem cell leukemia/T-cell acute lymphoblastic leukemia [T-ALL] 1) is an essential transcription factor in normal and malignant hematopoiesis. It is required for specification of the blood program during development, adult hematopoietic stem cell survival and quiescence, and terminal maturation of select blood lineages. Following ectopic expression, SCL contributes to oncogenesis in T-ALL. Remarkably, SCL’s activities are all mediated through nucleation of a core quaternary protein complex (SCL:E-protein:LMO1/2 [LIM domain only 1 or 2]:LDB1 [LIM domain-binding protein 1]) and dynamic recruitment of conserved combinatorial associations of additional regulators in a lineage- and stage-specific context. The finely tuned control of SCL’s regulatory functions (lineage priming, activation, and repression of gene expression programs) provides insight into fundamental developmental and transcriptional mechanisms, and highlights mechanistic parallels between normal and oncogenic processes. Importantly, recent discoveries are paving the way to the development of innovative therapeutic opportunities in SCL+ T-ALL.

Prediction and validation of protein-protein interactors from genome-wide DNA-binding data using a knowledge-based machine-learning approach

The ability to accurately predict the DNA targets and interacting cofactors of transcriptional regulators from genome-wide data can significantly advance our understanding of gene regulatory networks. NKX2-5 is a homeodomain transcription factor that sits high in the cardiac gene regulatory network and is essential for normal heart development. We previously identified genomic targets for NKX2-5 in mouse HL-1 atrial cardiomyocytes using DNA-adenine methyltransferase identification (DamID). Here, we apply machine learning algorithms and propose a knowledge-based feature selection method for predicting NKX2-5 protein : protein interactions based on motif grammar in genome-wide DNA-binding data. We assessed model performance using leave-one-out cross-validation and a completely independent DamID experiment performed with replicates. In addition to identifying previously described NKX2-5-interacting proteins, including GATA, HAND and TBX family members, a number of novel interactors were identified, with direct protein : protein interactions between NKX2-5 and retinoid X receptor (RXR), paired-related homeobox (PRRX) and Ikaros zinc fingers (IKZF) validated using the yeast two-hybrid assay. We also found that the interaction of RXRα with NKX2-5 mutations found in congenital heart disease (Q187H, R189G and R190H) was altered. These findings highlight an intuitive approach to accessing protein-protein interaction information of transcription factors in DNA-binding experiments.

Chromatin programming by developmentally regulated transcription factors: lessons from the study of haematopoietic stem cell specification and differentiation

Although the body plan of individuals is encoded in their genomes, each cell type expresses a different gene expression programme and therefore has access to only a subset of this information. Alterations to gene expression programmes are the underlying basis for the differentiation of multiple cell types and are driven by tissue-specific transcription factors (TFs) that interact with the epigenetic regulatory machinery to programme the chromatin landscape into transcriptionally active and inactive states. The haematopoietic system has long served as a paradigm for studying the molecular principles that regulate gene expression in development. In this review article, we summarize the current knowledge on the mechanism of action of TFs regulating haematopoietic stem cell specification and differentiation, and place this information into the context of general principles governing development.

Comprehensive Proteomic Analysis of Human Erythropoiesis

Mass spectrometry-based proteomics now enables the absolute quantification of thousands of proteins in individual cell types. We used this technology to analyze the dynamic proteome changes occurring during human erythropoiesis. We quantified the absolute expression of 6,130 proteins during erythroid differentiation from late burst-forming units-erythroid (BFU-Es) to orthochromatic erythroblasts. A modest correlation between mRNA and protein expression was observed. We identified several proteins with unexpected expression patterns in erythroid cells, highlighting a breakpoint in the erythroid differentiation process at the basophilic stage. We also quantified the distribution of proteins between reticulocytes and pyrenocytes after enucleation. These analyses identified proteins that are actively sorted either with the reticulocyte or the pyrenocyte. Our study provides the absolute quantification of protein expression during a complex cellular differentiation process in humans, and it establishes a framework for future studies of disordered erythropoiesis.

The proto-oncogenic protein TAL1 controls TGF-β1 signaling through interaction with SMAD3

TGF-β1 is involved in many aspects of tissue development and homeostasis including hematopoiesis. The TAL1 transcription factor is also an important player of this latter process and is expressed very early in the myeloid and erythroid lineages. We previously established a link between TGF-β1 signaling and TAL1 by showing that the cytokine was able to induce its proteolytic degradation by the ubiquitin proteasome pathway. In this manuscript we show that TAL1 interacts with SMAD3 that acts in the pathway downstream of TGF-β1 association with its receptor. TAL1 expression strengthens the positive or negative effect of SMAD3 on various genes. Both transcription factors activate the inhibitory SMAD7 factor through the E box motif present in its transcriptional promoter. DNA precipitation assays showed that TAL1 present in Jurkat or K562 cells binds to this SMAD binding element in a SMAD3 dependent manner. SMAD3 and TAL1 also inhibit several genes including ID1, hTERT and TGF-β1 itself. In this latter case TAL1 and SMAD3 can impair the positive effect exerted by E47. Our results indicate that TAL1 expression can modulate TGF-β1 signaling by interacting with SMAD3 and by increasing its transcriptional properties. They also suggest the existence of a negative feedback loop between TAL1 expression and TGF-β1 signaling.

Insight into GATA1 transcriptional activity through interrogation of cis elements disrupted in human erythroid disorders

Whole-exome sequencing has been incredibly successful in identifying causal genetic variants and has revealed a number of novel genes associated with blood and other diseases. One limitation of this approach is that it overlooks mutations in noncoding regulatory elements. Furthermore, the mechanisms by which mutations in transcriptionalcis-regulatory elements result in disease remain poorly understood. Here we used CRISPR/Cas9 genome editing to interrogate three such elements harboring mutations in human erythroid disorders, which in all cases are predicted to disrupt a canonical binding motif for the hematopoietic transcription factor GATA1. Deletions of as few as two to four nucleotides resulted in a substantial decrease (>80%) in target gene expression. Isolated deletions of the canonical GATA1 binding motif completely abrogated binding of the cofactor TAL1, which binds to a separate motif. Having verified the functionality of these three GATA1 motifs, we demonstrate strong evolutionary conservation of GATA1 motifs in regulatory elements proximal to other genes implicated in erythroid disorders, and show that targeted disruption of such elements results in altered gene expression. By modeling transcription factor binding patterns, we show that multiple transcription factors are associated with erythroid gene expression, and have created predictive maps modeling putative disruptions of their binding sites at key regulatory elements. Our study provides insight into GATA1 transcriptional activity and may prove a useful resource for investigating the pathogenicity of noncoding variants in human erythroid disorders.