Experimental validation of predicted mammalian erythroid cis-regulatory modules

Protocol for efficient CRISPR-Cas9-mediated fluorescent tag knockin in hard-to-transfect erythroid cell lines

Precise insertion of fluorescent tags by CRISPR-Cas9-mediated homologous recombination (HR) in mammalian genes is a powerful tool allowing to study gene function and protein gene products. Here, we present a protocol for efficient HR-mediated targeted insertion of fluorescent markers in the genome of hard-to-transfect erythroid cell lines MEL (mouse erythroleukemic) and MEDEP (mouse ES cell-derived erythroid progenitor line). We describe steps for plasmid construction, electroporation, amplification, and verification of genome editing. We then detail procedures for isolating positive clones and validating knockin clones. For complete details on the use and execution of this protocol, please refer to Deleuze et al.1.

Efficient genome editing in erythroid cells unveils novel MYB target genes and regulatory functions

Targeted genome editing holds great promise in biology. However, efficient genome modification, including gene knock-in (KI), remains an unattained goal in multiple cell types and loci due to poor transfection efficiencies and low target genes expression, impeding the positive selection of recombined cells. Here, we describe a genome editing approach to achieve efficient gene targeting using hard to transfect erythroid cell lines. We demonstrate robust fluorescent protein KI efficiency in low expressed transcription factor (TF) genes (e.g., Myb or Zeb1). We further show the ability to target two independent loci in individual cells, exemplified by MYB-GFP and NuMA-Cherry double KI, allowing multicolor labeling of regulatory factors at physiological endogenous levels. Our KI tagging approach allowed us to perform genome-wide TF analysis at increased signal-to-noise ratios, and highlighted previously unidentified MYB target genes and pathways. Overall, we establish a versatile CRISPR-Cas9-based platform, offering attractive opportunities for the dissection of the erythroid differentiation process.

Deciphering the molecular basis of tissue-specific gene expression in plants: Can synthetic biology help?

Gene expression differences between distinct cell types are orchestrated by specific sets of transcription factors and epigenetic regulators acting upon the genome. In plants, the mechanisms underlying tissue-specific gene activity remain largely unexplored. Although transcriptional and epigenetic profiling of individual organs, tissues, and more recently, of single cells can easily detect the molecular signatures of different biological samples, how these unique cell identities are established at the mechanistic level is only beginning to be decoded. Computational methods, including machine learning, used in combination with experimental approaches, enable the identification and validation of candidate cis-regulatory elements driving cell-specific expression. Synthetic biology shows great promise not only as a means of testing candidate DNA motifs but also for establishing the general rules of nature driving promoter architecture and for the rational design of genetic circuits in research and agriculture to confer tissue-specific expression to genes or molecular pathways of interest.

RefSeq Functional Elements as experimentally assayed nongenic reference standards and functional interactions in human and mouse

Eukaryotic genomes contain many nongenic elements that function in gene regulation, chromosome organization, recombination, repair, or replication, and mutation of those elements can affect genome function and cause disease. Although numerous epigenomic studies provide high coverage of gene regulatory regions, those data are not usually exposed in traditional genome annotation and can be difficult to access and interpret without field-specific expertise. The National Center for Biotechnology Information (NCBI) therefore provides RefSeq Functional Elements (RefSeqFEs), which represent experimentally validated human and mouse nongenic elements derived from the literature. The curated data set is comprised of richly annotated sequence records, descriptive records in the NCBI Gene database, reference genome feature annotation, and activity-based interactions between nongenic regions, target genes, and each other. The data set provides succinct functional details and transparent experimental evidence, leverages data from multiple experimental sources, is readily accessible and adaptable, and uses a flexible data model. The data have multiple uses for basic functional discovery, bioinformatics studies, genetic variant interpretation; as known positive controls for epigenomic data evaluation; and as reference standards for functional interactions. Comparisons to other gene regulatory data sets show that the RefSeqFE data set includes a wider range of feature types representing more areas of biology, but it is comparatively smaller and subject to data selection biases. RefSeqFEs thus provide an alternative and complementary resource for experimentally assayed functional elements, with future data set growth expected.

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

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

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

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

Regulatory Single-Nucleotide Variant Predictor Increases Predictive Performance of Functional Regulatory Variants

In silico methods for detecting functionally relevant genetic variants are important for identifying genetic markers of human inherited disease. Much research has focused on protein-coding variants since coding regions have well-defined physicochemical and functional properties. However, many bioinformatics tools are not applicable to variants outside coding regions. Here, we increase the classification performance of our regulatory single-nucleotide variant predictor (RSVP) for variants that cause regulatory abnormalities from an AUC of 0.90-0.97 by incorporating genomic regions identified by the ENCODE project into RSVP. RSVP is comparable to a recently published tool, Genome-Wide Annotation of Variants (GWAVA); both RSVP and GWAVA perform better on regulatory variants than a traditional variant predictor, combined annotation-dependent depletion (CADD). However, our method outperforms GWAVA on variants located at similar distances to the transcription start site as the positive set (AUC: 0.96) as compared with GWAVA (AUC: 0.71). Much of this disparity is due to RSVP's incorporation of features pertaining to the nearest gene (expression, GO terms, etc.), which are not included in GWAVA. Our findings hold out the promise of a framework for the assessment of all functional regulatory variants, providing a means to predict which rare or de novo variants are of pathogenic significance.

Insight into GATA1 transcriptional activity through interrogation of cis elements disrupted in human erythroid disorders

Whole-exome sequencing has been incredibly successful in identifying causal genetic variants and has revealed a number of novel genes associated with blood and other diseases. One limitation of this approach is that it overlooks mutations in noncoding regulatory elements. Furthermore, the mechanisms by which mutations in transcriptionalcis-regulatory elements result in disease remain poorly understood. Here we used CRISPR/Cas9 genome editing to interrogate three such elements harboring mutations in human erythroid disorders, which in all cases are predicted to disrupt a canonical binding motif for the hematopoietic transcription factor GATA1. Deletions of as few as two to four nucleotides resulted in a substantial decrease (>80%) in target gene expression. Isolated deletions of the canonical GATA1 binding motif completely abrogated binding of the cofactor TAL1, which binds to a separate motif. Having verified the functionality of these three GATA1 motifs, we demonstrate strong evolutionary conservation of GATA1 motifs in regulatory elements proximal to other genes implicated in erythroid disorders, and show that targeted disruption of such elements results in altered gene expression. By modeling transcription factor binding patterns, we show that multiple transcription factors are associated with erythroid gene expression, and have created predictive maps modeling putative disruptions of their binding sites at key regulatory elements. Our study provides insight into GATA1 transcriptional activity and may prove a useful resource for investigating the pathogenicity of noncoding variants in human erythroid disorders.

Accelerated Evolution of Enhancer Hotspots in the Mammal Ancestor

