Search

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. [ABSTRACT FROM AUTHOR]

Micronutrient selenium (Se) plays a key role in redox regulation through its incorporation into selenoproteins as the 21st amino acid selenocysteine (Sec). Because Se deficiency appears to be a cofactor in the anemia associated with chronic inflammatory diseases, we reasoned that selenoproteins may contribute to erythropoietic recovery from anemia, referred to as stress erythropoiesis. Here, we report that loss of selenoproteins through Se deficiency or by mutation of the Sec tRNA (tRNA[Sec]) gene (Trsp) severely impairs stress erythropoiesis at 2 stages. Early stress erythroid progenitors failed to expand and properly differentiate into burst-forming unit-erythroid cells, whereas late-stage erythroid progenitors exhibited a maturation defect that affected the transition of proerythroblasts to basophilic erythroblasts. These defects were, in part, a result of the loss of selenoprotein W (SelenoW), whose expression was reduced at both transcript and protein levels in Se-deficient erythroblasts. Mutation of SelenoW in the bone marrow cells significantly decreased the expansion of stress burst-forming unit-erythroid cell colonies, which recapitulated the phenotypes induced by Se deficiency or mutation of Trsp. Similarly, mutation of SelenoW in murine erythroblast (G1E) cell line led to defects in terminal differentiation. In addition to the erythroid defects, the spleens of Se-deficient mice contained fewer red pulp macrophages and exhibited impaired development of erythroblastic island macrophages, which make up the niche supporting erythroblast development. Taken together, these data reveal a critical role of selenoproteins in the expansion and development of stress erythroid progenitors, as well as the erythroid niche during acute anemia recovery. [ABSTRACT FROM AUTHOR]

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

The variation in local rates of mutations can affect both the evolution of genes and their function in normal and cancer cells. Deciphering the molecular determinants of this variation will be aided by the elucidation of distinct types of mutations, as they differ in regional preferences and in associations with genomic features. Chromatin organization contributes to regional variation in mutation rates, but its contribution differs among mutation types. In both germline and somatic mutations, base substitutions are more abundant in regions of closed chromatin, perhaps reflecting error accumulation late in replication. By contrast, a distinctive mutational state with very high levels of insertions and deletions (indels) and substitutions is enriched in regions of open chromatin. These associations indicate an intricate interplay between the nucleotide sequence of DNA and its dynamic packaging into chromatin, and have important implications for current biomedical research. This Review focuses on recent studies showing associations between chromatin state and mutation rates, including pairwise and multivariate investigations of germline and somatic (particularly cancer) mutations. [ABSTRACT FROM AUTHOR]

R5637 --> 661.7 --> Thousands of genome‐wide epigenetic (epigenomic) datasets determined in the past decade have generated maps of chromatin accessibility, histone modifications, and binding by transcription factors and structural proteins in a large number of cell types in humans and model species. Such maps have revealed many informative associations between groups of epigenetic features and different classes of cis‐regulatory elements (CREs), and they enable a prediction of the locations of candidate CREs (cCREs) in many genomes. The actuation (opening of chromatin) and epigenetic state (combination of features such as histone modifications) acquired at specific cCREs has been associated with specific gene activation. However, despite the large number and rich variety of available epigenomic datasets, it is difficult for researchers to effectively utilize all the data relevant to their projects. Furthermore, comprehensive rules and catalogs have not been established for which CREs regulate specific target genes and how. Integrative analysis can help meet these needs, and our VISION project provides a ValIdated Systematic IntegratiONof epigenomic data in hematopoiesis in human and mouse blood cells. After careful curation and normalization of epigenomic data from many sources, we systematically integrated these data to generate a comprehensive view of the regulatory landscape (epigenetic states) and predicted cCREs in differentiating hematopoietic cell types jointly in both mouse and human. Multiple machine learning approaches reveal target genes predicted to be regulated by each element in different cell types. These approaches include predicting the regulatory impact of each element on targets by regression models, finding linkages between distinctive clusters of elements and potential target genes from self‐organizing maps, and visual overlays of chromatin interaction maps. Directed mutation and expression analysis around selected genes supports the utility of the assignments. These catalogs of highly annotated regulatory elements are useful for many applications, including interpretation of the phenotypic impact of genetic variants in noncoding genomic regions. Resources from the VISION project are available at our website usevision.org. [ABSTRACT FROM AUTHOR]

Complex traits such as susceptibility to diseases are determined in part by variants at multiple genetic loci. Genome-wide association studies can identify these loci, but most phenotypeassociated variants lie distal to protein-coding regions and are likely involved in regulating gene expression. Understanding how these genetic variants affect complex traits depends on the ability to predict and test the function of the genomic elements harboring them. Community efforts such as the ENCODE Project provide a wealth of data about epigenetic features associated with gene regulation. These data enable the prediction of testable functions for many phenotype-associated variants. [ABSTRACT FROM AUTHOR]

Differential gene expression is the fundamental mechanism underlying animal development and cell differentiation. However, it is a challenge to identify comprehensively and accurately the DNA sequences that are required to regulate gene expression: namely, cis-regulatory modules (CRMs). Three major features, either singly or in combination, are used to predict CRMs: clusters of transcription factor binding site motifs, non-coding DNA that is under evolutionary constraint and biochemical marks associated with CRMs, such as histone modifications and protein occupancy. The validation rates for predictions indicate that identifying diagnostic biochemical marks is the most reliable method, and understanding is enhanced by the analysis of motifs and conservation patterns within those predicted CRMs. [ABSTRACT FROM AUTHOR]

