Chromosome looping at the human alpha-globin locus is mediated via the major upstream regulatory element (HS -40)

Ancestry of the major long-range regulatory site of the α-globin genes in the Portuguese population with the common 3.7 kb α-thalassemia deletion

BACKGROUND

The α-Major Regulatory Element (α-MRE), also known as HS-40, is located upstream of the α-globin gene cluster and has a crucial role in the long-range regulation of the α-globin gene expression. This enhancer is polymorphic and several haplotypes were identified in different populations, with haplotype D almost exclusively found in African populations. The purpose of this research was to identify the HS-40 haplotype associated with the 3.7 kb α-thalassemia deletion (-α3.7del) in the Portuguese population, and determine its ancestry and influence on patients' hematological phenotype.

METHODS AND RESULTS

We selected 111 Portuguese individuals previously analyzed by Gap-PCR to detect the presence of the -α3.7del: 50 without the -α3.7del, 34 heterozygous and 27 homozygous for the -α3.7del. The HS-40 region was amplified by PCR followed by Sanger sequencing. Four HS-40 haplotypes were found (A to D). The distribution of HS-40 haplotypes and genotypes are significantly different between individuals with and without the -α3.7del, being haplotype D and genotype AD the most prevalent in patients with this deletion in homozygosity. Furthermore, multiple correspondence analysis revealed that individuals without the -α3.7del are grouped with other European populations, while samples with the -α3.7del are separated from these and found more closely related to the African population.

CONCLUSION

This study revealed for the first time an association of the HS-40 haplotype D with the -α3.7del in the Portuguese population, and its likely African ancestry. These results may have clinical importance as in vitro analysis of haplotype D showed a decrease in its enhancer activity on α-globin gene.

Interplay between α-thalassemia and β-hemoglobinopathies: Translating genotype-phenotype relationships into therapies

α-Thalassemia represents one of the most important genetic modulators of β-hemoglobinopathies. During this last decade, the ongoing interest in characterizing genotype-phenotype relationships has yielded incredible insights into α-globin gene regulation and its impact on β-hemoglobinopathies. In this review, we provide a holistic update on α-globin gene expression stemming from DNA to RNA to protein, as well as epigenetic mechanisms that can impact gene expression and potentially influence phenotypic outcomes. Here, we highlight defined α-globin targeted strategies and rationalize the use of distinct molecular targets based on the restoration of balanced α/β-like globin chain synthesis. Considering the therapies that either increase β-globin synthesis or reactivate γ-globin gene expression, the modulation of α-globin chains as a disease modifier for β-hemoglobinopathies still remains largely uncharted in clinical studies.

Enhancer Function in the 3D Genome

In this review, we consider various aspects of enhancer functioning in the context of the 3D genome. Particular attention is paid to the mechanisms of enhancer-promoter communication and the significance of the spatial juxtaposition of enhancers and promoters in 3D nuclear space. A model of an activator chromatin compartment is substantiated, which provides the possibility of transferring activating factors from an enhancer to a promoter without establishing direct contact between these elements. The mechanisms of selective activation of individual promoters or promoter classes by enhancers are also discussed.

Functional analysis of three new alpha-thalassemia deletions involving MCS-R2 reveals the presence of an additional enhancer element in the 5' boundary region

We report three novel deletions involving the Multispecies Conserved Sequences (MCS) R2, also known as the Major Regulative Element (MRE), in patients showing the α-thalassemia phenotype. The three new rearrangements showed peculiar positions of the breakpoints. 1) The (αα)ES is a telomeric 110 kb deletion ending inside the MCS-R3 element. 2) The (αα)FG, 984 bp-long, ends 51 bp upstream to MCS-R2; both are associated with a severe α-thalassemia phenotype. 3) The (αα)CT, 5058 bp-long starts at position +93 of MCS-R2 and is the only one associated to a mild α-thalassemia phenotype. To understand the specific role of different segments of the MCS-R2 element and of its boundary regions we carried out transcriptional and expression analysis. Transcriptional analysis of patients' reticulocytes showed that (αα)ES was unable to produce α2-globin mRNA, while a high level of expression of the α2-globin genes (56%) was detected in (αα)CT deletion, characterized by the presence of the first 93 bp of MCS-R2. Expression analysis of constructs containing breakpoints and boundary regions of the deletions (αα)CT and (αα)FG, showed comparable activity both for MCS-R2 and the boundary region (-682/-8). Considering that the (αα)CT deletion, almost entirely removing MCS-R2, has a less severe phenotype than the (αα)FG α0thalassemia deletion, removing both MCS-R2 almost entirely and an upstream 679 bp, we infer for the first time that an enhancer element must exist in this region that helps to increase the expression of the α-globin genes. The genotype-phenotype relationship of other previously published MCS-R2 deletions strengthened our hypothesis.

Molecular Basis and Genetic Modifiers of Thalassemia

Thalassemia syndromes are common monogenic disorders and represent a significant health issue worldwide. In this review, the authors elaborate on fundamental genetic knowledge about thalassemias, including the structure and location of globin genes, the production of hemoglobin during development, the molecular lesions causing α-, β-, and other thalassemia syndromes, the genotype-phenotype correlation, and the genetic modifiers of these conditions. In addition, they briefly discuss the molecular techniques applied for diagnosis and innovative cell and gene therapy strategies to cure these conditions.

Segregation of α- and β-Globin Gene Cluster in Vertebrate Evolution: Chance or Necessity?

The review is devoted to the patterns of evolution of α- and β-globin gene domains. A hypothesis is presented according to which segregation of the ancestral cluster of α/β-globin genes in Amniota occurred due to the performance by α-globins and β-globins of non-canonical functions not related to oxygen transport.

