High-throughput localization of functional elements by quantitative chromatin profiling

Genome-wide mechanisms of nuclear receptor action

Nuclear receptors are involved in a myriad of physiological processes, responding to ligands and binding to DNA at sequence-specific cis-regulatory elements. This binding occurs in the context of chromatin, a critical factor in regulating eukaryotic transcription. Recent high-throughput assays have examined nuclear receptor action genome-wide, advancing our understanding of receptor binding to regulatory elements. Here, we discuss current knowledge of genome-wide response element occupancy by receptors and the function of transcription factor networks in regulating nuclear receptor action. We highlight emerging roles for the epigenome, chromatin remodeling, histone modification, histone variants and long-range chromosomal interactions in nuclear receptor binding and receptor-dependent gene regulation. These mechanisms contribute importantly to the action of nuclear receptors in health and disease.

Rapid high-throughput analysis of DNaseI hypersensitive sites using a modified Multiplex Ligation-dependent Probe Amplification approach

BACKGROUND

Mapping DNaseI hypersensitive sites is commonly used to identify regulatory regions in the genome. However, currently available methods are either time consuming and laborious, expensive or require large numbers of cells. We aimed to develop a quick and straightforward method for the analysis of DNaseI hypersensitive sites that overcomes these problems.

RESULTS

We have developed a modified Multiplex Ligation-dependent Probe Amplification (MLPA) approach for the identification and analysis of genomic regulatory regions. The utility of this approach was demonstrated by simultaneously analysing 20 loci from the ENCODE project for DNaseI hypersensitivity in a range of different cell lines. We were able to obtain reproducible results with as little as 5 x 10(4) cells per DNaseI treatment. Our results broadly matched those previously reported by the ENCODE project, and both technical and biological replicates showed high correlations, indicating the sensitivity and reproducibility of this method.

CONCLUSION

This new method will considerably facilitate the identification and analysis of DNaseI hypersensitive sites. Due to the multiplexing potential of MLPA (up to 50 loci can be examined) it is possible to analyse dozens of DNaseI hypersensitive sites in a single reaction. Furthermore, the high sensitivity of MLPA means that fewer than 10(5) cells per DNaseI treatment can be used, allowing the discovery and analysis of tissue specific regulatory regions without the need for pooling. This method is quick and easy and results can be obtained within 48 hours after harvesting of cells or tissues. As no special equipment is required, this method can be applied by any laboratory interested in the analysis of DNaseI hypersensitive regions.

The HBS1L-MYB intergenic interval associated with elevated HbF levels shows characteristics of a distal regulatory region in erythroid cells

HBS1L-MYB intergenic polymorphism (HMIP) on chromosome 6q23 is associated with elevated fetal hemoglobin levels and has pleiotropic effects on several hematologic parameters. To investigate potential regulatory activity in the region, we have measured sensitivity of the sequences to DNase I cleavage that identified 3 tissue-specific DNase I hypersensitive sites in the core intergenic interval. Chromatin immunoprecipitation with microarray (ChIP-chip) analysis showed strong histone acetylation in a defined interval of 65 kb corresponding to the core HBS1L-MYB intergenic region in primary human erythroid cells but not in non–MYB-expressing HeLa cells. ChIP-chip analysis also identified several potential cis-regulatory elements as strong GATA-1 signals that coincided with the DNase I hypersensitive sites present in MYB-expressing erythroid cells. We suggest that HMIP contains regulatory sequences that could be important in hematopoiesis by controlling MYB expression. This study provides the functional link between genetic association of HMIP with control of fetal hemoglobin and other hematologic parameters. We also present a large-scale analysis of histone acetylation as well as RNA polymerase II and GATA-1 interactions on chromosome 6q, and α and β globin gene loci. The data suggest that GATA-1 regulates numerous genes of various functions on chromosome 6q.

Genetic and epigenetic networks controlling T helper 1 cell differentiation

Significant progress has been made during the past years in our understanding of the mechanisms that control the differentiation of naïve CD4(+) T cells into effector T-cell subsets with distinct functional properties. Previous work allowed the identification of key molecules involved in regulating this highly complex process, such as cytokines and their receptors, signal transducers and transcription factors. More recently, the emphasis of research in this field has been to elucidate how the multiplicity of signals is integrated to shape a T helper subset-specific gene-expression program controlling differentiation and effector functions. In this review we will highlight advances that have been made in unravelling the genetic and epigenetic networks controlling differentiation of naïve CD4(+) T cells into interferon-gamma(IFN-gamma)-secreting T helper type 1 (Th1) cells.

A thermodynamic approach to PCR primer design

We developed a primer design method, Pythia, in which state of the art DNA binding affinity computations are directly integrated into the primer design process. We use chemical reaction equilibrium analysis to integrate multiple binding energy calculations into a conservative measure of polymerase chain reaction (PCR) efficiency, and a precomputed index on genomic sequences to evaluate primer specificity. We show that Pythia can design primers with success rates comparable with those of current methods, but yields much higher coverage in difficult genomic regions. For example, in RepeatMasked sequences in the human genome, Pythia achieved a median coverage of 89% as compared with a median coverage of 51% for Primer3. For parameter settings yielding sensitivities of 81%, our method has a recall of 97%, compared with the Primer3 recall of 48%. Because our primer design approach is based on the chemistry of DNA interactions, it has fewer and more physically meaningful parameters than current methods, and is therefore easier to adjust to specific experimental requirements. Our software is freely available at http://pythia.sourceforge.net.