Mammals have evolved remarkably different sensory, reproductive, metabolic, and skeletal systems. To explore the genetic basis for these differences, we developed a comparative genomics approach to scan whole-genome multiple sequence alignments to identify regions that evolved rapidly in an ancestral lineage but are conserved within extant species. This pattern suggests that ancestral changes in function were maintained in descendants. After applying this test to therian mammals, we identified 4,797 accelerated regions, many of which are noncoding and located near developmental transcription factors. We then used mouse transgenic reporter assays to test if noncoding accelerated regions are enhancers and to determine how therian-specific substitutions affect their activity in vivo. We discovered enhancers with expression specific to the therian version in brain regions involved in the hormonal control of milk ejection, uterine contractions, blood pressure, temperature, and visual processing. This work underscores the idea that changes in developmental gene expression are important for mammalian evolution, and it pinpoints candidate genes for unique aspects of mammalian biology.

Occupancy by key transcription factors is a more accurate predictor of enhancer activity than histone modifications or chromatin accessibility

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-015-0009-5) contains supplementary material, which is available to authorized users.

New insights on the transcriptional regulation of CD69 gene through a potent enhancer located in the conserved non-coding sequence 2

The CD69 type II C-type lectin is one of the earliest indicators of leukocyte activation acting in lymphocyte migration and cytokine secretion. CD69 expression in hematopoietic lineage undergoes rapid changes depending on the cell-lineage, the activation state or the localization of the cell where it is expressed, suggesting a complex and tightly controlled regulation. Here we provide new insights on the transcriptional regulation of CD69 gene in mammal species. Through in silico studies, we analyzed several regulatory features of the 4 upstream conserved non-coding sequences (CNS 1-4) previously described, confirming a major function of CNS2 in the transcriptional regulation of CD69. In addition, multiple transcription binding sites are identified in the CNS2 region by DNA cross-species conservation analysis. By functional approaches we defined a core region of 226bp located within CNS2 as the main enhancer element of CD69 transcription in the hematopoietic cells analyzed. By chromatin immunoprecipitation, binding of RUNX1 to the core-CNS2 was shown in a T cell line. In addition, we found an activating but not essential role of RUNX1 in CD69 gene transcription by site-directed mutagenesis and RNA silencing, probably through the interaction with this potent enhancer specifically in the hematopoietic lineage. In summary, in this study we contribute with new evidences to the landscape of the transcriptional regulation of the CD69 gene.

X-linked sideroblastic anemia due to ALAS2 intron 1 enhancer element GATA-binding site mutations

X-linked sideroblastic anemia (XLSA) is the most common form of congenital sideroblastic anemia. In affected males, it is uniformly associated with partial loss-of-function missense mutations in the erythroid-specific heme biosynthesis protein 5-aminolevulinate synthase 2 (ALAS2). Here, we report five families with XLSA owing to mutations in a GATA transcription factor binding site located in a transcriptional enhancer element in intron 1 of the ALAS2 gene. As such, this study defines a new class of mutations that should be evaluated in patients undergoing genetic testing for a suspected diagnosis of XLSA.

Nuclear RNA sequencing of the mouse erythroid cell transcriptome

In addition to protein coding genes a substantial proportion of mammalian genomes are transcribed. However, most transcriptome studies investigate steady-state mRNA levels, ignoring a considerable fraction of the transcribed genome. In addition, steady-state mRNA levels are influenced by both transcriptional and posttranscriptional mechanisms, and thus do not provide a clear picture of transcriptional output. Here, using deep sequencing of nuclear RNAs (nucRNA-Seq) in parallel with chromatin immunoprecipitation sequencing (ChIP-Seq) of active RNA polymerase II, we compared the nuclear transcriptome of mouse anemic spleen erythroid cells with polymerase occupancy on a genome-wide scale. We demonstrate that unspliced transcripts quantified by nucRNA-seq correlate with primary transcript frequencies measured by RNA FISH, but differ from steady-state mRNA levels measured by poly(A)-enriched RNA-seq. Highly expressed protein coding genes showed good correlation between RNAPII occupancy and transcriptional output; however, genome-wide we observed a poor correlation between transcriptional output and RNAPII association. This poor correlation is due to intergenic regions associated with RNAPII which correspond with transcription factor bound regulatory regions and a group of stable, nuclear-retained long non-coding transcripts. In conclusion, sequencing the nuclear transcriptome provides an opportunity to investigate the transcriptional landscape in a given cell type through quantification of unspliced primary transcripts and the identification of nuclear-retained long non-coding RNAs.

Genomic approaches towards finding cis-regulatory modules in animals

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.

What fraction of the human genome is functional?

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.

Dynamics of the epigenetic landscape during erythroid differentiation after GATA1 restoration

Interplays among lineage-specific nuclear proteins, chromatin modifying enzymes, and the basal transcription machinery govern cellular differentiation, but their dynamics of action and coordination with transcriptional control are not fully understood. Alterations in chromatin structure appear to establish a permissive state for gene activation at some loci, but they play an integral role in activation at other loci. To determine the predominant roles of chromatin states and factor occupancy in directing gene regulation during differentiation, we mapped chromatin accessibility, histone modifications, and nuclear factor occupancy genome-wide during mouse erythroid differentiation dependent on the master regulatory transcription factor GATA1. Notably, despite extensive changes in gene expression, the chromatin state profiles (proportions of a gene in a chromatin state dominated by activating or repressive histone modifications) and accessibility remain largely unchanged during GATA1-induced erythroid differentiation. In contrast, gene induction and repression are strongly associated with changes in patterns of transcription factor occupancy. Our results indicate that during erythroid differentiation, the broad features of chromatin states are established at the stage of lineage commitment, largely independently of GATA1. These determine permissiveness for expression, with subsequent induction or repression mediated by distinctive combinations of transcription factors.

Computational identification of transcriptional regulators in human endotoxemia

One of the great challenges in the post-genomic era is to decipher the underlying principles governing the dynamics of biological responses. As modulating gene expression levels is among the key regulatory responses of an organism to changes in its environment, identifying biologically relevant transcriptional regulators and their putative regulatory interactions with target genes is an essential step towards studying the complex dynamics of transcriptional regulation. We present an analysis that integrates various computational and biological aspects to explore the transcriptional regulation of systemic inflammatory responses through a human endotoxemia model. Given a high-dimensional transcriptional profiling dataset from human blood leukocytes, an elementary set of temporal dynamic responses which capture the essence of a pro-inflammatory phase, a counter-regulatory response and a dysregulation in leukocyte bioenergetics has been extracted. Upon identification of these expression patterns, fourteen inflammation-specific gene batteries that represent groups of hypothetically ‘coregulated’ genes are proposed. Subsequently, statistically significant cis-regulatory modules (CRMs) are identified and decomposed into a list of critical transcription factors (34) that are validated largely on primary literature. Finally, our analysis further allows for the construction of a dynamic representation of the temporal transcriptional regulatory program across the host, deciphering possible combinatorial interactions among factors under which they might be active. Although much remains to be explored, this study has computationally identified key transcription factors and proposed a putative time-dependent transcriptional regulatory program associated with critical transcriptional inflammatory responses. These results provide a solid foundation for future investigations to elucidate the underlying transcriptional regulatory mechanisms under the host inflammatory response. Also, the assumption that coexpressed genes that are functionally relevant are more likely to share some common transcriptional regulatory mechanism seems to be promising, making the proposed framework become essential in unravelling context-specific transcriptional regulatory interactions underlying diverse mammalian biological processes.