Recent Approaches for Manipulating Globin Gene Expression in Treating Hemoglobinopathies

Tissue oxygenation throughout life depends on the activity of hemoglobin (Hb) one of the hemeproteins that binds oxygen in the lungs and secures its delivery throughout the body. Hb is composed of four monomers encoded by eight different genes the expression of which is tightly regulated during development, resulting in the formation of distinct hemoglobin tetramers in each developmental stage. Mutations that alter hemoglobin structure or its regulated expression result in a large group of diseases typically referred to as hemoglobinopathies that are amongst the most common genetic defects worldwide. Unprecedented efforts in the last decades have partially unraveled the complex mechanisms that control globin gene expression throughout development. In addition, genome wide association studies have revealed protective genetic traits capable of ameliorating the clinical manifestations of severe hemoglobinopathies. This knowledge has fueled the exploration of innovative therapeutic approaches aimed at modifying the genome or the epigenome of the affected cells to either restore hemoglobin function or to mimic the effect of protective traits. Here we describe the key steps that control the switch in gene expression that concerns the different globin genes during development and highlight the latest efforts in altering globin regulation for therapeutic purposes.

Thalassemias: From gene to therapy

Thalassemias (α, β, γ, δ, δβ, and εγδβ) are the most common genetic disorders worldwide and constitute a heterogeneous group of hereditary diseases characterized by the deficient synthesis of one or more hemoglobin (Hb) chain(s). This leads to the accumulation of unstable non-thalassemic Hb chains, which precipitate and cause intramedullary destruction of erythroid precursors and premature lysis of red blood cells (RBC) in the peripheral blood. Non-thalassemic Hbs display high oxygen affinity and no cooperativity. Thalassemias result from many different genetic and molecular defects leading to either severe or clinically silent hematologic phenotypes. Thalassemias α and β are particularly diffused in the regions spanning from the Mediterranean basin through the Middle East, Indian subcontinent, Burma, Southeast Asia, Melanesia, and the Pacific Islands, whereas δβ-thalassemia is prevalent in some Mediterranean regions including Italy, Greece, and Turkey. Although in the world thalassemia and malaria areas overlap apparently, the RBC protection against malaria parasites is openly debated. Here, we provide an overview of the historical, geographic, genetic, structural, and molecular pathophysiological aspects of thalassemias. Moreover, attention has been paid to molecular and epigenetic pathways regulating globin gene expression and globin switching. Challenges of conventional standard treatments, including RBC transfusions and iron chelation therapy, splenectomy and hematopoietic stem cell transplantation from normal donors are reported. Finally, the progress made by rapidly evolving fields of gene therapy and gene editing strategies, already in pre-clinical and clinical evaluation, and future challenges as novel curative treatments for thalassemia are discussed.

The role of globins in cardiovascular physiology

Globin proteins exist in every cell type of the vasculature, from erythrocytes to endothelial cells, vascular smooth muscle cells, and peripheral nerve cells. Many globin subtypes are also expressed in muscle tissues (including cardiac and skeletal muscle), in other organ-specific cell types, and in cells of the central nervous system. The ability of each of these globins to interact with molecular oxygen (O₂) and nitric oxide (NO) is preserved across these contexts. Endothelial α-globin is an example of extra-erythrocytic globin expression. Other globins, including myoglobin, cytoglobin, and neuroglobin are observed in other vascular tissues. Myoglobin is observed primarily in skeletal muscle and smooth muscle cells surrounding the aorta or other large arteries. Cytoglobin is found in vascular smooth muscle but can also be expressed in non-vascular cell types, especially in oxidative stress conditions after ischemic insult. Neuroglobin was first observed in neuronal cells, and its expression appears to be restricted mainly to the central and peripheral nervous systems. Brain and central nervous system neurons expressing neuroglobin are positioned close to many arteries within the brain parenchyma and can control smooth muscle contraction and, thus, tissue perfusion and vascular reactivity. Overall, reactions between NO and globin heme-iron contribute to vascular homeostasis by regulating vasodilatory NO signals and scaveging reactive species in cells of the mammalian vascular system. Here, we discuss how globin proteins affect vascular physiology with a focus on NO biology, and offer perspectives for future study of these functions.

Co-Regulated Genes and Gene Clusters

There are many co-regulated genes in eukaryotic cells. The coordinated activation or repression of such genes occurs at specific stages of differentiation, or under the influence of external stimuli. As a rule, co-regulated genes are dispersed in the genome. However, there are also gene clusters, which contain paralogous genes that encode proteins with similar functions. In this aspect, they differ significantly from bacterial operons containing functionally linked genes that are not paralogs. In this review, we discuss the reasons for the existence of gene clusters in vertebrate cells and propose that clustering is necessary to ensure the possibility of selective activation of one of several similar genes.

The two waves in single-cell 3D genomics

For decades, biochemical methods for the analysis of genome structure and function provided cell-population-averaged data that allowed general principles and tendencies to be disclosed. Microscopy-based studies, which immanently involve single-cell analysis, did not provide sufficient spatial resolution to investigate the particularly small details of 3D genome folding. Nevertheless, these studies demonstrated that mutual positions of chromosome territories within cell nuclei and individual genomic loci within chromosomal territories can vary significantly in individual cells. The development of new technologies in biochemistry and the advent of super-resolution microscopy in the last decade have made possible the full-scale study of 3D genome organization in individual cells. Maps of the 3D genome build based on C-data and super-resolution microscopy are highly consistent and, therefore, biologically relevant. The internal structures of individual chromosomes, loci, and topologically associating domains (TADs) are resolved as well as cell-cycle dynamics. 3D modeling allows one to investigate the physical mechanisms underlying genome folding. Finally, joint profiling of genome topology and epigenetic features will allow 3D genomics to handle complex cell-to-cell heterogeneity. In this review, we summarize the present state of studies into 3D genome organization in individual cells, analyze the technical problems of single-cell studies, and outline perspectives of 3D genomics.