Histone acetylation characterizes chromatin presetting by NF1 and Oct1 and enhances glucocorticoid receptor binding to the MMTV promoter

Transcription from the mouse mammary tumor virus (MMTV) promoter is induced by the glucocorticoid receptor (GR). This switch was reconstituted in Xenopus oocytes. Previously, we showed that Nuclear Factor 1 (NF1) and Octamer Transcription Factor 1 (Oct1) bind constitutively to the MMTV promoter and thereby induce translational nucleosome positioning representing an intermediary, i.e. preset, state of nucleosome organization. Here we further characterize this NF1 and Oct1 induced preset chromatin in relation to the inactive and the hormone-activated state. The preset chromatin exhibits increased histone acetylation but does not cause dissociation of histone H1 as oppose to the hormone-activated state. Furthermore, upon hormone induction the preset MMTV chromatin displays an enhanced and prolonged GR binding capacity and transcription during an intrinsic and time-dependent silencing of the injected template. The silencing process correlates with a reduced histone acetylation. However, a histone deacetylase inhibitor, trichostatin A (TSA), does not counteract silencing in spite of its distinct stimulation of GR-DNA binding. The latter indicates the importance of histone acetylation to maintain DNA access for inducible factor binding. We discuss how constitutively bound factors such as NF1 and Oct1 may participate in the maintenance of tissue specificity of hormone responsive genes.

Glucocorticoid receptor dynamics and gene regulation

The glucocorticoid receptor regulates the expression of a large number of genes in mammalian cells. The interaction of this receptor with regulatory elements has been discovered to be highly dynamic, with occupancy states measured in seconds, rather than minutes or hours. This finding has led to a paradigm shift in our understanding of receptor function throughout the genome. The mechanisms involved in these rapid exchange events, as well as the implications for receptor function, are discussed.

Isolation of active regulatory elements from eukaryotic chromatin using FAIRE (Formaldehyde Assisted Isolation of Regulatory Elements)

The binding of sequence-specific regulatory factors and the recruitment of chromatin remodeling activities cause nucleosomes to be evicted from chromatin in eukaryotic cells. Traditionally, these active sites have been identified experimentally through their sensitivity to nucleases. Here we describe the details of a simple procedure for the genome-wide isolation of nucleosome-depleted DNA from human chromatin, termed FAIRE (Formaldehyde Assisted Isolation of Regulatory Elements). We also provide protocols for different methods of detecting FAIRE-enriched DNA, including use of PCR, DNA microarrays, and next-generation sequencing. FAIRE works on all eukaryotic chromatin tested to date. To perform FAIRE, chromatin is crosslinked with formaldehyde, sheared by sonication, and phenol-chloroform extracted. Most genomic DNA is crosslinked to nucleosomes and is sequestered to the interphase, whereas DNA recovered in the aqueous phase corresponds to nucleosome-depleted regions of the genome. The isolated regions are largely coincident with the location of DNaseI hypersensitive sites, transcriptional start sites, enhancers, insulators, and active promoters. Given its speed and simplicity, FAIRE has utility in establishing chromatin profiles of diverse cell types in health and disease, isolating DNA regulatory elements en masse for further characterization, and as a screening assay for the effects of small molecules on chromatin organization.

Dynamic access of the glucocorticoid receptor to response elements in chromatin

Transcriptional activation as a rate-limiting step of gene expression is often triggered by an environmental stimulus that is transmitted through a signaling cascade to specific transcription factors. Transcription factors must then find appropriate target genes in the context of chromatin. Subsequent modulation of local chromatin domains is now recognized as a major mechanism of gene regulation. The interactions of transcription factors with chromatin structures have recently been observed to be highly dynamic, with residence times measured in seconds. Thus, the concept of static, multi-protein complexes forming at regulatory elements in the genome has been replaced by a new paradigm that envisages rapid and continuous exchange events with the template. These highly dynamic interactions are a property of both DNA-protein and protein-protein interactions and are inherent to every stage of the transcriptional response. In this review we discuss the dynamics of a nuclear receptor, and its transcriptional response in the chromatin context.

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.

Failure of terminal erythroid differentiation in EKLF-deficient mice is associated with cell cycle perturbation and reduced expression of E2F2

Erythroid Krüppel-like factor (EKLF) is a Krüppel-like transcription factor identified as a transcriptional activator and chromatin modifier in erythroid cells. EKLF-deficient (Eklf(-/-)) mice die at day 14.5 of gestation from severe anemia. In this study, we demonstrate that early progenitor cells fail to undergo terminal erythroid differentiation in Eklf(-/-) embryos. To discover potential EKLF target genes responsible for the failure of erythropoiesis, transcriptional profiling was performed with RNA from wild-type and Eklf(-/-) early erythroid progenitor cells. These analyses identified significant perturbation of a network of genes involved in cell cycle regulation, with the critical regulator of the cell cycle, E2f2, at a hub. E2f2 mRNA and protein levels were markedly decreased in Eklf(-/-) early erythroid progenitor cells, which showed a delay in the G(1)-to-S-phase transition. Chromatin immunoprecipitation analysis demonstrated EKLF occupancy at the proximal E2f2 promoter in vivo. Consistent with the role of EKLF as a chromatin modifier, EKLF binding sites in the E2f2 promoter were located in a region of EKLF-dependent DNase I sensitivity in early erythroid progenitor cells. We propose a model in which EKLF-dependent activation and modification of the E2f2 locus is required for cell cycle progression preceding terminal erythroid differentiation.

