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

Developmental roles of 21 Drosophila transcription factors are determined by quantitative differences in binding to an overlapping set of thousands of genomic regions

Distinct developmental fates in Drosophila melanogaster are specified by quantitative differences in transcription factor occupancy on a common set of bound regions.

Background

We previously established that six sequence-specific transcription factors that initiate anterior/posterior patterning in Drosophila bind to overlapping sets of thousands of genomic regions in blastoderm embryos. While regions bound at high levels include known and probable functional targets, more poorly bound regions are preferentially associated with housekeeping genes and/or genes not transcribed in the blastoderm, and are frequently found in protein coding sequences or in less conserved non-coding DNA, suggesting that many are likely non-functional.

Results

Here we show that an additional 15 transcription factors that regulate other aspects of embryo patterning show a similar quantitative continuum of function and binding to thousands of genomic regions in vivo. Collectively, the 21 regulators show a surprisingly high overlap in the regions they bind given that they belong to 11 DNA binding domain families, specify distinct developmental fates, and can act via different cis-regulatory modules. We demonstrate, however, that quantitative differences in relative levels of binding to shared targets correlate with the known biological and transcriptional regulatory specificities of these factors.

Conclusions

It is likely that the overlap in binding of biochemically and functionally unrelated transcription factors arises from the high concentrations of these proteins in nuclei, which, coupled with their broad DNA binding specificities, directs them to regions of open chromatin. We suggest that most animal transcription factors will be found to show a similar broad overlapping pattern of binding in vivo, with specificity achieved by modulating the amount, rather than the identity, of bound factor.

Genome-wide views of chromatin structure

Eukaryotic genomes are packaged into a nucleoprotein complex known as chromatin, which affects most processes that occur on DNA. Along with genetic and biochemical studies of resident chromatin proteins and their modifying enzymes, mapping of chromatin structure in vivo is one of the main pillars in our understanding of how chromatin relates to cellular processes. In this review, we discuss the use of genomic technologies to characterize chromatin structure in vivo, with a focus on data from budding yeast and humans. The picture emerging from these studies is the detailed chromatin structure of a typical gene, where the typical behavior gives insight into the mechanisms and deep rules that establish chromatin structure. Important deviation from the archetype is also observed, usually as a consequence of unique regulatory mechanisms at special genomic loci. Chromatin structure shows substantial conservation from yeast to humans, but mammalian chromatin has additional layers of complexity that likely relate to the requirements of multicellularity such as the need to establish faithful gene regulatory mechanisms for cell differentiation.

Experimental advances in the characterization of metazoan gene regulatory networks

Gene regulatory networks (GRNs) play a vital role in metazoan development and function, and deregulation of these networks is often implicated in disease. GRNs depict the dynamic interactions between genomic and regulatory state components. The genomic components comprise genes and their associated cis-regulatory elements. The regulatory state components consist primarily of transcriptional complexes that bind the latter elements. With the availability of complete genome sequences, several approaches have recently been developed which promise to significantly enhance our ability to identify either the genomic or regulatory state components, or the interactions between these two. In this review, I provide an in-depth overview of these approaches and detail how each contributes to a more comprehensive understanding of GRN composition and function.

Endonuclease-sensitive regions of human spermatozoal chromatin are highly enriched in promoter and CTCF binding sequences

During the haploid phase of mammalian spermatogenesis, nucleosomal chromatin is ultimately repackaged by small, highly basic protamines to generate an extremely compact, toroidal chromatin architecture that is critical to normal spermatozoal function. In common with several species, however, the human spermatozoon retains a small proportion of its chromatin packaged in nucleosomes. As nucleosomal chromatin in spermatozoa is structurally more open than protamine-packaged chromatin, we considered it likely to be more accessible to exogenously applied endonucleases. Accordingly, we have used this premise to identify a population of endonuclease-sensitive DNA sequences in human and murine spermatozoa. Our results show unequivocally that, in contrast to the endonuclease-resistant sperm chromatin packaged by protamines, regions of increased endonuclease sensitivity are closely associated with gene regulatory regions, including many promoter sequences and sequences recognized by CCCTC-binding factor (CTCF). Similar differential packaging of promoters is observed in the spermatozoal chromatin of both mouse and man. These observations imply the existence of epigenetic marks that distinguish gene regulatory regions in male germ cells and prevent their repackaging by protamines during spermiogenesis. The ontology of genes under the control of endonuclease-sensitive regulatory regions implies a role for this phenomenon in subsequent embryonic development.

