Search

Determining the mechanisms by which genes are switched on and off during development and differentiation is a key aim of current biomedical research. Gene transcription has been widely observed to occur in a discontinuous fashion, with short bursts of activity interspersed with longer periods of inactivity. It is currently not known if or how this dynamic behaviour changes as mammalian cells differentiate. To investigate this, using a newly developed on-microscope analysis, we monitored mouse α-globin transcription in live cells throughout sequential stages of erythropoiesis. We find that changes in the overall levels of α-globin transcription are most closely associated with changes in the fraction of time a gene spends in the active transcriptional state. We identify differences in the patterns of transcriptional bursting throughout differentiation, with maximal transcriptional activity occurring in the mid-phase of differentiation. Early in differentiation, we observe increased fluctuation in the patterns of transcriptional activity whereas at the peak of gene expression, in early and intermediate erythroblasts, transcription appears to be relatively stable and efficient. Later during differentiation as α-globin expression declines, we again observed more variability in transcription within individual cells. We propose that the observed changes in transcriptional behaviour may reflect changes in the stability of enhancer-promoter interactions and the formation of active transcriptional compartments as gene expression is turned on and subsequently declines at sequential stages of differentiation.

Erythropoiesis requires a combination of ubiquitous and tissue-specific transcription factors (TFs). Here, through DNA affinity purification followed by mass spectrometry, we have identified the widely expressed protein MAZ (Myc-associated zinc finger) as a TF that binds to the promoter of the erythroid-specific human α-globin gene. Genome-wide mapping in primary human erythroid cells revealed that MAZ also occupies active promoters as well as GATA1-bound enhancer elements of key erythroid genes. Consistent with an important role during erythropoiesis, knockdown of MAZ reduces α-globin expression in K562 cells and impairs differentiation in primary human erythroid cells. Genetic variants in the MAZ locus are associated with changes in clinically important human erythroid traits. Taken together, these findings reveal the zinc-finger TF MAZ to be a previously unrecognized regulator of the erythroid differentiation program.

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., Globin loci harbor genes that are expressed embryonically and silenced postnatally. Here the authors show that zeta-globin silencing depends upon selective hypoacetylation of its TAD subdomain, which blocks its interaction with the alpha-globin super-enhancer, and zeta-globin can be reactivated by acetylation.

Mammalian genomes are subdivided into large (50-2000 kb) regions of chromatin referred to as Topologically Associating Domains (TADs or sub-TADs). Chromatin within an individual TAD contacts itself more frequently than with regions in surrounding TADs thereby directing enhancer-promoter interactions. In many cases, the borders of TADs are defined by convergently orientated boundary elements associated with CCCTC-binding factor (CTCF), which stabilises the cohesin complex on chromatin and prevents its translocation. This delimits chromatin loop extrusion which is thought to underlie the formation of TADs. However, not all CTCF-bound sites act as boundaries and, importantly, not all TADs are flanked by convergent CTCF sites. Here, we examined the CTCF binding sites within a ∼70 kb sub-TAD containing the duplicated mouse α-like globin genes and their five enhancers (5’-R1-R2-R3-Rm-R4-α1-α2-3’). The 5’ border of this sub-TAD is defined by a pair of CTCF sites. Surprisingly, we show that deletion of the CTCF binding sites within and downstream of the α-globin locus leaves the sub-TAD largely intact. The predominant 3’ border of the sub-TAD is defined by a steep reduction in contacts: this corresponds to the transcribed α2-globin gene rather than the CTCF sites at the 3’-end of the sub-TAD. Of interest, the almost identical α1- and α2-globin genes interact differently with the enhancers, resulting in preferential expression of the proximal α1-globin gene which behaves as a partial boundary between the enhancers and the distal α2-globin gene. Together, these observations provide direct evidence that actively transcribed genes can behave as boundary elements.Significance StatementMammalian genomes are complex, organised 3D structures, partitioned into Topologically Associating Domains (TADs): chromatin regions that preferentially self-interact. These chromatin interactions are thought to be driven by a mechanism that continuously extrudes chromatin loops, forming structures delimited by chromatin boundary elements and reflecting the activity of enhancers and promoters. Boundary elements bind architectural proteins such as CCCTC-binding factor (CTCF). Previously, an overlap between the functional roles of enhancers and promoters has been shown. However, whether there is overlap between enhancers/promoters and boundary elements is not known. Here, we show that actively transcribed genes can also behave as boundary elements, similar to CTCF boundaries. In both cases, multi-protein complexes bound to these regions may stall the process of chromatin loop extrusion.

