Genomic approaches that aid in the identification of transcription factor target genes

Chromatin Immunoprecipitation (ChIP) to Study DNA-Protein Interactions

Chromatin immunoprecipitation (ChIP) is a method used to examine the genomic localization of a target of interest (e.g., proteins, protein posttranslational modifications, or DNA elements). As ChIP provides a snapshot of in vivo DNA-protein interactions, it lends insight to the mechanisms of gene expression and genome regulation. This chapter provides a detailed protocol focused on native-ChIP (N-ChIP), a robust approach to profile stable DNA-protein interactions. We also describe best practices for ChIP , including defined controls to ensure specific and efficient target enrichment and methods for data normalization.

Chromatin immunoprecipitation assays: analyzing transcription factor binding and histone modifications in vivo

Studies in the past decade have shown that differential gene expression depends not only on the binding of specific transcription factors to discrete promoter elements but also on the epigenetic modification of the DNA as well as histones associated with the promoter. While techniques like electrophoretic mobility shift assays could detect and characterize the binding of specific transcription factors present in cell lysates to DNA sequences in in vitro binding conditions, they were not effective in assessing the binding in intact cells. Development of chromatin immunoprecipitation technique in the past decade enabled the analysis of the association of regulatory molecules with specific promoters or changes in histone modifications in vivo, without overexpressing any component. ChIP assays can provide a snapshot of how a regulatory transcription factor affects the expression of a single gene, or a variety of genes at the same time. Availability of high quality antibodies that recognizes histones modified in a specific fashion further expanded the use of ChIP assays to analyze even minute changes in histone modification and nucleosomes structure. This chapter outlines the general strategies and protocols used to carry out ChIP assays to study the differential recruitment of transcription factors as well as histone modifications.

ChIP on chip and ChIP-Seq assays: genome-wide analysis of transcription factor binding and histone modifications

Deregulation of transcriptional activity of many genes has been causatively linked to human diseases including cancer. Altered patterns of gene expression in normal and cancer cells are the result of inappropriate expression of transcription factors and chromatin modifying proteins. Chromatin immunoprecipitation assay is a well-established tool for investigating the interactions between regulatory proteins and DNA at distinct stages of gene activation. ChIP coupled with DNA microarrays, known as ChIP on chip, or sequencing of DNA associated with the factors (ChIP-Seq) allow us to determine the entire spectrum of in vivo DNA binding sites for a given protein. This has been of immense value because ChIP on chip assays and ChIP-Seq experiments can provide a snapshot of the transcriptional regulatory mechanisms on a genome-wide scale. This chapter outlines the general strategies used to carry out ChIP-chip assays to study the differential recruitment of regulatory molecules based on the studies conducted in our lab as well as other published protocols; these can be easily modified to a ChIP-Seq analysis.

Fundamentals of vitamin D hormone-regulated gene expression

Initial research focused upon several known genetic targets provided early insight into the mechanism of action of the vitamin D hormone (1,25-dihydroxyvitamin D3 (1,25(OH)2D3)). Recently, however, a series of technical advances involving the coupling of chromatin immunoprecipitation (ChIP) to unbiased methodologies that initially involved tiled DNA microarrays (ChIP-chip analysis) and now Next Generation DNA Sequencing techniques (ChIP-seq analysis) has opened new avenues of research into the mechanisms through which 1,25(OH)2D3 regulates gene expression. In this review, we summarize briefly the results of this early work and then focus on more recent studies in which ChIP-chip and ChIP-seq analyses have been used to explore the mechanisms of 1,25(OH)2D3 action on a genome-wide scale providing specific target genes as examples. The results of this work have advanced our understanding of the mechanisms involved at both genetic and epigenetic levels and have revealed a series of new principles through which the vitamin D hormone functions to control the expression of genes. This article is part of a Special Issue entitled '16th Vitamin D Workshop'.

Regulation of gene expression: probing DNA–protein interactionsin vivoandin vitro

Tremendous efforts have been put together over the last several years to complete the entire sequencing of the human genome. As we enter the proteomic era, when the major aim is understanding which gene encodes which protein, the time has also come to identify their precise function inside the astonishing signaling network required to accomplish all cellular functions. Understanding when, why and how a gene is expressed has now become a necessity toward identifying all the regulatory pathways that mediate cellular processes such as differentiation, migration, replication, DNA repair and apoptosis. Regulation of gene transcription is a process that is primarily under the influence of nuclear-located transcription factors. Consequently, identifying which protein activates or represses a specific gene is a prerequisite for understanding cell fate and function. The current state of, and recent advances in, transcriptional regulation approaches are reviewed here, with special emphasis on new technologies required when probing for DNA-protein interactions. This review explores different strategies aimed at identifying both the regulatory sequences of any given gene and the trans-acting regulatory factors that recognize these elements as their target sites in the nucleus. Ongoing developments in the fields of nanotechnology, RNA silencing and protein modeling toward the investigation of DNA-protein interactions and their relevance in the battle against cancer are discussed.

Systematic genetic analysis of transcription factors to map the fission yeast transcription-regulatory network

Mapping transcriptional-regulatory networks requires the identification of target genes, binding specificities and signalling pathways of transcription factors. However, the characterization of each transcription factor sufficiently for deciphering such networks remains laborious. The recent availability of overexpression and deletion strains for almost all of the transcription factor genes in the fission yeast Schizosaccharomyces pombe provides a valuable resource to better investigate transcription factors using systematic genetics. In the present paper, I review and discuss the utility of these strain collections combined with transcriptome profiling and genome-wide chromatin immunoprecipitation to identify the target genes of transcription factors.