Improved ChIP-chip analysis by a mixture model approach

Background

Microarray analysis of immunoprecipitated chromatin (ChIP-chip) has evolved from a novel technique to a standard approach for the systematic study of protein-DNA interactions. In ChIP-chip, sites of protein-DNA interactions are identified by signals from the hybridization of selected DNA to tiled oligomers and are graphically represented as peaks. Most existing methods were designed for the identification of relatively sparse peaks, in the presence of replicates.

Results

We propose a data normalization method and a statistical method for peak identification from ChIP-chip data based on a mixture model approach. In contrast to many existing methods, including methods that also employ mixture model approaches, our method is more flexible by imposing less restrictive assumptions and allowing a relatively large proportion of peak regions. In addition, our method does not require experimental replicates and is computationally efficient. We compared the performance of our method with several representative existing methods on three datasets, including a spike-in dataset. These comparisons demonstrate that our approach is more robust and has comparable or higher power than the other methods, especially in the context of abundant peak regions.

Conclusion

Our data normalization and peak detection methods have improved performance to detect peak regions in ChIP-chip data.

Transcriptional dysregulation in NIPBL and cohesin mutant human cells

Genome-wide studies using cells from patients with Cornelia de Lange Syndrome reveal a role for cohesin in regulating gene expression in human cells.

Cohesin regulates sister chromatid cohesion during the mitotic cell cycle with Nipped-B-Like (NIPBL) facilitating its loading and unloading. In addition to this canonical role, cohesin has also been demonstrated to play a critical role in regulation of gene expression in nondividing cells. Heterozygous mutations in the cohesin regulator NIPBL or cohesin structural components SMC1A and SMC3 result in the multisystem developmental disorder Cornelia de Lange Syndrome (CdLS). Genome-wide assessment of transcription in 16 mutant cell lines from severely affected CdLS probands has identified a unique profile of dysregulated gene expression that was validated in an additional 101 samples and correlates with phenotypic severity. This profile could serve as a diagnostic and classification tool. Cohesin binding analysis demonstrates a preference for intergenic regions suggesting a cis-regulatory function mimicking that of a boundary/insulator interacting protein. However, the binding sites are enriched within the promoter regions of the dysregulated genes and are significantly decreased in CdLS proband, indicating an alternative role of cohesin as a transcription factor.

Appropriate segregation of chromosomes to daughter cells depends upon proper cohesion of sister chromatids during mitosis. The multiprotein cohesin complex and its regulators are key factors in this process. Intriguingly, recent work has shown that the cohesin complex also has other cellular roles, including a role in regulating gene expression. Additionally, mutations in cohesin structural and regulatory components have been linked to human multisystem developmental disorders such as Cornelia de Lange Syndrome (CdLS), but the role cohesin is playing in the pathogenesis of this disorder is unknown. To define the role that cohesin plays in regulating gene expression in human cells, we analyzed gene expression and genome-wide cohesin binding patterns in cells from normal subjects and from CdLS probands with mutations in the cohesin regulator NIPBL or in the cohesin structural component SMC1A. We found a strikingly conserved pattern of gene dysregulation in these different cell lines that correlates with disease severity and a significant correlation between gene dysregulation and cohesin binding around misexpressed genes. The observed pattern of binding and misexpression is consistent with cohesin having a putative role as a boundary/insulator interacting protein or transcription factor, the activity of which is disrupted in CdLS probands.

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.

Global mapping of protein-DNA interactions in vivo by digital genomic footprinting

The orchestrated binding of transcriptional activators and repressors to specific DNA sequences in the context of chromatin defines the regulatory program of eukaryotic genomes. We developed a digital approach to assay regulatory protein occupancy on genomic DNA in vivo by dense mapping of individual DNase I cleavages from intact nuclei using massively parallel DNA sequencing. Analysis of > 23 million cleavages across the Saccharomyces cerevisiae genome revealed thousands of protected regulatory protein footprints, enabling de novo derivation of factor binding motifs as well as the identification of hundreds of novel binding sites for major regulators. We observed striking correspondence between nucleotide-level DNase I cleavage patterns and protein-DNA interactions determined by crystallography. The data also yielded a detailed view of larger chromatin features including positioned nucleosomes flanking factor binding regions. Digital genomic footprinting provides a powerful approach to delineate the cis-regulatory framework of any organism with an available genome sequence.

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.