A Small Key for a Heavy Door: Genetic Therapies for the Treatment of Hemoglobinopathies

Throughout the past decades, the search for a treatment for severe hemoglobinopathies has gained increased interest within the scientific community. The discovery that ɤ-globin expression from intact HBG alleles complements defective HBB alleles underlying β-thalassemia and sickle cell disease, has provided a promising opening for research directed at relieving ɤ-globin repression mechanisms and, thereby, improve clinical outcomes for patients. Various gene editing strategies aim to reverse the fetal-to-adult hemoglobin switch to up-regulate ɤ-globin expression through disabling either HBG repressor genes or repressor binding sites in the HBG promoter regions. In addition to these HBB mutation-independent strategies involving fetal hemoglobin (HbF) synthesis de-repression, the expanding genome editing toolkit is providing increased accuracy to HBB mutation-specific strategies encompassing adult hemoglobin (HbA) restoration for a personalized treatment of hemoglobinopathies. Moreover, besides genome editing, more conventional gene addition strategies continue under investigation to restore HbA expression. Together, this research makes hemoglobinopathies a fertile ground for testing various innovative genetic therapies with high translational potential. Indeed, the progressive understanding of the molecular clockwork underlying the hemoglobin switch together with the ongoing optimization of genome editing tools heightens the prospect for the development of effective and safe treatments for hemoglobinopathies. In this context, clinical genetics plays an equally crucial role by shedding light on the complexity of the disease and the role of ameliorating genetic modifiers. Here, we cover the most recent insights on the molecular mechanisms underlying hemoglobin biology and hemoglobinopathies while providing an overview of state-of-the-art gene editing platforms. Additionally, current genetic therapies under development, are equally discussed.

Analysis of rare thalassemia caused by HS-40 regulatory site deletion

ABSTRACT Objectives: To investigate the effect of HS-40 regulatory site deletion on α-globin gene expression and its clinical significance. Methods: Venous blood samples of subjects were analyzed using a hematology analyzer and high- performance liquid chromatography; fetal cord blood was analyzed by a capillary electrophoresis analyzer. Gap-polymerase chain reaction (PCR), reverse dot blot (RDB), and multiple-link-dependent probe amplification (MLPA) were used for genotyping of thalassemia. Results: The proband was POLR3 K, HS-40 heterozygous deletion; the proband's wife was -SEA/αα; the fetus was POLR3 K, HS-40 heterozygous deletion combined with -SEA deletion; all of them had microcytic hypochromic anemia. Fetal umbilical cord blood electrophoresis revealed a suspected Hb Bart's band to be 88.4%, and the fetus was, thus, diagnosed as Hb Bart's fetus. The red blood cell parameters of the sporadic case showed that he had microcytic hypochromic anemia. Hemoglobin (Hb) electrophoresis analysis showed Hb H to be 5.3%, leading to a diagnosis of Hb H disease. Gap-PCR and RDB identified the genotype to be -α3.7/αα, β^A/β^A. MLPA detected heterozygous deletion or -α3.7 deletion on one allele and deletion of the HS-40 regulatory site on the other allele. Conclusion: The deletion of HS-40 regulatory site reduced expression of α-globin. HS-40 heterozygous deletion manifested as mild anemia, which was of microcytic hypochromic type. When compounded with -α3.7/αα, it manifested as Hb H disease; and when compounded with -SEA/αɑ, it manifested as Hb Bart's fetus.

MZF1 regulates α-globin gene transcription via long-range interactions in erythroid differentiation

Precise spatiotemporal gene expression regulation is crucial for human erythropoiesis. However, dramatic changes in the chromatin structure and transcriptome involved in α-globin gene expression during erythropoiesis still not fully understand. To identify candidate regulators for α-globin gene regulation, we carried out an integrated approach by integrating publicly available transcriptomic and epigenomic data. We computed active enhancers by overlapping enriched regions marked with H3K4me1 and H3K27ac and correlated their activity with mRNA expression. Next, we cataloged potential transcription factors via de novo motif analysis. We highlighted the discovery of potential novel transcription factor MZF1 of the α-globin gene in erythroid differentiation. To validate the role of MZF1, we quantified the expression level of MZF1 and α-globin gene in HSPCs, early erythroid progenitors and late erythroid precursors cells. Both the mRNA and protein expression patterns of MZF1 were consistent with the α-globin gene. Also, the qPCR result showed that the expression of the α-globin gene was significantly increased by the MZF1 overexpression. To further investigate the role of MZF1 regulating α-globin gene transcriptional activity during erythroid differentiation, we performed ChIP-qPCR at the α-globin locus. Our results showed that MZF1 recruitment both at 4 upstream HS sites and α-globin gene promoter in erythroid precursor cells. To determine the importance of the MZF1 to enhancer-promoter interaction at the α-globin locus, we compared interaction frequency before and after knockdown of MZF1 by chromosome conformation capture (3C) assay. Upon MZF1 depletion, both the expression of the α-globin gene and all 3C signals were significantly decreased. Taken together, MZF1 plays an important role in regulating α-globin gene expression by binding to long-region enhancers and α-globin gene promoter and facilitates the organization of specific 3D chromatin architecture in erythroid differentiation.

