Genomic Binding Sites Research Articles

Long Non-Coding RNAs Control Hematopoietic Stem Cell Function

Hematopoietic stem cells (HSCs) possess unique gene expression programs that enforce their identity and regulate lineage commitment. Long non-coding RNAs (lncRNAs) have emerged as important regulators of gene expression and cell fate decisions, although their functions in HSCs are unclear. Here we profiled the transcriptome of purified HSCs by deep sequencing and identified 323 unannotated lncRNAs. Comparing their expression in differentiated lineages revealed 159 lncRNAs enriched in HSCs, some of which are likely HSC specific (LncHSCs). These lncRNA genes share epigenetic features with protein-coding genes, including regulated expression via DNA methylation, and knocking down two LncHSCs revealed distinct effects on HSC self-renewal and lineage commitment. We mapped the genomic binding sites of one of these candidates and found enrichment for key hematopoietic transcription factor binding sites, especially E2A. Together, these results demonstrate that lncRNAs play important roles in regulating HSCs, providing an additional layer to the genetic circuitry controlling HSC function.

Wanda Blenska

<h3>ABSTRACT</h3> Lineage-defining transcription factors (LTFs) play key roles in tumor cell growth, making them highly attractive, but currently “undruggable”, small cell lung cancer (SCLC) vulnerabilities. Delineating LTF genomic binding sites...

Locus-Specific Targeting to the X Chromosome Revealed by the RNA Interactome of CTCF

CTCF is a master regulator that plays important roles in genome architecture and gene expression. How CTCF is recruited in a locus-specific manner is notfully understood. Evidence from epigenetic processes, such as X chromosome inactivation (XCI), indicates that CTCF associates functionally with RNA. Using genome-wide approaches to investigate the relationship between its RNA interactome and epigenomic landscape, here we report that CTCF binds thousands of transcripts in mouse embryonic stem cells, many in close proximity to CTCF's genomic binding sites. CTCF is a specific and high-affinity RNA-binding protein (Kd< 1nM). During XCI, CTCF differentially binds the active and inactive X chromosomes and interacts directly with Tsix, Xite, and Xist RNAs. Tsix and Xite RNAs target CTCF to the X inactivation center, thereby inducing homologous X chromosome pairing. Our work elucidates one mechanism by which CTCF is recruited in a locus-specific manner and implicates CTCF-RNA interactions in long-range chromosomal interactions.

Identification of CAR/RXR&amp;lt;i&amp;gt;α&amp;lt;/i&amp;gt; Hetero-dimer Binding Sites in the Human Genome by a Modified Yeast One-Hybrid Assay

The constitutive androstane receptor (CAR) is a transcription factor that belongs to the nuclear receptor superfamily. CAR binds as a heterodimer with the retinoid X receptor α (RXRα) to CAR response elements (CAREs) and regulates the expression of various drug metabolizing enzymes and transporters. To identify CAR/RXRα binding sites in the human genome, we performed a modified yeast one-hybrid assay that enables rapid and efficient identification of genomic targets for DNA-binding proteins. DNA fragments were recovered from positive yeast colonies by PCR and sequenced. A motif enrichment analysis revealed that the most frequent motif was a direct repeat (DR) of RGKTCA-like core sequence spaced by 4 bp. Next, we predicted 149 putative CAR/RXRα binding sites from 414 unique clones, by searching for DRs, everted repeats (ERs) and inverted repeats (IRs) of the RGKTCA-like core motif. Based on gel mobility shift assays, the CAR/RXRα heterodimer could directly interact with the 108 predicted sequences, which included not only classical CAREs but also a wide variety of arrangements. Furthermore, we identified 17 regulatory polymorphisms on the CAR/RXRα-binding sites that may influence individual variation in the expression of CAR-regulated genes. These results provide insights into the molecular mechanisms underlying the physiological and pathological actions of CAR/RXRα het-erodimers.

Fast masking of repeated primer binding sites in eukaryotic genomes.

In this article we describe the working principle and a list of practical applications for GenomeMasker-a program that finds and masks all repeated DNA motifs in fully sequenced genomes. The GenomeMasker exhaustively finds and masks all repeated DNA motifs in studied genomes. The software is optimized for PCR primer design. The algorithm is designed for high-throughput work, allowing masking of large DNA regions, even entire eukaryotic genomes. Additionally, the software is able to predict all alternative PCR products from studied genomes for thousands of candidate PCR primer pairs. Practical applications of the GenomeMasker are shown for command-line version of the GenomeMasker, which can be downloaded from http://bioinfo.ut.ee/download/. Graphical Web interfaces with limited options are available at http://bioinfo.ut.ee/genometester/ and http://bioinfo.ut.ee/snpmasker/.