Genome-scale identification of Caenorhabditis elegans regulatory elements by tiling-array mapping of DNase I hypersensitive sites

Background

A major goal of post-genomics research is the integrated analysis of genes, regulatory elements and the chromatin architecture on a genome-wide scale. Mapping DNase I hypersensitive sites within the nuclear chromatin is a powerful and well-established method of identifying regulatory element candidates.

Results

Here, we report the first genome-wide analysis of DNase I hypersensitive sites (DHSs) in Caenorhabditis elegans. The data was obtained by hybridizing DNase I-treated and end-captured material from young adult worms to a high-resolution tiling microarray. The data show that C. elegans DHSs were significantly enriched within intergenic regions located 2 kb upstream and downstream of coding genes, and also that a considerable fraction of all DHSs mapped to intergenic positions distant to annotated coding genes. Annotated transcribed loci were generally depleted in DHSs relative to intergenic regions, but DHSs were nonetheless enriched in coding exons and UTRs, whereas introns were significantly depleted in DHSs. Many DHSs appeared to be associated with annotated non-coding RNAs and recently detected transcripts of unknown function. It has been reported that nematode highly conserved non-coding elements were associated with cis-regulatory elements, and we also found that DHSs, particularly distal intergenic DHSs, were significantly enriched in regions that were conserved between the C. elegans and C. briggsae genomes.

Conclusion

We describe the first genome-wide analysis of C. elegans DHSs, and show that the distribution of DHSs is strongly associated with functional elements in the genome.

Genome-wide analyses of recombination prone regions predict role of DNA structural motif in recombination

HapMap findings reveal surprisingly asymmetric distribution of recombinogenic regions. Short recombinogenic regions (hotspots) are interspersed between large relatively non-recombinogenic regions. This raises the interesting possibility of DNA sequence and/or other cis- elements as determinants of recombination. We hypothesized the involvement of non-canonical sequences that can result in local non-B DNA structures and tested this using the G-quadruplex DNA as a model. G-quadruplex or G4 DNA is a unique form of four-stranded non-B DNA structure that engages certain G-rich sequences, presence of such motifs has been noted within telomeres. In support of this hypothesis, genome-wide computational analyses presented here reveal enrichment of potential G4 (PG4) DNA forming sequences within 25618 human hotspots relative to 9290 coldspots (p<0.0001). Furthermore, co-occurrence of PG4 DNA within several short sequence elements that are associated with recombinogenic regions was found to be significantly more than randomly expected. Interestingly, analyses of more than 50 DNA binding factors revealed that co-occurrence of PG4 DNA with target DNA binding sites of transcription factors c-Rel, NF-kappa B (p50 and p65) and Evi-1 was significantly enriched in recombination-prone regions. These observations support involvement of G4 DNA in recombination, predicting a functional model that is consistent with duplex-strand separation induced by formation of G4 motifs in supercoiled DNA and/or when assisted by other cellular factors.

Novel insights into adipogenesis from omics data

Obesity, the excess accumulation of adipose tissue, is one of the most pressing health problems in both the Western world and in developing countries. Adipose tissue growth results from two processes: the increase in number of adipocytes (hyperplasia) that develop from precursor cells, and the growth of individual fat cells (hypertrophy) due to incorporation of triglycerides. Adipogenesis, the process of fat cell development, has been extensively studied using various cell and animal models. While these studies pointed out a number of key factors involved in adipogenesis, the list of molecular components is far from complete. The advance of high-throughput technologies has sparked many experimental studies aimed at the identification of novel molecular components regulating adipogenesis. This paper examines the results of recent studies on adipogenesis using high-throughput technologies. Specifically, it provides an overview of studies employing microarrays for gene expression profiling and studies using gel based and non-gel based proteomics as well as a chromatin immunoprecipitation followed by microarray analysis (ChIP-chip) or sequencing (ChIP-seq). Due to the maturity of the technology, the bulk of the available data was generated using microarrays. Therefore these data sets were not only reviewed but also underwent meta analysis. The review also shows that large-scale omics technologies in conjunction with sophisticated bioinformatics analyses can provide not only a list of novel players, but also a global view on biological processes and molecular networks. Finally, developing technologies and computational challenges associated with the data analyses are highlighted, and an outlook on the questions not previously addressed is provided.

