The active spatial organization of the beta-globin locus requires the transcription factor EKLF

LDB1 establishes multi-enhancer networks to regulate gene expression

How specific enhancer-promoter pairing is established is still mostly unclear. Besides the CTCF/cohesin machinery, only a few nuclear factors have been studied for a direct role in physically connecting regulatory elements. Here, we show via acute degradation experiments that LDB1 directly and broadly promotes enhancer-promoter loops. Most LDB1-mediated contacts, even those spanning hundreds of kb, can form in the absence of CTCF, cohesin, or YY1 as determined via the use of multiple degron systems. Moreover, an engineered LDB1-driven chromatin loop is cohesin independent. Cohesin-driven loop extrusion does not stall at LDB1 occupied sites but may aid the formation of a subset of LDB1 anchored loops. Leveraging the dynamic reorganization of nuclear architecture during the transition from mitosis to G1-phase, we establish a relationship between LDB1-dependent interactions in the context of TAD organization and gene activation. Lastly, Tri-C and Region Capture Micro-C reveal that LDB1 organizes multi-enhancer networks to activate transcription. This establishes LDB1 as a direct driver of regulatory network inter-connectivity.

The Long Non-Coding RNA MALAT1 Modulates NR4A1 Expression through a Downstream Regulatory Element in Specific Cancer Cell Types

Long non-coding RNAs (lncRNAs) have been shown to modulate gene expression and are involved in the initiation and progression of various cancer types. Despite the wealth of studies describing transcriptome changes upon lncRNA knockdown, there is limited information describing lncRNA-mediated effects on regulatory elements (REs) modulating gene expression. In this study, we investigated how the metastasis-associated lung adenocarcinoma transcript 1 (MALAT1) lncRNA regulates primary target genes using time-resolved MALAT1 knockdown followed by parallel RNA-seq and ATAC-seq assays. The results revealed that MALAT1 primarily regulates specific protein-coding genes and a substantial decrease in the accessibility downstream of the NR4A1 gene that was associated with a decreased NR4A1 expression. Moreover, the presence of an NR4A1-downstream RE was demonstrated by CRISPR-i assays to define a functional MALAT1/NR4A1 axis. By analyzing TCGA data, we identified a positive correlation between NR4A1 expression and NR4A1-downstream RE accessibility in breast cancer but not in pancreatic cancer. Accordingly, this regulatory mechanism was experimentally validated in breast cancer cells (MCF7) but not in pancreatic duct epithelial carcinoma (PANC1) cells. Therefore, our results demonstrated that MALAT1 is involved in a molecular mechanism that fine-tunes NR4A1 expression by modulating the accessibility of a downstream RE in a cell type-specific manner.

Hi-C2B: Optimized Detection of Chromosomal Contacts Within Synchronized Meiotic S. cerevisiae Cells

Hi-C, a genome-wide chromosome conformation capture assay, is a powerful tool used to study three-dimensional genome organization by converting physical pairwise interactions into counts of pairwise interactions. To study the many temporally regulated facets of meiotic recombination in S. cerevisiae, the Hi-C assay must be robust such that fine- and wide-scale comparisons between genetic datasets can be made. Here we describe an updated protocol for Hi-C (Hi-C2B) that generates reproducible libraries of interaction data with low noise and for a relatively low cost.

Erythroid Krüppel-Like Factor (KLF1): A Surprisingly Versatile Regulator of Erythroid Differentiation

Erythroid Krüppel-like factor (KLF1), first discovered in 1992, is an erythroid-restricted transcription factor (TF) that is essential for terminal differentiation of erythroid progenitors. At face value, KLF1 is a rather inconspicuous member of the 26-strong SP/KLF TF family. However, 30 years of research have revealed that KLF1 is a jack of all trades in the molecular control of erythropoiesis. Initially described as a one-trick pony required for high-level transcription of the adult HBB gene, we now know that it orchestrates the entire erythroid differentiation program. It does so not only as an activator but also as a repressor. In addition, KLF1 was the first TF shown to be directly involved in enhancer/promoter loop formation. KLF1 variants underlie a wide range of erythroid phenotypes in the human population, varying from very mild conditions such as hereditary persistence of fetal hemoglobin and the In(Lu) blood type in the case of haploinsufficiency, to much more serious non-spherocytic hemolytic anemias in the case of compound heterozygosity, to dominant congenital dyserythropoietic anemia type IV invariably caused by a de novo variant in a highly conserved amino acid in the KLF1 DNA-binding domain. In this chapter, we present an overview of the past and present of KLF1 research and discuss the significance of human KLF1 variants.

The Mediator complex regulates enhancer-promoter interactions

Enhancer-mediated gene activation generally requires physical proximity between enhancers and their target gene promoters. However, the molecular mechanisms by which interactions between enhancers and promoters are formed are not well understood. Here, we investigate the function of the Mediator complex in the regulation of enhancer-promoter interactions, by combining rapid protein depletion and high-resolution MNase-based chromosome conformation capture approaches. We show that depletion of Mediator leads to reduced enhancer-promoter interaction frequencies, which are associated with a strong decrease in gene expression. In addition, we find increased interactions between CTCF-binding sites upon Mediator depletion. These changes in chromatin architecture are associated with a redistribution of the Cohesin complex on chromatin and a reduction in Cohesin occupancy at enhancers. Together, our results indicate that the Mediator and Cohesin complexes contribute to enhancer-promoter interactions and provide insights into the molecular mechanisms by which communication between enhancers and promoters is regulated.

A novel human cellular model of CDA IV enables comprehensive analysis revealing the molecular basis of the disease phenotype