The article presents comment on the study that analyzes epigenetic data such as transcription factor occupancy, histone modifications, and chromatin interaction frequencies in a cell line model of hematopoietic stem/progenitor cells (HSPCs). Details of the study that consider an analogy of gene expression being a song or symphony played by musicians are mentioned. The result of the study indicates that the model offers basis for consolidating information about regulatory complexes is offered.

We have sequenced and fully annotated a 65,871-bp region of mouse Chromosome 17 including the Hba-ps4 α-globin pseudogene. Comparative sequence analysis with the functional α-globin loci at human Chromosome 16p13.3 and mouse Chromosome 11 shows that this segment of mouse Chromosome 17 contains a group of three α-like pseudogenes (Hba-psm-Hba-ps4-Hba-q3), similar to the duplicated sets found at the functional mouse cluster on Chromosome 11. In addition, exons 7 to 12 of the mLuc7L gene are present just downstream from the pseudogene cluster, indicating that this clone contains the region in which human 16p13.3 switches in synteny between mouse Chromosomes 11 and 17. Comparison of the sequences around the α-like clusters on the two mouse chromosomes reveals the presence of conserved tandem repeats. We propose that these repetitive elements have played a role in the fragmentation of the mouse α cluster during evolution. [ABSTRACT FROM AUTHOR]

Comparing the genomes of two different species allow the exploration of a host of intriguing evolutionary and genetic questions. [ABSTRACT FROM AUTHOR]

Deals with the importance of comparative analyses of genome sequences to the improvement of the health of individuals and society. Factor considered to be the ultimate genetic map; Limitation of the knowledge on the DNA sequence; Major principles of comparative genomics.

Presents computational approaches in analyzing interspecies genomic sequence alignments. Way of distinguishing functional noncoding sequences from neutral DNA; Identification of putative regulatory regions throughout the human genome; Characteristics of the mammalian genes.

Presents a study on the evolutionary change in the human genome. Alignment of the human and mouse genome sequences; Number of protein families; Characteristics of the genes.

Building on the pioneering efforts of Professor Huisman, several different databases of hemoglobin variants have been developed, each with progressively increased capacity for sophisticated queries and prompt updating. These resources are reviewed in the context of a larger plan for providing related resources on hemoglobins, benign and pathological variation in these proteins and the genes that encode them, and the regulation of the globin genes. [ABSTRACT FROM AUTHOR]

Provides information on an analysis of hemoglobins, focusing on the evolution of different patterns of gene expression. In-depth look at plant and animal hemoglobin; Illustration of the diversity of proposed functions and regulation of hemoglobins; Mechanisms of regulating hemoglobin and gene expression; Results of the analysis.

The article focuses on the implications of the genomewide association studies (GWAS) on the genome therapy for the treatment of disorders of hemoglobins. It mentions about the characterization of common single-nucleotide polymorphisms (SNPs) in the globin genes of monogenic diseases of hemoglobin. It discusses the variation in fetal hemoglobin (HbF) in adult which attenuates the symptoms with sickle cell disease or thalassemias.

Based on a massive transcription factor location analysis within a single cell type, in this issue Yan et al. find that the great majority of occupancies occur within dense clusters of up to 100 factors that almost invariably contain cohesins. Retention of cohesins at cluster sites during mitosis raises the possibility that they contribute to transcriptional memory during the cell cycle. [ABSTRACT FROM AUTHOR]

Presents views of evolution and regulation of vertebrate beta-like globin gene clusters from an orphaned gene in marsupials. Model for the independent evolution of beta-like globin gene clusters in birds and mammals; Physiological roles of hemoglobins; Alteration in the view of hemoglobin gene evolution in the warm-blooded vertebrates.

The article offers information on the methods used to identify the cis-regulatory sequences by looking at their function and the determinants of organismic development including tissue-specificity and level of gene expression. Topics discussed include functional identification of regulatory elements within active chromatin (FIREWACh), site-specific integration fluorescence–activated cell sorting followed by sequencing (SIF-seq), and chromatin accessibility and histone modifications.

What is the best way to identify regulatory DNA sequences such as enhancers, promoters, insulators and silencers? A new study shows that specific binding by a coactivator protein identifies enhancers that are invisible to common detection methods based on evolutionary constraint. [ABSTRACT FROM AUTHOR]

Background: Epigenetic modification of chromatin plays a pivotal role in regulating gene expression during cell differentiation. The scale and complexity of epigenetic data pose significant challenges for biologists to identify the regulatory events controlling cell differentiation. Results: To reduce the complexity, we developed a package, called Snapshot, for clustering and visualizing candidate cis-regulatory elements (cCREs) based on their epigenetic signals during cell differentiation. This package first introduces a binarized indexing strategy for clustering the cCREs. It then provides a series of easily interpretable figures for visualizing the signal and epigenetic state patterns of the cCREs clusters during the cell differentiation. It can also use different hierarchies of cell types to highlight the epigenetic history specific to any particular cell lineage. We demonstrate the utility of Snapshot using data from a consortium project for ValIdated Systematic IntegratiON (VISION) of epigenomic data in hematopoiesis. Conclusion: The package Snapshot can identify all distinct clusters of genomic locations with unique epigenetic signal patterns during cell differentiation. It outperforms other methods in terms of interpreting and reproducing the identified cCREs clusters. The package of Snapshot is available at GitHub: https://github.com/guanjue/Snapshot. [ABSTRACT FROM AUTHOR]