BEST1 expression in the retinal pigment epithelium is modulated by OTX family members

A number of genes preferentially expressed in the retinal pigment epithelium (RPE) are associated with retinal degenerative disease. One of these, BEST1, encodes bestrophin-1, a protein that when mutated causes Best macular dystrophy. As a model for RPE gene regulation, we have been studying the mechanisms that control BEST1 expression, and recently demonstrated that members of the MITF-TFE family modulate BEST1 transcription. The human BEST1 upstream region from -154 to +38 bp is sufficient to direct expression in the RPE, and positive-regulatory elements exist between -154 and -104 bp. Here, we show that the -154 to -104 bp region is necessary for RPE expression in transgenic mice and contains a predicted OTX-binding site (Site 1). Since another non-canonical OTX site (Site 2) is located nearby, we tested the function of these sites using BEST1 promoter/luciferase constructs by in vivo electroporation and found that mutation of both sites reduces promoter activity. Three OTX family proteins - OTX1, OTX2 and CRX - bound to both Sites 1 and 2 in vitro, and all of them increased BEST1 promoter activity. Surprisingly, we found that human and bovine RPE expressed not only OTX2 but also CRX, the CRX genomic region in bovine RPE was hypersensitive to DNase I, consistent with active transcription, and that both OTX2 and CRX bound to the BEST1 proximal promoter in vivo. These results demonstrate for the first time CRX expression in the RPE, and suggest that OTX2 and CRX may act as positive modulators of the BEST1 promoter in the RPE.

Chromatin profiling across the human tumour necrosis factor gene locus reveals a complex, cell type-specific landscape with novel regulatory elements

The TNF locus on chromosome 6p21 encodes a family of proteins with key roles in the immune response whose dysregulation leads to severe disease. Transcriptional regulation is important, with cell type and stimulus-specific enhancer complexes involving the proximal TNF promoter. We show how quantitative chromatin profiling across a 34 kb region spanning the TNF locus has allowed us to identify a number of novel DNase hypersensitive sites and characterize more distant regulatory elements. We demonstrate DNase hypersensitive sites corresponding to the lymphotoxin alpha (LTA) and tumour necrosis factor (TNF) promoter regions, a CpG island in exon 4 of lymphotoxin beta (LTB), the 3′ end of nuclear factor of kappa light polypeptide gene enhancer in B-cells inhibitor-like 1 (NFKBIL1) and 3.4 kb upstream of LTA. These sites co-localize to highly conserved DNA sequences and show evidence of cell type specificity when lymphoblastoid, Jurkat, U937, HeLa and HEK293T cell lines are analysed using Southern blotting. For Jurkat T cells, we define histone modifications across the locus. Peaks of acetylated histone H3 and H4, together with tri-methyl K4 of histone H3, correspond to hypersensitive sites, notably in exon 4 of LTB. We provide evidence of a functional role for an intergenic DNase I hypersensitive site distal to LTA in Jurkat cells based on reporter gene analysis, with evidence of recruitment of upstream stimulatory factors (USF) transcription factors.

Automated mapping of large-scale chromatin structure in ENCODE

Motivation: A recently developed DNaseI assay has given us our first genome-wide view of chromatin structure. In addition to cataloging DNaseI hypersensitive sites, these data allows us to more completely characterize overall features of chromatin accessibility. We employed a Bayesian hierarchical change-point model (CPM), a generalization of a hidden Markov Model (HMM), to characterize tiled microarray DNaseI sensitivity data available from the ENCODE project.

Results: Our analysis shows that the accessibility of chromatin to cleavage by DNaseI is well described by a four state model of local segments with each state described by a continuous mixture of Gaussian variables. The CPM produces a better fit to the observed data than the HMM. The large posterior probability for the four-state CPM suggests that the data falls naturally into four classes of regions, which we call major and minor DNaseI hypersensitive sites (DHSs), regions of intermediate sensitivity, and insensitive regions. These classes agree well with a model of chromatin in which local disruptions (DHSs) are concentrated within larger domains of intermediate sensitivity, the accessibility islands. The CPM assigns 92% of the bases within the ENCODE regions to the insensitive regions. The 5.8% of the bases that are in regions of intermediate sensitivity are clearly enriched in functional elements, including genes and activating histone modifications, while the remaining 2.2% of the bases in hypersensitive regions are very strongly enriched in these elements.

Availability: The CPM software is available upon request from the authors.

Contact:jstam@stamlab.org; noble@gs.washington.edu; Charles_Lawrence@brown.edu

Supplementary information:Supplementary data are available at Bioinformatics online. Source code is available at http://noble.gs.washington.edu/proj/segment.

Association study of the NEDD9 gene with the risk of developing Alzheimer's and Parkinson's disease