Red blood cell disorders can result in severe anemia. One such disease congenital dyserythropoietic anemia IV (CDA IV) is caused by the heterozygous mutation E325K in the transcription factor KLF1. However, studying the molecular basis of CDA IV is severely impeded by the paucity of suitable and adequate quantities of material from patients with anemia and the rarity of the disease. We, therefore, took a novel approach, creating a human cellular disease model system for CDA IV that accurately recapitulates the disease phenotype. Next, using comparative proteomics, we reveal extensive distortion of the proteome and a wide range of disordered biological processes in CDA IV erythroid cells. These include downregulated pathways the governing cell cycle, chromatin separation, DNA repair, cytokinesis, membrane trafficking, and global transcription, and upregulated networks governing mitochondrial biogenesis. The diversity of such pathways elucidates the spectrum of phenotypic abnormalities that occur with CDA IV and impairment to erythroid cell development and survival, collectively explaining the CDA IV disease phenotype. The data also reveal far more extensive involvement of KLF1 in previously assigned biological processes, along with novel roles in the regulation of intracellular processes not previously attributed to this transcription factor. Overall, the data demonstrate the power of such a model cellular system to unravel the molecular basis of disease and how studying the effects of a rare mutation can reveal fundamental biology.

3D genome organization and epigenetic regulation in autoimmune diseases

Three-dimensional (3D) genomics is an emerging field of research that investigates the relationship between gene regulatory function and the spatial structure of chromatin. Chromatin folding can be studied using chromosome conformation capture (3C) technology and 3C-based derivative sequencing technologies, including chromosome conformation capture-on-chip (4C), chromosome conformation capture carbon copy (5C), and high-throughput chromosome conformation capture (Hi-C), which allow scientists to capture 3D conformations from a single site to the entire genome. A comprehensive analysis of the relationships between various regulatory components and gene function also requires the integration of multi-omics data such as genomics, transcriptomics, and epigenomics. 3D genome folding is involved in immune cell differentiation, activation, and dysfunction and participates in a wide range of diseases, including autoimmune diseases. We describe hierarchical 3D chromatin organization in this review and conclude with characteristics of C-techniques and multi-omics applications of the 3D genome. In addition, we describe the relationship between 3D genome structure and the differentiation and maturation of immune cells and address how changes in chromosome folding contribute to autoimmune diseases.

Identification of a genomic DNA sequence that quantitatively modulates KLF1 transcription factor expression in differentiating human hematopoietic cells

The onset of erythropoiesis is under strict developmental control, with direct and indirect inputs influencing its derivation from the hematopoietic stem cell. A major regulator of this transition is KLF1/EKLF, a zinc finger transcription factor that plays a global role in all aspects of erythropoiesis. Here, we have identified a short, conserved enhancer element in KLF1 intron 1 that is important for establishing optimal levels of KLF1 in mouse and human cells. Chromatin accessibility of this site exhibits cell-type specificity and is under developmental control during the differentiation of human CD34+ cells towards the erythroid lineage. This site binds GATA1, SMAD1, TAL1, and ETV6. In vivo editing of this region in cell lines and primary cells reduces KLF1 expression quantitatively. However, we find that, similar to observations seen in pedigrees of families with KLF1 mutations, downstream effects are variable, suggesting that the global architecture of the site is buffered towards keeping the KLF1 genetic region in an active state. We propose that modification of intron 1 in both alleles is not equivalent to complete loss of function of one allele.

Hematopoietic/erythroid enhancers activate nearby target genes by extending histone H3K27ac and transcribing intergenic RNA

Enhancers activate gene transcription remotely, which requires tissue specific transcription factors binding to them. GATA1 and TAL1 are hematopoietic/erythroid-specific factors and often bind together to enhancers, activating target genes. Interestingly, we found that some hematopoietic/erythroid genes are transcribed in a GATA1-dependent but TAL1-independnet manner. They appear to have enhancers within a relatively short distance. In this study, we paired highly transcribed hematopoietic/erythroid genes with the nearest GATA1/TAL1-binding enhancers and analyzed these putative enhancer-gene pairs depending on distance between them. Enhancers located at various distances from genes in the pairs, which was not related to transcription level of the genes. However, genes with enhancers at short distances away tended to be transcriptionally unaffected by TAL1 depletion. Histone H3K27ac extended from the enhancers to target genes. The H3K27ac extension was maintained without TAL1, even though it disappeared owing to the loss of GATA1. Intergenic RNA was highly transcribed from the enhancers to nearby target genes, independent of TAL1. Taken together, TAL1-independent transcription of hematopoietic/erythroid genes appears to be promoted by enhancers present in a short distance. These enhancers are likely to activate nearby target genes by tracking the intervening regions.

Forced enhancer-promoter rewiring to alter gene expression in animal models

Transcriptional enhancers can be in physical proximity of their target genes via chromatin looping. The enhancer at the β-globin locus (locus control region [LCR]) contacts the fetal-type (HBG) and adult-type (HBB) β-globin genes during corresponding developmental stages. We have demonstrated previously that forcing proximity between the LCR and HBG genes in cultured adult-stage erythroid cells can activate HBG transcription. Activation of HBG expression in erythroid cells is of benefit to patients with sickle cell disease. Here, using the β-globin locus as a model, we provide proof of concept at the organismal level that forced enhancer rewiring might present a strategy to alter gene expression for therapeutic purposes. Hematopoietic stem and progenitor cells (HSPCs) from mice bearing human β-globin genes were transduced with lentiviral vectors expressing a synthetic transcription factor (ZF-Ldb1) that fosters LCR-HBG contacts. When engrafted into host animals, HSPCs gave rise to adult-type erythroid cells with elevated HBG expression. Vectors containing ZF-Ldb1 were optimized for activity in cultured human and rhesus macaque erythroid cells. Upon transplantation into rhesus macaques, erythroid cells from HSPCs expressing ZF-Ldb1 displayed elevated HBG production. These findings in two animal models suggest that forced redirection of gene-regulatory elements may be used to alter gene expression to treat disease.