Many evolutionary studies over the past decade have estimated αsel, the proportion of all nucleotides in the human genome that are subject to purifying selection because of their biological function. Most of these studies have estimated the nucleotide substitution rates from genome sequence alignments across many diverse mammals. Some αsel estimates will be affected by the heterogeneity of substitution rates in neutral sequence across the genome. Most will also be inaccurate if change in the functional sequence repertoire occurs rapidly relative to the separation of lineages that are being compared. Evidence gathered from both evolutionary and experimental analyses now indicate that rates of "turnover" of functional, predominantly noncoding, sequence are, indeed, high. They are sufficiently high that an estimated 50% of mouse constrained noncoding sequence is predicted not to be shared with rat, a closely related rodent. The rapidity of turnover results in, at least, a twofold underestimate of αsel by analyses that measure constraint across the eutherian phylogeny. Approaches that take account of turnover estimate that the steady-state value of αsel lies between 10% and 15%. Experimental studies corroborate the predicted rates of loss and gain of noncoding functional sites. These studies show the limitations inherent in the use of deep sequence conservation for identifying functional sequence. Experimental investigations focusing on lineage-specific, noncoding, and functional sequence are now essential if we are to appreciate the complete functional repertoire of the human genome. [ABSTRACT FROM AUTHOR]

Accessing and analyzing the exponentially expanding genomic sequence and functional data pose a challenge for biomedical researchers. Here we describe an interactive system, Galaxy, that combines the power of existing genome annotation databases with a simple Web portal to enable users to search remote resources, combine data from independent queries, and visualize the results. The heart of Galaxy is a flexible history system that stores the queries from each user; performs operations such as intersections, unions, and subtractions; and links to other computational tools. Galaxy can be accessed at http://g2.bx.psu.edu. [ABSTRACT FROM AUTHOR]

Describes the use of comparative sequence analysis in identifying functional coding and noncoding elements conserved throughout evolution in the U.S. Presence of conserved transcription factor-binding sites; Adjustment of the threshold levels of conservation; Alignment of the microbial genomes.

Chromatin immunoprecipitation followed by massively parallel, high throughput sequencing (ChIP-seq) is the method of choice for genomewide identification of DNA segments bound by specific transcription factors or in chromatin with particular histone modifications. However, the quality of ChIP-seq datasets varies widely, with a substantial fraction being of intermediate to poor quality. Thus, it is important to discern and control the factors that contribute to variation in ChIP-seq. In this study, we focused on sonication, a user-controlled variable, to produce sheared chromatin. We systematically varied the amount of shearing of fixed chromatin from a mouse erythroid cell line, carefully measuring the distribution of resultant fragment lengths prior to ChIP-seq. This systematic study was complemented with a retrospective analysis of additional experiments. We found that the level of sonication had a pronounced impact on the quality of ChIP-seq signals. Over-sonication consistently reduced quality, while the impact of under-sonication differed among transcription factors, with no impact on sites bound by CTCF but frequently leading to the loss of sites occupied by TAL1 or bound by POL2. The bound sites not observed in low-quality datasets were inferred to be a mix of both direct and indirect binding. We leveraged these findings to produce a set of CTCF ChIP-seq datasets in rare, primary hematopoietic progenitor cells. Our observation that the amount of chromatin sonication is a key variable in success of ChIP-seq experiments indicates that monitoring the level of sonication can improve ChIP-seq quality and reproducibility and facilitate ChIP-seq in rare cell types. [ABSTRACT FROM AUTHOR]

Summary Epigenetic modifications reflect key aspects of transcriptional regulation, and many epigenomic datasets have been generated under different biological contexts to provide insights into regulatory processes. However, the technical noise in epigenomic datasets and the many dimensions (features) examined make it challenging to effectively extract biologically meaningful inferences from these datasets. We developed a package that reduces noise while normalizing the epigenomic data by a novel normalization method, followed by integrative dimensional reduction by learning and assigning epigenetic states. This package, called S3V2-IDEAS, can be used to identify epigenetic states for multiple features, or identify discretized signal intensity levels and a master peak list across different cell types for a single feature. We illustrate the outputs and performance of S3V2-IDEAS using 137 epigenomics datasets from the VISION project that provides V al I dated S ystematic I ntegrati ON of epigenomic data in hematopoiesis. Availability and implementation S3V2-IDEAS pipeline is freely available as open source software released under an MIT license at: https://github.com/guanjue/S3V2%5fIDEAS%5fESMP. Supplementary information Supplementary data are available at Bioinformatics online. [ABSTRACT FROM AUTHOR]

Increasing fetal hemoglobin (HbF) levels in adult red blood cells provides clinical benefit to patients with sickle cell disease and some forms of b-thalassemia. To identify potentially druggable HbF regulators in adult human erythroid cells, we employed a protein kinase domain–focused CRISPR-Cas9–based genetic screen with a newly optimized single-guide RNA scaffold. The screen uncovered the heme-regulated inhibitor HRI (also known as EIF2AK1), an erythroid-specific kinase that controls protein translation, as an HbF repressor. HRI depletion markedly increased HbF production in a specific manner and reduced sickling in cultured erythroid cells. Diminished expression of the HbF repressor BCL11A accounted in large part for the effects of HRI depletion. Taken together, these results suggest HRI as a potential therapeutic target for hemoglobinopathies. [ABSTRACT FROM AUTHOR]