Isolation of intracellular protein--DNA complexes using HaloCHIP, an antibody-free alternative to chromatin immunoprecipitation

Mapping of protein binding sites within the genome has been significantly advanced by microarray and sequencing technologies, yet the method traditionally used to isolate protein-DNA complexes, chromatin immunoprecipitation, has remained dependent of the use of antibodies. Furthermore, cross-linking is commonly used to trap protein-DNA complexes and the challenge of using antibodies has come in recognition of the cross-linked epitopes, sometimes limiting the success of the approach. Here we present a method, HaloCHIP, which utilizes a HaloTag protein fusion and corresponding interaction resin, HaloLink, for capture of cross-linked protein-DNA complexes directly from a cellular lysate. This process alleviates the need for using an antibody, yields the DNA fragments bound to a particular protein of interest, and allows for a variety of downstream analyses such as PCR, qPCR, microarrays, and sequencing.

New methods for separating causes from effects in genomics data

Background

The discovery of molecular pathways is a challenging problem and its solution relies on the identification of causal molecular interactions in genomics data. Causal molecular interactions can be discovered using randomized experiments; however such experiments are often costly, infeasible, or unethical. Fortunately, algorithms that infer causal interactions from observational data have been in development for decades, predominantly in the quantitative sciences, and many of them have recently been applied to genomics data. While these algorithms can infer unoriented causal interactions between involved molecular variables (i.e., without specifying which one is the cause and which one is the effect), causally orienting all inferred molecular interactions was assumed to be an unsolvable problem until recently. In this work, we use transcription factor-target gene regulatory interactions in three different organisms to evaluate a new family of methods that, given observational data for just two causally related variables, can determine which one is the cause and which one is the effect.

Results

We have found that a particular family of causal orientation methods (IGCI Gaussian) is often able to accurately infer directionality of causal interactions, and that these methods usually outperform other causal orientation techniques. We also introduced a novel ensemble technique for causal orientation that combines decisions of individual causal orientation methods. The ensemble method was found to be more accurate than any best individual causal orientation method in the tested data.

Conclusions

This work represents a first step towards establishing context for practical use of causal orientation methods in the genomics domain. We have found that some causal orientation methodologies yield accurate predictions of causal orientation in genomics data, and we have improved on this capability with a novel ensemble method. Our results suggest that these methods have the potential to facilitate reconstruction of molecular pathways by minimizing the number of required randomized experiments to find causal directionality and by avoiding experiments that are infeasible and/or unethical.

An Integrated Statistical Approach to Compare Transcriptomics Data Across Experiments: A Case Study on the Identification of Candidate Target Genes of the Transcription Factor PPARα

An effective strategy to elucidate the signal transduction cascades activated by a transcription factor is to compare the transcriptional profiles of wild type and transcription factor knockout models. Many statistical tests have been proposed for analyzing gene expression data, but most tests are based on pair-wise comparisons. Since the analysis of microarrays involves the testing of multiple hypotheses within one study, it is generally accepted that one should control for false positives by the false discovery rate (FDR). However, it has been reported that this may be an inappropriate metric for comparing data across different experiments. Here we propose an approach that addresses the above mentioned problem by the simultaneous testing and integration of the three hypotheses (contrasts) using the cell means ANOVA model. These three contrasts test for the effect of a treatment in wild type, gene knockout, and globally over all experimental groups. We illustrate our approach on microarray experiments that focused on the identification of candidate target genes and biological processes governed by the fatty acid sensing transcription factor PPARα in liver. Compared to the often applied FDR based across experiment comparison, our approach identified a conservative but less noisy set of candidate genes with same sensitivity and specificity. However, our method had the advantage of properly adjusting for multiple testing while integrating data from two experiments, and was driven by biological inference. Taken together, in this study we present a simple, yet efficient strategy to compare differential expression of genes across experiments while controlling for multiple hypothesis testing.

Regulation of target gene expression by the vitamin D receptor - an update on mechanisms

Virtually all of the known biological actions of the hormonal ligand 1,25-dihydroxyvitamin D(3) (1,25(OH)(2)D(3)) are mediated by the vitamin D receptor (VDR). Following binding and activation by the ligand, the VDR localizes in the nucleus to the regulatory regions of target genes and recruits chromatin-active coregulatory complexes which, in turn, modulate transcriptional output. The failure of the VDR to function due to crippling mutations results in total hereditary resistance to 1,25(OH)(2)D(3) in both mice and humans. In this review, we summarize the structural and functional properties of the VDR and the role of 1,25(OH)(2)D(3) in receptor activation, and then describe the results of recent studies using genome-wide analyses that define the overarching principles through which the VDR modulates genes expression. We also focus on the recent analysis of a specific 1,25(OH)(2)D(3) regulated gene that provides confirmation of the principles identified through these genome-wide methodologies. Taken together, these studies suggest an unanticipated increase in the complexity of the molecular processes that govern gene regulation by hormones and other factors.

Assessing quality and completeness of human transcriptional regulatory pathways on a genome-wide scale

Background

Pathway databases are becoming increasingly important and almost omnipresent in most types of biological and translational research. However, little is known about the quality and completeness of pathways stored in these databases. The present study conducts a comprehensive assessment of transcriptional regulatory pathways in humans for seven well-studied transcription factors: MYC, NOTCH1, BCL6, TP53, AR, STAT1, and RELA. The employed benchmarking methodology first involves integrating genome-wide binding with functional gene expression data to derive direct targets of transcription factors. Then the lists of experimentally obtained direct targets are compared with relevant lists of transcriptional targets from 10 commonly used pathway databases.