EKLF/Klf1 regulates erythroid transcription by its pioneering activity and selective control of RNA Pol II pause-release

EKLF/Klf1 is a zinc-finger transcription activator essential for erythroid lineage commitment and terminal differentiation. Using ChIP-seq, we investigate EKLF DNA binding and transcription activation mechanisms during mouse embryonic erythropoiesis. We utilize the Nan/+ mouse that expresses the EKLF-E339D (Nan) variant mutated in its conserved zinc-finger region and address the mechanism of hypomorphic and neomorphic changes in downstream gene expression. First, we show that Nan-EKLF limits normal EKLF binding to a subset of its sites. Second, we find that ectopic binding of Nan-EKLF occurs largely at enhancers and activates transcription through pioneering activity. Third, we find that for a subset of ectopic targets, gene activation is achieved in Nan/+ only by Nan-EKLF binding to distal enhancers, leading to RNA polymerase II pause-release. These results have general applicability to understanding how a DNA binding variant factor confers dominant disruptive effects on downstream gene expression even in the presence of its normal counterpart.

Globin vector regulatory elements are active in early hematopoietic progenitor cells

The globin genes are archetypal tissue-specific genes that are silent in most tissues but for late-stage erythroblasts upon terminal erythroid differentiation. The transcriptional activation of the β-globin gene is under the control of proximal and distal regulatory elements located on chromosome 11p15.4, including the β-globin locus control region (LCR). The incorporation of selected LCR elements in lentiviral vectors encoding β and β-like globin genes has enabled successful genetic treatment of the β-thalassemias and sickle cell disease. However, recent occurrences of benign clonal expansions in thalassemic patients and myelodysplastic syndrome in patients with sickle cell disease call attention to the non-erythroid functions of these powerful vectors. Here we demonstrate that lentivirally encoded LCR elements, in particular HS1 and HS2, can be activated in early hematopoietic cells including hematopoietic stem cells and myeloid progenitors. This activity is position-dependent and results in the transcriptional activation of a nearby reporter gene in these progenitor cell populations. We further show that flanking a globin vector with an insulator can effectively restrain this non-erythroid activity without impairing therapeutic globin expression. Globin lentiviral vectors harboring powerful LCR HS elements may thus expose to the risk of trans-activating cancer-related genes, which can be mitigated by a suitable insulator.

Inflammatory Immune-Associated eRNA: Mechanisms, Functions and Therapeutic Prospects

The rapid development of multiple high-throughput sequencing technologies has made it possible to explore the critical roles and mechanisms of functional enhancers and enhancer RNAs (eRNAs). The inflammatory immune response, as a fundamental pathological process in infectious diseases, cancers and immune disorders, coordinates the balance between the internal and external environment of the organism. It has been shown that both active enhancers and intranuclear eRNAs are preferentially expressed over inflammation-related genes in response to inflammatory stimuli, suggesting that enhancer transcription events and their products influence the expression and function of inflammatory genes. Therefore, in this review, we summarize and discuss the relevant inflammatory roles and regulatory mechanisms of eRNAs in inflammatory immune cells, non-inflammatory immune cells, inflammatory immune diseases and tumors, and explore the potential therapeutic effects of enhancer inhibitors affecting eRNA production for diseases with inflammatory immune responses.

Hemogen /BRG1 cooperativity modulates promoter and enhancer activation during erythropoiesis

Hemogen, also known as EDAG, is a hematopoietic tissue-specific gene that regulates the proliferation and differentiation of hematopoietic cells. However, the mechanism underlying hemogen function in erythropoiesis is unknown. We found that depletion of hemogen in human CD34+ erythroid progenitor cells and HUDEP2 cells significantly reduced the expression of genes associated with heme and hemoglobin synthesis, supporting a positive role of hemogen in erythroid maturation. In human K562 cells, hemogen antagonized the occupancy of co-repressors NuRD complex and facilitated LDB1 complex-mediated chromatin looping. Hemogen recruited SWI/SNF complex ATPase BRG1 as a co-activator to regulate nucleosome accessibility and H3K27ac enrichment for promoter and enhancer activity. To ask if hemogen/BRG1 cooperativity is conserved in mammalian systems, we generated hemogen KO/KI mice and investigated hemogen/BRG1 function in murine erythropoiesis. Loss of hemogen in E12.5-E16.5 fetal liver cells impeded erythroid differentiation through reducing the production of mature erythroblasts. ChIP-seq in WT and hemogen KO animal revealed BRG1 is largely dependent on hemogen to regulate chromatin accessibility at erythroid gene promoters and enhancers. In summary, hemogen/BRG1 interaction in mammals is essential for fetal erythroid maturation and hemoglobin production through its active role in promoter and enhancer activity and chromatin organization.

Epigenomic analysis of KLF1 haploinsufficiency in primary human erythroblasts

Haploinsufficiency for the erythroid-specific transcription factor KLF1 is associated with hereditary persistence of fetal hemoglobin (HPFH). Increased HbF ameliorates the symptoms of β-hemoglobinopathies and downregulation of KLF1 activity has been proposed as a potential therapeutic strategy. However, the feasibility of this approach has been challenged by the observation that KLF1 haploinsufficient individuals with the same KLF1 variant, within the same family, display a wide range of HbF levels. This phenotypic variability is not readily explained by co-inheritance of known HbF-modulating variants in the HBB, HBS1L-MYB and/or BCL11A loci. We studied cultured erythroid progenitors obtained from Maltese individuals in which KLF1 p.K288X carriers display HbF levels ranging between 1.3 and 12.3% of total Hb. Using a combination of gene expression analysis, chromatin accessibility assays and promoter activity tests we find that variation in expression of the wildtype KLF1 allele may explain a significant part of the variability in HbF levels observed in KLF1 haploinsufficiency. Our results have general bearing on the variable penetrance of haploinsufficiency phenotypes and on conflicting interpretations of pathogenicity of variants in other transcriptional regulators such as EP300, GATA2 and RUNX1.