When needles look like hay: how to find tissue-specific enhancers in model organism genomes

A major prerequisite for the investigation of tissue-specific processes is the identification of cis-regulatory elements. No generally applicable technique is available to distinguish them from any other type of genomic non-coding sequence. Therefore, researchers often have to identify these elements by elaborate in vivo screens, testing individual regions until the right one is found. Here, based on many examples from the literature, we summarize how functional enhancers have been isolated from other elements in the genome and how they have been characterized in transgenic animals. Covering computational and experimental studies, we provide an overview of the global properties of cis-regulatory elements, like their specific interactions with promoters and target gene distances. We describe conserved non-coding elements (CNEs) and their internal structure, nucleotide composition, binding site clustering and overlap, with a special focus on developmental enhancers. Conflicting data and unresolved questions on the nature of these elements are highlighted. Our comprehensive overview of the experimental shortcuts that have been found in the different model organism communities and the new field of high-throughput assays should help during the preparation phase of a screen for enhancers. The review is accompanied by a list of general guidelines for such a project.

Assessing computational methods of cis-regulatory module prediction

Computational methods attempting to identify instances of cis-regulatory modules (CRMs) in the genome face a challenging problem of searching for potentially interacting transcription factor binding sites while knowledge of the specific interactions involved remains limited. Without a comprehensive comparison of their performance, the reliability and accuracy of these tools remains unclear. Faced with a large number of different tools that address this problem, we summarized and categorized them based on search strategy and input data requirements. Twelve representative methods were chosen and applied to predict CRMs from the Drosophila CRM database REDfly, and across the human ENCODE regions. Our results show that the optimal choice of method varies depending on species and composition of the sequences in question. When discriminating CRMs from non-coding regions, those methods considering evolutionary conservation have a stronger predictive power than methods designed to be run on a single genome. Different CRM representations and search strategies rely on different CRM properties, and different methods can complement one another. For example, some favour homotypical clusters of binding sites, while others perform best on short CRMs. Furthermore, most methods appear to be sensitive to the composition and structure of the genome to which they are applied. We analyze the principal features that distinguish the methods that performed well, identify weaknesses leading to poor performance, and provide a guide for users. We also propose key considerations for the development and evaluation of future CRM-prediction methods.

Transcriptional regulation involves multiple transcription factors binding to DNA sequences. A limited repertoire of transcription factors performs this complex regulatory step through various spatial and temporal interactions between themselves and their binding sites. These transcription factor binding interactions are clustered as distinct modules: cis-regulatory modules (CRMs). Computational methods attempting to identify instances of CRMs in the genome face a challenging problem because a majority of these interactions between transcription factors remain unknown. To investigate the reliability and accuracy of these methods, we chose twelve representative methods and applied them to predict CRMs on both the fly and human genomes. Our results show that the optimal choice of method varies depending on species and composition of the sequences in question. Different CRM representations and search strategies rely on different CRM properties, and different methods can complement one another. We provide a guide for users and key considerations for developers. We also expect that, along with new technology generating new types of genomic data, future CRM prediction methods will be able to reveal transcription binding interactions in three-dimensional space.

Comparative genomic analysis of the false killer whale (Pseudorca crassidens) LMBR1 locus

The sequencing and comparative genomic analysis of LMBR1 loci in mammals or other species, including human, would be very important in understanding evolutionary genetic changes underlying the evolution of limb development. In this regard, comparative genomic annotation of the false killer whale LMBR1 locus could shed new light on the evolution of limb development. We sequenced two false killer whale BAC clones, corresponding to 156 kb and 144 kb, respectively, harboring the tightly linked RNF32, LMBR1, and NOM1 genes. Our annotation of the false killer whale LMBR1 gene showed that it consists of 17 exons (1473 bp), in contrast to 18 exons (1596 bp) in human, and it displays 93.1% and 95.6% nucleotide and amino acid sequence similarity, respectively, compared with the human gene. In particular, we discovered that exon 10, deleted in the false killer whale LMBR1 gene, is present only in primates, and this fact strongly implies that exon 10 might be crucial in determining primate-specific limb development. ZRS and TFBS sequences have been well conserved across 11 species, suggesting that these regions could be involved in an important function of limb development and limb patterning. The neighboring gene RNF32 showed several lineage-conserved exons, such as exons 2 through 9 conserved in eutherian mammals, exons 3 through 9 conserved in mammals, and exons 5 through 9 conserved in vertebrates. The other neighboring gene, NOM1, had undergone a substitution (ATG→GTA) at the start codon, giving rise to a 36 bp shorter N-terminal sequence compared with the human sequence. Our comparative analysis of the false killer whale LMBR1 genomic locus provides important clues regarding the genetic regions that may play crucial roles in limb development and patterning.

A single cis element maintains repression of the key developmental regulator Gata2

In development, lineage-restricted transcription factors simultaneously promote differentiation while repressing alternative fates. Molecular dissection of this process has been challenging as transcription factor loci are regulated by many trans-acting factors functioning through dispersed cis elements. It is not understood whether these elements function collectively to confer transcriptional regulation, or individually to control specific aspects of activation or repression, such as initiation versus maintenance. Here, we have analyzed cis element regulation of the critical hematopoietic factor Gata2, which is expressed in early precursors and repressed as GATA-1 levels rise during terminal differentiation. We engineered mice lacking a single cis element −1.8 kb upstream of the Gata2 transcriptional start site. Although Gata2 is normally repressed in late-stage erythroblasts, the −1.8 kb mutation unexpectedly resulted in reactivated Gata2 transcription, blocked differentiation, and an aberrant lineage-specific gene expression pattern. Our findings demonstrate that the −1.8 kb site selectively maintains repression, confers a specific histone modification pattern and expels RNA Polymerase II from the locus. These studies reveal how an individual cis element establishes a normal developmental program via regulating specific steps in the mechanism by which a critical transcription factor is repressed.

Different cell types are formed and maintained by proteins called transcription factors that directly bind to specific DNA sequences to activate or repress gene expression. While numerous DNA sequences bound by transcription factors are established, many questions remain unanswered regarding how they function at specific sites located at distinct chromosomal regions. As a model to study this process, we examined the regulation of a gene controlling red blood cell development, Gata2, by the transcription factor GATA1. In the DNA sequence upstream of Gata2, there are several sites that GATA1 is known to bind to; however, it is unclear whether these binding sites work together or independently to control expression of Gata2. To study this, we engineered mice to specifically remove one of these GATA1-binding sites. We found that removal of this single site reactivated expression of Gata2 in a specific stage of red blood cell development where Gata2 is normally not expressed, caused a block in differentiation of these cells, and changed the histone modification pattern specifically in the region upstream of Gata2. This work supports a model in which individual transcription factor binding sites within regions of multiple binding sites can independently and distinctly regulate gene expression during development.