Results

The results of this study show that for the majority of pathway databases, the overlap between experimentally obtained target genes and targets reported in transcriptional regulatory pathway databases is surprisingly small and often is not statistically significant. The only exception is MetaCore pathway database which yields statistically significant intersection with experimental results in 84% cases. Additionally, we suggest that the lists of experimentally derived direct targets obtained in this study can be used to reveal new biological insight in transcriptional regulation and suggest novel putative therapeutic targets in cancer.

Conclusions

Our study opens a debate on validity of using many popular pathway databases to obtain transcriptional regulatory targets. We conclude that the choice of pathway databases should be informed by solid scientific evidence and rigorous empirical evaluation.

Reviewers

This article was reviewed by Prof. Wing Hung Wong, Dr. Thiago Motta Venancio (nominated by Dr. L Aravind), and Prof. Geoff J McLachlan.

Microarray-based identification of transcription factor target genes

Microarray analysis is widely used to identify transcriptional changes associated with genetic perturbation or signaling events. Here we describe its application in the identification of plant transcription factor target genes with emphasis on the design of suitable DNA constructs for controlling TF activity, the experimental setup, the statistical analysis of the microarray data, and the validation of target genes.

A functional analysis of the CREB signaling pathway using HaloCHIP-chip and high throughput reporter assays

Background

Regulation of gene expression is essential for normal development and cellular growth. Transcriptional events are tightly controlled both spatially and temporally by specific DNA-protein interactions. In this study we finely map the genome-wide targets of the CREB protein across all known and predicted human promoters, and characterize the functional consequences of a subset of these binding events using high-throughput reporter assays. To measure CREB binding, we used HaloCHIP, an antibody-free alternative to the ChIP method that utilizes the HaloTag fusion protein, and also high-throughput promoter-luciferase reporter assays, which provide rapid and quantitative screening of promoters for transcriptional activation or repression in living cells.

Results

In analysis of CREB genome-wide binding events using a comprehensive DNA microarray of human promoters, we observe for the first time that CREB has a strong preference for binding at bidirectional promoters and unlike unidirectional promoters, these binding events often occur downstream of transcription start sites. Comparison between HaloCHIP-chip and ChIP-chip data reveal this to be true for both methodologies, indicating it is not a bias of the technology chosen. Transcriptional data obtained from promoter-luciferase reporter arrays also show an unprecedented, high level of activation of CREB-bound promoters in the presence of the co-activator protein TORC1.

Conclusion

These data suggest for the first time that TORC1 provides directional information when CREB is bound at bidirectional promoters and possible pausing of the CREB protein after initial transcriptional activation. Also, this combined approach demonstrates the ability to more broadly characterize CREB protein-DNA interactions wherein not only DNA binding sites are discovered, but also the potential of the promoter sequence to respond to CREB is evaluated.

Deregulated direct targets of the hepatitis B virus (HBV) protein, HBx, identified through chromatin immunoprecipitation and expression microarray profiling

The hepatitis B-X (HBx) protein is strongly associated with hepatocellular carcinoma. It is implicated not to directly cause cancer but to play a role in hepatocellular carcinoma as a co-factor. The oncogenic potential of HBx primarily lies in its interaction with transcriptional regulators resulting in aberrant gene expression and deregulated cellular pathways. Utilizing ultraviolet irradiation to simulate a tumor-initiating event, we integrated chip-based chromatin immunoprecipitation (ChIP-chip) with expression microarray profiling and identified 184 gene targets directly deregulated by HBx. One-hundred forty-four transcription factors interacting with HBx were computationally inferred. We experimentally validated that HBx interacts with some of the predicted transcription factors (pTF) as well as the promoters of the deregulated target genes of these pTFs. Significantly, we demonstrated that the pTF interacts with the promoters of the deregulated HBx target genes and that deregulation by HBx of these HBx target genes carrying the pTF consensus sequences can be reversed using pTF small interfering RNAs. The roles of these deregulated direct HBx target genes and their relevance in cancer was inferred via querying against biogroup/cancer-related microarray databases using web-based NextBio(TM) software. Six pathways, including the Jak-STAT pathway, were predicted to be significantly deregulated when HBx binds indirectly to direct target gene promoters. In conclusion, this study represents the first ever demonstration of the utilization of ChIP-chip to identify deregulated direct gene targets from indirect protein-DNA binding as well as transcriptional factors directly interacting with HBx. Increased knowledge of the gene/transcriptional factor targets of HBx will enhance our understanding of the role of HBx in hepatocellular carcinogenesis and facilitate the design of better strategies in combating hepatitis B virus-associated hepatocellular carcinoma.

Lighting up individual DNA binding proteins with quantum dots

The ability to determine the precise loci and occupancy of DNA-binding proteins is instrumental to our understanding of cellular processes like gene expression and regulation. We propose a single-molecule approach for the direct visualization of proteins bound to their template DNA. Fluorescent quantum dots (QD) are used to label proteins bound to DNA, allowing multicolor, nanometer-resolution localization. Protein-DNA complexes are linearly extended and imaged to determine the precise location of the protein binding sites. The method is demonstrated by detecting individual QD-labeled T7-RNA polymerases on the T7 bacteriophage genome. This work demonstrates the potential of this approach to precisely read protein binding position or, alternatively, "write" such information on extended DNA with QDs via sequence-specific molecular recognition.

Chromatin immunoprecipitation assays: analyzing transcription factor binding and histone modifications in vivo