Epigenetic-Mediated Regulation of Gene Expression for Biological Control and Cancer: Cell and Tissue Structure, Function, and Phenotype

Epigenetic gene regulatory mechanisms play a central role in the biological control of cell and tissue structure, function, and phenotype. Identification of epigenetic dysregulation in cancer provides mechanistic into tumor initiation and progression and may prove valuable for a variety of clinical applications. We present an overview of epigenetically driven mechanisms that are obligatory for physiological regulation and parameters of epigenetic control that are modified in tumor cells. The interrelationship between nuclear structure and function is not mutually exclusive but synergistic. We explore concepts influencing the maintenance of chromatin structures, including phase separation, recognition signals, factors that mediate enhancer-promoter looping, and insulation and how these are altered during the cell cycle and in cancer. Understanding how these processes are altered in cancer provides a potential for advancing capabilities for the diagnosis and identification of novel therapeutic targets.

CTCF-mediated chromatin looping provides a topological framework for the formation of phase-separated transcriptional condensates

CTCF is crucial to the organization of mammalian genomes into loop structures. According to recent studies, the transcription apparatus is compartmentalized and concentrated at super-enhancers to form phase-separated condensates and drive the expression of cell-identity genes. However, it remains unclear whether and how transcriptional condensates are coupled to higher-order chromatin organization. Here, we show that CTCF is essential for RNA polymerase II (Pol II)-mediated chromatin interactions, which occur as hyperconnected spatial clusters at super-enhancers. We also demonstrate that CTCF clustering, unlike Pol II clustering, is independent of liquid-liquid phase-separation and resistant to perturbation of transcription. Interestingly, clusters of Pol II, BRD4, and MED1 were found to dissolve upon CTCF depletion, but were reinstated upon restoration of CTCF, suggesting a potent instructive function for CTCF in the formation of transcriptional condensates. Overall, we provide evidence suggesting that CTCF-mediated chromatin looping acts as an architectural prerequisite for the assembly of phase-separated transcriptional condensates.

The 3D Genome: From Structure to Function

The genome is the most functional part of a cell, and genomic contents are organized in a compact three-dimensional (3D) structure. The genome contains millions of nucleotide bases organized in its proper frame. Rapid development in genome sequencing and advanced microscopy techniques have enabled us to understand the 3D spatial organization of the genome. Chromosome capture methods using a ligation approach and the visualization tool of a 3D genome browser have facilitated detailed exploration of the genome. Topologically associated domains (TADs), lamin-associated domains, CCCTC-binding factor domains, cohesin, and chromatin structures are the prominent identified components that encode the 3D structure of the genome. Although TADs are the major contributors to 3D genome organization, they are absent in Arabidopsis. However, a few research groups have reported the presence of TAD-like structures in the plant kingdom.

Transcriptional enhancers and their communication with gene promoters

Transcriptional enhancers play a key role in the initiation and maintenance of gene expression programmes, particularly in metazoa. How these elements control their target genes in the right place and time is one of the most pertinent questions in functional genomics, with wide implications for most areas of biology. Here, we synthesise classic and recent evidence on the regulatory logic of enhancers, including the principles of enhancer organisation, factors that facilitate and delimit enhancer-promoter communication, and the joint effects of multiple enhancers. We show how modern approaches building on classic insights have begun to unravel the complexity of enhancer-promoter relationships, paving the way towards a quantitative understanding of gene control.

Seeing the forest through the trees: prioritising potentially functional interactions from Hi-C

Eukaryotic genomes are highly organised within the nucleus of a cell, allowing widely dispersed regulatory elements such as enhancers to interact with gene promoters through physical contacts in three-dimensional space. Recent chromosome conformation capture methodologies such as Hi-C have enabled the analysis of interacting regions of the genome providing a valuable insight into the three-dimensional organisation of the chromatin in the nucleus, including chromosome compartmentalisation and gene expression. Complicating the analysis of Hi-C data, however, is the massive amount of identified interactions, many of which do not directly drive gene function, thus hindering the identification of potentially biologically functional 3D interactions. In this review, we collate and examine the downstream analysis of Hi-C data with particular focus on methods that prioritise potentially functional interactions. We classify three groups of approaches: structural-based discovery methods, e.g. A/B compartments and topologically associated domains, detection of statistically significant chromatin interactions, and the use of epigenomic data integration to narrow down useful interaction information. Careful use of these three approaches is crucial to successfully identifying potentially functional interactions within the genome.

Chromatin Conformation in Development and Disease

Chromatin domains and loops are important elements of chromatin structure and dynamics, but much remains to be learned about their exact biological role and nature. Topological associated domains and functional loops are key to gene expression and hold the answer to many questions regarding developmental decisions and diseases. Here, we discuss new findings, which have linked chromatin conformation with development, differentiation and diseases and hypothesized on various models while integrating all recent findings on how chromatin architecture affects gene expression during development, evolution and disease.

Enhancers navigate the three-dimensional genome to direct cell fate decisions

The activity and selectivity of transcriptional enhancers determine gene expression patterns that enable a zygote to become a complex organism. How enhancers convey regulatory information is a central conundrum in biology. Here, we discuss recent progress provided by rapidly evolving technologies in understanding enhancer-promoter interactions in the context of overall nuclear genome organization.

Models behind the mystery of establishing enhancer-promoter interactions