Hypoxia regulates BMP4 expression in the murine spleen during the recovery from acute anemia

Background

Bone marrow erythropoiesis is primarily homeostatic, producing new erythrocytes at a constant rate. However at times of acute anemia, new erythrocytes must be rapidly produced much faster than bone marrow steady state erythropoiesis. At these times stress erythropoiesis predominates. Stress erythropoiesis occurs in the fetal liver during embryogenesis and in the adult spleen and liver. In adult mice, stress erythropoiesis utilizes a specialized population of stress erythroid progenitors that are resident in the spleen. In response to acute anemia, these progenitors rapidly expand and differentiate in response to three signals, BMP4, SCF and hypoxia. In absence of acute anemic stress, two of these signals, BMP4 and hypoxia, are not present and the pathway is not active. The initiating event in the activation of this pathway is the up-regulation of BMP4 expression in the spleen.

Methodology/Principal Findings

In this paper we analyze the regulation of BMP4 expression in the spleen by hypoxia. Using stromal cell lines, we establish a role for hypoxia transcription factor HIFs (Hypoxia Inducible Factors) in the transcription of BMP4. We identified putative Hypoxia Responsive Elements (HREs) in the BMP4 gene using bioinformatics. Analysis of these elements showed that in vivo, Hif2α binds two cis regulatory sites in the BMP4 gene, which regulate BMP4 expression during the recovery from acute anemia.

Conclusions and Significance

These data show that hypoxia plays a key role in initiating the BMP4 dependent stress erythropoiesis pathway by regulating BMP4 expression.

Primary sequence and epigenetic determinants of in vivo occupancy of genomic DNA by GATA1

DNA sequence motifs and epigenetic modifications contribute to specific binding by a transcription factor, but the extent to which each feature determines occupancy in vivo is poorly understood. We addressed this question in erythroid cells by identifying DNA segments occupied by GATA1 and measuring the level of trimethylation of histone H3 lysine 27 (H3K27me3) and monomethylation of H3 lysine 4 (H3K4me1) along a 66 Mb region of mouse chromosome 7. While 91% of the GATA1-occupied segments contain the consensus binding-site motif WGATAR, only ∼0.7% of DNA segments with such a motif are occupied. Using a discriminative motif enumeration method, we identified additional motifs predictive of occupancy given the presence of WGATAR. The specific motif variant AGATAA and occurrence of multiple WGATAR motifs are both strong discriminators. Combining motifs to pair a WGATAR motif with a binding site motif for GATA1, EKLF or SP1 improves discriminative power. Epigenetic modifications are also strong determinants, with the factor-bound segments highly enriched for H3K4me1 and depleted of H3K27me3. Combining primary sequence and epigenetic determinants captures 52% of the GATA1-occupied DNA segments and substantially increases the specificity, to one out of seven segments with the required motif combination and epigenetic signals being bound.

Annotating the regulatory genome

Determining the timing and molecular repertoire responsible for gene expression is fundamental to understanding a gene's function. Heritable differences in this character are increasingly regarded as explanatory for complex and common traits. For many known trait-predisposing genes, studies have sought to elucidate the associated logic behind gene regulation. However, there exist many challenges in deciphering these mechanisms. Among them, it is recognized that we have limited understanding of regulatory complexity, the current models of gene regulation have low specificity and any gene's regulatory logic is dependent on biological context. Addressing these limitations and defining the regulatory genome is an ongoing challenge for molecular biology. We discuss current efforts to define and annotate the regulatory genome by focusing on curation and text-mining activities. We further highlight the type of information and curation process for describing regulatory elements within the ORegAnno database ( www.oreganno.org ) and how the general standards for such information are changing.

Erythroid GATA1 function revealed by genome-wide analysis of transcription factor occupancy, histone modifications, and mRNA expression

The transcription factor GATA1 regulates an extensive program of gene activation and repression during erythroid development. However, the associated mechanisms, including the contributions of distal versus proximal cis-regulatory modules, co-occupancy with other transcription factors, and the effects of histone modifications, are poorly understood. We studied these problems genome-wide in a Gata1 knockout erythroblast cell line that undergoes GATA1-dependent terminal maturation, identifying 2616 GATA1-responsive genes and 15,360 GATA1-occupied DNA segments after restoration of GATA1. Virtually all occupied DNA segments have high levels of H3K4 monomethylation and low levels of H3K27me3 around the canonical GATA binding motif, regardless of whether the nearby gene is induced or repressed. Induced genes tend to be bound by GATA1 close to the transcription start site (most frequently in the first intron), have multiple GATA1-occupied segments that are also bound by TAL1, and show evolutionary constraint on the GATA1-binding site motif. In contrast, repressed genes are further away from GATA1-occupied segments, and a subset shows reduced TAL1 occupancy and increased H3K27me3 at the transcription start site. Our data expand the repertoire of GATA1 action in erythropoiesis by defining a new cohort of target genes and determining the spatial distribution of cis-regulatory modules throughout the genome. In addition, we begin to establish functional criteria and mechanisms that distinguish GATA1 activation from repression at specific target genes. More broadly, these studies illustrate how a "master regulator" transcription factor coordinates tissue differentiation through a panoply of DNA and protein interactions.

The role and regulation of friend of GATA-1 (FOG-1) during blood development in the zebrafish

The nuclear protein FOG-1 binds transcription factor GATA-1 to facilitate erythroid and megakaryocytic maturation. However, little is known about the function of FOG-1 during myeloid and lymphoid development or how FOG-1 expression is regulated in any tissue. We used in situ hybridization, gain- and loss-of-function studies in zebrafish to address these problems. Zebrafish FOG-1 is expressed in early hematopoietic cells, as well as heart, viscera, and paraspinal neurons, suggesting that it has multifaceted functions in organogenesis. We found that FOG-1 is dispensable for endoderm specification but is required for endoderm patterning affecting the expression of late-stage T-cell markers, independent of GATA-1. The suppression of FOG-1, in the presence of normal GATA-1 levels, induces severe anemia and thrombocytopenia and expands myeloid-progenitor cells, indicating that FOG-1 is required during erythroid/myeloid commitment. To functionally interrogate whether GATA-1 regulates FOG-1 in vivo, we used bioinformatics combined with transgenic assays. Thus, we identified 2 cis-regulatory elements that control the tissue-specific gene expression of FOG-1. One of these enhancers contains functional GATA-binding sites, indicating the potential for a regulatory loop in which GATA factors control the expression of their partner protein FOG-1.

LRF is an essential downstream target of GATA1 in erythroid development and regulates BIM-dependent apoptosis