The Role of Nuclear Matrix Attachment Regions in Plants

Regions of DNA that bind to the nuclear matrix, or nucleoskeleton, are known as Matrix Attachment Regions (MARs). MARs are thought to play an important role in higher-order structure and chromatin organization within the nucleus. MARs are also thought to act as boundaries of chromosomal domains that act to separate regions of gene-rich, decondensed euchromatin from highly repetitive, condensed heterochromatin. Herein I will present evidence that MARs do indeed act as domain boundaries and can prevent the spread of silencing into active genes. Many fundamental questions remain unanswered about how MARs function in the nucleus. New findings in epigenetics indicate that MARs may also play an important role in the organization of genes and the eventual transport of their mRNAs through the nuclear pore.

Genome-wide profiling of salt fractions maps physical properties of chromatin

We applied genome-wide profiling to successive salt-extracted fractions of micrococcal nuclease-treated Drosophila chromatin. Chromatin fractions extracted with 80 mM or 150 mM NaCl after digestion contain predominantly mononucleosomes and represent classical "active" chromatin. Profiles of these low-salt soluble fractions display phased nucleosomes over transcriptionally active genes that are locally depleted of histone H3.3 and correspond closely to profiles of histone H2Av (H2A.Z) and RNA polymerase II. This correspondence suggests that transcription can result in loss of H3.3+H2Av nucleosomes and generate low-salt soluble nucleosomes. Nearly quantitative recovery of chromatin is obtained with 600 mM NaCl; however, the remaining insoluble chromatin is enriched in actively transcribed regions. Salt-insoluble chromatin likely represents oligonucleosomes that are attached to large protein complexes. Both low-salt extracted and insoluble chromatin are rich in sequences that correspond to epigenetic regulatory elements genome-wide. The presence of active chromatin at both extremes of salt solubility suggests that these salt fractions capture bound and unbound intermediates in active processes, thus providing a simple, powerful strategy for mapping epigenome dynamics.

Genome-scale validation of deep-sequencing libraries

Chromatin immunoprecipitation followed by high-throughput (HTP) sequencing (ChIP-seq) is a powerful tool to establish protein-DNA interactions genome-wide. The primary limitation of its broad application at present is the often-limited access to sequencers. Here we report a protocol, Mab-seq, that generates genome-scale quality evaluations for nucleic acid libraries intended for deep-sequencing. We show how commercially available genomic microarrays can be used to maximize the efficiency of library creation and quickly generate reliable preliminary data on a chromosomal scale in advance of deep sequencing. We also exploit this technique to compare enriched regions identified using microarrays with those identified by sequencing, demonstrating that they agree on a core set of clearly identified enriched regions, while characterizing the additional enriched regions identifiable using HTP sequencing.

Chromatin domains in higher eukaryotes: insights from genome-wide mapping studies

In genomes of higher eukaryotes, adjacent genes often show coordinated regulation of their expression. Compartmentalization of multiple neighboring genes into a shared chromatin environment can facilitate this coordinated expression. New mapping techniques have begun to reveal that such multigene chromatin domains are a common feature of fly and mammalian genomes. Many different types of chromatin domains have been identified based on the genomic binding patterns of various proteins and histone modifications. In addition, maps of genome-nuclear lamina associations and of looping interactions between loci provide the first systematic views of the three-dimensional folding of interphase chromosomes. These genome-wide datasets uncover new architectural principles of eukaryotic genomes and indicate that multigene chromatin domains are prevalent and important regulatory units.

Efficient algorithms for the computational design of optimal tiling arrays

The representation of a genome by oligonucleotide probes is a prerequisite for the analysis of many of its basic properties, such as transcription factor binding sites, chromosomal breakpoints, gene expression of known genes and detection of novel genes, in particular those coding for small RNAs. An ideal representation would consist of a high density set of oligonucleotides with similar melting temperatures that do not cross-hybridize with other regions of the genome and are equidistantly spaced. The implementation of such design is typically called a tiling array or genome array. We formulate the minimal cost tiling path problem for the selection of oligonucleotides from a set of candidates. Computing the selection of probes requires multi-criterion optimization, which we cast into a shortest path problem. Standard algorithms running in linear time allow us to compute globally optimal tiling paths from millions of candidate oligonucleotides on a standard desktop computer for most problem variants. The solutions to this multi-criterion optimization are spatially adaptive to the problem instance. Our formulation incorporates experimental constraints with respect to specific regions of interest and trade offs between hybridization parameters, probe quality and tiling density easily. A web application is available at http://tileomatic.org.