The 3D Genome as a Target for Anticancer Therapy

The role of 3D genome organization in the precise regulation of gene expression is well established. Accordingly, the mechanistic connections between 3D genome alterations and disease development are becoming increasingly apparent. This opinion article provides a snapshot of our current understanding of the 3D genome alterations associated with cancers. We discuss potential connections of the 3D genome and cancer transcriptional addiction phenomenon as well as molecular mechanisms of action of 3D genome-disrupting drugs. Finally, we highlight issues and perspectives raised by the discovery of the first pharmaceutical strongly affecting 3D genome organization.

Characterization of two novel Alu element-mediated α-globin gene cluster deletions causing α0-thalassemia by targeted next-generation sequencing

α-thalassemia is an inherited blood disorder commonly caused by deletions or point mutations involving one or both α-globin genes. Recent studies shed new light on the critical role of upstream enhancers multi-species conserved sequences (MCSs) in the ordered regulation of α-globin gene expression. Herein, we reported two unrelated probands with deletions in α-globin genes and MCSs, respectively. The proband from Family A is a compound heterozygote carrying a known α⁺ mutation (-α^3.7) and a novel 60.2 kb deletion causing the absence of both α-globin genes. The proband from Family B, on the other hand, is a compound heterozygote with a known α⁰ mutation (--^SEA) and a novel deletion involving only upstream regulatory elements MCS-R1, R2 and R3, while the α-globin genes remain intact. Notably, both these two patients suffered varied extent of anemia, indicating that the loss of enhancer elements could equally lead to reduced synthesis of α-globin. Upon these observations, we then confirmed the exact breakpoints of these two novel deletions using a targeted next-generation sequencing (NGS) previously established by our group, which may enable further elucidation of the rearrangement mechanisms on these deletions and functional dissection of MCSs. Taken together, our study reports a reliable NGS-based molecular screening approach for accurate identification of copy number variations (CNVs) in the α-globin cluster and the genetic diagnosis of these two probands may help to extend the spectrum of α-thalassemia mutations in Chinese population.

Mouse model of the human serotonin transporter-linked polymorphic region

Genetic factors play a significant role in risk for mood and anxiety disorders. Polymorphisms in genes that regulate the brain monoamine systems, such as catabolic enzymes and transporters, are attractive candidates for being risk factors for emotional disorders given the weight of evidence implicating monoamines involvement in these conditions. Several common genetic variants have been identified in the human serotonin transporter (5-HTT) gene, including a repetitive sequence located in the promoter region of the locus called the serotonin transporter-linked polymorphic region (5-HTT-LPR). This polymorphism has been associated with a number of mental traits in both humans and primates, including depression, neuroticism, and harm avoidance. Some, but not all, studies found a link between the polymorphism and 5-HTT levels, leaving open the question of whether the polymorphism affects risk for mental traits via changes in 5-HTT expression. To investigate the impact of the polymorphism on gene expression, serotonin homeostasis, and behavioral traits, we set out to develop a mouse model of the human 5-HTT-LPR. Here we describe the creation and characterization of a set of mouse lines with single-copy human transgenes carrying the short and long 5-HTT-LPR variants. Although we were not able to detect differences in expression between the short and long variants, we encountered several technical issues concerning the design of our humanized mice that are likely to have influenced our findings. Our study serves as a cautionary note for future studies aimed at studying human transgene regulation in the context of the living mouse.

Heterogeneous Loop Model to Infer 3D Chromosome Structures from Hi-C

Adapting a well-established formalism in polymer physics, we develop a minimalist approach to infer three-dimensional folding of chromatin from Hi-C data. The three-dimensional chromosome structures generated from our heterogeneous loop model (HLM) are used to visualize chromosome organizations that can substantiate the measurements from fluorescence in situ hybridization, chromatin interaction analysis by paired-end tag sequencing, and RNA-seq signals. We demonstrate the utility of the HLM with several case studies. Specifically, the HLM-generated chromosome structures, which reproduce the spatial distribution of topologically associated domains from fluorescence in situ hybridization measurement, show the phase segregation between two types of topologically associated domains explicitly. We discuss the origin of cell-type-dependent gene-expression level by modeling the chromatin globules of α-globin and SOX2 gene loci for two different cell lines. We also use the HLM to discuss how the chromatin folding and gene-expression level of Pax6 loci, associated with mouse neural development, are modulated by interactions with two enhancers. Finally, HLM-generated structures of chromosome 19 of mouse embryonic stem cells, based on single-cell Hi-C data collected over each cell-cycle phase, visualize changes in chromosome conformation along the cell-cycle. Given a contact frequency map between chromatic loci supplied from Hi-C, HLM is a computationally efficient and versatile modeling tool to generate chromosome structures that can complement interpreting other experimental data.

Humanising the mouse genome piece by piece

To better understand human health and disease, researchers create a wide variety of mouse models that carry human DNA. With recent advances in genome engineering, the targeted replacement of mouse genomic regions with orthologous human sequences has become increasingly viable, ranging from finely tuned humanisation of individual nucleotides and amino acids to the incorporation of many megabases of human DNA. Here, we examine emerging technologies for targeted genomic humanisation, we review the spectrum of existing genomically humanised mouse models and the insights such models have provided, and consider the lessons learned for designing such models in the future.

Generation of transgenic mice has become routine in studying gene function and disease mechanisms, but often this is not enough to fully understand human biology. Here, the authors review the current state of the art of targeted genomic humanisation strategies and their advantages over classic approaches.