Enhancers and promoters are transcriptional regulatory elements whose facilitated interactions increase gene expression. Enhancer DNA sequences can be located far away from the promoter sequences that they regulate. Currently, the mechanism facilitating the establishment of enhancer-promoter interactions remains unclear. However, mutations causing errors in these interactions have been linked to cancer and disease, further conveying the need to understand the full mechanism. This review discusses multiple models that have been proposed to describe how enhancers go the distance to interact with promoters. Evidence supporting loop formation models is reviewed in addition to more complex hypotheses involving aspects of 3D chromatin organization and phase separation.

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.

The relationship between genome structure and function

Precise patterns of gene expression in metazoans are controlled by three classes of regulatory elements: promoters, enhancers and boundary elements. During differentiation and development, these elements form specific interactions in dynamic higher-order chromatin structures. However, the relationship between genome structure and its function in gene regulation is not completely understood. Here we review recent progress in this field and discuss whether genome structure plays an instructive role in regulating gene expression or is a reflection of the activity of the regulatory elements of the genome.

A computational framework for identifying the transcription factors involved in enhancer-promoter loop formation

The pairwise interaction between transcription factors (TFs) plays an important role in enhancer-promoter loop formation. Although thousands of TFs in the human genome have been found, only a few TF pairs have been demonstrated to be related to loop formation. It is still a challenge to determine which TF pairs could be involved in the enhancer-promoter regulation network. This work describes a computational framework to identify TF pairs in enhancer-promoter regulation. By integrating different levels of data derived from Promoter Capture Hi-C, chromatin immunoprecipitation sequencing (ChIP-seq) of histone marks, RNA-seq, protein-protein interaction (PPI), and TF motif, we identified 361 significant TF pairs and constructed a TF interaction network. From the network, we found several hub-TFs, which may have important roles in the regulation of long-range interactions. Our studies extended TF pairs identified in other experimental and computational approaches. These findings will help the further study of long-range interactions between enhancers and promoters.

Spatial Genome Re-organization between Fetal and Adult Hematopoietic Stem Cells

Fetal hematopoietic stem cells (HSCs) undergo a developmental switch to become adult HSCs with distinct functional properties. To better understand the molecular mechanisms underlying the developmental switch, we have conducted deep sequencing of the 3D genome, epigenome, and transcriptome of fetal and adult HSCs in mouse. We find that chromosomal compartments and topologically associating domains (TADs) are largely conserved between fetal and adult HSCs. However, there is a global trend of increased compartmentalization and TAD boundary strength in adult HSCs. In contrast, intra-TAD chromatin interactions are much more dynamic and wide-spread, involving over a thousand gene promoters and distal enhancers. These developmental-stage-specific enhancer-promoter interactions are mediated by different sets of transcription factors, such as TCF3 and MAFB in fetal HSCs, versus NR4A1 and GATA3 in adult HSCs. Loss-of-function studies of TCF3 confirm the role of TCF3 in mediating condition-specific enhancer-promoter interactions and gene regulation in fetal HSCs.

A developmental transition occurs between fetal and adult hematopoietic stem cells. How the 3D genome folding contributes to this transition is poorly understood. Chen et al. show global genome organization is largely conserved, but a large fraction of enhancer-promoter interactions is reorganized and regulate genes contributing to the phenotypic differences.

Transcription factors: building hubs in the 3D space

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

A Krüppel-like factor 1 (KLF1) Mutation Associated with Severe Congenital Dyserythropoietic Anemia Alters Its DNA-Binding Specificity

Krüppel-like factor 1 (KLF1/EKLF) is a transcription factor that globally activates genes involved in erythroid cell development. Various mutations are identified in the human KLF1 gene. The E325K mutation causes congenital dyserythropoietic anemia (CDA) type IV, characterized by severe anemia and non-erythroid-cell-related symptoms. The CDA mutation is in the second zinc finger of KLF1 at a position functionally involved in its interactions with DNA. The molecular parameters of how CDA-KLF1 exerts its biological effects have not been addressed. Here, using an in vitro selection strategy, we determined the preferred DNA-binding site for CDA-KLF1. Binding to the deduced consensus sequence is supported by in vitro gel shifts and by in vivo functional reporter gene studies. Two significant changes compared to wild-type (WT) binding are observed: G is selected as the middle nucleotide, and the 3' portion of the consensus sequence is more degenerate. As a consequence, CDA-KLF1 did not bind the WT consensus sequence. However, activation of ectopic sites is promoted. Continuous activation of WT target genes occurs if they fortuitously contain the novel CDA site nearby. Our findings provide a molecular understanding of how a single mutation in the KLF1 zinc finger exerts effects on erythroid physiology in CDA type IV.

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

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

Implications of CpG islands on chromosomal architectures and modes of global gene regulation

CpG islands (CGIs) have long been implicated in the regulation of vertebrate gene expression. However, the involvement of CGIs in chromosomal architectures and associated gene expression regulations has not yet been thoroughly explored. By combining large-scale integrative data analyses and experimental validations, we show that CGIs clearly reconcile two competing models explaining nuclear gene localizations. We first identify CGI-containing (CGI+) and CGI-less (CGI-) genes are non-randomly clustered within the genome, which reflects CGI-dependent spatial gene segregation in the nucleus and corresponding gene regulatory modes. Regardless of their transcriptional activities, CGI+ genes are mainly located at the nuclear center and encounter frequent long-range chromosomal interactions. Meanwhile, nuclear peripheral CGI- genes forming heterochromatin are activated and internalized into the nuclear center by local enhancer-promoter interactions. Our findings demonstrate the crucial implications of CGIs on chromosomal architectures and gene positioning, linking the critical importance of CGIs in determining distinct mechanisms of global gene regulation in three-dimensional space in the nucleus.