Summary: Hidden Markov models (HMMs) are flexible and widely used in scientific studies. Particularly in genomics and genetics, there are multiple distinct regimes in the genome within each of which the relationships among multivariate features are distinct. Examples include differential gene regulation depending on gene functions and experimental conditions, and varying combinatorial patterns of multiple transcription factors. We developed a software package called MRHMMs (Multivariate Regression Hidden Markov Models and the variantS) that accommodates a variety of HMMs that can be flexibly applied to many biological studies and beyond. MRHMMs supplements existing HMM software packages in two aspects. First, MRHMMs provides a diverse set of emission probability structures, including mixture of multivariate normal distributions and (logistic) regression models. Second, MRHMMs is computationally efficient for analyzing large data-sets generated in current genome-wide studies. Especially, the software is written in C for the speed advantage and further amenable to implement alternative models to meet users’ own purposes.Availability and implementation: http://sourceforge.net/projects/mrhmms/Contact: ghoshd@psu.edu or yuzhang@stat.psu.eduSupplementary information: Supplementary data are available at Bioinformatics online. [ABSTRACT FROM PUBLISHER]

Transcriptional enhancers have been extensively characterized, but cis -regulatory elements involved in acute gene repression have received less attention. Transcription factor GATA1 promotes erythroid differentiation by activating and repressing distinct gene sets. Here, we study the mechanism by which GATA1 silences the proliferative gene Kit during murine erythroid cell maturation and define stages from initial loss of activation to heterochromatinization. We find that GATA1 inactivates a potent upstream enhancer but concomitantly creates a discrete intronic regulatory region marked by H3K27ac, short noncoding RNAs, and de novo chromatin looping. This enhancer-like element forms transiently and serves to delay Kit silencing. The element is ultimately erased via the FOG1/NuRD deacetylase complex, as revealed by the study of a disease-associated GATA1 variant. Hence, regulatory sites can be self-limiting by dynamic co-factor usage. Genome-wide analyses across cell types and species uncover transiently active elements at numerous genes during repression, suggesting that modulation of silencing kinetics is widespread. [Display omitted] • Time-resolved analyses reveal transient enhancer elements during gene silencing • A GATA1-driven transient enhancer delays GATA1-mediated silencing of the Kit gene • FOG1/NuRD is required for the self-limitation of GATA1-occupied transient elements • Transient, potential silencing delay elements seem to be widespread Transcription factor GATA1 can activate and repress transcription. Vermunt et al. show that GATA1-induced silencing of the Kit gene is associated with GATA1-induced formation of a transient enhancer that delays Kit silencing. Co-factor usage determines the self-limiting nature of GATA1-occupied transient regulatory elements. [ABSTRACT FROM AUTHOR]

During mitosis, RNA polymerase II (Pol II) and many transcription factors dissociate from chromatin, and transcription ceases globally. Transcription is known to restart in bulk by telophase, but whether de novo transcription at the mitosis–G1 transition is in any way distinct from later in interphase remains unknown. We tracked Pol II occupancy genome-wide in mammalian cells progressing from mitosis through late G1. Unexpectedly, during the earliest rounds of transcription at the mitosis–G1 transition, ∼50% of active genes and distal enhancers exhibit a spike in transcription, exceeding levels observed later in G1 phase. Enhancer–promoter chromatin contacts are depleted during mitosis and restored rapidly upon G1 entry but do not spike. Of the chromatin-associated features examined, histone H3 Lys27 acetylation levels at individual loci in mitosis best predict the mitosis–G1 transcriptional spike. Single-molecule RNA imaging supports that the mitosis–G1 transcriptional spike can constitute the maximum transcriptional activity per DNA copy throughout the cell division cycle. The transcriptional spike occurs heterogeneously and propagates to cell-to-cell differences in mature mRNA expression. Our results raise the possibility that passage through the mitosis–G1 transition might predispose cells to diverge in gene expression states. [ABSTRACT FROM AUTHOR]

Summary Long non-coding (lnc) RNAs can regulate gene expression and protein functions. However, the proportion of lncRNAs with biological activities among the thousands expressed in mammalian cells is controversial. We studied Lockd (lncRNA downstream of Cdkn1b ), a 434-nt polyadenylated lncRNA originating 4 kb 3′ to the Cdkn1b gene. Deletion of the 25-kb Lockd locus reduced Cdkn1b transcription by approximately 70% in an erythroid cell line. In contrast, homozygous insertion of a polyadenylation cassette 80 bp downstream of the Lockd transcription start site reduced the entire lncRNA transcript level by >90% with no effect on Cdkn1b transcription. The Lockd promoter contains a DNase-hypersensitive site, binds numerous transcription factors, and physically associates with the Cdkn1b promoter in chromosomal conformation capture studies. Therefore, the Lockd gene positively regulates Cdkn1b transcription through an enhancer-like cis element, whereas the lncRNA itself is dispensable, which may be the case for other lncRNAs. [ABSTRACT FROM AUTHOR]