Nucleosome positioning in human HOX gene clusters

The distribution of nucleosomes along the genome is a significant aspect of chromatin structure and is thought to influence gene regulation through modulation of DNA accessibility. However, properties of nucleosome organization remain poorly understood, particularly in mammalian genomes. Toward this goal we used tiled microarrays to identify stable nucleosome positions along the HOX gene clusters in human cell lines. We show that nucleosome positions exhibit sequence properties and long-range organization that are different from those characterized in other organisms. Despite overall variability of internucleosome distances, specific loci contain regular nucleosomal arrays with 195-bp periodicity. Moreover, such arrays tend to occur preferentially toward the 3' ends of genes. Through comparison of different cell lines, we find that active transcription is correlated with increased positioning of nucleosomes, suggesting an unexpected role for transcription in the establishment of well-positioned nucleosomes.

SigWin-detector: a Grid-enabled workflow for discovering enriched windows of genomic features related to DNA sequences

Background

Chromosome location is often used as a scaffold to organize genomic information in both the living cell and molecular biological research. Thus, ever-increasing amounts of data about genomic features are stored in public databases and can be readily visualized by genome browsers. To perform in silico experimentation conveniently with this genomics data, biologists need tools to process and compare datasets routinely and explore the obtained results interactively. The complexity of such experimentation requires these tools to be based on an e-Science approach, hence generic, modular, and reusable. A virtual laboratory environment with workflows, workflow management systems, and Grid computation are therefore essential.

Findings

Here we apply an e-Science approach to develop SigWin-detector, a workflow-based tool that can detect significantly enriched windows of (genomic) features in a (DNA) sequence in a fast and reproducible way. For proof-of-principle, we utilize a biological use case to detect regions of increased and decreased gene expression (RIDGEs and anti-RIDGEs) in human transcriptome maps. We improved the original method for RIDGE detection by replacing the costly step of estimation by random sampling with a faster analytical formula for computing the distribution of the null hypothesis being tested and by developing a new algorithm for computing moving medians. SigWin-detector was developed using the WS-VLAM workflow management system and consists of several reusable modules that are linked together in a basic workflow. The configuration of this basic workflow can be adapted to satisfy the requirements of the specific in silico experiment.

Conclusion

As we show with the results from analyses in the biological use case on RIDGEs, SigWin-detector is an efficient and reusable Grid-based tool for discovering windows enriched for features of a particular type in any sequence of values. Thus, SigWin-detector provides the proof-of-principle for the modular e-Science based concept of integrative bioinformatics experimentation.

Maps of cis-Regulatory Nodes in Megabase Long Genome Segments are an Inevitable Intermediate Step Toward Whole Genome Functional Mapping

The availability of complete human and other metazoan genome sequences has greatly facilitated positioning and analysis of various genomic functional elements, with initial emphasis on coding sequences. However, complete functional maps of sequenced eukaryotic genomes should include also positions of all non-coding regulatory elements. Unfortunately, experimental data on genomic positions of a multitude of regulatory sequences, such as enhancers, silencers, insulators, transcription terminators, and replication origins are very limited, especially at the whole genome level. Since most genomic regulatory elements (e.g. enhancers) are generally gene-, tissue-, or cell-specific, the prediction of these elements by computational methods is difficult and often ambiguous. Therefore, the development of high-throughput experimental approaches for identifying and mapping genomic functional elements is highly desirable. At the same time, the creation of whole-genome map of hundreds of thousands of regulatory elements in several hundreds of tissue/cell types is presently far beyond our capabilities. A possible alternative for the whole genome approach is to concentrate efforts on individual genomic segments and then to integrate the data obtained into a whole genome functional map. Moreover, the maps of polygenic fragments with functional cis-regulatory elements would provide valuable data on complex regulatory systems, including their variability and evolution. Here, we reviewed experimental approaches to the realization of these ideas, including our own developments of experimental techniques for selection of cis-acting functionally active DNA fragments from large (megabase-sized) segments of mammalian genomes.

Genome-wide analysis of alternative promoters of human genes using a custom promoter tiling array

Background