Alzheimer's disease (AD) and Parkinson's disease (PD), the two most common neurodegenerative disorders in the elderly, have been hypothesized to share genetic determinants. Recently, Li et al. proposed that a variant in the NEDD9 gene may be one of these common genetic factors. We attempted to confirm this initial observation by conducting an equivalent analysis in terms of pathologies and sample size. We genotyped the NEDD9 rs760678 SNP in three independent AD case-control studies (n = 3176) and two independent PD case-control studies (n = 1855). However, we failed to detect an association of this SNP with the risk of developing AD or PD, in any of these populations. In conclusion, these data indicate that the rs760678 SNP of the NEDD9 gene is at best a weak genetic determinant of AD or PD.

Long-range regulation of cytokine gene expression

In metazoans, transcription is regulated by promoters and additional elements, which may be located far from their target gene(s). Moreover, genes (including those encoding cytokines and cytokine receptors) are commonly clustered in the genome, providing the opportunity for the shared, competitive, or sequential use of regulatory elements. New techniques, discussed here, are generating an avalanche of high-resolution genome-wide data through which candidate regulatory elements have been identified in specific cell types (including T cells), their functions inferred, and their physical interactions in three-dimensional space demonstrated. As a result, a nearly comprehensive list of regulatory elements in the Th2 cytokine locus, a growing list of elements in the interferon-gamma gene locus, and maps of their three-dimensional interactions are now available, though much remains to be learned about the molecular mechanisms at play, the dynamics of these interactions, and their functional importance.

Interaction of the glucocorticoid receptor with the chromatin landscape

The generality and spectrum of chromatin-remodeling requirements for nuclear receptor function are unknown. We have characterized glucocorticoid receptor (GR) binding events and chromatin structural transitions across GR-induced or -repressed genes. This analysis reveals that GR binding invariably occurs at nuclease-accessible sites (DHS). A remarkable diversity of mechanisms, however, render these sites available for GR binding. Accessibility of the GR binding sites is either constitutive or hormone inducible. Within each category, some DHS sites require the Brg1-containing Swi/Snf complex, but others are Brg1 independent, implicating a different remodeling complex. The H2A.Z histone variant is highly enriched at both inducible and constitutive DHS sites and is subject to exchange during hormone activation. The DHS profile is highly cell specific, implicating cell-selective organization of the chromatin landscape as a critical determinant of tissue-selective receptor function. Furthermore, the widespread requirement for chromatin remodeling supports the recent hypothesis that the rapid exchange of receptor proteins occurs during nucleosome reorganization.

Detection and characterization of silencers and enhancer-blockers in the greater CFTR locus

Silencers and enhancer-blockers (EBs) are cis-acting, negative regulatory elements (NREs) that control interactions between promoters and enhancers. Although relatively uncharacterized in terms of biological mechanisms, these elements are likely to be abundant in the genome. We developed an experimental strategy to identify silencers and EBs using transient transfection assays. A known insulator and EB from the chicken beta-globin locus, cHS4, served as a control element for these assays. We examined 47 sequences from a 1.8-Mb region of human chromosome 7 for silencer and EB activities. The majority of functional elements displayed directional and promoter-specific activities. A limited number of sequences acted in a dual manner, as both silencers and EBs. We examined genomic data, epigenetic modifications, and sequence motifs within these regions. Strong silencer elements contained a novel CT-rich motif, often in multiple copies. Deletion of the motif from three regions caused a measurable loss of silencing ability in these sequences. Moreover, five duplicate occurrences of this motif were identified in the cHS4 insulator. These motifs provided an explanation for an uncharacterized silencing activity we measured in the insulator element. Overall, we identified 15 novel NREs, which contribute new insights into the prevalence and composition of sequences that negatively regulate gene expression.

Genome-wide analysis of protein-DNA interactions

The human genome is predominantly composed of nonprotein-coding sequences whose function remains largely undefined. A significant portion of the noncoding DNA is believed to serve as transcriptional regulatory elements that control gene expression in specific cell types at appropriate developmental stages. Identifying these regulatory sequences and determining the mechanisms by which they act present a great challenge in the postgenomic era. Previous investigations using genetic, molecular, and biochemical approaches have uncovered a large number of proteins involved in regulating transcription. Knowledge of the genomic locations of DNA binding for these proteins in the nucleus should define the identity and nature of the transcriptional regulatory sequences and reveal the gene regulatory networks in cells. Chromatin immunoprecipitation (ChIP) is a common method for detecting interactions between a protein and a DNA sequence in vivo. In recent years, this method has been combined with DNA microarrays and other high-throughput technologies to enable genome-wide identification of DNA-binding sites for various nuclear proteins. Here, we review recent advances in ChIP-based methods for genome-wide detection of protein-DNA interactions, and discuss their significance in enhancing our knowledge of the gene regulatory networks and epigenetic mechanisms in cells.

Oct4 dependence of chromatin structure within the extended Nanog locus in ES cells

Embryonic stem (ES) cells offer insight into early developmental fate decisions, and their controlled differentiation may yield vast regenerative potential. The molecular determinants supporting ES cell self-renewal are incompletely understood. The homeodomain proteins Nanog and Oct4 are essential for mouse ES cell self-renewal. Using a high-throughput approach, we discovered DNaseI hypersensitive sites and potential regulatory elements along a 160-kb region of the genome that includes GDF3, Dppa3, and Nanog. We analyzed gene expression, chromatin occupancy, and higher-order chromatin structure throughout this gene locus and found that expression of the reprogramming factor Oct4 is required to maintain its integrity.