Erythropoiesis requires a combination of ubiquitous and tissue-specific transcription factors. Here, through DNA affinity purification followed by mass spectrometry, we have identified the widely expressed protein MAZ (Myc-associated zinc finger) as a transcription factor that binds to the promoter of the erythroid-specific human α-globin gene. Genome-wide mapping in primary human erythroid cells revealed that MAZ also occupies active promoters as well as GATA1-bound enhancer elements of key erythroid genes. Consistent with an important role during erythropoiesis, knockdown of MAZ in primary human erythroid cells impairs erythroid differentiation, and genetic variants in the MAZ locus are associated with clinically important human erythroid traits. Taken together, these findings reveal the Zinc-finger transcription factor MAZ to be a previously unrecognised regulator of the erythroid differentiation program.

A substantial amount of organismal complexity is thought to be encoded by enhancers which specify the location, timing, and levels of gene expression. In mammals there are more enhancers than promoters which are distributed both between and within genes. Here we show that activated, intragenic enhancers frequently act as alternative tissue-specific promoters producing a class of abundant, spliced, multiexonic poly(A) + RNAs (meRNAs) which reflect the host gene's structure. meRNAs make a substantial and unanticipated contribution to the complexity of the transcriptome, appearing as alternative isoforms of the host gene. The low protein-coding potential of meRNAs suggests that many meRNAs may be byproducts of enhancer activation or underlie as-yet-unidentified RNA-encoded functions. Distinguishing between meRNAs and mRNAs will transform our interpretation of dynamic changes in transcription both at the level of individual genes and of the genome as a whole. © 2012 Elsevier Inc..

To understand how mammalian genes are regulated from their natural chromosomal environment, we have analysed the molecular events occurring throughout a 150 kb chromatin segment containing the alpha globin gene locus as it changes from a poised, silent state in erythroid progenitors, to the fully activated state in late, erythroid cells. Active transcription requires the late recruitment of general transcription factors, mediator and Pol II not only to the promoter but also to its remote regulatory elements. Natural mutants of the alpha cluster show that whereas recruitment of the pre-initiation complex to the upstream elements occurs independently, recruitment to the promoter is largely dependent on the regulatory elements. An improved, quantitative chromosome conformation capture analysis demonstrates that this recruitment is associated with a conformational change, in vivo, apposing the promoter with its remote regulators, consistent with a chromosome looping mechanism. These findings point to a general mechanism by which a gene can be held in a poised state until the appropriate stage for expression, coordinating the level and timing of gene expression during terminal differentiation.

© 2016 Nature America, Inc. All rights reserved. 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.

We have characterized a newly identified 16.6 kb deletion which removes a significant proportion of the human alpha-globin cluster including the psizeta1, alpha(D), psialpha1 and alpha2-globin genes but leaves the duplicated alpha1 gene intact. This complicated rearrangement results from a combination of slippage and strand switching at sites of microhomology during replication. Functional analysis shows that expression of the remaining alpha1 gene is increased, rather than down-regulated by this deletion. This could be related to its proximity to the remote upstream alpha-globin regulatory elements or reduced competition for these elements in the absence of the dominant alpha2-globin gene. The finding of a very mild phenotype associated with such an extensive deletion in the alpha-globin cluster implies that much of the DNA removed by the deletion is likely to be functionally unimportant. These findings suggest that other than the upstream regulatory elements and promoter proximal elements there are unlikely to be additional positive cis-acting sequences in the alpha-globin cluster.