Independent lines of evidence suggested that a large fraction of human genes possess multiple promoters driving gene expression from distinct transcription start sites. Understanding which promoter is employed in which cellular context is required to unravel gene regulatory networks within the cell.

Results

We have developed a custom microarray platform that tiles roughly 35,000 alternative putative promoters from nearly 7,000 genes in the human genome. To demonstrate the utility of this array platform, we have analyzed the patterns of promoter usage in 17β-estradiol (E2)-treated and untreated MCF7 cells and show widespread usage of alternative promoters. Most intriguingly, we show that the downstream promoter in E2-sensitive multiple promoter genes tends to be very close to the 3'-terminus of the gene, suggesting exotic mechanisms of expression regulation in these genes.

Conclusion

The usage of alternative promoters greatly multiplies the transcriptional complexity available within the human genome. The fact that many of these promoters are incapable of driving the synthesis of a meaningful protein-encoding transcript further complicates the story.

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.

High-resolution mapping and characterization of open chromatin across the genome

Mapping DNase I hypersensitive (HS) sites is an accurate method of identifying the location of genetic regulatory elements, including promoters, enhancers, silencers, insulators, and locus control regions. We employed high-throughput sequencing and whole-genome tiled array strategies to identify DNase I HS sites within human primary CD4+ T cells. Combining these two technologies, we have created a comprehensive and accurate genome-wide open chromatin map. Surprisingly, only 16%-21% of the identified 94,925 DNase I HS sites are found in promoters or first exons of known genes, but nearly half of the most open sites are in these regions. In conjunction with expression, motif, and chromatin immunoprecipitation data, we find evidence of cell-type-specific characteristics, including the ability to identify transcription start sites and locations of different chromatin marks utilized in these cells. In addition, and unexpectedly, our analyses have uncovered detailed features of nucleosome structure.

Model-based analysis of two-color arrays (MA2C)

A novel normalization method based on the GC content of probes is developed for two-color tiling arrays. The proposed method, together with robust estimates of the model parameters, is shown to perform superbly on published data sets. A robust algorithm for detecting peak regions is also formulated and shown to perform well compared to other approaches. The tools have been implemented as a stand-alone Java program called MA2C, which can display various plots of statistical analysis for quality control.

Mapping DNase-I hypersensitive sites on human isochores

Mapping DNase-I hypersensitive sites (HS) was used in the past to identify regulatory elements of specific genes. More recently, thousands of HS were identified in the human genome by using high-throughput methods. These approaches showed a general enrichment of HS near or within known genes, within CpG islands, within human-mouse conserved regions and in GC-rich regions of the genome. Here we show that HS: (i) are characterized by a much higher GC level (approximately 56%) than the average GC level of the human genome (approximately 41%); (ii) are overwhelmingly located in the GC-richest compartment of the genome, which is predominantly associated with an open chromatin structure; (iii) and are slightly more and slightly less frequent than genes, respectively, in the gene-rich and in the gene-poor isochore families.

Genomics and the immune system

While the hereditary information encoded in the Watson-Crick base pairing of genomes is largely static within a given individual, access to this information is controlled by dynamic mechanisms. The human genome is pervasively transcribed, but the roles played by the majority of the non-protein-coding genome sequences are still largely unknown. In this review we focus on insights to gene transcriptional regulation by placing special emphasis on genome-wide approaches, and on how non-coding RNAs, which derive from global transcription of the genome, in turn control gene expression. We review recent progress in the field with highlights on the immune system.

Interferon regulatory factors are transcriptional regulators of adipogenesis

We have sought to identify transcriptional pathways in adipogenesis using an integrated experimental and computational approach. Here, we employ high-throughput DNase hypersensitivity analysis to find regions of altered chromatin structure surrounding key adipocyte genes. Regions that display differentiation-dependent changes in hypersensitivity were used to predict binding sites for proteins involved in adipogenesis. A high-scoring example was a binding motif for interferon regulatory factor (IRF) family members. Expression of all nine mammalian IRF mRNAs is regulated during adipogenesis, and several bind to the identified motifs in a differentiation-dependent manner. Furthermore, several IRF proteins repress differentiation. This analysis suggests an important role for IRF proteins in adipocyte biology and demonstrates the utility of this approach in identifying cis- and trans-acting factors not previously suspected to participate in adipogenesis.

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.

Analysis of transcription, chromatin dynamics and epigenetic changes in neural genes