Studies in the past decade have shown that differential gene expression depends not only on the binding of specific transcription factors to discrete promoter elements but also on the epigenetic modification of the DNA as well as histones associated with the promoter. While techniques like electrophoretic mobility shift assays could detect and characterize the binding of specific transcription factors present in cell lysates to DNA sequences in in vitro binding conditions, they were not effective in assessing the binding in intact cells. Development of chromatin immunoprecipitation technique in the past decade enabled the analysis of the association of regulatory molecules with specific promoters or changes in histone modifications in vivo, without overexpressing any component. ChIP assays can provide a snapshot of how a regulatory transcription factor affects the expression of a single gene or a variety of genes at the same time. Availability of high-quality antibodies that recognizes histones modified in a specific fashion further expanded the use of ChIP assays to analyze even minute changes in histone modification and nucleosomes structure. This chapter outlines the general strategies and protocols used to carry out ChIP assays to study the differential recruitment of transcription factors as well as histone modifications.

Chromatin immunoprecipitation (ChIP) for analysis of histone modifications and chromatin-associated proteins

Disruption of epigenetic regulators of transcription is a central mechanism of oncogenesis. Many of the advances in the understanding of these mechanisms are attributable to the successful development of chromatin immunoprecipitation (ChIP) for in vivo detection of histone modifications as well as chromatin binding regulatory proteins. This is a powerful technique for analyzing histone modifications as well as binding sites for proteins that bind either directly or indirectly to DNA. Here we present two ChIP protocols. The first is particularly useful for identifying histone modifications or binding at specific, known genomic sites. The second, employing serial analysis of gene expression, is particularly powerful for the discovery of previously unidentified sites of modification or binding.

ChIP on chip assays: genome-wide analysis of transcription factor binding and histone modifications

Deregulation of transcriptional activity of many genes has been causatively linked to human diseases including cancer. Altered patterns of gene expression in normal and cancer cells are the result of inappropriate expression of transcription factors and chromatin-modifying proteins. Chromatin immunoprecipitation assay is a well-established tool for investigating the interactions between regulatory proteins and DNA at distinct stages of gene activation. ChIP coupled with DNA microarrays, known as ChIP on chip, allow us to determine the entire spectrum of in vivo DNA-binding sites for a given protein. This has been of immense value because ChIP on chip assays can provide a snapshot of the transcriptional regulatory mechanisms on a genome-wide scale. This article outlines the general strategies used to carry out ChIP-chip assays to study the differential recruitment of regulatory molecules based on the studies conducted in our lab as well as other published protocols.

Identification and characterization of CCAAT/Enhancer Binding proteindelta (C/EBPdelta) target genes in G0 growth arrested mammary epithelial cells

BACKGROUND

CCAAT/Enhancer Binding Proteindelta (C/EBPdelta) is a member of the highly conserved C/EBP family of leucine zipper (bZIP) proteins. C/EBPdelta is highly expressed in G0 growth arrested mammary epithelial cells (MECs) and "loss of function" alterations in C/EBPdelta have been associated with impaired contact inhibition, increased genomic instability and increased cell migration. Reduced C/EBPdelta expression has also been reported in breast cancer and acute myeloid leukemia (AML). C/EBPdelta functions as a transcriptional activator, however, only a limited number of C/EBPdelta target genes have been reported. As a result, the role of C/EBPdelta in growth control and the potential mechanisms by which "loss of function" alterations in C/EBPdelta contribute to tumorigenesis are poorly understood. The goals of the present study were to identify C/EBPdelta target genes using Chromatin Immunoprecipitation coupled with a CpG Island (HCG12K) Array gene chip ("ChIP-chip") assay and to assess the expression and potential functional roles of C/EBPdelta target genes in growth control.

RESULTS

ChIP-chip assays identified approximately 100 C/EBPdelta target gene loci which were classified by gene ontology (GO) into cell adhesion, cell cycle regulation, apoptosis, signal transduction, intermediary metabolism, gene transcription, DNA repair and solute transport categories. Conventional ChIP assays validated the ChIP-chip results and demonstrated that 14/14 C/EBPdelta target loci were bound by C/EBPdelta in G0 growth arrested MCF-12A MECs. Gene-specific RT-PCR analysis also demonstrated C/EBPdelta-inducible expression of 14/14 C/EBPdelta target genes in G0 growth arrested MCF-12A MECs. Finally, expression of endogenous C/EBPdelta and selected C/EBPdelta target genes was also demonstrated in contact-inhibited G0 growth arrested nontransformed human MCF-10A MECs and in mouse HC11 MECs. The results demonstrate consistent activation and downstream function of C/EBPdelta in growth arrested human and murine MECs.

CONCLUSION

C/EBPdelta target genes were identified by a global gene array approach and classified into functional categories that are consistent with biological contexts in which C/EBPdelta is induced, such as contact-mediated G0 growth arrest, apoptosis, metabolism and inflammation. The identification and validation of C/EBPdelta target genes provides new insights into the mechanistic role of C/EBPdelta in mammary epithelial cell biology and sheds new light on the potential impact of "loss of function" alterations in C/EBPdelta in tumorigenesis.

Identification of direct regulatory targets of the transcription factor Sox10 based on function and conservation

Background

Sox10, a member of the Sry-related HMG-Box gene family, is a critical transcription factor for several important cell lineages, most notably the neural crest stem cells and the derivative peripheral glial cells and melanocytes. Thus far, only a handful of direct target genes are known for this transcription factor limiting our understanding of the biological network it governs.

Results

