Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes

Knowledge of locations and activities of cis-regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our Validated Systematic Integration (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state Regulatory Potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbored distinctive transcription factor binding motifs that were similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we showed that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

Haemoglobin Bart's Hydrops Fetalis: Charting the Past and Envisioning the Future

Haemoglobin Bart's hydrops fetalis syndrome (BHFS) represents the most severe form of α-thalassaemia, arising from deletion of the duplicated α-globin genes from both alleles. The absence of α-globin leads to the formation of non-functional haemoglobin Bart's (γ4) or haemoglobin H (HbH: β4) resulting in severe anaemia, tissue hypoxia, and, in some cases, variable congenital or neurocognitive abnormalities. BHFS is the most common cause of hydrops fetalis in Southeast Asia; however, owing to global migration, the burden of this condition is increasing worldwide. With the availability of intensive perinatal care and intrauterine transfusions, an increasing number of patients survive with this condition. The current approach to long-term management of survivors involves regular blood transfusions and iron chelation, a task made challenging by the need for intensified transfusions to suppress the production of non-functional HbH-containing erythrocytes. While our knowledge of outcomes of this condition is evolving, it seems, in comparison to individuals with transfusion-dependent β-thalassaemia, those with BHFS may face an elevated risk of complications arising from chronic anaemia and hypoxia, ongoing haemolysis, iron overload, and from their respective treatments. Although stem cell transplantation remains a viable option for a select few, it is not without potential side effects. Looking ahead, potential advancements in the form of genetic engineering and innovative therapeutic approaches, such as the reactivation of embryonic α-like globin gene expression, hold promise for furthering the treatment of this condition. Prevention remains a crucial aspect of care, particularly in areas with high prevalence or limited resources.

Super-enhancers include classical enhancers and facilitators to fully activate gene expression

Super-enhancers are compound regulatory elements that control expression of key cell identity genes. They recruit high levels of tissue-specific transcription factors and co-activators such as the Mediator complex and contact target gene promoters with high frequency. Most super-enhancers contain multiple constituent regulatory elements, but it is unclear whether these elements have distinct roles in activating target gene expression. Here, by rebuilding the endogenous multipartite α-globin super-enhancer, we show that it contains bioinformatically equivalent but functionally distinct element types: classical enhancers and facilitator elements. Facilitators have no intrinsic enhancer activity, yet in their absence, classical enhancers are unable to fully upregulate their target genes. Without facilitators, classical enhancers exhibit reduced Mediator recruitment, enhancer RNA transcription, and enhancer-promoter interactions. Facilitators are interchangeable but display functional hierarchy based on their position within a multipartite enhancer. Facilitators thus play an important role in potentiating the activity of classical enhancers and ensuring robust activation of target genes.

RNA polymerase II pausing temporally coordinates cell cycle progression and erythroid differentiation

Controlled release of promoter-proximal paused RNA polymerase II (RNA Pol II) is crucial for gene regulation. However, studying RNA Pol II pausing is challenging, as pause-release factors are almost all essential. In this study, we identified heterozygous loss-of-function mutations in SUPT5H, which encodes SPT5, in individuals with β-thalassemia. During erythropoiesis in healthy human cells, cell cycle genes were highly paused as cells transition from progenitors to precursors. When the pathogenic mutations were recapitulated by SUPT5H editing, RNA Pol II pause release was globally disrupted, and as cells began transitioning from progenitors to precursors, differentiation was delayed, accompanied by a transient lag in erythroid-specific gene expression and cell cycle kinetics. Despite this delay, cells terminally differentiate, and cell cycle phase distributions normalize. Therefore, hindering pause release perturbs proliferation and differentiation dynamics at a key transition during erythropoiesis, identifying a role for RNA Pol II pausing in temporally coordinating the cell cycle and erythroid differentiation.

Ancient genomic linkage couples metabolism with erythroid development

Generation of mature cells from progenitors requires tight coupling of differentiation and metabolism. During erythropoiesis, erythroblasts are required to massively upregulate globin synthesis then clear extraneous material and enucleate to produce erythrocytes1-3. Nprl3 has remained in synteny with the α-globin genes for >500 million years4, and harbours the majority of the α-globin enhancers5. Nprl3 is a highly conserved inhibitor of mTORC1, which controls cellular metabolism. However, whether Nprl3 itself serves an erythroid role is unknown. Here, we show that Nprl3 is a key regulator of erythroid metabolism. Using Nprl3-deficient fetal liver and adult competitive bone marrow - fetal liver chimeras, we show that NprI3 is required for sufficient erythropoiesis. Loss of Nprl3 elevates mTORC1 signalling, suppresses autophagy and disrupts erythroblast glycolysis and redox control. Human CD34+ progenitors lacking NPRL3 produce fewer enucleated cells and demonstrate dysregulated mTORC1 signalling in response to nutrient availability and erythropoietin. Finally, we show that the α-globin enhancers upregulate NprI3 expression, and that this activity is necessary for optimal erythropoiesis. Therefore, the anciently conserved linkage of NprI3, α-globin and their associated enhancers has enabled coupling of metabolic and developmental control in erythroid cells. This may enable erythropoiesis to adapt to fluctuating nutritional and environmental conditions.

Direct correction of haemoglobin E β-thalassaemia using base editors

Haemoglobin E (HbE) β-thalassaemia causes approximately 50% of all severe thalassaemia worldwide; equating to around 30,000 births per year. HbE β-thalassaemia is due to a point mutation in codon 26 of the human HBB gene on one allele (GAG; glutamatic acid → AAG; lysine, E26K), and any mutation causing severe β-thalassaemia on the other. When inherited together in compound heterozygosity these mutations can cause a severe thalassaemic phenotype. However, if only one allele is mutated individuals are carriers for the respective mutation and have an asymptomatic phenotype (β-thalassaemia trait). Here we describe a base editing strategy which corrects the HbE mutation either to wildtype (WT) or a normal variant haemoglobin (E26G) known as Hb Aubenas and thereby recreates the asymptomatic trait phenotype. We have achieved editing efficiencies in excess of 90% in primary human CD34 + cells. We demonstrate editing of long-term repopulating haematopoietic stem cells (LT-HSCs) using serial xenotransplantation in NSG mice. We have profiled the off-target effects using a combination of circularization for in vitro reporting of cleavage effects by sequencing (CIRCLE-seq) and deep targeted capture and have developed machine-learning based methods to predict functional effects of candidate off-target mutations.

RNA Polymerase II pausing temporally coordinates cell cycle progression and erythroid differentiation

The controlled release of promoter-proximal paused RNA polymerase II (Pol II) into productive elongation is a major step in gene regulation. However, functional analysis of Pol II pausing is difficult because factors that regulate pause release are almost all essential. In this study, we identified heterozygous loss-of-function mutations in SUPT5H , which encodes SPT5, in individuals with β-thalassemia unlinked to HBB mutations. During erythropoiesis in healthy human cells, cell cycle genes were highly paused at the transition from progenitors to precursors. When the pathogenic mutations were recapitulated by SUPT5H editing, Pol II pause release was globally disrupted, and the transition from progenitors to precursors was delayed, marked by a transient lag in erythroid-specific gene expression and cell cycle kinetics. Despite this delay, cells terminally differentiate, and cell cycle phase distributions normalize. Therefore, hindering pause release perturbs proliferation and differentiation dynamics at a key transition during erythropoiesis, revealing a role for Pol II pausing in the temporal coordination between the cell cycle and differentiation.

Scalable in vitro production of defined mouse erythroblasts

Mouse embryonic stem cells (mESCs) can be manipulated in vitro to recapitulate the process of erythropoiesis, during which multipotent cells undergo lineage specification, differentiation and maturation to produce erythroid cells. Although useful for identifying specific progenitors and precursors, this system has not been fully exploited as a source of cells to analyse erythropoiesis. Here, we establish a protocol in which characterised erythroblasts can be isolated in a scalable manner from differentiated embryoid bodies (EBs). Using transcriptional and epigenetic analysis, we demonstrate that this system faithfully recapitulates normal primitive erythropoiesis and fully reproduces the effects of natural and engineered mutations seen in primary cells obtained from mouse models. We anticipate this system to be of great value in reducing the time and costs of generating and maintaining mouse lines in a number of research scenarios.

Recapitulation of erythropoiesis in congenital dyserythropoietic anaemia type I (CDA-I) identifies defects in differentiation and nucleolar abnormalities

The investigation of inherited disorders of erythropoiesis has elucidated many of the principles underlying the production of normal red blood cells and how this is perturbed in human disease. Congenital Dyserythropoietic Anaemia type 1 (CDA-I) is a rare form of anaemia caused by mutations in two genes of unknown function: CDAN1 and CDIN1 (previously called C15orf41), whilst in some cases, the underlying genetic abnormality is completely unknown. Consequently, the pathways affected in CDA-I remain to be discovered. To enable detailed analysis of this rare disorder we have validated a culture system which recapitulates all of the cardinal haematological features of CDA-I, including the formation of the pathognomonic 'spongy' heterochromatin seen by electron microscopy. Using a variety of cell and molecular biological approaches we discovered that erythroid cells in this condition show a delay during terminal erythroid differentiation, associated with increased proliferation and widespread changes in chromatin accessibility. We also show that the proteins encoded by CDAN1 and CDIN1 are enriched in nucleoli which are structurally and functionally abnormal in CDA-I. Together these findings provide important pointers to the pathways affected in CDA-I which for the first time can now be pursued in the tractable culture system utilised here.

Reactivation of a developmentally silenced embryonic globin gene

The α- and β-globin loci harbor developmentally expressed genes, which are silenced throughout post-natal life. Reactivation of these genes may offer therapeutic approaches for the hemoglobinopathies, the most common single gene disorders. Here, we address mechanisms regulating the embryonically expressed α-like globin, termed ζ-globin. We show that in embryonic erythroid cells, the ζ-gene lies within a ~65 kb sub-TAD (topologically associating domain) of open, acetylated chromatin and interacts with the α-globin super-enhancer. By contrast, in adult erythroid cells, the ζ-gene is packaged within a small (~10 kb) sub-domain of hypoacetylated, facultative heterochromatin within the acetylated sub-TAD and that it no longer interacts with its enhancers. The ζ-gene can be partially re-activated by acetylation and inhibition of histone de-acetylases. In addition to suggesting therapies for severe α-thalassemia, these findings illustrate the general principles by which reactivation of developmental genes may rescue abnormalities arising from mutations in their adult paralogues.

The mouse alpha-globin cluster: a paradigm for studying genome regulation and organization

The mammalian globin gene clusters provide a paradigm for studying the relationship between genome structure and function. As blood stem cells undergo lineage specification and differentiation to form red blood cells, the chromatin structure and expression of the α-globin cluster change. The gradual activation of the α-globin genes in well-defined cell populations has enabled investigation of the structural and functional roles of its enhancers, promoters and boundary elements. Recent studies of gene regulatory processes involving these elements at the mouse α-globin cluster have brought new insights into the general principles underlying the three-dimensional structure of the genome and its relationship to gene expression throughout time.

Enhancers predominantly regulate gene expression during differentiation via transcription initiation

Gene transcription occurs via a cycle of linked events, including initiation, promoter-proximal pausing, and elongation of RNA polymerase II (Pol II). A key question is how transcriptional enhancers influence these events to control gene expression. Here, we present an approach that evaluates the level and change in promoter-proximal transcription (initiation and pausing) in the context of differential gene expression, genome-wide. This combinatorial approach shows that in primary cells, control of gene expression during differentiation is achieved predominantly via changes in transcription initiation rather than via release of Pol II pausing. Using genetically engineered mouse models, deleted for functionally validated enhancers of the α- and β-globin loci, we confirm that these elements regulate Pol II recruitment and/or initiation to modulate gene expression. Together, our data show that gene expression during differentiation is regulated predominantly at the level of initiation and that enhancers are key effectors of this process.

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.

Loss of Extreme Long-Range Enhancers in Human Neural Crest Drives a Craniofacial Disorder

Non-coding mutations at the far end of a large gene desert surrounding the SOX9 gene result in a human craniofacial disorder called Pierre Robin sequence (PRS). Leveraging a human stem cell differentiation model, we identify two clusters of enhancers within the PRS-associated region that regulate SOX9 expression during a restricted window of facial progenitor development at distances up to 1.45 Mb. Enhancers within the 1.45 Mb cluster exhibit highly synergistic activity that is dependent on the Coordinator motif. Using mouse models, we demonstrate that PRS phenotypic specificity arises from the convergence of two mechanisms: confinement of Sox9 dosage perturbation to developing facial structures through context-specific enhancer activity and heightened sensitivity of the lower jaw to Sox9 expression reduction. Overall, we characterize the longest-range human enhancers involved in congenital malformations, directly demonstrate that PRS is an enhanceropathy, and illustrate how small changes in gene expression can lead to morphological variation.

Genetic and functional insights into CDA-I prevalence and pathogenesis

BACKGROUND

Congenital dyserythropoietic anaemia type I (CDA-I) is a hereditary anaemia caused by biallelic mutations in the widely expressed genes CDAN1 and C15orf41. Little is understood about either protein and it is unclear in which cellular pathways they participate.

METHODS

Genetic analysis of a cohort of patients with CDA-I identifies novel pathogenic variants in both known causative genes. We analyse the mutation distribution and the predicted structural positioning of amino acids affected in Codanin-1, the protein encoded by CDAN1. Using western blotting, immunoprecipitation and immunofluorescence, we determine the effect of particular mutations on both proteins and interrogate protein interaction, stability and subcellular localisation.

RESULTS

We identify six novel CDAN1 mutations and one novel mutation in C15orf41 and uncover evidence of further genetic heterogeneity in CDA-I. Additionally, population genetics suggests that CDA-I is more common than currently predicted. Mutations are enriched in six clusters in Codanin-1 and tend to affect buried residues. Many missense and in-frame mutations do not destabilise the entire protein. Rather C15orf41 relies on Codanin-1 for stability and both proteins, which are enriched in the nucleolus, interact to form an obligate complex in cells.

CONCLUSION

Stability and interaction data suggest that C15orf41 may be the key determinant of CDA-I and offer insight into the mechanism underlying this disease. Both proteins share a common pathway likely to be present in a wide variety of cell types; however, nucleolar enrichment may provide a clue as to the erythroid specific nature of CDA-I. The surprisingly high predicted incidence of CDA-I suggests that better ascertainment would lead to improved patient care.

Dynamics of the 4D genome during in vivo lineage specification and differentiation

Mammalian gene expression patterns are controlled by regulatory elements, which interact within topologically associating domains (TADs). The relationship between activation of regulatory elements, formation of structural chromatin interactions and gene expression during development is unclear. Here, we present Tiled-C, a low-input chromosome conformation capture (3C) technique. We use this approach to study chromatin architecture at high spatial and temporal resolution through in vivo mouse erythroid differentiation. Integrated analysis of chromatin accessibility and single-cell expression data shows that regulatory elements gradually become accessible within pre-existing TADs during early differentiation. This is followed by structural re-organization within the TAD and formation of specific contacts between enhancers and promoters. Our high-resolution data show that these enhancer-promoter interactions are not established prior to gene expression, but formed gradually during differentiation, concomitant with progressive upregulation of gene activity. Together, these results provide new insight into the close, interdependent relationship between chromatin architecture and gene regulation during development.

An integrative view of the regulatory and transcriptional landscapes in mouse hematopoiesis

Thousands of epigenomic data sets have been generated in the past decade, but it is difficult for researchers to effectively use all the data relevant to their projects. Systematic integrative analysis can help meet this need, and the VISION project was established for validated systematic integration of epigenomic data in hematopoiesis. Here, we systematically integrated extensive data recording epigenetic features and transcriptomes from many sources, including individual laboratories and consortia, to produce a comprehensive view of the regulatory landscape of differentiating hematopoietic cell types in mouse. By using IDEAS as our integrative and discriminative epigenome annotation system, we identified and assigned epigenetic states simultaneously along chromosomes and across cell types, precisely and comprehensively. Combining nuclease accessibility and epigenetic states produced a set of more than 200,000 candidate cis-regulatory elements (cCREs) that efficiently capture enhancers and promoters. The transitions in epigenetic states of these cCREs across cell types provided insights into mechanisms of regulation, including decreases in numbers of active cCREs during differentiation of most lineages, transitions from poised to active or inactive states, and shifts in nuclease accessibility of CTCF-bound elements. Regression modeling of epigenetic states at cCREs and gene expression produced a versatile resource to improve selection of cCREs potentially regulating target genes. These resources are available from our VISION website to aid research in genomics and hematopoiesis.

ATR-16 syndrome: mechanisms linking monosomy to phenotype

Background

Deletions removing 100s–1000s kb of DNA, and variable numbers of poorly characterised genes, are often found in patients with a wide range of developmental abnormalities. In such cases, understanding the contribution of the deletion to an individual’s clinical phenotype is challenging.

Methods

Here, as an example of this common phenomenon, we analysed 41 patients with simple deletions of ~177 to ~2000 kb affecting one allele of the well-characterised, gene dense, distal region of chromosome 16 (16p13.3), referred to as ATR-16 syndrome. We characterised deletion extents and screened for genetic background effects, telomere position effect and compensatory upregulation of hemizygous genes.

Results

We find the risk of developmental and neurological abnormalities arises from much smaller distal chromosome 16 deletions (~400 kb) than previously reported. Beyond this, the severity of ATR-16 syndrome increases with deletion size, but there is no evidence that critical regions determine the developmental abnormalities associated with this disorder. Surprisingly, we find no evidence of telomere position effect or compensatory upregulation of hemizygous genes; however, genetic background effects substantially modify phenotypic abnormalities.

Conclusions

Using ATR-16 as a general model of disorders caused by CNVs, we show the degree to which individuals with contiguous gene syndromes are affected is not simply related to the number of genes deleted but depends on their genetic background. We also show there is no critical region defining the degree of phenotypic abnormalities in ATR-16 syndrome and this has important implications for genetic counselling.

Systematic integration of GATA transcription factors and epigenomes via IDEAS paints the regulatory landscape of hematopoietic cells

Members of the GATA family of transcription factors play key roles in the differentiation of specific cell lineages by regulating the expression of target genes. Three GATA factors play distinct roles in hematopoietic differentiation. In order to better understand how these GATA factors function to regulate genes throughout the genome, we are studying the epigenomic and transcriptional landscapes of hematopoietic cells in a model-driven, integrative fashion. We have formed the collaborative multi-lab VISION project to conduct ValIdated Systematic IntegratiON of epigenomic data in mouse and human hematopoiesis. The epigenomic data included nuclease accessibility in chromatin, CTCF occupancy, and histone H3 modifications for 20 cell types covering hematopoietic stem cells, multilineage progenitor cells, and mature cells across the blood cell lineages of mouse. The analysis used the Integrative and Discriminative Epigenome Annotation System (IDEAS), which learns all common combinations of features (epigenetic states) simultaneously in two dimensions-along chromosomes and across cell types. The result is a segmentation that effectively paints the regulatory landscape in readily interpretable views, revealing constitutively active or silent loci as well as the loci specifically induced or repressed in each stage and lineage. Nuclease accessible DNA segments in active chromatin states were designated candidate cis-regulatory elements in each cell type, providing one of the most comprehensive registries of candidate hematopoietic regulatory elements to date. Applications of VISION resources are illustrated for the regulation of genes encoding GATA1, GATA2, GATA3, and Ikaros. VISION resources are freely available from our website http://usevision.org.

A revised model for promoter competition based on multi-way chromatin interactions at the α-globin locus

Specific communication between gene promoters and enhancers is critical for accurate regulation of gene expression. However, it remains unclear how specific interactions between multiple regulatory elements contained within a single chromatin domain are coordinated. Recent technological advances which can detect multi-way chromatin interactions at single alleles can provide insights into how multiple regulatory elements cooperate or compete for transcriptional activation. Here, we use such an approach to investigate how interactions of the α-globin enhancers are distributed between multiple promoters in a mouse model in which the α-globin domain is extended to include several additional genes. Our data show that gene promoters do not form mutually exclusive interactions with enhancers, but all interact simultaneously in a single complex. These findings suggest that promoters do not structurally compete for interactions with enhancers, but form a regulatory hub structure, which is consistent with recent models of transcriptional activation occurring in non-membrane bound nuclear compartments.

Potential new approaches to the management of the Hb Bart's hydrops fetalis syndrome: the most severe form of α-thalassemia

The α-thalassemia trait, associated with deletions removing both α-globin genes from 1 chromosome (genotype ζ αα/ζ--), is common throughout Southeast Asia. Consequently, many pregnancies in couples of Southeast Asian origin carry a 1 in 4 risk of producing a fetus inheriting no functional α-globin genes (ζ--/ζ--), leading to hemoglobin (Hb) Bart's hydrops fetalis syndrome (BHFS). Expression of the embryonic α-globin genes (ζ-globin) is normally limited to the early stages of primitive erythropoiesis, and so when the ζ-globin genes are silenced, at ∼6 weeks of gestation, there should be no α-like globin chains to pair with the fetal γ-globin chains of Hb, which consequently form nonfunctional tetramers (γ₄) known as Hb Bart's. When deletions leave the ζ-globin gene intact, a low level of ζ-globin gene expression continues in definitive erythroid cells, producing small amounts of Hb Portland (ζ₂γ₂), a functional form of Hb that allows the fetus to survive up to the second or third trimester. Untreated, all affected individuals die at these stages of development. Prevention is therefore of paramount importance. With improvements in early diagnosis, intrauterine transfusion, and advanced perinatal care, there are now a small number of individuals with BHFS who have survived, with variable outcomes. A deeper understanding of the mechanism underlying the switch from ζ- to α-globin expression could enable persistence or reactivation of embryonic globin synthesis in definitive cells, thereby providing new therapeutic options for such patients.

Single-allele chromatin interactions identify regulatory hubs in dynamic compartmentalized domains

The promoters of mammalian genes are commonly regulated by multiple distal enhancers, which physically interact within discrete chromatin domains. How such domains form and how the regulatory elements within them interact in single cells is not understood. To address this we developed Tri-C, a new Chromosome Conformation Capture (3C) approach to identify concurrent chromatin interactions at individual alleles. Analysis by Tri-C reveals heterogeneous patterns of single-allele interactions between CTCF boundary elements, indicating that the formation of chromatin domains likely results from a dynamic process. Within these domains, we observe specific higher-order structures involving simultaneous interactions between multiple enhancers and promoters. Such regulatory hubs provide a structural basis for understanding how multiple cis-regulatory elements act together to establish robust regulation of gene expression.

A tissue-specific self-interacting chromatin domain forms independently of enhancer-promoter interactions

Self-interacting chromatin domains encompass genes and their cis-regulatory elements; however, the three-dimensional form a domain takes, whether this relies on enhancer-promoter interactions, and the processes necessary to mediate the formation and maintenance of such domains, remain unclear. To examine these questions, here we use a combination of high-resolution chromosome conformation capture, a non-denaturing form of fluorescence in situ hybridisation and super-resolution imaging to study a 70 kb domain encompassing the mouse α-globin regulatory locus. We show that this region forms an erythroid-specific, decompacted, self-interacting domain, delimited by frequently apposed CTCF/cohesin binding sites early in terminal erythroid differentiation, and does not require transcriptional elongation for maintenance of the domain structure. Formation of this domain does not rely on interactions between the α-globin genes and their major enhancers, suggesting a transcription-independent mechanism for establishment of the domain. However, absence of the major enhancers does alter internal domain interactions. Formation of a loop domain therefore appears to be a mechanistic process that occurs irrespective of the specific interactions within.

Robust CRISPR/Cas9 Genome Editing of the HUDEP-2 Erythroid Precursor Line Using Plasmids and Single-Stranded Oligonucleotide Donors

The study of cellular processes and gene regulation in terminal erythroid development has been greatly facilitated by the generation of an immortalised erythroid cell line derived from Human Umbilical Derived Erythroid Precursors, termed HUDEP-2 cells. The ability to efficiently genome edit HUDEP-2 cells and make clonal lines hugely expands their utility as the insertion of clinically relevant mutations allows study of potentially every genetic disease affecting red blood cell development. Additionally, insertion of sequences encoding short protein tags such as Strep, FLAG and Myc permits study of protein behaviour in the normal and disease state. This approach is useful to augment the analysis of patient cells as large cell numbers are obtainable with the additional benefit that the need for specific antibodies may be circumvented. This approach is likely to lead to insights into disease mechanisms and provide reagents to allow drug discovery. HUDEP-2 cells provide a favourable alternative to the existing immortalised erythroleukemia lines as their karyotype is much less abnormal. These cells also provide sufficient material for a broad range of analyses as it is possible to generate in vitro-differentiated erythroblasts in numbers 4-7 fold higher than starting cell numbers within 9-12 days of culture. Here we describe an efficient, robust and reproducible plasmid-based methodology to introduce short (<20 bp) DNA sequences into the genome of HUDEP-2 cells using the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR associated protein 9 Cas9 system combined with single-stranded oligodeoxynucleotide (ssODN) donors. This protocol produces genetically modified lines in ~30 days and could also be used to generate knock-out and knock-in mutations.

Molecular Basis and Genetic Modifiers of Thalassemia

Thalassemia is a disorder of hemoglobin characterized by reduced or absent production of one of the globin chains in human red blood cells with relative excess of the other. Impaired synthesis of β-globin results in β-thalassemia, whereas defective synthesis of α-globin leads to α-thalassemia. Despite being a monogenic disorder, thalassemia exhibits remarkable clinical heterogeneity that is directly related to the intracellular imbalance between α- and β-like globin chains. Novel insights into the genetic modifiers have contributed to the understanding of the correlation between genotype and phenotype and are being explored as therapeutic pathways to cure this life-limiting disease.

Robust detection of chromosomal interactions from small numbers of cells using low-input Capture-C

Chromosome conformation capture (3C) techniques are crucial to understanding tissue-specific regulation of gene expression, but current methods generally require large numbers of cells. This hampers the investigation of chromatin architecture in rare cell populations. We present a new low-input Capture-C approach that can generate high-quality 3C interaction profiles from 10 000-20 000 cells, depending on the resolution used for analysis. We also present a PCR-free, sequencing-free 3C technique based on NanoString technology called C-String. By comparing C-String and Capture-C interaction profiles we show that the latter are not skewed by PCR amplification. Furthermore, we demonstrate that chromatin interactions detected by Capture-C do not depend on the degree of cross-linking by performing experiments with varying formaldehyde concentrations.

Phenotypic and molecular characterization of a serum-free miniature erythroid differentiation system suitable for high-throughput screening and single-cell assays

In vitro erythroid differentiation systems are used to study the mechanisms underlying normal and abnormal erythropoiesis and to test the effects of various extracellular factors on erythropoiesis. The use of serum or conditioned medium in liquid cultures and the seeding of cultures with heterogeneous peripheral blood mononuclear cells confound the reproducibility of these systems. Newer erythroid differentiation culture systems have overcome some of these limitations by using a fully defined, serum-free medium and initiating cultures using purified CD34⁺ cells. Although widely used in bulk cultures, these protocols have not been rigorously tested in high-throughput or single-cell assays. Here, we describe a serum-free erythroid differentiation system suitable for small-scale and single-cell experiments. This system generates large numbers of terminally differentiated erythroid cells of very high purity. Here we have adapted this culture system to a 96-well format and have developed a protocol to grow erythroid colonies from single erythroid progenitors in minute culture volumes.

How to Tackle Challenging ChIP-Seq, with Long-Range Cross-Linking, Using ATRX as an Example

Chromatin immunoprecipitation coupled with high-throughput, next-generation DNA sequencing (ChIP-seq) has enabled researchers to establish the genome-wide patterns of chromatin modifications and binding of chromatin-associated proteins. Well-established protocols produce robust ChIP-seq data for many proteins by sequencing the DNA obtained following immunoprecipitation of fragmented chromatin using a wide range of specific antibodies. In general, the quality of these data mainly depends on the specificity and avidity of the antibody used. However, even using optimal antibodies, ChIP-seq can become more challenging when the protein associates with chromatin via protein-protein interactions rather than directly binding DNA. An example of such a protein is the alpha-thalassaemia mental retardation X-linked (ATRX) protein; a chromatin remodeler that associates with the histone chaperone DAXX, in the deposition of the replication-independent histone variant H3.3 and plays an important role in maintaining chromatin integrity. Inherited mutations of ATRX cause syndromal mental retardation (ATR-X Syndrome) whereas acquired mutations are associated with myelodysplasia, acute myeloid leukemia (ATMDS syndrome), and a range of solid tumors. Therefore, high quality ChIP-seq data have been needed to analyze the genome-wide distribution of ATRX, to advance our understanding of its normal role and to comprehend how mutations contribute to human disease. Here, we describe an optimized ChIP-seq protocol for ATRX which can also be used to produce high quality data sets for other challenging proteins which are indirectly associated with DNA and complement the ChIP-seq toolkit for genome-wide analyses of histone chaperon complexes and associated chromatin remodelers. Although not a focus of this chapter, we will also provide some insight for the analysis of the large dataset generated by ChIP-seq. Even though this protocol has been fully optimized for ATRX, it should also provide guidance for efficient ChIP-seq analysis, using the appropriate antibodies, for other proteins interacting indirectly with DNA.

Between form and function: the complexity of genome folding

It has been known for over a century that chromatin is not randomly distributed within the nucleus. However, the question of how DNA is folded and the influence of such folding on nuclear processes remain topics of intensive current research. A longstanding, unanswered question is whether nuclear organization is simply a reflection of nuclear processes such as transcription and replication, or whether chromatin is folded by independent mechanisms and this per se encodes function? Evidence is emerging that both may be true. Here, using the α-globin gene cluster as an illustrative model, we provide an overview of the most recent insights into the layers of genome organization across different scales and how this relates to gene activity.

DNA methylation of intragenic CpG islands depends on their transcriptional activity during differentiation and disease

The human genome contains ∼30,000 CpG islands (CGIs). While CGIs associated with promoters nearly always remain unmethylated, many of the ∼9,000 CGIs lying within gene bodies become methylated during development and differentiation. Both promoter and intragenic CGIs may also become abnormally methylated as a result of genome rearrangements and in malignancy. The epigenetic mechanisms by which some CGIs become methylated but others, in the same cell, remain unmethylated in these situations are poorly understood. Analyzing specific loci and using a genome-wide analysis, we show that transcription running across CGIs, associated with specific chromatin modifications, is required for DNA methyltransferase 3B (DNMT3B)-mediated DNA methylation of many naturally occurring intragenic CGIs. Importantly, we also show that a subgroup of intragenic CGIs is not sensitive to this process of transcription-mediated methylation and that this correlates with their individual intrinsic capacity to initiate transcription in vivo. We propose a general model of how transcription could act as a primary determinant of the patterns of CGI methylation in normal development and differentiation, and in human disease.

Tissue-specific CTCF-cohesin-mediated chromatin architecture delimits enhancer interactions and function in vivo

The genome is organized via CTCF-cohesin-binding sites, which partition chromosomes into 1-5 megabase (Mb) topologically associated domains (TADs), and further into smaller sub-domains (sub-TADs). Here we examined in vivo an ∼80 kb sub-TAD, containing the mouse α-globin gene cluster, lying within a ∼1 Mb TAD. We find that the sub-TAD is flanked by predominantly convergent CTCF-cohesin sites that are ubiquitously bound by CTCF but only interact during erythropoiesis, defining a self-interacting erythroid compartment. Whereas the α-globin regulatory elements normally act solely on promoters downstream of the enhancers, removal of a conserved upstream CTCF-cohesin boundary extends the sub-TAD to adjacent upstream CTCF-cohesin-binding sites. The α-globin enhancers now interact with the flanking chromatin, upregulating expression of genes within this extended sub-TAD. Rather than acting solely as a barrier to chromatin modification, CTCF-cohesin boundaries in this sub-TAD delimit the region of chromatin to which enhancers have access and within which they interact with receptive promoters.

The chromatin remodelling factor ATRX suppresses R-loops in transcribed telomeric repeats

ATRX is a chromatin remodelling factor found at a wide range of tandemly repeated sequences including telomeres (TTAGGG)_n ATRX mutations are found in nearly all tumours that maintain their telomeres via the alternative lengthening of telomere (ALT) pathway, and ATRX is known to suppress this pathway. Here, we show that recruitment of ATRX to telomeric repeats depends on repeat number, orientation and, critically, on repeat transcription. Importantly, the transcribed telomeric repeats form RNA-DNA hybrids (R-loops) whose abundance correlates with the recruitment of ATRX Here, we show loss of ATRX is also associated with increased R-loop formation. Our data suggest that the presence of ATRX at telomeres may have a central role in suppressing deleterious DNA secondary structures that form at transcribed telomeric repeats, and this may account for the increased DNA damage, stalling of replication and homology-directed repair previously observed upon loss of ATRX function.

An international registry of survivors with Hb Bart's hydrops fetalis syndrome

Hemoglobin (Hb) Bart's hydrops fetalis syndrome (BHFS) resulting from α0-thalassemia is considered a universally fatal disorder. However, over the last 3 decades, improvements in intrauterine interventions and perinatal intensive care have resulted in increasing numbers of BHFS survivors. We have initiated an international registry containing information on 69 patients, of which 31 are previously unpublished. In this perspective, we analyze the available clinical information to document the natural history of BHFS. In the future, once we have accrued sufficient cases, we aim to build on this study and provide information to allow counseling of at-risk couples. To date, 39 patients have survived beyond the age of 5 years, 18 of whom are now older than 10 years. Based on the available cases, we find evidence to suggest that intrauterine therapy provides benefits during the perinatal and neonatal period; however, it may not provide additional benefits to long-term growth and neurodevelopmental outcomes. Growth retardation is a major adverse long-term outcome among BHFS patients with ∼40% being severely affected in terms of weight and ∼50% in terms of height. There is also an increased risk of neurodevelopmental delay as we find 20% (11/55) of BHFS survivors suffer from a serious delay of ≥6 months. Most patients in the registry require lifelong transfusion and often have associated congenital abnormalities and comorbidities. This perspective is a first step in gathering information to allow provision of informed counseling on the predicted outcomes of affected babies.

Functional characterisation of cis-regulatory elements governing dynamic Eomes expression in the early mouse embryo

The T-box transcription factor (TF) Eomes is a key regulator of cell fate decisions during early mouse development. The cis-acting regulatory elements that direct expression in the anterior visceral endoderm (AVE), primitive streak (PS) and definitive endoderm (DE) have yet to be defined. Here, we identified three gene-proximal enhancer-like sequences (PSE_a, PSE_b and VPE) that faithfully activate tissue-specific expression in transgenic embryos. However, targeted deletion experiments demonstrate that PSE_a and PSE_b are dispensable, and only VPE is required for optimal Eomes expression in vivo. Embryos lacking this enhancer display variably penetrant defects in anterior-posterior axis orientation and DE formation. Chromosome conformation capture experiments reveal VPE-promoter interactions in embryonic stem cells (ESCs), prior to gene activation. The locus resides in a large (500 kb) pre-formed compartment in ESCs and activation during DE differentiation occurs in the absence of 3D structural changes. ATAC-seq analysis reveals that VPE, PSE_a and four additional putative enhancers display increased chromatin accessibility in DE that is associated with Smad2/3 binding coincident with transcriptional activation. By contrast, activation of the Eomes target genes Foxa2 and Lhx1 is associated with higher order chromatin reorganisation. Thus, diverse regulatory mechanisms govern activation of lineage specifying TFs during early development.

Summary: Expression of the mouse T-box factor Eomes is controlled by a key gene-proximal enhancer-like element, with changes in chromatin accessibility influencing its activity in definitive endoderm.

How best to identify chromosomal interactions: a comparison of approaches

Chromosome conformation capture (3C) methods are central to understanding the link between nuclear structure and function, and the physical interactions between distal regulatory elements and promoters. However, no one method is appropriate to address all biological questions, as each variant differs markedly in resolution, reproducibility, throughput and biases. A thorough appreciation of the strengths and weaknesses of each technique is critical when choosing the correct method for a specific application or for gauging how best to interpret different sources of data. In addition, the analysis method must be carefully considered, as this choice can profoundly affect the output. In this Review, we describe and compare the different available 3C-based approaches, with a focus on the analysis of mammalian genomes.

Genetic dissection of the α-globin super-enhancer in vivo

Many genes determining cell identity are regulated by clusters of Mediator-bound enhancer elements collectively referred to as super-enhancers. These super-enhancers have been proposed to manifest higher-order properties important in development and disease. Here we report a comprehensive functional dissection of one of the strongest putative super-enhancers in erythroid cells. By generating a series of mouse models, deleting each of the five regulatory elements of the α-globin super-enhancer individually and in informative combinations, we demonstrate that each constituent enhancer seems to act independently and in an additive fashion with respect to hematological phenotype, gene expression, chromatin structure and chromosome conformation, without clear evidence of synergistic or higher-order effects. Our study highlights the importance of functional genetic analyses for the identification of new concepts in transcriptional regulation.

MicroRNAs of the miR-290-295 Family Maintain Bivalency in Mouse Embryonic Stem Cells

Numerous developmentally regulated genes in mouse embryonic stem cells (ESCs) are marked by both active (H3K4me3)- and polycomb group (PcG)-mediated repressive (H3K27me3) histone modifications. This bivalent state is thought to be important for transcriptional poising, but the mechanisms that regulate bivalent genes and the bivalent state remain incompletely understood. Examining the contribution of microRNAs (miRNAs) to the regulation of bivalent genes, we found that the miRNA biogenesis enzyme DICER was required for the binding of the PRC2 core components EZH2 and SUZ12, and for the presence of the PRC2-mediated histone modification H3K27me3 at many bivalent genes. Genes that lost bivalency were preferentially upregulated at the mRNA and protein levels. Finally, reconstituting Dicer-deficient ESCs with ESC miRNAs restored bivalent gene repression and PRC2 binding at formerly bivalent genes. Therefore, miRNAs regulate bivalent genes and the bivalent state itself.

Transcription factors LRF and BCL11A independently repress expression of fetal hemoglobin

Genes encoding human β-type globin undergo a developmental switch from embryonic to fetal to adult-type expression. Mutations in the adult form cause inherited hemoglobinopathies or globin disorders, including sickle cell disease and thalassemia. Some experimental results have suggested that these diseases could be treated by induction of fetal-type hemoglobin (HbF). However, the mechanisms that repress HbF in adults remain unclear. We found that the LRF/ZBTB7A transcription factor occupies fetal γ-globin genes and maintains the nucleosome density necessary for γ-globin gene silencing in adults, and that LRF confers its repressive activity through a NuRD repressor complex independent of the fetal globin repressor BCL11A. Our study may provide additional opportunities for therapeutic targeting in the treatment of hemoglobinopathies.

Understanding α-globin gene regulation and implications for the treatment of β-thalassemia

Over the past three decades, a vast amount of new information has been uncovered describing how the globin genes are regulated. This knowledge has provided significant insights into the general understanding of the regulation of human genes. It is now known that molecular defects within and around the α- and β-globin genes, as well as in the distant regulatory elements, can cause thalassemia. Unbalanced production of globin chains owing to defective synthesis of one, and the continued unopposed synthesis of another, is the central causative factor in the cellular pathology and pathophysiology of thalassemia. A large body of clinical, genetic, and experimental evidence suggests that altering globin chain imbalance by reducing the production of α-globin synthesis ameliorates the disease severity in patients with β-thalassemia. With the development of new genetic-based therapeutic tools that have a potential to decrease the expression of a selected gene in a tissue-specific manner, the possibility of decreasing expression of the α-globin gene to improve the clinical severity of β-thalassemia could become a reality.

Multiplexed analysis of chromosome conformation at vastly improved sensitivity

Since methods for analysing chromosome conformation in mammalian cells are either low resolution or low throughput and are technically challenging they are not widely used outside of specialised laboratories. We have re-designed the Capture-C method producing a new approach, called next generation (NG) Capture-C. This produces unprecedented levels of sensitivity and reproducibility and can be used to analyse many genetic loci and samples simultaneously. Importantly, high-resolution data can be produced on as few as 100,000 cells and SNPs can be used to generate allele specific tracks. The method is straightforward to perform and should therefore greatly facilitate the task of linking SNPs identified by genome wide association studies with the genes they influence. The complete and detailed protocol presented here, with new publicly available tools for library design and data analysis, will allow most laboratories to analyse chromatin conformation at levels of sensitivity and throughput that were previously impossible.

Suppression of the alternative lengthening of telomere pathway by the chromatin remodelling factor ATRX

Fifteen per cent of cancers maintain telomere length independently of telomerase by the homologous recombination (HR)-associated alternative lengthening of telomeres (ALT) pathway. A unifying feature of these tumours are mutations in ATRX. Here we show that expression of ectopic ATRX triggers a suppression of the pathway and telomere shortening. Importantly ATRX-mediated ALT suppression is dependent on the histone chaperone DAXX. Re-expression of ATRX is associated with a reduction in replication fork stalling, a known trigger for HR and loss of MRN from telomeres. A G-quadruplex stabilizer partially reverses the effect of ATRX, inferring ATRX may normally facilitate replication through these sequences that, if they persist, promote ALT. We propose that defective telomere chromatinization through loss of ATRX promotes the persistence of aberrant DNA secondary structures, which in turn present a barrier to DNA replication, leading to replication fork stalling, collapse, HR and subsequent recombination-mediated telomere synthesis in ALT cancers.

ATRX, a chromatin remodelling factor, is mutated in cancers that maintain telomere length by alternative lengthening of telomeres (ALT). Here, the authors show that ectopic expression of ATRX triggers telomere shortening, ALT suppression and reduced replication fork stalling.