To address the mechanism by which the human globin genes are activated during erythropoiesis, we have used a tiled microarray to analyze the pattern of transcription factor binding and associated histone modifications across the telomeric region of human chromosome 16 in primary erythroid and nonerythroid cells. This 220-kb region includes the α globin genes and 9 widely expressed genes flanking the α globin locus. This un-biased, comprehensive analysis of transcription factor binding and histone modifications (acetylation and methylation) described here not only identified all known cis-acting regulatory elements in the human α globin cluster but also demonstrated that there are no additional erythroid-specific regulatory elements in the 220-kb region tested. In addition, the pattern of histone modification distinguished promoter elements from potential enhancer elements across this region. Finally, comparison of the human and mouse orthologous regions in a unique mouse model, with both regions coexpressed in the same animal, showed significant differences that may explain how these 2 clusters are regulated differently in vivo.

Natural deletions of the region upstream of the human alpha-globin gene cluster, together with expression studies in cell lines and transgenic mice, identified a single element (HS -40) as necessary and perhaps sufficient for high-level expression of the alpha-globin genes. A similar element occupies the corresponding position upstream of the mouse (m) alpha-globin genes (mHS -26) and was thought to have similar functional properties. We knocked out mHS -26 by homologous recombination and observed the surprising result that instead of the expected severe alpha-thalassemia phenotype, the mice had a mild disease. Transcription levels of the mouse genes were reduced by about 50%, but homozygotes were healthy, with normal hemoglobin levels and only mild decreases in mean corpuscular volume and mean corpuscular hemoglobin. These results may indicate differences in the regulation of the alpha-globin clusters in mice and humans or that additional cis-acting elements remain to be characterized in one or both clusters.

To date, all of the chromosomal deletions that cause -thalassemia remove the structural  genes and/or their regulatory element (HS –40). A unique deletion occurs in a single family that juxtaposes a region that normally lies approximately 18-kilobase downstream of the human  cluster, next to a structurally normal -globin gene, and silences its expression. During development, the CpG island associated with the -globin promoter in the rearranged chromosome becomes densely methylated and insensitive to endonucleases, demonstrating that the normal chromatin structure around the -globin gene is perturbed by this mutation and that the gene is inactivated by a negative chromosomal position effect. These findings highlight the importance of the chromosomal environment in regulating globin gene expression.

We have examined the pattern of human globin gene switching in transgenic mice containing three different γ and β gene constructs (HS2GγAγδβ, HS2Aγβneo, and HS2Aγenβ) and compared the results with previously described transgenics (HS2Aγβ, HS2GγAγ-117δβ, and LCRεGγAγδβ). Developmental regulation was observed in all cases with identical patterns in lines bearing the same construct. Three different patterns of switching were observed: LCRεGγAγδβ and HS2Aγβneo mice switched rapidly, HS2GγAγδβ and HS2GγAγ-117δβ at an intermediate rate, and HS2Aγβ and HS2Aγenβ mice showed delayed switching, with a plateau in late fetal-early neonatal life and readily detectable levels of γ mRNA in adults. No difference was observed in the time of switching of the HS2GγAγδβ mice compared with those with the Aγ-117 hereditary persistence of fetal hemoglobin mutation, but adult levels of γ mRNA were significantly higher (≈5%) in lines carrying the mutation than in those without (≈1%). Reversion to the rapid switch of the LCRεGγAγδβ mice was observed in three lines with the HS2Aγβ neo construct in which expression of the tk-neo gene was approximately equal to that of the globin genes. The inclusion of the Aγ enhancer in HS2Aγβ mice did not alter the pattern of switching, or reduce the relatively high levels of γ mRNA in these lines. However, unlike other HS2 mice, the combination of HS2 and the Aγ enhancer resulted in copy number-dependent expression in HS2Aγenβ lines, with intrauterine death at ≈12.5 days gestation at high copy numbers. These results demonstrate that numerous elements throughout the β globin gene cluster interact to produce the correct pattern of developmental regulation of these genes. Furthermore, extinction of γ gene expression in adult life is not completely autonomous and is incomplete when HS2 is the only LCR element present.