We describe identification of multiple direct regulatory target genes of Sox10 through a procedure based on function and conservation. By combining RNA interference technique and DNA microarray technology, we have identified a set of genes that show significant down-regulation upon introduction of Sox10 specific siRNA into Schwannoma cells. Subsequent comparative genomics analyses led to potential binding sites for Sox10 protein conserved across several mammalian species within the genomic region proximal to these genes. Multiple sites belonging to 4 different genes (proteolipid protein, Sox10, extracellular superoxide dismutase, and pleiotrophin) were shown to directly interact with Sox10 by chromatin immunoprecipitation assay. We further confirmed the direct regulation through the identified cis-element for one of the genes, extracellular superoxide dismutase, using electrophoretic mobility shift assay and reporter assay.

Conclusion

In sum, the process of combining differential expression profiling and comparative genomics successfully led to further defining the role of Sox10, a critical transcription factor for the development of peripheral glia. Our strategy utilizing relatively accessible techniques and tools should be applicable to studying the function of other transcription factors.

Polycomb complex 2 is required for E-cadherin repression by the Snail1 transcription factor

The transcriptional factor Snail1 is a repressor of E-cadherin (CDH1) gene expression essential for triggering epithelial-mesenchymal transition. Snail1 represses CDH1, directly binding its promoter and inducing the synthesis of the Zeb1 repressor. In this article, we show that repression of CDH1 by Snail1, but not by Zeb1, is dependent on the activity of Polycomb repressive complex 2 (PRC2). Embryonic stem (ES) cells null for Suz12, one of the components of PRC2, show higher levels of Cdh1 mRNA than control ES cells. In tumor cells, interference of PRC2 activity prevents the ability of Snail1 to downregulate CDH1 and partially derepresses CDH1. Chromatin immunoprecipitation assays demonstrated that Snail1 increases the binding of Suz12 to the CDH1 promoter and the trimethylation of lysine 27 in histone H3. Moreover, Snail1 interacts with Suz12 and Ezh2, as shown by coimmunoprecipitation experiments. In conclusion, these results demonstrate that Snail1 recruits PRC2 to the CDH1 promoter and requires the activity of this complex to repress E-cadherin expression.

Chromatin immunoprecipitation assays: molecular analysis of chromatin modification and gene regulation

Gene expression pattern in cancer cells differ significantly from their normal counter parts, owing to mutations in oncogenes and tumor suppressor genes, their downstream targets, or owing to increased proliferation, and altered apoptotic potential. Various microarray based techniques have been widely utilized to study the differential expression of genes in cancer in recent years. Along with this, attempts have been made to study the transcriptional regulatory mechanisms and chromatin modifications facilitating such differential gene expression. One of the widely used assays for this purpose is the chromatin immunoprecipitation (ChIP) assay, which enables the analysis of the association of regulatory molecules with specific promoters or changes in histone modifications in vivo, without overexpressing any component. This has been of immense value, because ChIP assays can provide a snapshot of the regulatory mechanisms involved in the expression of a single gene, or a variety of genes at the same time. This review article outlines the general strategies and protocols used to carry out ChIP assays to study the differential recruitment of transcription factors, based on the experience in studying E2F1 and histone modifications as well as other published protocols. In addition, the use of ChIP assays to carry out global analysis of transcription factor recruitment is also addressed.

Transcriptional control of erythropoiesis: emerging mechanisms and principles

Transcriptional networks orchestrate fundamental biological processes, including hematopoiesis, in which hematopoietic stem cells progressively differentiate into specific progenitors cells, which in turn give rise to the diverse blood cell types. Whereas transcription factors recruit coregulators to chromatin, leading to targeted chromatin modification and recruitment of the transcriptional machinery, many questions remain unanswered regarding the underlying molecular mechanisms. Furthermore, how diverse cell type-specific transcription factors function cooperatively or antagonistically in distinct cellular contexts is poorly understood, especially since genes in higher eukaryotes commonly encompass broad chromosomal regions (100 kb and more) and are littered with dispersed regulatory sequences. In this article, we describe an important set of transcription factors and coregulators that control erythropoiesis and highlight emerging transcriptional mechanisms and principles. It is not our intent to comprehensively survey all factors implicated in the transcriptional control of erythropoiesis, but rather to underscore specific mechanisms, which have potential to be broadly relevant to transcriptional control in diverse systems.

A novel non-overlapping bi-clustering algorithm for network generation using living cell array data

MOTIVATION

The living cell array quantifies the contribution of activated transcription factors upon the expression levels of their target genes. The direct manipulation of the regulatory mechanisms offers enormous possibilities for deciphering the machinery that activates and controls gene expression. We propose a novel bi-clustering algorithm for generating non-overlapping clusters of reporter genes and conditions and demonstrate how this information can be interpreted in order to assist in the construction of transcription factor interaction networks.

The transcriptional repressor NIPP1 is an essential player in EZH2-mediated gene silencing

EZH2 is a Polycomb group (PcG) protein that promotes the late-stage development of cancer by silencing a specific set of genes, at least in part through trimethylation of associated histone H3 on Lys 27 (H3K27). Nuclear inhibitor of protein phosphatase-1 (NIPP1) is a ubiquitously expressed transcriptional repressor that has binding sites for the EZH2 interactor EED. Here, we examine the contribution of NIPP1 to EZH2-mediated gene silencing. Studies on NIPP1-deficient cells disclose a widespread and essential role of NIPP1 in the trimethylation of H3K27 by EZH2, not only in the onset of this trimethylation during embryonic development, but also in the maintenance of this repressive mark in proliferating cells. Consistent with this notion, EZH2 and NIPP1 silence a common set of genes, as revealed by gene-expression profiling, and NIPP1 is associated with established Polycomb target genes and with genomic regions that are enriched in Polycomb targets. Furthermore, most NIPP1 target genes are trimethylated on H3K27 and the knockdown of either NIPP1 or EZH2 is often associated with a loss of this modification. Our data reveal that NIPP1 is required for the global trimethylation of H3K27 and is implicated in gene silencing by EZH2.