Identification and Characterization of Human Genomic Binding Sites for Retinoic Acid Receptor/Retinoid X Receptor Heterodimers

All-trans retinoic acid (ATRA) triggers a wide range of effects on vertebrate development by regulating cell proliferation, differentiation, and apoptosis. ATRA activates retinoic acid receptors (RARs) which heterodimerize with retinoid X receptors (RXRs). RAR/RXR heterodimers function as ATRA-dependent transcriptional regulators by binding to retinoic acid response elements (RAREs). To identify RAR/RXR heterodimer-binding sites in the human genome, we performed a modified yeast one-hybrid assays and identified 193 RAR/RXR heterodimer-binding fragments in the human genome. The putative target genes included genes involved in development process and cell differentiation. Gel mobility shift assays indicated that 160 putative RAREs could directly interact with the RAR/RXR heterodimer. Moreover, 19 functional regulatory single nucleotide polymorphisms (rSNPs) on the RAR/RXR-binding sequences were identified by analyzing the difference in the DNA-binding affinities. These results provide insights into the molecular mechanisms underlying the physiological and pathological actions of RAR/RXR heterodimers.

A variably occupied CTCF binding site in the ultrabithorax gene in the Drosophila bithorax complex.

Although the majority of genomic binding sites for the insulator protein CCCTC-binding factor (CTCF) are constitutively occupied, a subset show variable occupancy. Such variable sites provide an opportunity to assess context-specific CTCF functions in gene regulation. Here, we have identified a variably occupied CTCF site in the Drosophila Ultrabithorax (Ubx) gene. This site is occupied in tissues where Ubx is active (third thoracic leg imaginal disc) but is not bound in tissues where the Ubx gene is repressed (first thoracic leg imaginal disc). Using chromatin conformation capture, we show that this site preferentially interacts with the Ubx promoter region in the active state. The site lies close to Ubx enhancer elements and is also close to the locations of several gypsy transposon insertions that disrupt Ubx expression, leading to the bx mutant phenotype. gypsy insertions carry the Su(Hw)-dependent gypsy insulator and were found to affect both CTCF binding at the variable site and the chromatin topology. This suggests that insertion of the gypsy insulator in this region interferes with CTCF function and supports a model for the normal function of the variable CTCF site as a chromatin loop facilitator, promoting interaction between Ubx enhancers and the Ubx transcription start site.

Vitamin D receptor and RXR in the post-genomic era.

Following the elucidation of the human genome and components of the epigenome, it is timely to revisit what is known of vitamin D receptor (VDR) function. Early transcriptomic studies using microarray approaches focused on the protein coding mRNA that were regulated by the VDR, usually following treatment with ligand. These studies quickly established the approximate size and surprising diversity of the VDR transcriptome, revealing it to be highly heterogenous and cell type and time dependent. Investigators also considered VDR regulation of non-protein coding RNA and again, cell and time dependency was observed. Attempts to integrate mRNA and miRNA regulation patterns are beginning to reveal patterns of co-regulation and interaction that allow for greater control of mRNA expression, and the capacity to govern more complex cellular events. Alternative splicing in the trasncriptome has emerged as a critical process in transcriptional control and there is evidence of the VDR interacting with components of the splicesome. ChIP-Seq approaches have proved to be pivotal to reveal the diversity of the VDR binding choices across cell types and following treatment, and have revealed that the majority of these are non-canonical in nature. The underlying causes driving the diversity of VDR binding choices remain enigmatic. Finally, genetic variation has emerged as important to impact the transcription factor affinity towards genomic binding sites, and recently the impact of this on VDR function has begun to be considered.

MBD2 and MBD3: elusive functions and mechanisms.