GATA-1-dependent transcription is essential for erythroid differentiation and maturation. Suppression of programmed cell death is also thought to be critical for this process; however, the link between these two features of erythropoiesis has remained elusive. Here, we show that the POZ-Krüppel family transcription factor, LRF (also known as Zbtb7a/Pokemon), is a direct target of GATA1 and plays an essential antiapoptotic role during terminal erythroid differentiation. We find that loss of Lrf leads to lethal anemia in embryos, due to increased apoptosis of late-stage erythroblasts. This programmed cell death is Arf and p53 independent and is instead mediated by upregulation of the proapoptotic factor Bim. We identify Lrf as a direct repressor of Bim transcription. In strong support of this mechanism, genetic Bim loss delays the lethality of Lrf-deficient embryos and rescues their anemia phenotype. Thus, our data define a key transcriptional cascade for effective erythropoiesis, whereby GATA-1 suppresses BIM-mediated apoptosis via LRF.

Transcription factor site dependencies in human, mouse and rat genomes

BACKGROUND

It is known that transcription factors frequently act together to regulate gene expression in eukaryotes. In this paper we describe a computational analysis of transcription factor site dependencies in human, mouse and rat genomes.

RESULTS

Our approach for quantifying tendencies of transcription factor binding sites to co-occur is based on a binding site scoring function which incorporates dependencies between positions, the use of information about the structural class of each transcription factor (major/minor groove binder), and also considered the possible implications of varying GC content of the sequences. Significant tendencies (dependencies) have been detected by non-parametric statistical methodology (permutation tests). Evaluation of obtained results has been performed in several ways: reports from literature (many of the significant dependencies between transcription factors have previously been confirmed experimentally); dependencies between transcription factors are not biased due to similarities in their DNA-binding sites; the number of dependent transcription factors that belong to the same functional and structural class is significantly higher than would be expected by chance; supporting evidence from GO clustering of targeting genes. Based on dependencies between two transcription factor binding sites (second-order dependencies), it is possible to construct higher-order dependencies (networks). Moreover results about transcription factor binding sites dependencies can be used for prediction of groups of dependent transcription factors on a given promoter sequence. Our results, as well as a scanning tool for predicting groups of dependent transcription factors binding sites are available on the Internet.

CONCLUSION

We show that the computational analysis of transcription factor site dependencies is a valuable complement to experimental approaches for discovering transcription regulatory interactions and networks. Scanning promoter sequences with dependent groups of transcription factor binding sites improve the quality of transcription factor predictions.

Chromatin architecture and transcription factor binding regulate expression of erythrocyte membrane protein genes

Erythrocyte membrane protein genes serve as excellent models of complex gene locus structure and function, but their study has been complicated by both their large size and their complexity. To begin to understand the intricate interplay of transcription, dynamic chromatin architecture, transcription factor binding, and genomic organization in regulation of erythrocyte membrane protein genes, we performed chromatin immunoprecipitation (ChIP) coupled with microarray analysis and ChIP coupled with massively parallel DNA sequencing in both erythroid and nonerythroid cells. Unexpectedly, most regions of GATA-1 and NF-E2 binding were remote from gene promoters and transcriptional start sites, located primarily in introns. Cooccupancy with FOG-1, SCL, and MTA-2 was found at all regions of GATA-1 binding, with cooccupancy of SCL and MTA-2 also found at regions of NF-E2 binding. Cooccupancy of GATA-1 and NF-E2 was found frequently. A common signature of histone H3 trimethylation at lysine 4, GATA-1, NF-E2, FOG-1, SCL, and MTA-2 binding and consensus GATA-1-E-box binding motifs located 34 to 90 bp away from NF-E2 binding motifs was found frequently in erythroid cell-expressed genes. These results provide insights into our understanding of membrane protein gene regulation in erythropoiesis and the regulation of complex genetic loci in erythroid and nonerythroid cells and identify numerous candidate regions for mutations associated with membrane-linked hemolytic anemia.

GATA-2 reinforces megakaryocyte development in the absence of GATA-1

GATA-2 is an essential transcription factor that regulates multiple aspects of hematopoiesis. Dysregulation of GATA-2 is a hallmark of acute megakaryoblastic leukemia in children with Down syndrome, a malignancy that is defined by the combination of trisomy 21 and a GATA1 mutation. Here, we show that GATA-2 is required for normal megakaryocyte development as well as aberrant megakaryopoiesis in Gata1 mutant cells. Furthermore, we demonstrate that GATA-2 indirectly controls cell cycle progression in GATA-1-deficient megakaryocytes. Genome-wide microarray analysis and chromatin immunoprecipitation studies revealed that GATA-2 regulates a wide set of genes, including cell cycle regulators and megakaryocyte-specific genes. Surprisingly, GATA-2 also negatively regulates the expression of crucial myeloid transcription factors, such as Sfpi1 and Cebpa. In the absence of GATA-1, GATA-2 prevents induction of a latent myeloid gene expression program. Thus, GATA-2 contributes to cell cycle progression and the maintenance of megakaryocyte identity of GATA-1-deficient cells, including GATA-1s-expressing fetal megakaryocyte progenitors. Moreover, our data reveal that overexpression of GATA-2 facilitates aberrant megakaryopoiesis.

Graded repression of PU.1/Sfpi1 gene transcription by GATA factors regulates hematopoietic cell fate

GATA-1 and PU.1 are essential hematopoietic transcription factors that control erythromegakaryocytic and myelolymphoid differentiation, respectively. These proteins antagonize each other through direct physical interaction to repress alternate lineage programs. We used immortalized Gata1(-) erythromegakaryocytic progenitor cells to study how PU.1/Sfpi1 expression is regulated by GATA-1 and GATA-2, a related factor that is normally expressed at earlier stages of hematopoiesis. Both GATA factors bind the PU.1/Sfpi1 gene at 2 highly conserved regions. In the absence of GATA-1, GATA-2 binding is associated with an undifferentiated state, intermediate level PU.1/Sfpi1 expression, and low-level expression of its downstream myeloid target genes. Restoration of GATA-1 function induces erythromegakaryocytic differentiation. Concomitantly, GATA-1 replaces GATA-2 at the PU.1/Sfpi1 locus and PU.1/Sfpi1 expression is extinguished. In contrast, when GATA-1 is not present, shRNA knockdown of GATA-2 increases PU.1/Sfpi1 expression by 3-fold and reprograms the cells to become macrophages. Our findings indicate that GATA factors act sequentially to regulate lineage determination during hematopoiesis, in part by exerting variable repressive effects at the PU.1/Sfpi1 locus.

EKLF/KLF1 controls cell cycle entry via direct regulation of E2f2