Using genomics to study how chromatin influences gene expression

A postgenome challenge is to understand how the code in DNA is converted into the biological processes underlying various cell fates. By establishing the appropriate technical tools, we are moving from an era in which such questions have been asked by studying individual genes to one in which large domains, whole chromosomes, and the entire human genome can be investigated. These developments will allow us to study in parallel the transcriptional program and components of the epigenetic program (nuclear position, timing of replication, chromatin structure and modification, DNA methylation) to determine the hierarchy and order of events required to switch genes on and off during differentiation and development.

Identification of tissue-specific cis-regulatory modules based on interactions between transcription factors

BACKGROUND

Evolutionary conservation has been used successfully to help identify cis-acting DNA regions that are important in regulating tissue-specific gene expression. Motivated by increasing evidence that some DNA regulatory regions are not evolutionary conserved, we have developed an approach for cis-regulatory region identification that does not rely upon evolutionary sequence conservation.

RESULTS

The conservation-independent approach is based on an empirical potential energy between interacting transcription factors (TFs). In this analysis, the potential energy is defined as a function of the number of TF interactions in a genomic region and the strength of the interactions. By identifying sets of interacting TFs, the analysis locates regions enriched with the binding sites of these interacting TFs. We applied this approach to 30 human tissues and identified 6232 putative cis-regulatory modules (CRMs) regulating 2130 tissue-specific genes. Interestingly, some genes appear to be regulated by different CRMs in different tissues. Known regulatory regions are highly enriched in our predicted CRMs. In addition, DNase I hypersensitive sites, which tend to be associated with active regulatory regions, significantly overlap with the predicted CRMs, but not with more conserved regions. We also find that conserved and non-conserved CRMs regulate distinct gene groups. Conserved CRMs control more essential genes and genes involved in fundamental cellular activities such as transcription. In contrast, non-conserved CRMs, in general, regulate more non-essential genes, such as genes related to neural activity.

CONCLUSION

These results demonstrate that identifying relevant sets of binding motifs can help in the mapping of DNA regulatory regions, and suggest that non-conserved CRMs play an important role in gene regulation.

Detection of DNA structural motifs in functional genomic elements

The completion of the human genome project has fueled the search for regulatory elements by a variety of different approaches. Many successful analyses have focused on examining primary DNA sequence and/or chromatin structure. However, it has been difficult to detect common sequence motifs within the feature of chromatin structure most closely associated with regulatory elements, DNase I hypersensitive sites (DHSs). Considering just the nucleotide sequence and/or the chromatin structure of regulatory elements may neglect a critical feature of what is recognized by the regulatory machinery--DNA structure. We introduce a new computational method to detect common DNA structural motifs in a large collection of DHSs that are found in the ENCODE regions of the human genome. We show that DHSs have common DNA structural motifs that show no apparent sequence consensus. One such structural motif is much more highly enriched in experimentally identified DHSs that are in CpG islands and near transcription start sites (TSSs), compared to DHSs not in CpG islands and farther from TSSs, suggesting that DNA structural motifs may participate in the formation of functional regulatory elements. We propose that studies of the conservation of DNA structure, independent of sequence conservation, will provide new information about the link between the nucleotide sequence of a DNA molecule and its experimentally demonstrated function.

Tracking chromatin states using controlled DNase I treatment and real-time PCR

A novel approach to DNase I-sensitivity analysis was applied to examining genes of the spermatogenic pathway, reflective of the substantial morphological and genomic changes that occur during this program of differentiation. A new real-time PCR-based strategy that considers the nuances of response to nuclease treatment was used to assess the nuclease susceptibility through differentiation. Data analysis was automated with the K-Lab PCR algorithm, facilitating the rapid analysis of multiple samples while eliminating the subjectivity usually associated with C_t analyses. The utility of this assay and analytical paradigm as applied to nuclease-sensitivity mapping is presented.

Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project

We report the generation and analysis of functional data from multiple, diverse experiments performed on a targeted 1% of the human genome as part of the pilot phase of the ENCODE Project. These data have been further integrated and augmented by a number of evolutionary and computational analyses. Together, our results advance the collective knowledge about human genome function in several major areas. First, our studies provide convincing evidence that the genome is pervasively transcribed, such that the majority of its bases can be found in primary transcripts, including non-protein-coding transcripts, and those that extensively overlap one another. Second, systematic examination of transcriptional regulation has yielded new understanding about transcription start sites, including their relationship to specific regulatory sequences and features of chromatin accessibility and histone modification. Third, a more sophisticated view of chromatin structure has emerged, including its inter-relationship with DNA replication and transcriptional regulation. Finally, integration of these new sources of information, in particular with respect to mammalian evolution based on inter- and intra-species sequence comparisons, has yielded new mechanistic and evolutionary insights concerning the functional landscape of the human genome. Together, these studies are defining a path for pursuit of a more comprehensive characterization of human genome function.

Comprehensive epigenetic profiling identifies multiple distal regulatory elements directing transcription of the gene encoding interferon-gamma