CpG binding protein (CFP1) occupies open chromatin regions of active genes, including enhancers and non-CpG islands

BACKGROUND

The mechanism by which protein complexes interact to regulate the deposition of post-translational modifications of histones remains poorly understood. This is particularly important at regulatory regions, such as CpG islands (CGIs), which are known to recruit Trithorax (TrxG) and Polycomb group proteins. The CxxC zinc finger protein 1 (CFP1, also known as CGBP) is a subunit of the TrxG SET1 protein complex, a major catalyst of trimethylation of H3K4 (H3K4me3).

RESULTS

Here, we used ChIP followed by high-throughput sequencing (ChIP-seq) to analyse genomic occupancy of CFP1 in two human haematopoietic cell types. We demonstrate that CFP1 occupies CGIs associated with active transcription start sites (TSSs), and is mutually exclusive with H3K27 trimethylation (H3K27me3), a marker of polycomb repressive complex 2. Strikingly, rather than being restricted to active CGI TSSs, CFP1 also occupies a substantial fraction of active non-CGI TSSs and enhancers of transcribed genes. However, relative to other TrxG subunits, CFP1 was specialised to TSSs. Finally, we found enrichment of CpG-containing DNA motifs in CFP1 peaks at CGI promoters.

CONCLUSIONS

We found that CFP1 is not solely recruited to CpG islands as it was originally defined, but also other regions including non-CpG island promoters and enhancers.

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.

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.

Computational construction of 3D chromatin ensembles and prediction of functional interactions of alpha-globin locus from 5C data

Conformation capture technologies measure frequencies of interactions between chromatin regions. However, understanding gene-regulation require knowledge of detailed spatial structures of heterogeneous chromatin in cells. Here we describe the nC-SAC (n-Constrained-Self Avoiding Chromatin) method that transforms experimental interaction frequencies into 3D ensembles of chromatin chains. nC-SAC first distinguishes specific from non-specific interaction frequencies, then generates 3D chromatin ensembles using identified specific interactions as spatial constraints. Application to α-globin locus shows that these constraints (∼20%) drive the formation of ∼99% all experimentally captured interactions, in which ∼30% additional to the imposed constraints is found to be specific. Many novel specific spatial contacts not captured by experiments are also predicted. A subset, of which independent ChIA-PET data are available, is validated to be RNAPII-, CTCF-, and RAD21-mediated. Their positioning in the architectural context of imposed specific interactions from nC-SAC is highly important. Our results also suggest the presence of a many-body structural unit involving α-globin gene, its enhancers, and POL3RK gene for regulating the expression of α-globin in silent cells.

The free-energy cost of interaction between DNA loops

From the viewpoint of thermodynamics, the formation of DNA loops and the interaction between them, which are all non-equilibrium processes, result in the change of free energy, affecting gene expression and further cell-to-cell variability as observed experimentally. However, how these processes dissipate free energy remains largely unclear. Here, by analyzing a mechanic model that maps three fundamental topologies of two interacting DNA loops into a 4-state model of gene transcription, we first show that a longer DNA loop needs more mean free energy consumption. Then, independent of the type of interacting two DNA loops (nested, side-by-side or alternating), the promotion between them always consumes less mean free energy whereas the suppression dissipates more mean free energy. More interestingly, we find that in contrast to the mechanism of direct looping between promoter and enhancer, the facilitated-tracking mechanism dissipates less mean free energy but enhances the mean mRNA expression, justifying the facilitated-tracking hypothesis, a long-standing debate in biology. Based on minimal energy principle, we thus speculate that organisms would utilize the mechanisms of loop-loop promotion and facilitated tracking to survive in complex environments. Our studies provide insights into the understanding of gene expression regulation mechanism from the view of energy consumption.

Activation of the alpha-globin gene expression correlates with dramatic upregulation of nearby non-globin genes and changes in local and large-scale chromatin spatial structure

BACKGROUND

In homeotherms, the alpha-globin gene clusters are located within permanently open genome regions enriched in housekeeping genes. Terminal erythroid differentiation results in dramatic upregulation of alpha-globin genes making their expression comparable to the rRNA transcriptional output. Little is known about the influence of the erythroid-specific alpha-globin gene transcription outburst on adjacent, widely expressed genes and large-scale chromatin organization. Here, we have analyzed the total transcription output, the overall chromatin contact profile, and CTCF binding within the 2.7 Mb segment of chicken chromosome 14 harboring the alpha-globin gene cluster in cultured lymphoid cells and cultured erythroid cells before and after induction of terminal erythroid differentiation.

RESULTS

We found that, similarly to mammalian genome, the chicken genomes is organized in TADs and compartments. Full activation of the alpha-globin gene transcription in differentiated erythroid cells is correlated with upregulation of several adjacent housekeeping genes and the emergence of abundant intergenic transcription. An extended chromosome region encompassing the alpha-globin cluster becomes significantly decompacted in differentiated erythroid cells, and depleted in CTCF binding and CTCF-anchored chromatin loops, while the sub-TAD harboring alpha-globin gene cluster and the upstream major regulatory element (MRE) becomes highly enriched with chromatin interactions as compared to lymphoid and proliferating erythroid cells. The alpha-globin gene domain and the neighboring loci reside within the A-like chromatin compartment in both lymphoid and erythroid cells and become further segregated from the upstream gene desert upon terminal erythroid differentiation.

CONCLUSIONS

Our findings demonstrate that the effects of tissue-specific transcription activation are not restricted to the host genomic locus but affect the overall chromatin structure and transcriptional output of the encompassing topologically associating domain.

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.

The Functionality and Evolution of Eukaryotic Transcriptional Enhancers