The ways in which gene transcription is investigated have undergone radical change since the turn of the millennium. Piece-meal approaches focussed upon model genes have increasingly been complemented by genome-wide approaches that allow interrogation of multiple cohorts of genes or even entire genomes. This sea change has been founded upon the increasing availability of whole genome sequences and the attendant evolution of microarray based discovery platforms. Collectively, these approaches are being used to build a global and dynamic perspective of transcription factor occupancy, co-factor recruitment and epigenetic signature. As yet, few of these approaches have been applied to the study of neuronal gene transcription, but this is set to change. Here, I review these key developments and point to their potential application to the study of transcriptional and epigenetic changes in neurons in health and disease.

Computational analyses of eukaryotic promoters

Computational analysis of eukaryotic promoters is one of the most difficult problems in computational genomics and is essential for understanding gene expression profiles and reverse-engineering gene regulation network circuits. Here I give a basic introduction of the problem and recent update on both experimental and computational approaches. More details may be found in the extended references. This review is based on a summer lecture given at Max Planck Institute at Berlin in 2005.

Determinants of cell- and gene-specific transcriptional regulation by the glucocorticoid receptor

The glucocorticoid receptor (GR) associates with glucocorticoid response elements (GREs) and regulates selective gene transcription in a cell-specific manner. Native GREs are typically thought to be composite elements that recruit GR as well as other regulatory factors into functional complexes. We assessed whether GR occupancy is commonly a limiting determinant of GRE function as well as the extent to which core GR binding sequences and GRE architecture are conserved at functional loci. We surveyed 100-kb regions surrounding each of 548 known or potentially glucocorticoid-responsive genes in A549 human lung cells for GR-occupied GREs. We found that GR was bound in A549 cells predominately near genes responsive to glucocorticoids in those cells and not at genes regulated by GR in other cells. The GREs were positionally conserved at each responsive gene but across the set of responsive genes were distributed equally upstream and downstream of the transcription start sites, with 63% of them >10 kb from those sites. Strikingly, although the core GR binding sequences across the set of GREs varied extensively around a consensus, the precise sequence at an individual GRE was conserved across four mammalian species. Similarly, sequences flanking the core GR binding sites also varied among GREs but were conserved at individual GREs. We conclude that GR occupancy is a primary determinant of glucocorticoid responsiveness in A549 cells and that core GR binding sequences as well as GRE architecture likely harbor gene-specific regulatory information.

The glucocorticoid receptor (GR) regulates a myriad of physiological functions, such as cell differentiation and metabolism, achieved through modulating transcription in a cell- and gene-specific manner. However, the determinants that specify cell- and gene-specific GR transcriptional regulation are not well established. We describe three properties that contribute to this specificity: (1) GR occupancy at genomic glucocorticoid response elements (GREs) appears to be a primary determinant of glucocorticoid responsiveness; (2) the DNA sequences bound by GR vary widely around a consensus, but the precise sequences of individual GREs are highly conserved, suggesting a role for these sequences in gene-specific GR transcriptional regulation; and (3) native chromosomal GREs were generally found to be composite elements, comprised of multiple factor binding sites that were highly variable in composition, but as with the GR binding sequences, highly conserved at individual GREs. In addition, we discovered that most GREs were positioned far from their GR target genes and that they were equally distributed upstream and downstream of the target genes. These findings, which may be applicable to other regulatory factors, provide fundamental insights for understanding cell- and gene-specific transcriptional regulation.

Optimal design of oligonucleotide microarrays for measurement of DNA copy-number

Copy-number variants (CNVs) occur frequently within the human genome, and may be associated with many human phenotypes. If disease association studies of CNVs are to be performed routinely, it is essential that the copy-number status be accurately genotyped. We systematically assessed the dynamic range response of an oligonucleotide microarray platform to accurately predict copy-number in a set of seven patients who had previously been shown to carry between 1 and 6 copies of an approximately 4 Mb region of 15q12.2-q13.1. We identify probe uniqueness, probe length, uniformity of probe melting temperature, overlap with SNPs and common repeats (particularly Alu elements) and guanine homopolymer content as parameters that significantly affect probe performance. Further, we prove the influence of these criteria on array performance by using these parameters to prospectively filter data from a second array design covering an independent genomic region and observing significant improvements in data quality. The informed selection of probes which have superior performance characteristics allows the prospective design of oligonucleotide arrays which show increased sensitivity and specificity compared with current designs. Although based on the analysis of data from comparative genomic hybridization experiments, we anticipate that our results are relevant to the design of improved oligonucleotide arrays for high-throughput copy-number genotyping of complex regions of the human genome.