Differentiation of erythroid cells requires precise control over the cell cycle to regulate the balance between cell proliferation and differentiation. The zinc finger transcription factor, erythroid Krüppel-like factor (EKLF/KLF1), is essential for proper erythroid cell differentiation and regulates many erythroid genes. Here we show that loss of EKLF leads to aberrant entry into S-phase of the cell cycle during both primitive and definitive erythropoiesis. This cell cycle defect was associated with a significant reduction in the expression levels of E2f2 and E2f4, key factors necessary for the induction of S-phase gene expression and erythropoiesis. We found and validated novel intronic enhancers in both the E2f2 and E2f4 genes, which contain conserved CACC, GATA, and E-BOX elements. The E2f2 enhancer was occupied by EKLF in vivo. Furthermore, we were able to partially restore cell cycle dynamics in EKLF(-/-) fetal liver upon additional genetic depletion of Rb, establishing a genetic causal link between reduced E2f2 and the EKLF cell cycle defect. Finally, we propose direct regulation of the E2f2 enhancer is a generic mechanism by which many KLFs regulate proliferation and differentiation.

Expression of the leukemia oncogene Lmo2 is controlled by an array of tissue-specific elements dispersed over 100 kb and bound by Tal1/Lmo2, Ets, and Gata factors

The Lmo2 gene encodes a transcriptional cofactor critical for the development of hematopoietic stem cells. Ectopic LMO2 expression causes leukemia in T-cell acute lymphoblastic leukemia (T-ALL) patients and severe combined immunodeficiency patients undergoing retroviral gene therapy. Tightly controlled Lmo2 expression is therefore essential, yet no comprehensive analysis of Lmo2 regulation has been published so far. By comparative genomics, we identified 17 highly conserved noncoding elements, 9 of which revealed specific acetylation marks in chromatin-immunoprecipitation and microarray (ChIP-chip) assays performed across 250 kb of the Lmo2 locus in 11 cell types covering different stages of hematopoietic differentiation. All candidate regulatory regions were tested in transgenic mice. An extended LMO2 proximal promoter fragment displayed strong endothelial activity, while the distal promoter showed weak forebrain activity. Eight of the 15 distal candidate elements functioned as enhancers, which together recapitulated the full expression pattern of Lmo2, directing expression to endothelium, hematopoietic cells, tail, and forebrain. Interestingly, distinct combinations of specific distal regulatory elements were required to extend endothelial activity of the LMO2 promoter to yolk sac or fetal liver hematopoietic cells. Finally, Sfpi1/Pu.1, Fli1, Gata2, Tal1/Scl, and Lmo2 were shown to bind to and transactivate Lmo2 hematopoietic enhancers, thus identifying key upstream regulators and positioning Lmo2 within hematopoietic regulatory networks.

A new cis-acting regulatory element driving gene expression in the zebrafish pineal gland

MOTIVATION

The identification of functional cis-acting DNA regulatory elements is a crucial step towards understanding gene regulation. Ab initio motif detection algorithms have been extensively used in search of regulatory elements. Yet, their success in providing experimentally validated regulatory elements in vertebrates has been limited.

RESULTS

Here we report in silico identification and in vivo validation of regulatory elements that determine enhanced gene expression in the pineal gland of zebrafish. Microarray data enabled detection of genes that exhibit high expression in the pineal gland. The promoter regions of these genes were computationally analyzed in order to identify overrepresented motifs. The highest ranking motif identified is a CRX/OTX binding site, known to govern expression in the pineal gland and retina. The second highest ranking motif was not reported before; we experimentally validated its function in vivo by mutational analysis. The methodology presented here may be applicable as a general scheme for finding regulatory elements that contribute to tissue-specific gene expression.

Assaying the regulatory potential of mammalian conserved non-coding sequences in human cells

The fraction of experimentally active conserved non-coding sequences within any given cell type is low, so classical assays are unlikely to expose their potential.

Background

Conserved non-coding sequences in the human genome are approximately tenfold more abundant than known genes, and have been hypothesized to mark the locations of cis-regulatory elements. However, the global contribution of conserved non-coding sequences to the transcriptional regulation of human genes is currently unknown. Deeply conserved elements shared between humans and teleost fish predominantly flank genes active during morphogenesis and are enriched for positive transcriptional regulatory elements. However, such deeply conserved elements account for <1% of the conserved non-coding sequences in the human genome, which are predominantly mammalian.

Results

We explored the regulatory potential of a large sample of these 'common' conserved non-coding sequences using a variety of classic assays, including chromatin remodeling, and enhancer/repressor and promoter activity. When tested across diverse human model cell types, we find that the fraction of experimentally active conserved non-coding sequences within any given cell type is low (approximately 5%), and that this proportion increases only modestly when considered collectively across cell types.

Conclusions

The results suggest that classic assays of cis-regulatory potential are unlikely to expose the functional potential of the substantial majority of mammalian conserved non-coding sequences in the human genome.

Transcriptional enhancement by GATA1-occupied DNA segments is strongly associated with evolutionary constraint on the binding site motif

Tissue development and function are exquisitely dependent on proper regulation of gene expression, but it remains controversial whether the genomic signals controlling this process are subject to strong selective constraint. While some studies show that highly constrained noncoding regions act to enhance transcription, other studies show that DNA segments with biochemical signatures of regulatory regions, such as occupancy by a transcription factor, are seemingly unconstrained across mammalian evolution. To test the possible correlation of selective constraint with enhancer activity, we used chromatin immunoprecipitation as an approach unbiased by either evolutionary constraint or prior knowledge of regulatory activity to identify DNA segments within a 66-Mb region of mouse chromosome 7 that are occupied by the erythroid transcription factor GATA1. DNA segments bound by GATA1 were identified by hybridization to high-density tiling arrays, validated by quantitative PCR, and tested for gene regulatory activity in erythroid cells. Whereas almost all of the occupied segments contain canonical WGATAR binding site motifs for GATA1, in only 45% of the cases is the motif deeply preserved (found at the orthologous position in placental mammals or more distant species). However, GATA1-bound segments with high enhancer activity tend to be the ones with an evolutionarily preserved WGATAR motif, and this relationship was confirmed by a loss-of-function assay. Thus, GATA1 binding sites that regulate gene expression during erythroid maturation are under strong selective constraint, while nonconstrained binding may have only a limited or indirect role in regulation.

Evolution of the mammalian transcription factor binding repertoire via transposable elements

Identification of lineage-specific innovations in genomic control elements is critical for understanding transcriptional regulatory networks and phenotypic heterogeneity. We analyzed, from an evolutionary perspective, the binding regions of seven mammalian transcription factors (ESR1, TP53, MYC, RELA, POU5F1, SOX2, and CTCF) identified on a genome-wide scale by different chromatin immunoprecipitation approaches and found that only a minority of sites appear to be conserved at the sequence level. Instead, we uncovered a pervasive association with genomic repeats by showing that a large fraction of the bona fide binding sites for five of the seven transcription factors (ESR1, TP53, POU5F1, SOX2, and CTCF) are embedded in distinctive families of transposable elements. Using the age of the repeats, we established that these repeat-associated binding sites (RABS) have been associated with significant regulatory expansions throughout the mammalian phylogeny. We validated the functional significance of these RABS by showing that they are over-represented in proximity of regulated genes and that the binding motifs within these repeats have undergone evolutionary selection. Our results demonstrate that transcriptional regulatory networks are highly dynamic in eukaryotic genomes and that transposable elements play an important role in expanding the repertoire of binding sites.