Enhancers regulate precise spatial and temporal patterns of gene expression in eukaryotes and, moreover, evolutionary changes in these modular cis-regulatory elements may represent the predominant genetic basis for phenotypic evolution. Here, we review approaches to identify and functionally analyze enhancers and their transcription factor binding sites, including assay for transposable-accessible chromatin-sequencing (ATAC-Seq) and clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9, respectively. We also explore enhancer functionality, including how transcription factor binding sites combine to regulate transcription, as well as research on shadow and super enhancers, and how enhancers can act over great distances and even in trans. Finally, we discuss recent theoretical and empirical data on how transcription factor binding sites and enhancers evolve. This includes how the function of enhancers is maintained despite the turnover of transcription factor binding sites as well as reviewing studies where mutations in enhancers have been shown to underlie morphological change.

The regulation of human globin promoters by CCAAT box elements and the recruitment of NF-Y

CCAAT boxes are motifs found within the proximal promoter of many genes, including the human globin genes. The highly conserved nature of CCAAT box motifs within the promoter region of both α-like and β-like globin genes emphasises the functional importance of the CCAAT sequence in globin gene regulation. Mutations within the β-globin CCAAT box result in β-thalassaemia, while mutations within the distal γ-globin CCAAT box cause the Hereditary Persistence of Foetal Haemoglobin, a benign condition which results in continued γ-globin expression during adult life. Understanding the transcriptional regulation of the globin genes is of particular interest, as reactivating the foetal γ-globin gene alleviates the symptoms of β-thalassaemia and sickle cell anaemia. NF-Y is considered to be the primary activating transcription factor which binds to globin CCAAT box motifs. Here we review recruitment of NF-Y to globin CCAAT boxes and the role NF-Y plays in regulating globin gene expression. This article is part of a Special Issue entitled: Nuclear Factor Y in Development and Disease, edited by Prof. Roberto Mantovani.

The Hierarchy of Transcriptional Activation: From Enhancer to Promoter

Regulatory elements (enhancers) that are remote from promoters play a critical role in the spatial, temporal, and physiological control of gene expression. Studies on specific loci, together with genome-wide approaches, suggest that there may be many common mechanisms involved in enhancer-promoter communication. Here, we discuss the multiprotein complexes that are recruited to enhancers and the hierarchy of events taking place between regulatory elements and promoters.

Effect of Interaction between Chromatin Loops on Cell-to-Cell Variability in Gene Expression

According to recent experimental evidence, the interaction between chromatin loops, which can be characterized by three factors-connection pattern, distance between regulatory elements, and communication form, play an important role in determining the level of cell-to-cell variability in gene expression. These quantitative experiments call for a corresponding modeling effect that addresses the question of how changes in these factors affect variability at the expression level in a systematic rather than case-by-case fashion. Here we make such an effort, based on a mechanic model that maps three fundamental patterns for two interacting DNA loops into a 4-state model of stochastic transcription. We first show that in contrast to side-by-side loops, nested loops enhance mRNA expression and reduce expression noise whereas alternating loops have just opposite effects. Then, we compare effects of facilitated tracking and direct looping on gene expression. We find that the former performs better than the latter in controlling mean expression and in tuning expression noise, but this control or tuning is distance-dependent, remarkable for moderate loop lengths, and there is a limit loop length such that the difference in effect between two communication forms almost disappears. Our analysis and results justify the facilitated chromatin-looping hypothesis.

Differential contribution of cis-regulatory elements to higher order chromatin structure and expression of the CFTR locus

Higher order chromatin structure establishes domains that organize the genome and coordinate gene expression. However, the molecular mechanisms controlling transcription of individual loci within a topological domain (TAD) are not fully understood. The cystic fibrosis transmembrane conductance regulator (CFTR) gene provides a paradigm for investigating these mechanisms.CFTR occupies a TAD bordered by CTCF/cohesin binding sites within which are cell-type-selective cis-regulatory elements for the locus. We showed previously that intronic and extragenic enhancers, when occupied by specific transcription factors, are recruited to the CFTR promoter by a looping mechanism to drive gene expression. Here we use a combination of CRISPR/Cas9 editing of cis-regulatory elements and siRNA-mediated depletion of architectural proteins to determine the relative contribution of structural elements and enhancers to the higher order structure and expression of the CFTR locus. We found the boundaries of the CFTRTAD are conserved among diverse cell types and are dependent on CTCF and cohesin complex. Removal of an upstream CTCF-binding insulator alters the interaction profile, but has little effect on CFTR expression. Within the TAD, intronic enhancers recruit cell-type selective transcription factors and deletion of a pivotal enhancer element dramatically decreases CFTR expression, but has minor effect on its 3D structure.

Three Dimensional Organization of Genome Might Have Guided the Dynamics of Gene Order Evolution in Eukaryotes

In eukaryotes, genes are nonrandomly organized into short gene-dense regions or "gene-clusters" interspersed by long gene-poor regions. How these gene-clusters have evolved is not entirely clear. Gene duplication may not account for all the gene-clusters since the genes in most of the clusters do not exhibit significant sequence similarity. In this study, using genome-wide data sets from budding yeast, fruit-fly, and human, we show that: 1) long-range evolutionary repositioning of genes strongly associate with their spatial proximity in the nucleus; 2) presence of evolutionary DNA break-points at involved loci hints at their susceptibility to undergo long-range genomic rearrangements; and 3) correlated epigenetic and transcriptional states of engaged genes highlight the underlying evolutionary constraints. The significance of observation 1, 2, and 3 are particularly stronger for the instances of inferred evolutionary gain, as compared with loss, of linear gene-clustering. These observations suggest that the long-range genomic rearrangements guided through 3D genome organization might have contributed to the evolution of gene order. We further hypothesize that the evolution of linear gene-clusters in eukaryotic genomes might have been mediated through spatial interactions among distant loci in order to optimize co-ordinated regulation of genes. We model this hypothesis through a heuristic model of gene-order evolution.