Although mutations causing monogenic disorders most frequently lie within the affected gene, sequence variation in complex disorders is more commonly found in noncoding regions. Furthermore, recent genome- wide studies have shown that common DNA sequence variants in noncoding regions are associated with "normal" variation in gene expression resulting in cell-specific and/or allele-specific differences. The mechanism by which such sequence variation causes changes in gene expression is largely unknown. We have addressed this by studying natural variation in the binding of key transcription factors (TFs) in the well-defined, purified cell system of erythropoiesis. We have shown that common polymorphisms frequently directly perturb the binding sites of key TFs, and detailed analysis shows how this causes considerable (~10-fold) changes in expression from a single allele in a tissue-specific manner. We also show how a SNP, located at some distance from the recognized TF binding site, may affect the recruitment of a large multiprotein complex and alter the associated chromatin modification of the variant regulatory element. This study illustrates the principles by which common sequence variation may cause changes in tissue-specific gene expression, and suggests that such variation may underlie an individual's propensity to develop complex human genetic diseases.

We have combined the circular chromosome conformation capture protocol with high-throughput, genome-wide sequence analysis to characterize the cis -acting regulatory network at a single locus. In contrast to methods which identify large interacting regions (10–1000 kb), the 4C approach provides a comprehensive, high-resolution analysis of a specific locus with the aim of defining, in detail, the cis -regulatory elements controlling a single gene or gene cluster. Using the human α-globin locus as a model, we detected all known local and long-range interactions with this gene cluster. In addition, we identified two interactions with genes located 300 kb (NME4) and 625 kb (FAM173a) from the α-globin cluster.

Somatic cell hybrids generated from transgenic mouse cells have been used to examine the developmental regulation of human gamma-to-beta-globin gene switching. In hybrids between mouse erythroleukemia (MEL) cells and transgenic erythroblasts taken at various stages of development, there was regulated expression of the human fetal gamma and adult beta genes, reproducing the in vivo pattern prior to fusion. Hybrids formed from embryonic blood cells produced predominantly gamma mRNA, whereas beta gene expression was observed in adult hybrids and a complete range of intermediate patterns was found in fetal liver hybrids. The adult environment of the MEL cells, therefore, did not appear to influence selective transcription from this gene complex. Irradiation of the embryonic erythroid cells prior to fusion resulted in hybrids containing only small fragments of donor chromosomes, but the pattern of gene expression did not differ from that of unirradiated hybrids. This finding suggests that continued expression of trans-acting factors from the donor erythroblasts is not necessary for continued expression of the human gamma gene in MEL cells. These results contrast with the lack of developmental regulation of the cluster after transfection of naked DNA into MEL cells and suggest that epigenetic processes established during normal development result in the gene cluster adopting a developmental stage-specific, stable conformation which is maintained through multiple rounds of replication and transcription in the MEL cell hybrids. On prolonged culture, hybrids that initially expressed only the gamma transgene switched to beta gene expression. The time period of switching, from approximately 10 to > 40 weeks, was similar to that seen previously in human fetal erythroblast x MEL cell hybrids but in this case bore no relationship to the time of in vivo switching. It seems unlikely, therefore, that switching in these hybrids is regulated by a developmental clock.

C16orf35 is a conserved and widely expressed gene lying adjacent to the human α-globin cluster in all vertebrate species. In-depth sequence analysis shows that C16orf35 (now called NPRL3) is an orthologue of the yeast gene Npr3 (nitrogen permease regulator 3) and, furthermore, is a paralogue of its protein partner Npr2. The yeast Npr2/3 dimeric protein complex senses amino acid starvation and appropriately adjusts cell metabolism via the TOR pathway. Here we have analysed a mouse model in which expression of Nprl3 has been abolished using homologous recombination. The predominant effect on RNA expression appears to involve genes that regulate protein synthesis and cell cycle, consistent with perturbation of the mTOR pathway. Embryos homozygous for this mutation die towards the end of gestation with a range of cardiovascular defects, including outflow tract abnormalities and ventriculoseptal defects consistent with previous observations, showing that perturbation of the mTOR pathway may affect development of the myocardium. NPRL3 is a candidate gene for harbouring mutations in individuals with developmental abnormalities of the cardiovascular system. © 2012 Springer Science+Business Media, LLC.