Erythroid development and differentiation from multiprogenitor cells into red blood cells requires precise transcriptional regulation. Key erythroid transcription factors, GATA1 and TAL1, cooperate, along with other proteins, to regulate many aspects of this process. How GATA1 and TAL1 are juxtaposed along the DNA and their cognate DNA binding site across the mouse genome remains unclear. We applied high-resolution ChIP-exo (chromatin immunoprecipitation followed by 5'-to-3' exonuclease treatment and then massively parallel DNA sequencing) to GATA1 and TAL1 to study their positional organization across the mouse genome during GATA1-dependent maturation. Two complementary methods, MultiGPS and peak pairing, were used to determine high-confidence binding locations by ChIP-exo. We identified ~10,000 GATA1 and ~15,000 TAL1 locations, which were essentially confirmed by ChIP-seq (chromatin immunoprecipitation followed by massively parallel DNA sequencing). Of these, 4,000 locations were bound by both GATA1 and TAL1. About three-quarters of them were tightly linked to a partial E-box located 7 or 8 bp upstream of a WGATAA motif. Both TAL1 and GATA1 generated distinct characteristic ChIP-exo peaks around WGATAA motifs that reflect their positional arrangement within a complex. We show that TAL1 and GATA1 form a precisely organized complex at a compound motif consisting of a TG 7 or 8 bp upstream of a WGATAA motif across thousands of genomic locations. [ABSTRACT FROM AUTHOR]

Erythroid development and differentiation from multiprogenitor cells into red blood cells requires precise transcriptional regulation. Key erythroid transcription factors, GATA1 and TAL1, cooperate, along with other proteins, to regulate many aspects of this process. How GATA1 and TAL1 are juxtaposed along the DNA and their cognate DNA binding site across the mouse genome remains unclear. We applied high-resolution ChIP-exo (chromatin immunoprecipitation followed by 5'-to-3' exonuclease treatment and then massively parallel DNA sequencing) to GATA1 and TAL1 to study their positional organization across the mouse genome during GATA1-dependent maturation. Two complementary methods, MultiGPS and peak pairing, were used to determine high-confidence binding locations by ChIP-exo. We identified ~ 10,000 GATA1 and ~ 15,000 TAL1 locations, which were essentially confirmed by ChIP-seq (chromatin immunoprecipitation followed by massively parallel DNA sequencing). Of these, ~4,000 locations were bound by both GATA1 and TAL1. About three-quarters of them were tightly linked to a partial E-box located 7 or 8 bp upstream of a WGATAA motif. Both TAL1 and GATA1 generated distinct characteristic ChIP-exo peaks around WGATAA motifs that reflect their positional arrangement within a complex. We show that TAL1 and GATA1 form a precisely organized complex at a compound motif consisting of a TG 7 or 8 bp upstream of a WGATAA motif across thousands of genomic locations. [ABSTRACT FROM AUTHOR]

Background: Regulated gene expression controls organismal development, and variation in regulatory patterns has been implicated in complex traits. Thus accurate prediction of enhancers is important for further understanding of these processes. Genome-wide measurement of epigenetic features, such as histone modifications and occupancy by transcription factors, is improving enhancer predictions, but the contribution of these features to prediction accuracy is not known. Given the importance of the hematopoietic transcription factor TAL1 for erythroid gene activation, we predicted candidate enhancers based on genomic occupancy by TAL1 and measured their activity. Contributions of multiple features to enhancer prediction were evaluated based on the results of these and other studies. Results: TAL1-bound DNA segments were active enhancers at a high rate both in transient transfections of cultured cells (39 of 79, or 56%) and transgenic mice (43 of 66, or 65%). The level of binding signal for TAL1 or GATA1 did not help distinguish TAL1-bound DNA segments as active versus inactive enhancers, nor did the density of regulation-related histone modifications. A meta-analysis of results from this and other studies (273 tested predicted enhancers) showed that the presence of TAL1, GATA1, EP300, SMAD1, H3K4 methylation, H3K27ac, and CAGE tags at DNase hypersensitive sites gave the most accurate predictors of enhancer activity, with a success rate over 80% and a median threefold increase in activity. Chromatin accessibility assays and the histone modifications H3K4me1 and H3K27ac were sensitive for finding enhancers, but they have high false positive rates unless transcription factor occupancy is also included. Conclusions: Occupancy by key transcription factors such as TAL1, GATA1, SMAD1, and EP300, along with evidence of transcription, improves the accuracy of enhancer predictions based on epigenetic features. [ABSTRACT FROM AUTHOR]

Inhibitors of bromodomain and extraterminal motif proteins (BETs) are being evaluated for the treatment of cancer and other diseases, yet much remains to be learned about how BET proteins function during normal physiology. We used genomic and genetic approaches to examine BET function in a hematopoietic maturation system driven by GATA1, an acetylated transcription factor previously shown to interact with BETs. We found that BRD2, BRD3, and BRD4 were variably recruited to GATA1-regulated genes, with BRD3 binding the greatest number of G ATA1-occupied sites. Pharmacologic BET inhibition impaired GATA1-mediated transcriptional activation, but not repression, genome-wide. Mechanistically, BETs promoted chromatin occupancy of GATA1 and sub sequently supported transcriptional activation. Using a combination of CRISPRC as 9 -mediated genomic engineering and shRNA approaches , we observed that depletion of either BRD2 or BRD4 alone bluntedery throid gene activation. Surprisingly, depletion of BRD3 only affectedery throid transcription in the context of BRD2 deficiency. Consistent with functional overlap among BET proteins, forced BRD3 expression substantially rescued defects caused by BRD2 deficiency. These results suggest that pharmacologic BET inhibition should be interpreted in the context of distinct steps in transcriptional activation and overlapping functions among BET family members. [ABSTRACT FROM AUTHOR]