Deoxyribonucleic acid methylation is a long known epigenetic mark involved in many biological processes and the ‘readers’ of this mark belong to several distinct protein families that ‘read’ and ‘translate’ the methylation mark into a function. Methyl-CpG binding domain proteins belong to one of these families that are associated with transcriptional activation/repression, regulation of chromatin structure, pluripotency, development, and differentiation. Discovered decades ago, the systematic determination of the genomic binding sites of these readers and their epigenome make-up at a genome-wide level revealed the tip of the functional iceberg. This review focuses on two members of the methyl binding proteins, namely MBD2 and MBD3 that reside in very similar complexes, yet appear to have very different biological roles. We provide a comprehensive comparison of their genome-wide binding features and emerging roles in gene regulation.

In situ dissection of RNA functional subunits by domain-specific chromatin isolation by RNA purification (dChIRP).

Here we describe domain-specific chromatin isolation by RNA purification (dChIRP), a technique for dissecting the functional domains of a target RNA in situ. For an RNA of interest, dChIRP can identify domain-level intramolecular and intermolecular RNA-RNA, RNA-protein, and RNA-DNA interactions and maps the RNA's genomic binding sites with higher precision than domain-agnostic methods. We illustrate how this technique has been applied to the roX1 lncRNA to resolve its domain-level architecture, discover its protein- and chromatin-interacting domains, and map its occupancy on the X chromosome.

Cell context dependent p53 genome-wide binding patterns and enrichment at repeats.

The p53 ability to elicit stress specific and cell type specific responses is well recognized, but how that specificity is established remains to be defined. Whether upon activation p53 binds to its genomic targets in a cell type and stress type dependent manner is still an open question. Here we show that the p53 binding to the human genome is selective and cell context-dependent. We mapped the genomic binding sites for the endogenous wild type p53 protein in the human cancer cell line HCT116 and compared them to those we previously determined in the normal cell line IMR90. We report distinct p53 genome-wide binding landscapes in two different cell lines, analyzed under the same treatment and experimental conditions, using the same ChIP-seq approach. This is evidence for cell context dependent p53 genomic binding. The observed differences affect the p53 binding sites distribution with respect to major genomic and epigenomic elements (promoter regions, CpG islands and repeats). We correlated the high-confidence p53 ChIP-seq peaks positions with the annotated human repeats (UCSC Human Genome Browser) and observed both common and cell line specific trends. In HCT116, the p53 binding was specifically enriched at LINE repeats, compared to IMR90 cells. The p53 genome-wide binding patterns in HCT116 and IMR90 likely reflect the different epigenetic landscapes in these two cell lines, resulting from cancer-associated changes (accumulated in HCT116) superimposed on tissue specific differences (HCT116 has epithelial, while IMR90 has mesenchymal origin). Our data support the model for p53 binding to the human genome in a highly selective manner, mobilizing distinct sets of genes, contributing to distinct pathways.

Two new insulator proteins, Pita and ZIPIC, target CP190 to chromatin.

Insulators are multiprotein–DNA complexes that regulate the nuclear architecture. The Drosophila CP190 protein is a cofactor for the DNA-binding insulator proteins Su(Hw), CTCF, and BEAF-32. The fact that CP190 has been found at genomic sites devoid of either of the known insulator factors has until now been unexplained. We have identified two DNA-binding zinc-finger proteins, Pita, and a new factor named ZIPIC, that interact with CP190 in vivo and in vitro at specific interaction domains. Genomic binding sites for these proteins are clustered with CP190 as well as with CTCF and BEAF-32. Model binding sites for Pita or ZIPIC demonstrate a partial enhancer-blocking activity and protect gene expression from PRE-mediated silencing. The function of the CTCF-bound MCP insulator sequence requires binding of Pita. These results identify two new insulator proteins and emphasize the unifying function of CP190, which can be recruited by many DNA-binding insulator proteins.

Dynamic shifts in occupancy by TAL1 are guided by GATA factors and drive large-scale reprogramming of gene expression during hematopoiesis