Identification of multiple MAPK-mediated transcription factors regulated by tobacco smoke in airway epithelial cells

Activation and regulation of transcription factors (TFs) are the major mechanisms regulating changes in gene expression upon environmental exposure. Tobacco smoke (TS) is a complex mixture of chemicals, each of which could act through different signal cascades, leading to the regulation of distinct TFs and alterations in subsequent gene expression. We proposed that TS exposure affects inflammatory gene expression at the transcriptional level by modulating the DNA binding activities of TFs. To investigate transcriptional regulation upon TS exposure, a protein/DNA array was applied to screen TFs that are affected by TS exposure. This array-based screening allowed us to simultaneously detect 244 different TFs. Our results indicated that multiple TFs were rapidly activated upon TS exposure. DNA-binding activity of differentially expressed TFs was confirmed by EMSA. Our results showed that at least 20 TFs displayed more than twofold expressional changes after smoke treatment. Ten smoke-induced TFs, including NF-kappaB, VDR, ISRE, and RSRFC4, were involved in MAPK signaling pathways. The NF-kappaB family, which is involved in inflammation-induced gene activation, was selected for further study to characterize TS exposure-induced transcriptional activation. Western blot analysis and immunofluorescence microscopy indicated that TS exposure induced phosphorylation of IkappaB and translocation of NF-kappaB p65/p50 heterodimers into the nucleus. This activity was abrogated by the MAPK inhibitors PD98059 and U0126. Our results confirmed that activation of MAPK signaling pathways by TS exposure increased transcriptional activity of NF-kappaB. These data provide a potential mechanism for TS-induced inflammatory gene expression.

Characteristics of oligonucleotide tiling arrays measured by hybridizing full-length cDNA clones: causes of signal variation and false positive signals

An assessment of the hybridization characteristics of oligonucleotide tiling arrays was carried out using 162 full-length sequenced cDNA clones in spike-in experiments. The properties of array probes that influence signal intensity were investigated, and their capability in the detection of the cDNA exons was evaluated. The signal intensities detected in exonic and nonexonic genomic regions were examined by focusing on the features of probe sequences that raise or lower the level of intensity and on the causes of false positive signals found in nonexonic regions. The effectiveness of measures used in published protocols to improve the separation between signal and background intensity distributions, including the use of replicates and threshold parameterization of signal intensity, was assessed. Sensitivity and specificity in the detection of exons were measured using various sets of threshold parameters, and the effects of each parameter on the detection efficiency and the rate of false positives were evaluated. It was also demonstrated that hybridization of full-length cDNA clones is an excellent method to investigate the characteristics of oligonucleotide tiling arrays.

Locating mammalian transcription factor binding sites: a survey of computational and experimental techniques

Fields such as genomics and systems biology are built on the synergism between computational and experimental techniques. This type of synergism is especially important in accomplishing goals like identifying all functional transcription factor binding sites in vertebrate genomes. Precise detection of these elements is a prerequisite to deciphering the complex regulatory networks that direct tissue specific and lineage specific patterns of gene expression. This review summarizes approaches for in silico, in vitro, and in vivo identification of transcription factor binding sites. A variety of techniques useful for localized- and high-throughput analyses are discussed here, with emphasis on aspects of data generation and verification.

New p63 targets in keratinocytes identified by a genome-wide approach

p63 is a developmentally regulated transcription factor related to p53. It is involved in the development of ectodermal tissues, including limb, skin and in general, multilayered epithelia. The DeltaNp63alpha isoform is thought to play a 'master' role in the asymmetric division of epithelial cells. It is also involved in the pathogenesis of several human diseases, phenotypically characterized by ectodermal dysplasia. Our understanding of transcriptional networks controlled by p63 is limited, owing to the low number of bona fide targets. To screen for new targets, we employed chromatin immunoprecipitation from keratinocytes (KCs) coupled to the microarray technology, using both CpG islands and promoter arrays. The former revealed 96 loci, the latter yielded 85 additional genes. We tested 40 of these targets in several functional assays, including: (i) in vivo binding by p63 in primary KCs; (ii) expression analysis in differentiating HaCaT cells and in cells overexpressing DeltaNp63alpha; (iii) promoter transactivation and (iv) immunostaining in normal tissues, confirming their regulation by p63. We discovered several new specific targets whose functional categorization links p63 to cell growth and differentiation.

Strategies to unravel molecular codes essential for the development of meso-diencephalic dopaminergic neurons

Understanding the development of neuronal systems has become an important asset in the attempt to solve complex questions about neuropathology as found in Parkinson's disease, schizophrenia and other complex neuronal diseases. The development of anatomical and functional divergent structures in the brain is achieved by a combination of early anatomical patterning and highly coordinated neuronal migration and differentiation events. Fundamental to the existence of divergent structures in the brain is the early region-specific molecular programming. Neuronal progenitors located along the neural tube can still adapt many different identities. Their exact position in the developing brain, however, determines early molecular specification by region-specific signalling molecules. These signals determine time and region-specific expression of early regulatory genes, leading to neuronal differentiation. Here, we focus on a well-described neuronal group, the meso-diencephalic dopaminergic neurons, of which heterogeneity based on anatomical position could account for the difference in vulnerability of specific subgroups as observed in Parkinson's disease. The knowledge of their molecular coding helps us to understand how the meso-diencephalic dopaminergic system is built and could provide clues that unravel mechanisms associated with the neuropathology in complex diseases such as Parkinson's disease.