Germline CATA1 mutations that result in the production of an amino-truncated protein termed GATA1s (where s indicates short) cause congenital hypoplastic anemia. In patients with trisomy 21, similar somatic GATA1s-producing mutations promote transient myeloproliferative disease and acute megakaryoblastic leukemia. Here, we demonstrate that induced pluripotent stem cells (iPSCs) from patients with GATA1-truncating mutations exhibit impaired erythroid potential, but enhanced megakaryopoiesis and myelopoiesis, recapitulatingthe major phenotypes of the associated diseases. Similarly, in developmentally arrested GATA1-deficient murine megakaryocyte-erythroid progenitors derived from murine embryonic stem cells (ESCs), expression of GATA1s promoted megakaryopoiesis, but not erythropoiesis. Transcriptome analysis revealed a selective deficiency in the ability of GATA1s to activate erythroid-specific genes within populations of hematopoietic progenitors. Although its DNA-binding domain was intact, chromatin immunoprecipitation studies showed that GATA1s binding at specific erythroid regulatory regions was impaired, while binding at many nonerythroid sites, including megakaryocytic and myeloid target genes, was normal. Together, these observations indicate that lineage-specific GATA1 cofactor associations are essential for normal chromatin occupancy and provide mechanistic insights into how GATAIs mutations cause human disease. More broadly, our studies underscore the value of ESCs and iPSCs to recapitulate and study disease phenotypes. [ABSTRACT FROM AUTHOR]

Background: Because species-specific gene expression is driven by species-specific regulation, understanding the relationship between sequence and function of the regulatory regions in different species will help elucidate how differences among species arise. Despite active experimental and computational research, relationships among sequence, conservation, and function are still poorly understood. Results: We compared transcription factor occupied segments (TFos) for 116 human and 35 mouse TFs in 546 human and 125 mouse cell types and tissues from the Human and the Mouse ENCODE projects. We based the map between human and mouse TFos on a one-to-one nucleotide cross-species mapper, bnMapper, that utilizes whole genome alignments (WGA). Our analysis shows that TFos are under evolutionary constraint, but a substantial portion (25.1% of mouse and 25.85% of human on average) of the TFos does not have a homologous sequence on the other species; this portion varies among cell types and TFs. Furthermore, 47.67% and 57.01% of the homologous TFos sequence shows binding activity on the other species for human and mouse respectively. However, 79.87% and 69.22% is repurposed such that it binds the same TF in different cells or different TFs in the same cells. Remarkably, within the set of repurposed TFos, the corresponding genome regions in the other species are preferred locations of novel TFos. These events suggest exaptation of some functional regulatory sequences into new function. Despite TFos repurposing, we did not find substantial changes in their predicted target genes, suggesting that CRMs buffer evolutionary events allowing little or no change in the TFos - target gene associations. Thus, the small portion of TFos with strictly conserved occupancy underestimates the degree of conservation of regulatory interactions. Conclusion: We mapped regulatory sequences from an extensive number of TFs and cell types between human and mouse using WGA. A comparative analysis of this correspondence unveiled the extent of the shared regulatory sequence across TFs and cell types under study. Importantly, a large part of the shared regulatory sequence is repurposed on the other species. This sequence, fueled by turnover events, provides a strong case for exaptation in regulatory elements. [ABSTRACT FROM AUTHOR]

Eukaryotic chromosomes replicate in a temporal order known as the replication-timing program. In mammals, replication timing is cell-type-specific with at least half the genome switching replication timing during development, primarily in units of 400-800 kilobases ('replication domains'), whose positions are preserved in different cell types, conserved between species, and appear to confine long-range effects of chromosome rearrangements. Early and late replication correlate, respectively, with open and closed three-dimensional chromatin compartments identified by high-resolution chromosome conformation capture (Hi-C), and, to a lesser extent, late replication correlates with lamina-associated domains (LADs). Recent Hi-C mapping has unveiled substructure within chromatin compartments called topologically associating domains (TADs) that are largely conserved in their positions between cell types and are similar in size to replication domains. However, TADs can be further sub-stratified into smaller domains, challenging the significance of structures at any particular scale. Moreover, attempts to reconcile TADs and LADs to replication-timing data have not revealed a common, underlying domain structure. Here we localize boundaries of replication domains to the early-replicating border of replication-timing transitions and map their positions in 18 human and 13 mouse cell types. We demonstrate that, collectively, replication domain boundaries share a near one-to-one correlation with TAD boundaries, whereas within a cell type, adjacent TADs that replicate at similar times obscure replication domain boundaries, largely accounting for the previously reported lack of alignment. Moreover, cell-type-specific replication timing of TADs partitions the genome into two large-scale sub-nuclear compartments revealing that replication-timing transitions are indistinguishable from late-replicating regions in chromatin composition and lamina association and accounting for the reduced correlation of replication timing to LADs and heterochromatin. Our results reconcile cell-type-specific sub-nuclear compartmentalization and replication timing with developmentally stable structural domains and offer a unified model for large-scale chromosome structure and function. [ABSTRACT FROM AUTHOR]