We used mouse ENCODE data along with complementary data from other laboratories to study the dynamics of occupancy and the role in gene regulation of the transcription factor TAL1, a critical regulator of hematopoiesis, at multiple stages of hematopoietic differentiation. We combined ChIP-seq and RNA-seq data in six mouse cell types representing a progression from multilineage precursors to differentiated erythroblasts and megakaryocytes. We found that sites of occupancy shift dramatically during commitment to the erythroid lineage, vary further during terminal maturation, and are strongly associated with changes in gene expression. In multilineage progenitors, the likely target genes are enriched for hematopoietic growth and functions associated with the mature cells of specific daughter lineages (such as megakaryocytes). In contrast, target genes in erythroblasts are specifically enriched for red cell functions. Furthermore, shifts in TAL1 occupancy during erythroid differentiation are associated with gene repression (dissociation) and induction (co-occupancy with GATA1). Based on both enrichment for transcription factor binding site motifs and co-occupancy determined by ChIP-seq, recruitment by GATA transcription factors appears to be a stronger determinant of TAL1 binding to chromatin than the canonical E-box binding site motif. Studies of additional proteins lead to the model that TAL1 regulates expression after being directed to a distinct subset of genomic binding sites in each cell type via its association with different complexes containing master regulators such as GATA2, ERG, and RUNX1 in multilineage cells and the lineage-specific master regulator GATA1 in erythroblasts.

Identification of Genomic Binding Sites for Candida glabrata Pdr1 Transcription Factor in Wild-Type and ρ 0 Cells

The fungal pathogen Candida glabrata is an emerging cause of candidiasis in part owing to its robust ability to acquire tolerance to the major clinical antifungal drug fluconazole. Similar to the related species Candida albicans, C. glabrata most typically gains azole tolerance via transcriptional induction of a suite of resistance genes, including a locus encoding an ABCG-type ATP-binding cassette (ABC) transporter that is referred to as CDR1 in Candida species. In C. glabrata, CDR1 expression is controlled primarily by the activity of a transcriptional activator protein called Pdr1. Strains exhibiting reduced azole susceptibility often contain substitution mutations in PDR1 that in turn lead to elevated mRNA levels of target genes with associated azole resistance. Pdr1 activity is also induced upon loss of the mitochondrial genome status and upon challenge by azole drugs. While extensive analyses of the transcriptional effects of Pdr1 have identified a number of genes that are regulated by this factor, we cannot yet separate direct from indirect target genes. Here we used chromatin immunoprecipitation (ChIP) coupled with high-throughput sequencing (ChIP-seq) to identify the promoters and associated genes directly regulated by Pdr1. These genes include many that are shared with the yeast Saccharomyces cerevisiae but others that are unique to C. glabrata, including the ABC transporter-encoding locus YBT1, genes involved in DNA repair, and several others. These data provide the outline for understanding the primary response genes involved in production of Pdr1-dependent azole resistance in C. glabrata.

The Long Noncoding RNAs NEAT1 and MALAT1 Bind Active Chromatin Sites

Mechanistic roles for many lncRNAs are poorly understood, in part because their direct interactions with genomic loci and proteins are difficult to assess. Using a method to purify endogenous RNAs and their associated factors, we mapped the genomic binding sites for two highly expressed human lncRNAs, NEAT1 and MALAT1. We show that NEAT1 and MALAT1 localize to hundreds of genomic sites in human cells, primarily over active genes. NEAT1 and MALAT1 exhibit colocalization to many of these loci, but display distinct gene body binding patterns at these sites, suggesting independent but complementary functions for these RNAs. We also identified numerous proteins enriched by both lncRNAs, supporting complementary binding and function, in addition to unique associated proteins. Transcriptional inhibition or stimulation alters localization of NEAT1 on active chromatin sites, implying that underlying DNA sequence does not target NEAT1 to chromatin, and that localization responds to cues involved in the transcription process.

Promoter-proximal transcription factor binding is transcriptionally active when coupled with nucleosome repositioning in immediate vicinity.

Previous studies have analyzed patterns of transcription, transcription factor (TF) binding or mapped nucleosome occupancy across the genome. These suggest that the three aspects are genetically connected but the cause and effect relationships are still unknown. For example, physiologic TF binding studies involve many TFs, consequently, it is difficult to assign nucleosome reorganization to the binding site occupancy of any particular TF. Therefore, several aspects remain unclear: does TF binding influence nucleosome (re)organizations locally or impact the chromatin landscape at a more global level; are all or only a fraction of TF binding a result of reorganization in nucleosome occupancy and do all TF binding and associated changes in nucleosome occupancy result in altered gene expression? With these in mind, following characterization of two states (before and after induction of a single TF of choice) we determined: (i) genomic binding sites of the TF, (ii) promoter nucleosome occupancy and (iii) transcriptome profiles. Results demonstrated that promoter-proximal TF binding influenced expression of the target gene when it was coupled to nucleosome repositioning at or close to its binding site in most cases. In contrast, only in few cases change in target gene expression was found when TF binding occurred without local nucleosome reorganization.