Predicting the three-dimensional folding of cis-regulatory regions in mammalian genomes using bioinformatic data and polymer models

The three-dimensional (3D) organization of chromosomes can be probed using methods like Capture-C. However, it is unclear how such population-level data relate to the organization within a single cell, and the mechanisms leading to the observed interactions are still largely obscure. We present a polymer modeling scheme based on the assumption that chromosome architecture is maintained by protein bridges, which form chromatin loops. To test the model, we perform FISH experiments and compare with Capture-C data. Starting merely from the locations of protein binding sites, our model accurately predicts the experimentally observed chromatin interactions, revealing a population of 3D conformations.

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.

α-Globin as a molecular target in the treatment of β-thalassemia

The thalassemias, together with sickle cell anemia and its variants, are the world's most common form of inherited anemia, and in economically undeveloped countries, they still account for tens of thousands of premature deaths every year. In developed countries, treatment of thalassemia is also still far from ideal, requiring lifelong transfusion or allogeneic bone marrow transplantation. Clinical and molecular genetic studies over the course of the last 50 years have demonstrated how coinheritance of modifier genes, which alter the balance of α-like and β-like globin gene expression, may transform severe, transfusion-dependent thalassemia into relatively mild forms of anemia. Most attention has been paid to pathways that increase γ-globin expression, and hence the production of fetal hemoglobin. Here we review the evidence that reduction of α-globin expression may provide an equally plausible approach to ameliorating clinically severe forms of β-thalassemia, and in particular, the very common subgroup of patients with hemoglobin E β-thalassemia that makes up approximately half of all patients born each year with severe β-thalassemia.

Evolution of α- and β-globin genes and their regulatory systems in light of the hypothesis of domain organization of the genome

The α- and β-globin gene domains are a traditional model for study of the domain organization of the eucaryotic genome because these genes encode hemoglobin, a physiologically important protein. The α-globin and β-globin gene domains are organized in completely different ways, while the expression of globin genes is tightly coordinated, which makes it extremely interesting to study the origin of these genes and the evolution of their regulatory systems. In this review, the organization of the α- and β-globin gene domains and their genomic environment in different taxonomic groups are comparatively analyzed. A new hypothesis of possible evolutionary pathways for segregated α- and β-globin gene domains of warm-blooded animals is proposed.

Uncovering enhancer functions using the α-globin locus

Over the last three decades, studies of the α- and β-globin genes clusters have led to elucidation of the general principles of mammalian gene regulation, such as RNA stability, termination of transcription, and, more importantly, the identification of remote regulatory elements. More recently, detailed studies of α-globin regulation, using both mouse and human loci, allowed the dissection of the sequential order in which transcription factors are recruited to the locus during lineage specification. These studies demonstrated the importance of the remote regulatory elements in the recruitment of RNA polymerase II (PolII) together with their role in the generation of intrachromosomal loops within the locus and the removal of polycomb complexes during differentiation. The multiple roles attributed to remote regulatory elements that have emerged from these studies will be discussed.

The architecture of gene expression: integrating dispersed cis-regulatory modules into coherent regulatory domains

Specificity and precision of expression are essential for the genes that regulate developmental processes. The specialized cis-acting modules, such as enhancers, that define gene expression patterns can be distributed across large regions, raising questions about the nature of the mechanisms that underline their action. Recent data has exposed the structural 3D context in which these long-range enhancers are operating. Here, we present how these studies shed new light on principles driving long-distance regulatory relationships. We discuss the molecular mechanisms that enable and accompany the action of long-range acting elements and the integration of multiple distributed regulatory inputs into the coherent and specific regulatory programs that are key to embryonic development.

Differential regulation of the α-globin locus by Krüppel-like Factor 3 in erythroid and non-erythroid cells

Background

Krüppel-like Factor 3 (KLF3) is a broadly expressed zinc-finger transcriptional repressor with diverse biological roles. During erythropoiesis, KLF3 acts as a feedback repressor of a set of genes that are activated by Krüppel-like Factor 1 (KLF1). Noting that KLF1 binds α-globin gene regulatory sequences during erythroid maturation, we sought to determine whether KLF3 also interacts with the α-globin locus to regulate transcription.

Results

We found that expression of a human transgenic α-globin reporter gene is markedly up-regulated in fetal and adult erythroid cells of Klf3^−/− mice. Inspection of the mouse and human α-globin promoters revealed a number of canonical KLF-binding sites, and indeed, KLF3 was shown to bind to these regions both in vitro and in vivo. Despite these observations, we did not detect an increase in endogenous murine α-globin expression in Klf3^−/− erythroid tissue. However, examination of murine embryonic fibroblasts lacking KLF3 revealed significant de-repression of α-globin gene expression. This suggests that KLF3 may contribute to the silencing of the α-globin locus in non-erythroid tissue. Moreover, ChIP-Seq analysis of murine fibroblasts demonstrated that across the locus, KLF3 does not occupy the promoter regions of the α-globin genes in these cells, but rather, binds to upstream, DNase hypersensitive regulatory regions.

Conclusions