Unlike the well defined T helper type 2 cytokine locus, little is known about the regulatory elements that govern the expression of Ifng, which encodes the 'signature' T helper type 1 cytokine interferon-gamma. Here our evolutionary analysis showed that the mouse Ifng locus diverged from the ancestral locus as a result of structural rearrangements producing deletion of the neighboring gene encoding interleukin 26 and disrupting synteny 57 kilobases upstream of Ifng. Proximal to that disruption, we identified by high-resolution mapping many regions with CD4+ T cell subset-specific epigenetic modifications. A subset of those regions represented enhancers, whereas others blocked the activity of upstream enhancers or insulated Ifng from neighboring chromatin. Our findings suggest that proper expression of Ifng is maintained through the collective action of multiple distal regulatory elements present in a region of about 100 kilobases flanking Ifng.

Real-time PCR mapping of DNaseI-hypersensitive sites using a novel ligation-mediated amplification technique

Mapping sites within the genome that are hypersensitive to digestion with DNaseI is an important method for identifying DNA elements that regulate transcription. The standard approach to locating these DNaseI-hypersensitive sites (DHSs) has been to use Southern blotting techniques, although we, and others, have recently published alternative methods using a range of technologies including high-throughput sequencing and genomic array tiling paths. In this article, we describe a novel protocol to use real-time PCR to map DHS. Advantages of the technique reported here include the small cell numbers required for each analysis, rapid, relatively low-cost experiments with minimal need for specialist equipment. Presented examples include comparative DHS mapping of known TAL1/SCL regulatory elements between human embryonic stem cells and K562 cells.

Mapping and characterization of DNase I hypersensitive sites in Arabidopsis chromatin

Recent genome-wide analyses of yeast and human chromatin revealed the widespread prevalence of DNase I hypersensitive sites (DNase I HSs) at gene regulatory regions with possible roles in eukaryotic gene regulation. The presence of DNase I HSs in plants has been described for only a few genes, and we analyzed the chromatin structure of an 80 kb genomic region containing 30 variably expressed genes by DNase I sensitivity assay at 500 bp resolution in Arabidopsis. Distinct DNase I HSs were found at the 5' and/or 3' ends of most genes irrespective of their expression levels. Further analysis of well-characterized genes showed that the DNase I HSs occurred near cis-regulatory elements in the promoters of these genes. Upon transcriptional activation of a heat-inducible gene, the DNase I HS was extended into the vicinity of a cis-element and adjacent TATA element in the promoter. Concomitant with this change in DNase I HS, histones were acetylated, removed from the promoter, and a transcription activator bound to this cis-element. These results suggest that the DNase I HSs participate in the transcriptional regulation of Arabidopsis genes by enhancing the access of chromatin remodeling factors and/or transcription factors to their target sites as seen in yeast and human chromatin.

HS-48 alone has no enhancement role on the expression of human α-globin gene cluster

To investigate the in vivo function of the newly defined DNase I hypersensitive site HS-48 on the whole human alpha-globin gene cluster, the region containing all the other known 5 hypersensitive sites HS-4 to HS-40 was deleted from a 117 kb bacterial artificial chromosome clone bearing the whole human alpha-globin gene cluster. Transgenic mice were generated from this construct. The RNase protection assays showed that with HS-48 left and all the other 5 hypersensitive sites deleted, the expression of human alpha-like globin genes was completely silenced in embryonic, fetal and adult stages in all tissues. This finding indicates that HS-48 alone has no enhancer activity on the expression of human alpha-like globin genes, and that the region of HS-4 to HS-40 already contains all the upstream cis-elements needed for regulating human alpha-like globin genes.

Association analysis for large-scale gene set data

High-throughput experiments in biology often produce sets of genes of potential interests. Some of those gene sets might be of considerable size. Therefore, computer-assisted analysis is necessary for the biological interpretation of the gene sets, and for creating working hypotheses, which can be tested experimentally. One obvious way to analyze gene set data is to associate the genes with a particular biological feature, for example, a given pathway. Statistical analysis could be used to evaluate if a gene set is truly associated with a feature. Over the past few years many tools that perform such analysis have been created. In this chapter, using WebGestalt as an example, it will be explained in detail how to associate gene sets with functional annotations, pathways, publication records, and protein domains.

Genome-wide prediction of conserved and nonconserved enhancers by histone acetylation patterns

Comparative genomic studies have been useful in identifying transcriptional regulatory elements in higher eukaryotic genomes, but many important regulatory elements cannot be detected by such analyses due to evolutionary variations and alignment tool limitations. Therefore, in this study we exploit the highly conserved nature of epigenetic modifications to identify potential transcriptional enhancers. By using a high-resolution genome-wide mapping technique, which combines the chromatin immunoprecipitation and serial analysis of gene expression assays, we have recently determined the distribution of lysine 9/14-diacetylated histone H3 in human T cells. We showed the existence of 46,813 regions with clusters of histone acetylation, termed histone acetylation islands, some of which correspond to known transcriptional regulatory elements. In the present study, we find that 4679 sequences conserved between human and pufferfish coincide with histone acetylation islands, and random sampling shows that 33% (13/39) of these can function as transcriptional enhancers in human Jurkat T cells. In addition, by comparing the human histone acetylation island sequences with mouse genome sequences, we find that despite the conservation of many of these regions between these species, 21,855 of these sequences are not conserved. Furthermore, we demonstrate that about 50% (26/51) of these nonconserved sequences have enhancer activity in Jurkat cells, and that many of the orthologous mouse sequences also have enhancer activity in addition to conserved epigenetic modification patterns in mouse T-cell chromatin. Therefore, by combining epigenetic modification and sequence data, we have established a novel genome-wide method for identifying regulatory elements not discernable by comparative genomics alone.