Runx1 and RORγt Cooperate to Upregulate IL-22 Expression in Th Cells through Its Distal Enhancer

IL-22 is a cytokine that plays a pivotal role in regulating tissue homeostasis at barrier surfaces and is produced by activated CD4+ Th cells. Currently, the molecular mechanisms regulating Il22 gene expression are still unclear. In this study, we have identified a crucial cis-regulatory element located 32 kb upstream of the mouse Il22 promoter, termed conserved noncoding sequence (CNS)-32. We demonstrated that CNS-32 acts as an enhancer in reporter assays and contains binding motifs for Runt-related transcription factor (Runx)1 and retinoic acid-related orphan receptor γt (RORγt). Mutation of these motifs significantly abrogated the reporter activity, suggesting a role for both factors in the control of enhancer-mediated Il22 expression. Runx1 and RORγt occupancy and elevated histone H4 acetylation at CNS-32 were evident, as naive T cells differentiated into IL-22-producing Th22 cells. Overexpression of Runx1 promoted IL-22 production by inducing RORγt and IL-23 receptor, all critical to Th22 cell induction. Although Runx1 alone enhanced IL-22 production in Th22 cells, it was further enhanced in the presence of RORγt. Conversely, short hairpin RNA-mediated knockdown of core-binding factor β, a cofactor essential for Runx1 activity, was effective in limiting IL-22 production. Collectively, our results suggest that IL-22 production is controlled by a regulatory circuit in which Runx1 induces RORγt and then partners with RORγt to direct Il22 expression through their targeting of the Il22 enhancer.

The mouse KLF1 Nan variant impairs nuclear condensation and erythroid maturation

Krüppel-like factor 1 (KLF1) is an essential transcription factor for erythroid development, as demonstrated by Klf1 knockout mice which die around E14 due to severe anemia. In humans, >140 KLF1 variants, causing different erythroid phenotypes, have been described. The KLF1 Nan variant, a single amino acid substitution (p.E339D) in the DNA binding domain, causes hemolytic anemia and is dominant over wildtype KLF1. Here we describe the effects of the KLF1 Nan variant during fetal development. We show that Nan embryos have defects in erythroid maturation. RNA-sequencing of the KLF1 Nan fetal liver cells revealed that Exportin 7 (Xpo7) was among the 782 deregulated genes. This nuclear exportin is implicated in terminal erythroid differentiation; in particular it is involved in nuclear condensation. Indeed, KLF1 Nan fetal liver cells had larger nuclei and reduced chromatin condensation. Knockdown of XPO7 in wildtype erythroid cells caused a similar phenotype. We propose that reduced expression of XPO7 is partially responsible for the erythroid defects observed in KLF1 Nan erythroid cells.

Chromosome territories and the global regulation of the genome

Spatial positioning is a fundamental principle governing nuclear processes. Chromatin is organized as a hierarchy from nucleosomes to Mbp chromatin domains (CD) or topologically associating domains (TADs) to higher level compartments culminating in chromosome territories (CT). Microscopic and sequencing techniques have substantiated chromatin organization as a critical factor regulating gene expression. For example, enhancers loop back to interact with their target genes almost exclusively within TADs, distally located coregulated genes reposition into common transcription factories upon activation, and Mbp CDs exhibit dynamic motion and configurational changes in vivo. A longstanding question in the nucleus field is whether an interactive nuclear matrix provides a direct link between structure and function. The findings of nonrandom radial positioning of CT within the nucleus suggest the possibility of preferential interaction patterns among populations of CT. Sequential labeling up to 10 CT followed by application of computer imaging and geometric graph mining algorithms revealed cell-type specific interchromosomal networks (ICN) of CT that are altered during the cell cycle, differentiation, and cancer progression. It is proposed that the ICN correlate with the global level of genome regulation. These approaches also demonstrated that the large scale 3-D topology of CT is specific for each CT. The cell-type specific proximity of certain chromosomal regions in normal cells may explain the propensity of distinct translocations in cancer subtypes. Understanding how genes are dysregulated upon disruption of the normal "wiring" of the nucleus by translocations, deletions, and amplifications that are hallmarks of cancer, should enable more targeted therapeutic strategies.

Genetic disarray follows mutant KLF1-E325K expression in a congenital dyserythropoietic anemia patient

Congenital dyserythropoietic anemia type IV is caused by a heterozygous mutation, Glu325Lys (E325K), in the KLF1 transcription factor. Molecular characteristics of this disease have not been clarified, partly due to its rarity. We expanded erythroid cells from a patient's peripheral blood and analyzed its global expression pattern. We find that a large number of erythroid pathways are disrupted, particularly those related to membrane transport, globin regulation, and iron utilization. The altered genetics lead to significant deficits in differentiation. Glu325 is within the KLF1 zinc finger domain at an amino acid critical for site specific DNA binding. The change to Lys is predicted to significantly alter the target site recognition sequence, both by subverting normal recognition and by enabling interaction with novel sites. Consistent with this, we find high level ectopic expression of genes not normally present in the red cell. These altered properties explain patients' clinical and phenotypic features, and elucidate the dominant character of the mutation.

Transvection-like interchromosomal interaction is not observed at the transcriptional level when tested in the Rosa26 locus in mouse