These findings reveal that the occupancy profile of KLF3 at the α-globin locus differs in erythroid and non-erythroid cells. In erythroid cells, KLF3 primarily binds to the promoters of the adult α-globin genes, but appears dispensable for normal transcriptional regulation. In non-erythroid cells, KLF3 distinctly binds to the HS-12 and HS-26 elements and plays a non-redundant, albeit modest, role in the silencing of α-globin expression.

The clustering of CpG islands may constitute an important determinant of the 3D organization of interphase chromosomes

We used the 4C-Seq technique to characterize the genome-wide patterns of spatial contacts of several CpG islands located on chromosome 14 in cultured chicken lymphoid and erythroid cells. We observed a clear tendency for the spatial clustering of CpG islands present on the same and different chromosomes, regardless of the presence or absence of promoters within these CpG islands. Accordingly, we observed preferential spatial contacts between Sp1 binding motifs and other GC-rich genomic elements, including the DNA sequence motifs capable of forming G-quadruplexes. However, an anchor placed in a gene/CpG island-poor area formed spatial contacts with other gene/CpG island-poor areas on chromosome 14 and other chromosomes. These results corroborate the two-compartment model of the spatial organization of interphase chromosomes and suggest that the clustering of CpG islands constitutes an important determinant of the 3D organization of the eukaryotic genome in the cell nucleus. Using the ChIP-Seq technique, we mapped the genome-wide CTCF deposition sites in the chicken lymphoid and erythroid cells that were used for the 4C analysis. We observed a good correlation between the density of CTCF deposition sites and the level of 4C signals for the anchors located in CpG islands but not for an anchor located in a gene desert. It is thus possible that CTCF contributes to the clustering of CpG islands observed in our experiments.

Analysis of hundreds of cis-regulatory landscapes at high resolution in a single, high-throughput experiment

Gene expression during development and differentiation is regulated in a cell- and stage-specific manner by complex networks of intergenic and intragenic cis-regulatory elements whose numbers and representation in the genome far exceed those of structural genes. Using chromosome conformation capture, it is now possible to analyze in detail the interaction between enhancers, silencers, boundary elements and promoters at individual loci, but these techniques are not readily scalable. Here we present a high-throughput approach (Capture-C) to analyze cis interactions, interrogating hundreds of specific interactions at high resolution in a single experiment. We show how this approach will facilitate detailed, genome-wide analysis to elucidate the general principles by which cis-acting sequences control gene expression. In addition, we show how Capture-C will expedite identification of the target genes and functional effects of SNPs that are associated with complex diseases, which most frequently lie in intergenic cis-acting regulatory elements.

Chromatin signatures at transcriptional start sites separate two equally populated yet distinct classes of intergenic long noncoding RNAs

BACKGROUND

Mammalian transcriptomes contain thousands of long noncoding RNAs (lncRNAs). Some lncRNAs originate from intragenic enhancers which, when active, behave as alternative promoters producing transcripts that are processed using the canonical signals of their host gene. We have followed up this observation by analyzing intergenic lncRNAs to determine the extent to which they might also originate from intergenic enhancers.

RESULTS

We integrated high-resolution maps of transcriptional initiation and transcription to annotate a conservative set of intergenic lncRNAs expressed in mouse erythroblasts. We subclassified intergenic lncRNAs according to chromatin status at transcriptional initiation regions, defined by relative levels of histone H3K4 mono- and trimethylation. These transcripts are almost evenly divided between those arising from enhancer-associated (elncRNA) or promoter-associated (plncRNA) elements. These two classes of 5' capped and polyadenylated RNA transcripts are indistinguishable with regard to their length, number of exons or transcriptional orientation relative to their closest neighboring gene. Nevertheless, elncRNAs are more tissue-restricted, less highly expressed and less well conserved during evolution. Of considerable interest, we found that expression of elncRNAs, but not plncRNAs, is associated with enhanced expression of neighboring protein-coding genes during erythropoiesis.

CONCLUSIONS

We have determined globally the sites of initiation of intergenic lncRNAs in erythroid cells, allowing us to distinguish two similarly abundant classes of transcripts. Different correlations between the levels of elncRNAs, plncRNAs and expression of neighboring genes suggest that functional lncRNAs from the two classes may play contrasting roles in regulating the transcript abundance of local or distal loci.

Combinatorial effects of multiple enhancer variants in linkage disequilibrium dictate levels of gene expression to confer susceptibility to common traits

DNA variants (SNPs) that predispose to common traits often localize within noncoding regulatory elements such as enhancers. Moreover, loci identified by genome-wide association studies (GWAS) often contain multiple SNPs in linkage disequilibrium (LD), any of which may be causal. Thus, determining the effect of these multiple variant SNPs on target transcript levels has been a major challenge. Here, we provide evidence that for six common autoimmune disorders (rheumatoid arthritis, Crohn's disease, celiac disease, multiple sclerosis, lupus, and ulcerative colitis), the GWAS association arises from multiple polymorphisms in LD that map to clusters of enhancer elements active in the same cell type. This finding suggests a “multiple enhancer variant” hypothesis for common traits, where several variants in LD impact multiple enhancers and cooperatively affect gene expression. Using a novel method to delineate enhancer–gene interactions, we show that multiple enhancer variants within a given locus typically target the same gene. Using available data from HapMap and B lymphoblasts as a model system, we provide evidence at numerous loci that multiple enhancer variants cooperatively contribute to altered expression of their gene targets. The effects on target transcript levels tend to be modest and can be either gain- or loss-of-function. Additionally, the genes associated with multiple enhancer variants encode proteins that are often functionally related and enriched in common pathways. Overall, the multiple enhancer variant hypothesis offers a new paradigm by which noncoding variants can confer susceptibility to common traits.