Mammals express thousands of long noncoding (lnc) RNAs, a few of which are known to function in tissue development. However, the entire repertoire of lncRNAs in most tissues and species is not defined. Indeed, most lncRNAs are not conserved, raising questions about function. We used RNA sequencing to identify 1109 polyadenylated lncRNAs expressed in erythroblasts, megakaryocytes, and megakaryocyte-erythroid precursors of mice, and 594 in erythroblasts of humans. More than half of these lncRNAs were unannotated, emphasizing the opportunity for new discovery through studies of specialized cell types. Analysis of the mouse erythro-megakaryocytic polyadenylated lncRNA transcriptome indicates that ~75% arise from promoters and 25% from enhancers, many of which are regulated by key transcription factors including GATA1 and TAL1. Erythroid lncRNA expression is largely conserved among 8 different mouse strains, yet only 15% of mouse lncRNAs are expressed in humans and vice versa, reflecting dramatic species-specificity. RNA interference assays of 21 abundant erythroid-specific murine lncRNAs in primary mouse erythroid precursors identified 7 whose knockdown inhibited terminal erythroid maturation. At least 6 of these 7 functional lncRNAs have no detectable expression in human erythroblasts, suggesting that lack of conservation between mammalian species does not predict lack of function. [ABSTRACT FROM AUTHOR]

Identification of cell type-specific enhancers is important for understanding the regulation of programs controlling cellular development and differentiation. Enhancers are typically marked by the co-transcriptional activator protein p300 or by groups of cell-expressed transcription factors. We hypothesized that a unique set of enhancers regulates gene expression in human erythroid cells, a highly specialized cell type evolved to provide adequate amounts of oxygen throughout the body. Using chromatin immunoprecipitation followed by massively parallel sequencing, genome-wide maps of candidate enhancers were constructed for p300 and four transcription factors, GATA1, NF-E2, KLF1, and SCL, using primary human erythroid cells. These data were combined with gene expression analyses, and candidate enhancers were identified. Consistent with their predicted function as candidate enhancers, there was statistically significant enrichment of p300 and combinations of co-localizing erythroid transcription factors within 1-50 kb of the transcriptional start site (TSS) of genes highly expressed in erythroid cells. Candidate enhancers were also enriched near genes with known erythroid cell function or phenotype. Candidate enhancers exhibited moderate conservation with mouse and minimal conservation with nonplacental vertebrates. Candidate enhancers were mapped to a set of erythroid-associated, biologically relevant, SNPs from the genome-wide association studies (GWAS) catalogue of NHGRI, National Institutes of Health. Fourteen candidate enhancers, representing 10 genetic loci, mapped to sites associated with biologically relevant erythroid traits. Fragments from these loci directed statistically significant expression in reporter gene assays. Identification of enhancers in human erythroid cells will allow a better understanding of erythroid cell development, differentiation, structure, and function and provide insights into inherited and acquired hematologic disease. [ABSTRACT FROM AUTHOR]

Patients with Down syndrome (trisomy 21, T21) have hematologic abnormalities throughout life. Newborns frequently exhibit abnormal blood counts and a clonal preleukemia. Human T21 fetal livers contain expanded erythro-megakaryocytic precursors with enhanced proliferative capacity. The impact of T21 on the earliest stages of embryonic hematopoiesis is unknown and nearly impossible to examine in human subjects. We modeled T21 yolk sac hematopoiesis using human induced pluripotent stem cells (iPSCs). Blood progenitor populations generated from T21 iPSCs were present at normal frequency and proliferated normally. However, their developmental potential was altered with enhanced erythropoiesis and reduced myelopoiesis, but normal megakaryocyte production. These abnormalities overlap with those of T21 fetal livers, but also reflect important differences. Our studies show that T21 confers distinct developmental stage- and species-specific hematopoietic defects. More generally, we illustrate how iPSCs can provide insight into early stages of normal and pathological human development. [ABSTRACT FROM AUTHOR]

Many software tools for comparative analysis of genomic sequence data have been released in recent decades. Despite this, it remains challenging to determine evolutionary relationships in gene clusters due to their complex histories involving duplications, deletions, inversions, and conversions. One concept describing these relationships is orthology. Orthologs derive from a common ancestor by speciation, in contrast to paralogs, which derive from duplication. Discriminating orthologs from paralogs is a necessary step in most multispecies sequence analyses, but doing so accurately is impeded by the occurrence of gene conversion events. We propose a refined method of orthology assignment based on two paradigms for interpreting its definition: by genomic context or by sequence content. X-orthology (based on context) traces orthology resulting from speciation and duplication only, while N-orthology (based on content) includes the influence of conversion events. We developed a computational method for automatically mapping both types of orthology on a per-nucleotide basis in gene cluster regions studied by comparative sequencing, and we make this mapping accessible by visualizing the output. All of these steps are incorporated into our newly extended CHAP 2 package. We evaluate our method using both simulated data and real gene clusters (including the well-characterized α-globin and β-globin clusters). We also illustrate use of CHAP 2 by analyzing four more loci: CCL (chemokine ligand), IFN (interferon), CYP2abf (part of cytochrome P450 family 2), and KIR (killer cell immunoglobulin-like receptors). These new methods facilitate and extend our understanding of evolution at these and other loci by adding automated accurate evolutionary inference to the biologist's toolkit. The CHAP 2 package is freely available from http://www.bx.psu.edu/miller_lab. [ABSTRACT FROM AUTHOR]