Evolutionary conservation of zinc finger transcription factor binding sites in promoters of genes co-expressed with WT1 in prostate cancer

BACKGROUND

Gene expression analyses have led to a better understanding of growth control of prostate cancer cells. We and others have identified the presence of several zinc finger transcription factors in the neoplastic prostate, suggesting a potential role for these genes in the regulation of the prostate cancer transcriptome. One of the transcription factors (TFs) identified in the prostate cancer epithelial cells was the Wilms tumor gene (WT1). To rapidly identify coordinately expressed prostate cancer growth control genes that may be regulated by WT1, we used an in silico approach.

RESULTS

Evolutionary conserved transcription factor binding sites (TFBS) recognized by WT1, EGR1, SP1, SP2, AP2 and GATA1 were identified in the promoters of 24 differentially expressed prostate cancer genes from eight mammalian species. To test the relationship between sequence conservation and function, chromatin of LNCaP prostate cancer and kidney 293 cells were tested for TF binding using chromatin immunoprecipitation (ChIP). Multiple putative TFBS in gene promoters of placental mammals were found to be shared with those in human gene promoters and some were conserved between genomes that diverged about 170 million years ago (i.e., primates and marsupials), therefore implicating these sites as candidate binding sites. Among those genes coordinately expressed with WT1 was the kallikrein-related peptidase 3 (KLK3) gene commonly known as the prostate specific antigen (PSA) gene. This analysis located several potential WT1 TFBS in the PSA gene promoter and led to the rapid identification of a novel putative binding site confirmed in vivo by ChIP. Conversely for two prostate growth control genes, androgen receptor (AR) and vascular endothelial growth factor (VEGF), known to be transcriptionally regulated by WT1, regulatory sequence conservation was observed and TF binding in vivo was confirmed by ChIP.

CONCLUSION

Overall, this targeted approach rapidly identified important candidate WT1-binding elements in genes coordinately expressed with WT1 in prostate cancer cells, thus enabling a more focused functional analysis of the most likely target genes in prostate cancer progression. Identifying these genes will help to better understand how gene regulation is altered in these tumor cells.

PAZAR: a framework for collection and dissemination of cis-regulatory sequence annotation

PAZAR is an open-access and open-source database of transcription factor and regulatory sequence annotation with associated web interface and programming tools for data submission and extraction.

PAZAR is an open-access and open-source database of transcription factor and regulatory sequence annotation with associated web interface and programming tools for data submission and extraction. Curated boutique data collections can be maintained and disseminated through the unified schema of the mall-like PAZAR repository. The Pleiades Promoter Project collection of brain-linked regulatory sequences is introduced to demonstrate the depth of annotation possible within PAZAR. PAZAR, located at , is open for business.

Exchange of GATA factors mediates transitions in looped chromatin organization at a developmentally regulated gene locus

Enhancers can regulate designate promoters over long distances by forming chromatin loops. Whether chromatin loops are lost or reconfigured during gene repression is largely unexplored. We examined the chromosome conformation of the Kit gene that is expressed during early erythropoiesis but is downregulated upon cell maturation. Kit expression is controlled by sequential occupancy of two GATA family transcription factors. In immature cells, a distal enhancer bound by GATA-2 is in physical proximity with the active Kit promoter. Upon cell maturation, GATA-1 displaces GATA-2 and triggers a loss of the enhancer/promoter interaction. Moreover, GATA-1 reciprocally increases the proximity in nuclear space among distinct downstream GATA elements. GATA-1-induced transitions in chromatin conformation are not simply the consequence of transcription inhibition and require the cofactor FOG-1. This work shows that a GATA factor exchange reconfigures higher-order chromatin organization, and suggests that de novo chromatin loop formation is employed by nuclear factors to specify repressive outcomes.

A GATA-1-regulated microRNA locus essential for erythropoiesis

MicroRNAs (miRNAs) control tissue development, but their mechanism of regulation is not well understood. We used a gene complementation strategy combined with microarray screening to identify miRNAs involved in the formation of erythroid (red blood) cells. Two conserved miRNAs, miR 144 and miR 451, emerged as direct targets of the critical hematopoietic transcription factor GATA-1. In vivo, GATA-1 binds a distal upstream regulatory element to activate RNA polymerase II-mediated transcription of a single common precursor RNA (pri-miRNA) encoding both mature miRNAs. Zebrafish embryos depleted of miR 451 by using antisense morpholinos form erythroid precursors, but their development into mature circulating red blood cells is strongly and specifically impaired. These results reveal a miRNA locus that is required for erythropoiesis and uncover a new regulatory axis through which GATA-1 controls this process.

A single SNP in an evolutionary conserved region within intron 86 of the HERC2 gene determines human blue-brown eye color

We have previously demonstrated that haplotypes of three single nucleotide polymorphisms (SNPs) within the first intron of the OCA2 gene are extremely strongly associated with variation in human eye color. In the present work, we describe additional fine association mapping of eye color SNPs in the intergenic region upstream of OCA2 and within the neighboring HERC2 (hect domain and RLD2) gene. We screened an additional 92 SNPs in 300-3000 European individuals and found that a single SNP in intron 86 of HERC2, rs12913832, predicted eye color significantly better (ordinal logistic regression R(2) = 0.68, association LOD = 444) than our previous best OCA2 haplotype. Comparison of sequence alignments of multiple species showed that this SNP lies in the center of a short highly conserved sequence and that the blue-eye-associated allele (frequency 78%) breaks up this conserved sequence, part of which forms a consensus binding site for the helicase-like transcription factor (HLTF). We were also able to demonstrate the OCA2 R419Q, rs1800407, coding SNP acts as a penetrance modifier of this new HERC2 SNP for eye color, and somewhat independently, of melanoma risk. We conclude that the conserved region around rs12913832 represents a regulatory region controlling constitutive expression of OCA2 and that the C allele at rs12913832 leads to decreased expression of OCA2, particularly within iris melanocytes, which we postulate to be the ultimate cause of blue eye color.

ORegAnno: an open-access community-driven resource for regulatory annotation

ORegAnno is an open-source, open-access database and literature curation system for community-based annotation of experimentally identified DNA regulatory regions, transcription factor binding sites and regulatory variants. The current release comprises 30 145 records curated from 922 publications and describing regulatory sequences for over 3853 genes and 465 transcription factors from 19 species. A new feature called the 'publication queue' allows users to input relevant papers from scientific literature as targets for annotation. The queue contains 4438 gene regulation papers entered by experts and another 54 351 identified by text-mining methods. Users can enter or 'check out' papers from the queue for manual curation using a series of user-friendly annotation pages. A typical record entry consists of species, sequence type, sequence, target gene, binding factor, experimental outcome and one or more lines of experimental evidence. An evidence ontology was developed to describe and categorize these experiments. Records are cross-referenced to Ensembl or Entrez gene identifiers, PubMed and dbSNP and can be visualized in the Ensembl or UCSC genome browsers. All data are freely available through search pages, XML data dumps or web services at: http://www.oreganno.org.