Integrated analysis of experimental data sets reveals many novel promoters in 1% of the human genome

The regulation of transcriptional initiation in the human genome is a critical component of global gene regulation, but a complete catalog of human promoters currently does not exist. In order to identify regulatory regions, we developed four computational methods to integrate 129 sets of ENCODE-wide chromatin immunoprecipitation data. They collectively predicted 1393 regions. Roughly 47% of the regions were unique to one method, as each method makes different assumptions about the data. Overall, predicted regions tend to localize to highly conserved, DNase I hypersensitive, and actively transcribed regions in the genome. Interestingly, a significant portion of the regions overlaps with annotated 3'-UTRs, suggesting that some of them might regulate anti-sense transcription. The majority of the predicted regions are >2 kb away from the 5'-ends of previously annotated human cDNAs and hence are novel. These novel regions may regulate unannotated transcripts or may represent new alternative transcription start sites of known genes. We tested 163 such regions for promoter activity in four cell lines using transient transfection assays, and 25% of them showed transcriptional activity above background in at least one cell line. We also performed 5'-RACE experiments on 62 novel regions, and 76% of the regions were associated with the 5'-ends of at least two RACE products. Our results suggest that there are at least 35% more functional promoters in the human genome than currently annotated.

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.

Finding cis-regulatory elements using comparative genomics: some lessons from ENCODE data

Identification of functional genomic regions using interspecies comparison will be most effective when the full span of relationships between genomic function and evolutionary constraint are utilized. We find that sets of putative transcriptional regulatory sequences, defined by ENCODE experimental data, have a wide span of evolutionary histories, ranging from stringent constraint shown by deep phylogenetic comparisons to recent selection on lineage-specific elements. This diversity of evolutionary histories can be captured, at least in part, by the suite of available comparative genomics tools, especially after correction for regional differences in the neutral substitution rate. Putative transcriptional regulatory regions show alignability in different clades, and the genes associated with them are enriched for distinct functions. Some of the putative regulatory regions show evidence for recent selection, including a primate-specific, distal promoter that may play a novel role in regulation.

Positive Cross-Regulatory Loop Ties GATA-3 to Estrogen Receptor α Expression in Breast Cancer

The transcription factor GATA-3 is required for normal mammary gland development, and its expression is highly correlated with estrogen receptor alpha (ER alpha) in human breast tumors. However, the functional role of GATA-3 in ER alpha-positive breast cancers is yet to be established. Here, we show that GATA-3 is required for estradiol stimulation of cell cycle progression in breast cancer cells. The role of GATA-3 in estradiol signaling requires the direct positive regulation of the expression of the ER alpha gene itself by GATA-3. GATA-3 binds to two cis-regulatory elements located within the ER alpha gene, and this is required for RNA polymerase II recruitment to ER alpha promoters. Reciprocally, ER alpha directly stimulates the transcription of the GATA-3 gene, indicating that these two factors are involved in a positive cross-regulatory loop. Moreover, GATA-3 and ER alpha regulate their own expression in breast cancer cells. Hence, this transcriptional coregulatory mechanism accounts for the robust coexpression of GATA-3 and ER alpha in human breast cancers. In addition, these results highlight the crucial role of GATA-3 for the response of ER alpha-positive breast cancers to estradiol. Moreover, they identify GATA-3 as a critical component of the master cell-type-specific transcriptional network including ER alpha and FoxA1 that dictates the phenotype of hormone-dependent breast cancer.

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.

Methods for identification of epigenetic elements in mammalian long multigenic genome sequences

Epigenetic elements of the genome, i.e. elements that determine stably inherited changes in gene expression without changes in the genomic DNA sequence, are essential tools of genetic regulation in higher eukaryotes. The complete sequencing of the human and other genomes allowed studies to be started on positioning of these elements within long multigenic regions of the genome, which is a prerequisite for a comprehensive functional annotation of genomes. This mini-review considers some recent experimental approaches to the high-throughput identification and mapping of epigenetic elements of mammalian genomes, including the mapping of methylated CpG sites, open and closed chromatin regions, and DNase I hypersensitivity sites.