Long-range associations between enhancers and their target gene promoters have been shown to play critical roles in executing genome function. Recent variations of chromosome capture technology have revealed a comprehensive view of intra- and interchromosomal contacts between specific genomic sites. The locus control region of the β-globin genes (β-LCR) is a super-enhancer that is capable of activating all of the β-like globin genes within the locus in cis through physical interaction by forming DNA loops. CTCF helps to mediate loop formation between LCR-HS5 and 3’HS1 in the human β-globin locus, in this way thought to contribute to the formation of a “chromatin hub”. The β-globin locus is also in close physical proximity to other erythrocyte-specific genes located long distances away on the same chromosome. In this case, erythrocyte-specific genes gather together at a shared “transcription factory” for co-transcription. Theoretically, enhancers could also activate target gene promoters at the identical loci, yet on different chromosomes in trans, a phenomenon originally described as transvection in Drosophilla. Although close physical proximity has been reported for the β-LCR and the β-like globin genes when integrated at the mouse homologous loci in trans, their structural and functional interactions were found to be rare, possibly because of a lack of suitable regulatory elements that might facilitate such trans interactions. Therefore, we re-evaluated presumptive transvection-like enhancer-promoter communication by introducing CTCF binding sites and erythrocyte-specific transcription units into both LCR-enhancer and β-promoter alleles, each inserted into the mouse ROSA26 locus on separate chromosomes. Following cross-mating of mice to place the two mutant loci at the identical chromosomal position and into active chromation in trans, their transcriptional output was evaluated. The results demonstrate that there was no significant functional association between the LCR and the β-globin gene in trans even in this idealized experimental context.

The interdependence of gene-regulatory elements and the 3D genome

Imaging studies, high-resolution chromatin conformation maps, and genome-wide occupancy data of architectural proteins have revealed that genome topology is tightly intertwined with gene expression. Cross-talk between gene-regulatory elements is often organized within insulated neighborhoods, and regulatory cues that induce transcriptional changes can reshape chromatin folding patterns and gene positioning within the nucleus. The cause-consequence relationship of genome architecture and gene expression is intricate, and its molecular mechanisms are under intense investigation. Here, we review the interdependency of transcription and genome organization with emphasis on enhancer-promoter contacts in gene regulation.

Temporal autoregulation during human PU.1 locus SubTAD formation

Epigenetic control of gene expression occurs within discrete spatial chromosomal units called topologically associating domains (TADs), but the exact spatial requirements of most genes are unknown; this is of particular interest for genes involved in cancer. We therefore applied high-resolution chromosomal conformation capture sequencing to map the three-dimensional (3D) organization of the human locus encoding the key myeloid transcription factor PU.1 in healthy monocytes and acute myeloid leukemia (AML) cells. We identified a dynamic ∼75-kb unit (SubTAD) as the genomic region in which spatial interactions between PU.1 gene regulatory elements occur during myeloid differentiation and are interrupted in AML. Within this SubTAD, proper initiation of the spatial chromosomal interactions requires PU.1 autoregulation and recruitment of the chromatin-adaptor protein LDB1 (LIM domain-binding protein 1). However, once these spatial interactions have occurred, LDB1 stabilizes them independently of PU.1 autoregulation. Thus, our data support that PU.1 autoregulates its expression in a "hit-and-run" manner by initiating stable chromosomal loops that result in a transcriptionally active chromatin architecture.

Adriamycin inhibits embryonic development in zebrafish through downregulation of Kruppel-like factor4

Adriamycin is an effective anticancer drug used in a wide range of cancers. Anticancer drugs modulate oncogenes and nodal regulatory molecules that affect cell differentiation and organismal development. In this study, we explore the effect of adriamycin on Kruppel-like factor4 (Klf4), an essential pluripotent factor by choosing zebrafish embryos as a model system. Klf4 is involved in the regulation of cellular growth, proliferation, and differentiation. In zebrafish embryogenesis, Klf4 is a major regulator of differentiation of polster in the anterior mesendoderm region of cells into hatching gland cells. The importance of this study is to check the effect of adriamycin on embryonic development. We found, adriamycin dose dependently altered the gene expression level of Klf4 that occurs in parallel to its detrimental effect on hatching. Supportively, cathepsin L and cyclase-associated protein1 are the other two markers of hatching that are altered along with Klf4.

Robust erythroid differentiation system for rhesus hematopoietic progenitor cells allowing preclinical screening of genetic treatment strategies for the hemoglobinopathies

BACKGROUND AIMS

γ-globin expression can be induced by various gene modification strategies, which could be beneficial for hemoglobin (Hb) disorders. To translate promising ideas into clinics, large animal models have proven valuable to evaluate safety and efficacy of the approaches; however, in vitro erythroid differentiation methods have not been established to determine whether they can be modeled in nonhuman primates.

METHODS

We optimized erythroid differentiation culture to produce high-level adult Hb from rhesus hematopoietic progenitor cells by using low (LC) or high cytokine concentration (HC) protocols with or without feeder cells. In addition, we established rhesus globin protein analysis using reverse-phase high performance liquid chromatography and mass spectrometry.

RESULTS

Robust adult Hb production at protein levels was observed in the LC protocol when feeder cells were used, whereas the HC protocol resulted in higher baseline fetal Hb levels (P < 0.01). We then compared lentiviral transduction of rhesus cells between serum-containing LC media and serum-free StemSpan-based differentiation media, revealing 100-fold more efficient transduction in serum-free differentiation media (P < 0.01). Finally, rhesus CD34⁺ cells were transduced with lentiviral vectors encoding artificial zinc finger proteins (ZF-Ldb1), which can reactivate γ-globin expression via tethering the transcriptional co-regulator Ldb1 to γ-globin promoters, and were differentiated in the optimized erythroid differentiation method. This resulted in marked increases of γ-globin levels compared with control groups (P < 0.01).

DISCUSSION

In conclusion, we developed an efficient rhesus erythroid differentiation protocol from hematopoietic progenitor cells with low fetal and high adult Hb production. Further studies are warranted to optimize gene modification and transplantation of rhesus hematopoietic progenitor cells.

KFL1 Gene Variants in α-Thalassemia Individuals with Increased Fetal Hemoglobin in a Chinese Population