Runx1-mediated hematopoietic stem-cell emergence is controlled by a Gata/Ets/SCL-regulated enhancer

The transcription factor Runx1/AML1 is an important regulator of hematopoiesis and is critically required for the generation of the first definitive hematopoietic stem cells (HSCs) in the major vasculature of the mouse embryo. As a pivotal factor in HSC ontogeny, its transcriptional regulation is of high interest but is largely undefined. In this study, we used a combination of comparative genomics and chromatin analysis to identify a highly conserved 531-bp enhancer located at position + 23.5 in the first intron of the 224-kb mouse Runx1 gene. We show that this enhancer contributes to the early hematopoietic expression of Runx1. Transcription factor binding in vivo and analysis of the mutated enhancer in transient transgenic mouse embryos implicate Gata2 and Ets proteins as critical factors for its function. We also show that the SCL/Lmo2/Ldb-1 complex is recruited to the enhancer in vivo. Importantly, transplantation experiments demonstrate that the intronic Runx1 enhancer targets all definitive HSCs in the mouse embryo, suggesting that it functions as a crucial cis-regulatory element that integrates the Gata, Ets, and SCL transcriptional networks to initiate HSC generation.

Computational identification and functional validation of regulatory motifs in cartilage-expressed genes

Chondrocyte gene regulation is important for the generation and maintenance of cartilage tissues. Several regulatory factors have been identified that play a role in chondrogenesis, including the positive transacting factors of the SOX family such as SOX9, SOX5, and SOX6, as well as negative transacting factors such as C/EBP and delta EF1. However, a complete understanding of the intricate regulatory network that governs the tissue-specific expression of cartilage genes is not yet available. We have taken a computational approach to identify cis-regulatory, transcription factor (TF) binding motifs in a set of cartilage characteristic genes to better define the transcriptional regulatory networks that regulate chondrogenesis. Our computational methods have identified several TFs, whose binding profiles are available in the TRANSFAC database, as important to chondrogenesis. In addition, a cartilage-specific SOX-binding profile was constructed and used to identify both known, and novel, functional paired SOX-binding motifs in chondrocyte genes. Using DNA pattern-recognition algorithms, we have also identified cis-regulatory elements for unknown TFs. We have validated our computational predictions through mutational analyses in cell transfection experiments. One novel regulatory motif, N1, found at high frequency in the COL2A1 promoter, was found to bind to chondrocyte nuclear proteins. Mutational analyses suggest that this motif binds a repressive factor that regulates basal levels of the COL2A1 promoter.

An unmethylated 3' promoter-proximal region is required for efficient transcription initiation

The promoter regions of approximately 40% of genes in the human genome are embedded in CpG islands, CpG-rich regions that frequently extend on the order of one kb 3′ of the transcription start site (TSS) region. CpGs 3′ of the TSS of actively transcribed CpG island promoters typically remain methylation-free, indicating that maintaining promoter-proximal CpGs in an unmethylated state may be important for efficient transcription. Here we utilize recombinase-mediated cassette exchange to introduce a Moloney Murine Leukemia Virus (MoMuLV)-based reporter, in vitro methylated 1 kb downstream of the TSS, into a defined genomic site. In a subset of clones, methylation spreads to within ∼320 bp of the TSS, yielding a dramatic decrease in transcript level, even though the promoter/TSS region remains unmethylated. Chromatin immunoprecipitation analyses reveal that such promoter-proximal methylation results in loss of RNA polymerase II and TATA-box-binding protein (TBP) binding in the promoter region, suggesting that repression occurs at the level of transcription initiation. While DNA methylation-dependent trimethylation of H3 lysine (K)9 is confined to the intragenic methylated region, the promoter and downstream regions are hypo-acetylated on H3K9/K14. Furthermore, DNase I hypersensitivity and methylase-based single promoter analysis (M-SPA) experiments reveal that a nucleosome is positioned over the unmethylated TATA-box in these clones, indicating that dense DNA methylation downstream of the promoter region is sufficient to alter the chromatin structure of an unmethylated promoter. Based on these observations, we propose that a DNA methylation-free region extending several hundred bases downstream of the TSS may be a prerequisite for efficient transcription initiation. This model provides a biochemical explanation for the typical positioning of TSSs well upstream of the 3′ end of the CpG islands in which they are embedded.

Genes, the functional units of heredity, are made up of DNA, which is packaged inside the nuclei of eukaryotic cells in association with a number of proteins in a structure called chromatin. In order for transcription, the process of transferring genetic information from DNA to RNA, to take place, chromatin must be decondensed to allow the transcription machinery to bind the genes that are to be transcribed. In mammals, promoters, the starting position of genes, are frequently embedded in “CpG islands,” regions with a relatively high density of the CpG dinucleotide. Paradoxically, while cytosines in the context of the CpG dinucleotide are generally methylated, CpGs flanking the start sites of genes typically remain methylation-free. As CpG methylation is associated with condensed chromatin, it is generally believed that promoter regions must remain free of methylation to allow for binding of the transcription machinery. Here, using a novel method for introducing methylated DNA into a defined genomic site, we demonstrate that DNA methylation in the promoter-proximal region of a gene is sufficient to block transcription via the generation of a chromatin structure that inhibits binding of the transcription machinery. Thus, methylation may inhibit transcription even when present outside the promoter region.

ESPERR: learning strong and weak signals in genomic sequence alignments to identify functional elements

Genomic sequence signals - such as base composition, presence of particular motifs, or evolutionary constraint - have been used effectively to identify functional elements. However, approaches based only on specific signals known to correlate with function can be quite limiting. When training data are available, application of computational learning algorithms to multispecies alignments has the potential to capture broader and more informative sequence and evolutionary patterns that better characterize a class of elements. However, effective exploitation of patterns in multispecies alignments is impeded by the vast number of possible alignment columns and by a limited understanding of which particular strings of columns may characterize a given class. We have developed a computational method, called ESPERR (evolutionary and sequence pattern extraction through reduced representations), which uses training examples to learn encodings of multispecies alignments into reduced forms tailored for the prediction of chosen classes of functional elements. ESPERR produces a greatly improved Regulatory Potential score, which can discriminate regulatory regions from neutral sites with excellent accuracy ( approximately 94%). This score captures strong signals (GC content and conservation), as well as subtler signals (with small contributions from many different alignment patterns) that characterize the regulatory elements in our training set. ESPERR is also effective for predicting other classes of functional elements, as we show for DNaseI hypersensitive sites and highly conserved regions with developmental enhancer activity. Our software, training data, and genome-wide predictions are available from our Web site (http://www.bx.psu.edu/projects/esperr).