Summary: Tissue-specific transcription patterns are preserved throughout cell divisions to maintain lineage fidelity. We investigated whether transcription factor GATA1 plays a role in transmitting hematopoietic gene expression programs through mitosis when transcription is transiently silenced. Live-cell imaging revealed that a fraction of GATA1 is retained focally within mitotic chromatin. ChIP-seq of highly purified mitotic cells uncovered that key hematopoietic regulatory genes are occupied by GATA1 in mitosis. The GATA1 coregulators FOG1 and TAL1 dissociate from mitotic chromatin, suggesting that GATA1 functions as platform for their postmitotic recruitment. Mitotic GATA1 target genes tend to reactivate more rapidly upon entry into G1 than genes from which GATA1 dissociates. Mitosis-specific destruction of GATA1 delays reactivation selectively of genes that retain GATA1 during mitosis. These studies suggest a requirement of mitotic “bookmarking” by GATA1 for the faithful propagation of cell-type-specific transcription programs through cell division. [ABSTRACT FROM AUTHOR]

Many software tools for comparative analysis of genomic sequence data have been released in recent decades. Despite this, it remains challenging to determine evolutionary relationships in gene clusters due to their complex histories involving duplications, deletions, inversions, and conversions. One concept describing these relationships is orthology. Orthologs derive from a common ancestor by speciation, in contrast to paralogs, which derive from duplication. Discriminating orthologs from paralogs is a necessary step in most multispecies sequence analyses, but doing so accurately is impeded by the occurrence of gene conversion events. We propose a refined method of orthology assignment based on two paradigms for interpreting its definition: by genomic context or by sequence content. X-orthology (based on context) traces orthology resulting from speciation and duplication only, while N-orthology (based on content) includes the influence of conversion events. We developed a computational method for automatically mapping both types of orthology on a per-nucleotide basis in gene cluster regions studied by comparative sequencing, and we make this mapping accessible by visualizing the output. All of these steps are incorporated into our newly extended CHAP 2 package. We evaluate our method using both simulated data and real gene clusters (including the well-characterized α-globin and β-globin clusters). We also illustrate use of CHAP 2 by analyzing four more loci: CCL (chemokine ligand), IFN (interferon), CYP2abf (part of cytochrome P450 family 2), and KIR (killer cell immunoglobulin-like receptors). These new methods facilitate and extend our understanding of evolution at these and other loci by adding automated accurate evolutionary inference to the biologist's toolkit. The CHAP 2 package is freely available from http://www.bx.psu.edu/miller_lab. [ABSTRACT FROM PUBLISHER]

This article describes a set of alignments of 28 vertebrate genome sequences that is provided by the UCSC Genome Browser. The alignments can be viewed on the Human Genome Browser (March 2006 assembly) at http://genome.ucsc.edu, downloaded in bulk by anonymous FTP from http://hgdownload.cse.ucsc.edu/goldenPath/hg18/multiz28way, or analyzed with the Galaxy server at http://g2.bx.psu.edu. This article illustrates the power of this resource for exploring vertebrate and mammalian evolution, using three examples. First, we present several vignettes involving insertions and deletions within protein-coding regions, including a look at some human-specific indels. Then we study the extent to which start codons and stop codons in the human sequence are conserved in other species, showing that start codons are in general more poorly conserved than stop codons. Finally, an investigation of the phylogenetic depth of conservation for several classes of functional elements in the human genome reveals striking differences in the rates and modes of decay in alignability. Each functional class has a distinctive period of stringent constraint, followed by decays that allow (for the case of regulatory regions) or reject (for coding regions and ultraconserved elements) insertions and deletions. [ABSTRACT FROM AUTHOR]

Summary: The transcription factor GATA-1 is required for terminal erythroid maturation and functions as an activator or repressor depending on gene context. Yet its in vivo site selectivity and ability to distinguish between activated versus repressed genes remain incompletely understood. In this study, we performed GATA-1 ChIP-seq in erythroid cells and compared it to GATA-1-induced gene expression changes. Bound and differentially expressed genes contain a greater number of GATA-binding motifs, a higher frequency of palindromic GATA sites, and closer occupancy to the transcriptional start site versus nondifferentially expressed genes. Moreover, we show that the transcription factor Zbtb7a occupies GATA-1-bound regions of some direct GATA-1 target genes, that the presence of SCL/TAL1 helps distinguish transcriptional activation versus repression, and that polycomb repressive complex 2 (PRC2) is involved in epigenetic silencing of a subset of GATA-1-repressed genes. These data provide insights into GATA-1-mediated gene regulation in vivo. [Copyright &y& Elsevier]

Hemoglobin production during erythropoiesis is mechanistically coupled to the acquisition and metabolism of iron. We discovered that iron regulates the expression of α-hemoglobin-stabilizing protein (AHSP), a molecular chaperone that binds and stabilizes free α-globin during hemoglobin synthesis. In primates, the 3′-untranslated region (UTR) of AHSP mRNA contains a nucleotide sequence resembling iron responsive elements (IREs), stem-loop structures that regulate gene expression post-transcriptionally by binding iron regulatory proteins (IRPs). The AHSP IRE-like stem-loop deviates from classical consensus sequences and binds IRPs poorly in electrophoretic mobility shift assays. However, in cytoplasmic extracts, AHSP mRNA co-immunoprecipitates with IRPs in a fashion that is dependent on the stem-loop structure and inhibited by iron. Moreover, this interaction enhances AHSP mRNA stability in erythroid and heterologous cells. Our findings demonstrate that IRPs can regulate mRNA expression through non-canonical IREs and extend the repertoire of known iron-regulated genes. In addition, we illustrate a new mechanism through which hemoglobin may be modulated according to iron status. [ABSTRACT FROM AUTHOR]