Adaptively inferring human transcriptional subnetworks

Although the human genome has been sequenced, progress in understanding gene regulation in humans has been particularly slow. Many computational approaches developed for lower eukaryotes to identify cis-regulatory elements and their associated target genes often do not generalize to mammals, largely due to the degenerate and interactive nature of such elements. Motivated by the switch-like behavior of transcriptional responses, we present a systematic approach that allows adaptive determination of active transcriptional subnetworks (cis-motif combinations, the direct target genes and physiological processes regulated by the corresponding transcription factors) from microarray data in mammals, with accuracy similar to that achieved in lower eukaryotes. Our analysis uncovered several new subnetworks active in human liver and in cell-cycle regulation, with similar functional characteristics as the known ones. We present biochemical evidence for our predictions, and show that the recently discovered G2/M-specific E2F pathway is wider than previously thought; in particular, E2F directly activates certain mitotic genes involved in hepatocellular carcinomas. Additionally, we demonstrate that this method can predict subnetworks in a condition-specific manner, as well as regulatory crosstalk across multiple tissues. Our approach allows systematic understanding of how phenotypic complexity is regulated at the transcription level in mammals and offers marked advantage in systems where little or no prior knowledge of transcriptional regulation is available.

Exploring Transcription Factor Binding Properties of Several Non-coding DNA Sequence Elements in the Human NF-IL6 Gene

We examined several DNA segments upstream of the transcription start site of the human NF-IL6 gene to evaluate the predictions of two computational models developed to identify potential regulatory elements in the non-coding regions of genes. One model, comparative genomics, is based on the hypothesis that functional regulatory sequences can be localized in alignments of genomic DNA from several species. The other model is based on the hypothesis that protein-binding sites in genomic DNA may include sequence elements that occur frequently in proximal promoters of genes. The segments selected for DNA binding and functional evaluations included: (1) two conserved regions identified in multi-species sequence alignments; (2) a region containing several localized hits with 9-mers that ranked highly in studies of proximal promoters of human genes; and (3) two regions that were either GC-rich and/or contained tracts of G. The assays were done under nearly identical experimental conditions, using a cell line (U937) representing human monocytes/macrophages. The experiments also aimed at evaluating what effect, if any, cellular stimulation could have on the interactions of nuclear proteins with naturally occurring GC-rich elements in a human genomic DNA. In DNA binding assays, several complexes were formed with the conserved regions identified in multi-species sequence alignment. Furthermore, these regions were active in functional assays. The region containing several matches with 9-mers derived from proximal promoters of human genes was not conserved but formed several complexes with nuclear proteins including Sp1, Egr-1, and an unidentified protein. In addition, this region was active in functional assays and responded to cellular stimulations. Overall, the results of the assays suggest an important role for the sequence context of genomic DNA in protein binding and selection.

Chromatin immunoprecipitation identifies photoreceptor transcription factor targets in mouse models of retinal degeneration: new findings and challenges

The transcription factors, Otx2, Crx, Nrl, and Nr2e3, expressed by retinal photoreceptor cells are essential for photoreceptor gene expression, development, and maintenance. Malfunction of any of these factors due to genetic mutations causes photoreceptor disease. Protein-protein interaction studies suggest that these factors may form a regulatory network centered on Crx. To understand how these factors regulate photoreceptor gene transcription in vivo, we have employed chromatin immunoprecipitation (ChIP) assays to assess the ability of these proteins to bind to regulatory sequences of photoreceptor genes in the retina of wild-type and mutant mice with photoreceptor degeneration. This paper summarizes the advantages and limitations of ChIP, using examples from our studies to demonstrate how this technique has contributed to our understanding of the regulation of photoreceptor gene expression. We report that Crx, Otx2, Nrl, and Nr2e3 co-occupy the promoter/enhancer, but not the region 3' of selected Crx target genes, in a retina-specific fashion. We identified Crx-dependent (Nr2e3) and Crx-independent (Otx2 and Nrl) target binding using Crx knockout mice (Crx-/-), suggesting that individual factors may use distinct mechanism(s) for binding and regulating target genes. Consistent with ChIP results, we also found that Otx2, a close family member of Crx, can activate the promoter of rod and cone genes in HEK293 cells, implicating Otx2 in regulating photoreceptor gene expression. These findings provide important information for understanding how photoreceptor transcription factors regulate photoreceptor gene expression and the mechanisms by which mutations in these factors cause transcriptional dysregulation and photoreceptor degeneration.

Design optimization methods for genomic DNA tiling arrays

A recent development in microarray research entails the unbiased coverage, or tiling, of genomic DNA for the large-scale identification of transcribed sequences and regulatory elements. A central issue in designing tiling arrays is that of arriving at a single-copy tile path, as significant sequence cross-hybridization can result from the presence of non-unique probes on the array. Due to the fragmentation of genomic DNA caused by the widespread distribution of repetitive elements, the problem of obtaining adequate sequence coverage increases with the sizes of subsequence tiles that are to be included in the design. This becomes increasingly problematic when considering complex eukaryotic genomes that contain many thousands of interspersed repeats. The general problem of sequence tiling can be framed as finding an optimal partitioning of non-repetitive subsequences over a prescribed range of tile sizes, on a DNA sequence comprising repetitive and non-repetitive regions. Exact solutions to the tiling problem become computationally infeasible when applied to large genomes, but successive optimizations are developed that allow their practical implementation. These include an efficient method for determining the degree of similarity of many oligonucleotide sequences over large genomes, and two algorithms for finding an optimal tile path composed of longer sequence tiles. The first algorithm, a dynamic programming approach, finds an optimal tiling in linear time and space; the second applies a heuristic search to reduce the space complexity to a constant requirement. A Web resource has also been developed, accessible at http://tiling.gersteinlab.org, to generate optimal tile paths from user-provided DNA sequences.