A cell-type-specific transcriptional network required for estrogen regulation of cyclin D1 and cell cycle progression in breast cancer

Estrogen stimulates the proliferation of the most common type of human breast cancer that expresses estrogen receptor alpha (ERalpha) through the activation of the cyclin D1 (CCND1) oncogene. However, our knowledge of ERalpha transcriptional mechanisms remains limited. Hence, it is still elusive why ERalpha ectopically expressed in ER-negative breast cancer cells (BCC) is functional on ectopic reporter constructs but lacks activity on many endogenous target genes, including CCND1. Here, we show that estradiol (E2) stimulation of CCND1 expression in BCC depends on a novel cell-type-specific enhancer downstream from the CCND1 coding region, which is the primary ERalpha recruitment site in estrogen-responsive cells. The pioneer factor FoxA1 is specifically required for the active chromatin state of this enhancer and as such is crucial for both CCND1 expression and subsequent cell cycle progression. Interestingly, even in BCC, CCND1 levels and proliferation are tightly controlled by E2 through the establishment of a negative feedforward loop involving the induction of NFIC, a putative tumor suppressor capable of directly repressing CCND1 transcription. Taken together, our results reveal an estrogen-regulated combinatorial network including cell-specific cis- and trans-regulators of CCND1 expression where ERalpha collaborates with other transcription factors associated with the ER-positive breast cancer phenotype, including FoxA1 and NFIC.

Unraveling transcription regulatory networks by protein-DNA and protein-protein interaction mapping

Metazoan genomes contain thousands of protein-coding and noncoding RNA genes, most of which are differentially expressed, i.e., at different locations or at different times during development, function, or pathology of the organism. Differential gene expression is achieved in part by the action of regulatory transcription factors (TFs) that bind to cis-regulatory elements that are often located in or near their target genes. Each TF likely regulates many targets in the context of intricate transcription regulatory networks. Up to 10% of a genome may encode TFs, but only a handful of these have been studied in detail. Here, I will discuss the different steps involved in the mapping and analysis of transcription regulatory networks, including the identification of network nodes (TFs and their target sequences) and edges (TF-TF dimers and TF-DNA target interactions), integration with other data types, and network properties and emerging principles that provide insights into differential gene expression.

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.

Identifying gene regulatory elements by genomic microarray mapping of DNaseI hypersensitive sites

The identification of cis-regulatory elements is central to understanding gene transcription. Hypersensitivity of cis-regulatory elements to digestion with DNaseI remains the gold-standard approach to locating such elements. Traditional methods used to identify DNaseI hypersensitive sites are cumbersome and can only be applied to short stretches of DNA at defined locations. Here we report the development of a novel genomic array-based approach to DNaseI hypersensitive site mapping (ADHM) that permits precise, large-scale identification of such sites from as few as 5 million cells. Using ADHM we identified all previously recognized hematopoietic regulatory elements across 200 kb of the mouse T-cell acute lymphocytic leukemia-1 (Tal1) locus, and, in addition, identified two novel elements within the locus, which show transcriptional regulatory activity. We further validated the ADHM protocol by mapping the DNaseI hypersensitive sites across 250 kb of the human TAL1 locus in CD34+ primary stem/progenitor cells and K562 cells and by mapping the previously known DNaseI hypersensitive sites across 240 kb of the human alpha-globin locus in K562 cells. ADHM provides a powerful approach to identifying DNaseI hypersensitive sites across large genomic regions.

Genome-scale mapping of DNase I sensitivity in vivo using tiling DNA microarrays

Localized accessibility of critical DNA sequences to the regulatory machinery is a key requirement for regulation of human genes. Here we describe a high-resolution, genome-scale approach for quantifying chromatin accessibility by measuring DNase I sensitivity as a continuous function of genome position using tiling DNA microarrays (DNase-array). We demonstrate this approach across 1% ( approximately 30 Mb) of the human genome, wherein we localized 2,690 classical DNase I hypersensitive sites with high sensitivity and specificity, and also mapped larger-scale patterns of chromatin architecture. DNase I hypersensitive sites exhibit marked aggregation around transcriptional start sites (TSSs), though the majority mark nonpromoter functional elements. We also developed a computational approach for visualizing higher-order features of chromatin structure. This revealed that human chromatin organization is dominated by large (100-500 kb) 'superclusters' of DNase I hypersensitive sites, which encompass both gene-rich and gene-poor regions. DNase-array is a powerful and straightforward approach for systematic exposition of the cis-regulatory architecture of complex genomes.

DNase-chip: a high-resolution method to identify DNase I hypersensitive sites using tiled microarrays