BPeaks: a bioinformatics tool to detect transcription factor binding sites from ChIPseq data in yeasts and other organisms with small genomes.

Peak calling is a critical step in ChIPseq data analysis. Choosing the correct algorithm as well as optimized parameters for a specific biological system is an essential task. In this article, we present an original peak-calling method (bPeaks) specifically designed to detect transcription factor (TF) binding sites in small eukaryotic genomes, such as in yeasts. As TF interactions with DNA are strong and generate high binding signals, bPeaks uses simple parameters to compare the sequences (reads) obtained from the immunoprecipitation (IP) with those from the control DNA (input). Because yeasts have small genomes (<20 Mb), our program has the advantage of using ChIPseq information at the single nucleotide level and can explore, in a reasonable computational time, results obtained with different sets of parameter values. Graphical outputs and text files are provided to rapidly assess the relevance of the detected peaks. Taking advantage of the simple promoter structure in yeasts, additional functions were implemented in bPeaks to automatically assign the peaks to promoter regions and retrieve peak coordinates on the DNA sequence for further predictions of regulatory motifs, enriched in the list of peaks. Applications of the bPeaks program to three different ChIPseq datasets from Saccharomyces cerevisiae, Candida albicans and Candida glabrata are presented. Each time, bPeaks allowed us to correctly predict the DNA binding sequence of the studied TF and provided relevant lists of peaks. The bioinformatics tool bPeaks is freely distributed to academic users. Supplementary data, together with detailed tutorials, are available online: http://bpeaks.gene-networks.net.

Regulatory crosstalk between lineage-survival oncogenes KLF5, GATA4 and GATA6 cooperatively promotes gastric cancer development

ObjectiveGastric cancer (GC) is a deadly malignancy for which new therapeutic strategies are needed. Three transcription factors, KLF5, GATA4 and GATA6, have been previously reported to exhibit genomic amplification in...

The active enhancer network operated by liganded RXR supports angiogenic activity in macrophages.

RXR signaling is predicted to have a major impact in macrophages, but neither the biological consequence nor the genomic basis of its ligand activation is known. Comprehensive genome-wide studies were carried out to map liganded RXR-mediated transcriptional changes, active binding sites, and cistromic interactions in the context of the macrophage genome architecture. The macrophage RXR cistrome has 5200 genomic binding sites, which are not impacted by ligand. Active enhancers are characterized by PU.1 binding, an increase of enhancer RNA, and P300 recruitment. Using these features, 387 liganded RXR-bound enhancers were linked to 226 genes, which predominantly reside in CTCF/cohesin-limited functional domains. These findings were molecularly validated using chromosome conformation capture (3C) and 3C combined with sequencing (3C-seq), and we show that selected long-range enhancers communicate with promoters via stable or RXR-induced loops and that some of the enhancers interact with each other, forming an interchromosomal network. A set of angiogenic genes, including Vegfa, has liganded RXR-controlled enhancers and provides the macrophage with a novel inducible program.

Biophysical fitness landscapes for transcription factor binding sites.

Phenotypic states and evolutionary trajectories available to cell populations are ultimately dictated by complex interactions among DNA, RNA, proteins, and other molecular species. Here we study how evolution of gene regulation in a single-cell eukaryote S. cerevisiae is affected by interactions between transcription factors (TFs) and their cognate DNA sites. Our study is informed by a comprehensive collection of genomic binding sites and high-throughput in vitro measurements of TF-DNA binding interactions. Using an evolutionary model for monomorphic populations evolving on a fitness landscape, we infer fitness as a function of TF-DNA binding to show that the shape of the inferred fitness functions is in broad agreement with a simple functional form inspired by a thermodynamic model of two-state TF-DNA binding. However, the effective parameters of the model are not always consistent with physical values, indicating selection pressures beyond the biophysical constraints imposed by TF-DNA interactions. We find little statistical support for the fitness landscape in which each position in the binding site evolves independently, indicating that epistasis is common in the evolution of gene regulation. Finally, by correlating TF-DNA binding energies with biological properties of the sites or the genes they regulate, we are able to rule out several scenarios of site-specific selection, under which binding sites of the same TF would experience different selection pressures depending on their position in the genome. These findings support the existence of universal fitness landscapes which shape evolution of all sites for a given TF, and whose properties are determined in part by the physics of protein-DNA interactions.