Krüppel-like factor 1 (KLF1) is a pleiotropic erythroid transcription factor that is a regulator of definitive erythropoiesis. The aim of this study was to detect KLF1 gene variants in α-thalassemia (α-thal) carriers with an increased Hb F level in a Chinese population, and determine the changes of hematological parameters as a result of interactions between KLF1 gene mutations and α-thal. Subjects with α-thal and Hb F levels of ≥1.0% were selected for further investigation. Direct sequencing was used to detect KLF1 gene mutations. Hematological parameters of subjects with α-thal and concomitant KLF1 gene mutations and those with α-thal alone were compared. The KLF1 gene variants were detected in 46 of 275 (16.7%) individuals with α-thal and Hb F levels of ≥1.0%. The detection rate of KLF1 gene mutations rose correspondingly when the Hb F level increased. For α⁰-thal carriers, significantly lower mean corpuscular volume (MCV) and mean corpuscular hemoglobin (Hb) (MCH) values were observed in KLF1 gene mutation-positive carriers than that in KLF1 gene mutation-free carriers; conversely, significantly higher Hb A₂ and Hb F levels were observed in the former condition rather than in the latter condition. The results of this study indicate that KLF1 gene variants are common in Chinese subjects with α-thal and increased Hb F levels, and KLF1 gene mutations decreased the red blood cell (RBC) indices in α-thal carriers as that in normal adults.

Enhancer and super-enhancer: Positive regulators in gene transcription

Enhancer is a positive regulator for spatiotemporal development in eukaryotes. As a cluster, super-enhancer is closely related to cell identity- and fate-determined processes. Both of them function tightly depending on their targeted transcription factors, cofactors, and genes through distal genomic interactions. They have been recognized as critical components and played positive roles in transcriptional regulatory network or factory. Recent advances of next-generation sequencing have dramatically expanded our ability and knowledge to interrogate the molecular mechanism of enhancer and super-enhancer for transcription. Here, we review the history, importance, advances and challenges on enhancer and super-enhancer field. This will benefit our understanding of their function mechanism for transcription underlying precise gene expression.

Trac-looping measures genome structure and chromatin accessibility

Long-range chromatin interactions play critical roles in genome organization and regulation of transcription. We now report transposase-mediated analysis of chromatin looping (Trac-looping) for simultaneous detection of multiscale genome-wide chromatin interactions among regulatory elements and chromatin accessibility. With this technique, a bivalent oligonucleotide linker is inserted between two interacting regions such that the chromatin interactions are captured without prior chromatin fragmentation and proximity-based ligation. Application of Trac-looping to human CD4⁺ T cells revealed substantial reorganization of enhancer-promoter interactions associated with changes in gene expression after T cell receptor stimulation.

Genetic and Epigenetic Mechanisms of β-Globin Gene Switching

Vertebrates have multiple forms of hemoglobin that differ in the composition of their polypeptide chains. During ontogenesis, the composition of these subunits changes. Genes encoding different α- and β-polypeptide chains are located in two multigene clusters on different chromosomes. Each cluster contains several genes that are expressed at different stages of ontogenesis. The phenomenon of stage-specific transcription of globin genes is referred to as globin gene switching. Mechanisms of expression switching, stage-specific activation, and repression of transcription of α- and β-globin genes are of interest from both theoretical and practical points of view. Alteration of balanced expression of globin genes, which usually occurs due to damage to adult β-globin genes, leads to development of severe diseases - hemoglobinopathies. In most cases, reactivation of the fetal hemoglobin gene in patients with β-thalassemia and sickle cell disease can reduce negative consequences of irreversible alterations of expression of the β-globin genes. This review focuses on the current state of research on genetic and epigenetic mechanisms underlying stage-specific switching of β-globin genes.

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.

miR-326 regulates HbF synthesis by targeting EKLF in human erythroid cells

Haploinsufficiency of erythroid Krüppel-like factor (EKLF/KLF1) has been shown recently to ameliorate the clinical severity of β-thalassemia by increased expression levels of fetal hemoglobin (HbF). The underlying mechanisms for role of EKLF in regulating HbF are of great interest but remain incompletely understood. In this study, we used a combination of in silico, in vitro, and in vivo approaches to identify microRNAs (miRs) involved in EKLF regulation and to validate the role of miR-326 in HbF modification. We found that miR-326 suppresses EKLF expression directly by targeting its 3' untranslated region. miR-326 overexpression in K562 cells or CD34⁺ hematopoietic progenitor cells resulted in reduced EKLF protein levels and was associated with elevated expression of γ-globin, whereas inhibition of physiological miR-326 levels increased EKLF and thus reduced γ-globin expression. Moreover, miR-326 expression is positively correlated with HbF levels in β-thalassemia patients. Our results suggest that miR-326 plays a key role in regulating EKLF expression and in modifying the HbF level, which may provide a new strategy for activating HbF in individuals with β-thalassemia or sickle cell disease.

Microfluidics-Based Chromosome Conformation Capture (3C) Technology for Examining Chromatin Organization with a Low Quantity of Cells

Detecting three-dimensional (3D) genome organization in the form of physical interactions between various genomic loci is of great importance for understanding transcriptional regulations and cellular fate. Chromosome Conformation Capture (3C) method is the gold standard for examining chromatin organization, but usually requires a large number of cells (>107). This hinders studies of scarce tissue samples from animals and patients using the method. Here we developed a microfluidics-based approach for examining chromosome conformation by 3C technology. Critical 3C steps, such as digestion and religation of BAC DNA and cross-linked chromatin, were implemented on a microfluidic chip using a low quantity of cells (<104). Using this technology, we analyzed the chromatin looping interactions in the human β-globin. We envision that our method will provide a powerful tool for low-input analysis of chromosome conformation and epigenetic regulations.