Mapping DNase I hypersensitive sites is an accurate method of identifying the location of gene regulatory elements, including promoters, enhancers, silencers and locus control regions. Although Southern blots are the traditional method of identifying DNase I hypersensitive sites, the conventional manual method is not readily scalable to studying large chromosomal regions, much less the entire genome. Here we describe DNase-chip, an approach that can rapidly identify DNase I hypersensitive sites for any region of interest, or potentially for the entire genome, by using tiled microarrays. We used DNase-chip to identify DNase I hypersensitive sites accurately from a representative 1% of the human genome in both primary and immortalized cell types. We found that although most DNase I hypersensitive sites were present in both cell types studied, some of them were cell-type specific. This method can be applied globally or in a targeted fashion to any tissue from any species with a sequenced genome.

The emerging science of epigenomics

One of the most exciting frontiers in both epigenetics and genome sciences is the new field of epigenomics. This new discipline promises novel insights into the genome because of its potential to detect quantitative alterations, multiplex modifications and regulatory sequences outside of genes. A number of new epigenomic strategies are emerging to exploit microarray formats with varying substrate choice, pre-processing and data analysis. These approaches are designed to detect large numbers of variations in DNA methylation and chromatin modification. Many groups are joining forces toward developing an organized Human Epigenome Project to exploit these new technologies to better understand the basis of normal development and human disease.

Genome-wide computational prediction of transcriptional regulatory modules reveals new insights into human gene expression

The identification of regulatory regions is one of the most important and challenging problems toward the functional annotation of the human genome. In higher eukaryotes, transcription-factor (TF) binding sites are often organized in clusters called cis-regulatory modules (CRM). While the prediction of individual TF-binding sites is a notoriously difficult problem, CRM prediction has proven to be somewhat more reliable. Starting from a set of predicted binding sites for more than 200 TF families documented in Transfac, we describe an algorithm relying on the principle that CRMs generally contain several phylogenetically conserved binding sites for a few different TFs. The method allows the prediction of more than 118,000 CRMs within the human genome. A subset of these is shown to be bound in vivo by TFs using ChIP-chip. Their analysis reveals, among other things, that CRM density varies widely across the genome, with CRM-rich regions often being located near genes encoding transcription factors involved in development. Predicted CRMs show a surprising enrichment near the 3' end of genes and in regions far from genes. We document the tendency for certain TFs to bind modules located in specific regions with respect to their target genes and identify TFs likely to be involved in tissue-specific regulation. The set of predicted CRMs, which is made available as a public database called PReMod (http://genomequebec.mcgill.ca/PReMod), will help analyze regulatory mechanisms in specific biological systems.

AUTOMATED VALIDATION OF POLYMERASE CHAIN REACTION AMPLICON MELTING CURVES

The polymerase chain reaction (PCR) is a fundamental tool of molecular biology. Quantitative PCR is the gold-standard methodology for determination of DNA copy numbers, quantitating transcription, and numerous other applications. A major barrier to large-scale application of PCR for quantitative genomic analyses is the current requirement for manual validation of individual PCRs to ensure generation of a single product. This typically requires visual inspection either of gel electrophoreses or temperature dissociation ("melting") curves of individual PCRs — a time-consuming and costly process.

Here we describe a robust computational solution to this fundamental problem. Using a training set of 10 080 reactions comprising multiple quantitative PCRs from each of 1728 unique human genomic amplicons, we developed a support vector machine classifier capable of discriminating single-product PCRs with better than 99% accuracy. This approach has broad utility, and eliminates a major bottleneck to widespread application of PCR for high-throughput genomic applications.

A reliable method to display authentic DNase I hypersensitive sites at long-ranges in single-copy genes from large genomes

The study of eukaryotic gene transcription depends on methods to discover distal cis-acting control sequences. Comparative bioinformatics is one powerful strategy to reveal these domains, but still requires conventional wet-bench techniques to elucidate their specificity and function. The DNase I hypersensitivity assay (DHA) is also a method to identify regulatory domains, but can also suggest their function. Technically however, the classical DHA is constrained to mapping gene loci in small increments of approximately 20 kb. This limitation hinders efficient and comprehensive analysis of distal gene regions. Here, we report an improved method termed mega-DHA that extends the range of existing DHAs to facilitate assaying intervals that approach 100 kb. We demonstrate its feasibility for efficient analysis of single-copy genes within a large and complex genome by assaying 230 kb of the human ADAMTS14-perforin-paladin gene cluster in four experiments. The results identify distinct networks of regulatory domains specific to expression of perforin and its two neighboring genes.

Automated validation of polymerase chain reactions using amplicon melting curves

PCR, the polymerase chain reaction, is a fundamental tool of molecular biology. Quantitative PCR is the gold-standard methodology for determination of DNA copy numbers, quantitating transcription, and numerous other applications. A major barrier to large-scale application of PCR for quantitative genomic analyses is the current requirement for manual validation of individual PCR reactions to ensure generation of a single product. This typically requires visual inspection either of gel electrophoreses or temperature dissociation ("melting") curves of individual PCR reactions - a time-consuming and costly process. Here we describe a robust computational solution to this fundamental problem. Using a training set of 10,080 reactions comprising multiple quantitative PCR reactions from each of 1,728 unique human genomic amplicons, we developed a support vector machine classifier capable of discriminating single-product PCR reactions with better than 99% accuracy. This approach has broad utility, and eliminates a major bottleneck to widespread application of PCR for high-throughput genomic applications.