Previous studies in the mouse have shown that high levels of α-globin gene expression in late erythropoiesis depend on long-range, physical interactions between remote upstream regulatory elements and the globin promoters. Using quantitative chromosome conformation capture (q3C), we have now analyzed all interactions between 4 such elements lying 10 to 50 kb upstream of the human α cluster and their interactions with the α-globin promoter. All of these elements interact with the α-globin gene in an erythroid-specific manner. These results were confirmed in a mouse model of human α globin expression in which the human cluster replaces the mouse cluster in situ (humanized mouse). We have also shown that expression and all of the long-range interactions depend largely on just one of these elements; removal of the previously characterized major regulatory element (called HS −40) results in loss of all the interactions and α-globin expression. Reinsertion of this element at an ectopic location restores both expression and the intralocus interactions. In contrast to other more complex systems involving multiple upstream elements and promoters, analysis of the human α-globin cluster during erythropoiesis provides a simple and tractable model to understand the mechanisms underlying long-range gene regulation.

Abstract 3381 It is well established that the level of gene expression can vary significantly between normal individuals, and that the majority of this variation is due to naturally occurring genomic variability caused by single nucleotide polymorphisms (SNPs). Therefore, identifying functional cis-regulatory polymorphisms and understanding how they influence gene expression is an important new task in many areas of medical research, including molecular hematology. We have previously shown that an entirely new form of alpha-thalassemia is caused by a gain of function regulatory SNP in an unremarkable non-coding region in the alpha-globin cluster. This SNP creates a novel, functional GATA site, which recruits a tissue-specific transcription factor (TF) complex. This creates a new promoter-like element, which interferes with activation of the globin genes (De Gobbi et al. Science 2006,312:1215–1257). Here, to investigate the extent and the impact of this class of regulatory SNP, using ChIP-Seq we characterized differences in the occupancy of Scl/TAL-1 (a tissue-specific TF critical for erythroid maturation) in the erythroblasts of two individuals from the same ethnic background (Caucasian 1, C1, and Caucasian 2, C2). Sequence reads from two biological replicates of each individual were merged and aligned to the human reference genome (NCBI36/hg18) and a total of 2936 Scl/TAL-1 bound regions were identified. Using two de novo motif finding algorithms (MEME and DREME), we identified GATA (WGATAR) and E-box (CAGMTG) sites as the preferred sequences associated with in vivo binding of Scl/TAL-1. In addition, other motifs were enriched at the Scl/TAL-1 targets; among these were binding sites for known TFs (Sp1/Klf, RUNX1 and NFE2). To identify differentially bound regions between C1 and C2, a two-class paired-test, Rank Product analysis (500 permutations, FDR Since it has been previously shown that the function of active transcriptional elements can be predicted on the basis of chromatin signatures (e.g. enhancers marked by H3K4me1 and promoters marked by H3K4me3), to further characterize the Scl/TAL-1 differentially occupied sites, we asked which chromatin signatures are associated with these regions. H3K4me1 and H3K4me3 ChIP-Seq experiments, together with analyses of publicly available data sets, showed that the most of the SNPs responsible for variation in the recruitment of Scl/TAL-1 (23/25) lie in DNA sequences that have chromatin signatures predictive of enhancer elements, suggesting a potential long-range function in modulating gene expression. Finally, Scl/TAL-1 ChIP-Seq analysis of erythroblasts of a third individual from a different ethnic background (African-Caribbean, A3) revealed more distinctive targets including a well known regulatory SNP at the promoter of the DARC gene (encoding the Duffy blood group), which alters a GATA binding motif in A3, conferring the malaria-resistant Duffy-null phenotype. Given the exponential growth in genome-wide association studies by which numerous SNPs are being either associated with hematological parameters, or implicated in the etiology of hematologic disorders, this study elucidates molecular mechanisms which might account for phenotypic diversity and highlights the importance of carrying out functional characterization of non-coding polymorphisms found to be associated with disease risk. Disclosures: No relevant conflicts of interest to declare.