A genomic approach to the identification and characterization of HOXA13 functional binding elements

HOX proteins are important transcriptional regulators in mammalian embryonic development and are dysregulated in human cancers. However, there are few known direct HOX target genes and their mechanisms of regulation are incompletely understood. To isolate and characterize gene segments through which HOX proteins regulate transcription we used cesium chloride centrifugation-based chromatin purification and immunoprecipitation (ChIP). From NIH 3T3-derived HOXA13-FLAG expressing cells, 33% of randomly selected, ChIP clones were reproducibly enriched. Hox-enriched fragments (HEFs) were more AT-rich compared with cloned fragments that failed reproducible ChIP. All HEFs augmented transcription of a heterologous promoter upon coexpression with HOXA13. One HEF was from intron 2 of Enpp2, a gene highly upregulated in these cells and has been implicated in cell motility. Using Enpp2 as a candidate direct target, we identified three additional HEFs upstream of the transcription start site. HOXA13 upregulated transcription from an Enpp2 promoter construct containing these sites, and each site was necessary for full HOXA13-induced expression. Lastly, given that HOX proteins have been demonstrated to interact with histone deacetylases and/or CBP, we explored whether histone acetylation changed at Enpp2 upon HOXA13-induced activation. No change in the general histone acetylation state was observed. Our results support models in which occupation of multiple HOX binding sites is associated with highly activated genes.

Chromatin domain activation via GATA-1 utilization of a small subset of dispersed GATA motifs within a broad chromosomal region

Cis elements that mediate transcription factor binding are abundant within genomes, but the rules governing occupancy of such motifs in chromatin are not understood. The transcription factor GATA-1 that regulates red blood cell development binds with high affinity to GATA motifs, and initial studies suggest that these motifs are often unavailable for occupancy in chromatin. Whereas GATA-2 regulates the differentiation of all blood cell lineages via GATA motif binding, the specificity of GATA-2 chromatin occupancy has not been studied. We found that conditionally active GATA-1 (ER-GATA-1) and GATA-2 occupy only a small subset of the conserved GATA motifs within the murine beta-globin locus. Kinetic analyses in GATA-1-null cells indicated that ER-GATA-1 preferentially occupied GATA motifs at the locus control region (LCR), in which chromatin accessibility is largely GATA-1-independent. Subsequently, ER-GATA-1 increased promoter accessibility and occupied the betamajor promoter. ER-GATA-1 increased erythroid Krüppel-like factor and SWI/SNF chromatin remodeling complex occupancy at restricted LCR sites. These studies revealed three phases of beta-globin locus activation: GATA-1-independent establishment of specific chromatin structure features, GATA-1-dependent LCR complex assembly, and GATA-1-dependent promoter complex assembly. The differential utilization of dispersed GATA motifs therefore establishes spatial/temporal regulation and underlies the multistep activation mechanism.

Chromatin immunoprecipitation-based screen to identify functional genomic binding sites for sequence-specific transactivators

In various human diseases, altered gene expression patterns are often the result of deregulated gene-specific transcription factor activity. To further understand disease on a molecular basis, the comprehensive analysis of transcription factor signaling networks is required. We developed an experimental approach, combining chromatin immunoprecipitation (ChIP) with a yeast-based assay, to screen the genome for transcription factor binding sites that link to transcriptionally regulated target genes. We used the tumor suppressor p53 to demonstrate the effectiveness of the method. Using primary and immortalized, nontransformed cultures of human mammary epithelial cells, we isolated over 100 genomic DNA fragments that contain novel p53 binding sites. This approach led to the identification and validation of novel p53 target genes involved in diverse signaling pathways, including growth factor signaling, protein kinase/phosphatase signaling, and RNA binding. Our results yield a more complete understanding of p53-regulated signaling pathways, and this approach could be applied to any number of transcription factors to further elucidate complex transcriptional networks.

Issues in the analysis of oligonucleotide tiling microarrays for transcript mapping

Traditional microarrays use probes complementary to known genes to quantitate the differential gene expression between two or more conditions. Genomic tiling microarray experiments differ in that probes that span a genomic region at regular intervals are used to detect the presence or absence of transcription. This difference means the same sets of biases and the methods for addressing them are unlikely to be relevant to both types of experiment. We introduce the informatics challenges arising in the analysis of tiling microarray experiments as open problems to the scientific community and present initial approaches for the analysis of this nascent technology.

Glucocorticoid receptor-dependent gene regulatory networks

While the molecular mechanisms of glucocorticoid regulation of transcription have been studied in detail, the global networks regulated by the glucocorticoid receptor (GR) remain unknown. To address this question, we performed an orthogonal analysis to identify direct targets of the GR. First, we analyzed the expression profile of mouse livers in the presence or absence of exogenous glucocorticoid, resulting in over 1,300 differentially expressed genes. We then executed genome-wide location analysis on chromatin from the same livers, identifying more than 300 promoters that are bound by the GR. Intersecting the two lists yielded 53 genes whose expression is functionally dependent upon the ligand-bound GR. Further network and sequence analysis of the functional targets enabled us to suggest interactions between the GR and other transcription factors at specific target genes. Together, our results further our understanding of the GR and its targets, and provide the basis for more targeted glucocorticoid therapies.