Search

Long noncoding RNAs (lncRNAs) are involved in the regulation of cell cycle, although only a few have been functionally characterized. By combining RNA sequencing and ChIP sequencing of cell cycle synchronized HaCaT cells we have previously identified lncRNAs highly enriched for cell cycle functions. Based on a cyclic expression profile and an overall high correlation to histone 3 lysine 4 trimethylation (H3K4me3) and RNA polymerase II (Pol II) signals, the lncRNA SNHG26 was identified as a top candidate. In the present study we report that downregulation of SNHG26 affects mitochondrial stress, proliferation, cell cycle phase distribution, and gene expression in cis- and in trans, and that this effect is reversed by upregulation of SNHG26. We also find that the effect on cell cycle phase distribution is cell type specific and stable over time. Results indicate an oncogenic role of SNHG26, possibly by affecting cell cycle progression through the regulation of downstream MYC-responsive genes.

The long non-coding RNA (lncRNA) EPB41L4A-AS1 is aberrantly expressed in various cancers and has been reported to be involved in metabolic reprogramming and as a repressor of the Warburg effect. Although the biological relevance of EPB41L4A-AS1 is evident, its functional role seems to vary depending on cell type and state of disease. By combining RNA sequencing and ChIP sequencing of cell cycle synchronized HaCaT cells we previously identified EPB41L4A-AS1 to be one of 59 lncRNAs with potential cell cycle functions. Here, we demonstrate that EPB41L4A-AS1 exists as bright foci and regulates gene expression in the nucleus in both cis and trans. Specifically, we find that EPB41L4A-AS1 positively regulates its sense overlapping gene EPB41L4A and influences expression of hundreds of other genes, including genes involved in cell proliferation. Finally, we show that EPB41L4A-AS1 affects cell cycle phase distribution, though these effects vary between cell types.

Acquiring a specific transcriptomic signature of the human and mouse cardiomyocyte (CM) will greatly increase our understanding of their biology and associated diseases that remain the most deadly across the world. In this study, using comprehensive transcriptomic mining of 91 cell types over 877 samples from bulk RNA-sequencing, single cell RNA-sequencing, and microarray techniques, we describe a unique 118-gene signature of human and mouse primary CMs. Once we had access to this CM-specific gene signature, we investigated the spatial heterogeneity of CMs throughout the heart tissue. Moreover, we compared the CM-specific gene signature to that of CMs derived from 10 differentiation protocols, and we identified the protocols that generate cells most similar to primary CMs. Finally, we looked at the specific differences between primary and differentiated CMs and found that differentiated cells underexpress genes related to CM development and maturity. The differentiated cells conversely overexpressed cell cycle-related genes, resulting in the progenitor features that remain in differentiated CMs compared to primary adult CMs. The presence of histone post translational modification H3K27ac from ChIP sequencing data sets were used to confirm transcriptomic findings. To the best of our knowledge, this is the most comprehensive study to date that unravels the unique transcriptomic signature of primary and differentiated CMs. This study provides important insights into our understanding of CM biology and the molecular mechanisms that make them such a unique cell type. Moreover, the specific transcriptomic signature of CMs could be used in developmental studies, stem cell therapy, regenerative medicine, and drug screening assays.

Most currently available transcriptome data of T-cell acute lymphoblastic leukemia (T-ALL) are based on polyA[+] RNA sequencing methods thus lacking non-polyadenylated transcripts. Here, we present the data of polyA[+] and total RNA sequencing in the context of in vitro TLX1 knockdown in ALL-SIL cells and a primary T-ALL cohort. We extended this dataset with ATAC sequencing and H3K4me1 and H3K4me3 ChIP sequencing data to map putative gene regulatory regions. In this data descriptor, we present a detailed report of how the data were generated and which bioinformatics analyses were performed. Through several technical validations, we showed that our sequencing data are of high quality and that our in vitro TLX1 knockdown was successful. We also validated the quality of the ATAC and ChIP sequencing data and showed that ATAC and H3K4me3 ChIP peaks are enriched at transcription start sites. We believe that this comprehensive set of sequencing data can be reused by others to further unravel the complex biology of T-ALL in general and TLX1 in particular.

Chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq) is a key technique for mapping the distribution and relative abundance of histone posttranslational modifications (PTMs) and chromatin-associated factors across genomes. There is a perceived challenge regarding the ability to quantitatively plot ChIP-seq data, and as such, approaches making use of exogenous additives, or “spike-ins” have recently been developed. Relying on the fact that the IP step of ChIP-seq is a competitive binding reaction, we present a quantitative framework for ChIP-seq analysis that circumvents the need to modify standard sample preparation pipelines with spike-in reagents. We also introduce a visualization technique that, when paired with our formal developments, produces a much more rich characterization of sequencing data.

A major challenge in evaluating quantitative ChIP-seq analyses, such as peak calling and differential binding, is a lack of reliable ground truth data. We present Tulip, a toolkit for rapidly simulating ChIP-seq data using statistical models of the experimental steps. Tulip may be used for a range of applications, including power analysis for experimental design, benchmarking of analysis tools, and modeling effects of processes such as replication on ChIP-seq signals.

Motivation Functional genomics experiments generate genomewide signal profiles that are dense information sources for annotating the regulatory elements. These profiles measure epigenetic activity at the nucleotide resolution and they exhibit distinctive patterns as they fluctuate along the genome. Most notable of these patterns are the valley patterns that are prevalently observed in assays such as ChIP Sequencing and bisulfite sequencing. The genomic positions of valleys pinpoint locations of cis-regulatory elements such as enhancers and insulators. Systematic identification of the valleys provides novel information for delineating the annotation of regulatory elements. Nevertheless, the valleys are not reported by majority of the analysis pipelines. Results We describe EpiSAFARI, a computational method for sensitive detection of valleys from diverse types of epigenetic profiles. EpiSAFARI employs a novel smoothing method for decreasing noise in signal profiles and accounts for technical factors such as sparse signals, mappability and nucleotide content. In performance comparisons, EpiSAFARI performs favorably in terms of accuracy. The histone modification valleys detected by EpiSAFARI exhibit high conservation, transcription factor binding and they are enriched in nascent transcription. In addition, the large clusters of histone valleys are found to be enriched at the promoters of the developmentally associated genes. Differential histone valleys exhibit concordance with differential DNase signal at cell line specific valleys. DNA methylation valleys exhibit elevated conservation and high transcription factor binding. Specifically, we observed enriched binding of transcription factors associated with chromatin structure around methyl-valleys. Availability and implementation EpiSAFARI is publicly available at https://github.com/harmancilab/EpiSAFARI. Supplementary information Supplementary data are available at Bioinformatics online.

Various studies have focussed on understanding the repertoire and biological function of the post-translational modifications that occur on testis-specific histone variants like TH2B, Transition Proteins etc. In our attempt to decipher the unique functions of histone variant TH2B, we discovered a new modification Serine 12 phosphorylation on TH2B (TH2BS12P) in spermatocytes. Our present study is aimed at understanding the function of the TH2BS12P modification in the context of processes that occur during meiotic prophase I. Immunofluorescence studies revealed that TH2BS12P histone mark is enriched in the unsynapsed axes of the sex body and is associated with XY body axes associated proteins like Scp3, γH2AX, pATM, ATR etc. We also observe that TH2BS12P is associated with DSB initiator Spo11 and with several recombination related proteins like pATM, ATR, Rad51, γH2AX etcin vivo. This modification was also found to associate with transcription and recombination related histone marks like H3K4me3 and H3K36me3 in the context of mononucleosomes. Genome-wide occupancy studies as determined by ChIP sequencing experiments revealed that TH2BS12P is localised to subset of recombination hotspots, but majorly associated with H3K4me3 containing genomic regions like gene promoters. Mass spectrometry analysis of proteins that bind to TH2BS12P containing mononucleosomes revealed many proteins linked with the functions of pericentric heterochromatin, transcription and recombination related pathways. We propose that TH2BS12P modification could act alone or in concert with other histone marks for recruitment of appropriate transcription or recombination protein machinery at specific genomic loci. This is the first report documenting the role of a post-translational modification of a germ cell specific histone variant in meiotic prophase I related events.

Identification of transcription factor binding sites (TFBSs) and cis-regulatory motifs (motifs for short) from genomics datasets, provides a powerful view of the rules governing the interactions between TFs and DNA. Existing motif prediction methods however, are limited by high false positive rates in TFBSs identification, contributions from non-sequence-specific binding, and complex and indirect binding mechanisms. High throughput next-generation sequencing data provides unprecedented opportunities to overcome these difficulties, as it provides multiple whole-genome scale measurements of TF binding information. Uncovering this information brings new computational and modeling challenges in high-dimensional data mining and heterogeneous data integration. To improve TFBS identification and novel motifs prediction accuracy in the human genome, we developed an advanced computational technique based on deep learning (DL) and high-performance computing, named DESSO. DESSO utilizes deep neural network and binomial distribution to optimize the motif prediction. Our results showed that DESSO outperformed existing tools in predicting distinct motifs from the 690 in vivo ENCODE ChIP-Sequencing (ChIP-Seq) datasets for 161 human TFs in 91 cell lines. We also found that protein-protein interactions (PPIs) are prevalent among human TFs, and a total of 61 potential tethering binding were identified among the 100 TFs in the K562 cell line. To further expand DESSO’s deep-learning capabilities, we included DNA shape features and found that (i) shape information has a strong predictive power for TF-DNA binding specificity; and (ii) it aided in identification of the shape motifs recognized by human TFs which in turn contributed to the interpretation of TF-DNA binding in the absence of sequence recognition. DESSO and the analyses it enabled will continue to improve our understanding of how gene expression is controlled by TFs and the complexities of DNA binding. The source code and the predicted motifs and TFBSs from the 690 ENCODE TF ChIP-Seq datasets are freely available at the DESSO web server: http://bmbl.sdstate.edu/DESSO.

With decreasing cost of next-generation sequencing (NGS), we are observing a rapid rise in the volume of ‘big data’ in academic research, healthcare and drug discovery sectors. The present bottleneck for extracting value from these ‘big data’ sets is data processing and analysis. Considering this, there is still a lack of reliable, automated and easy to use tools that will allow experimentalists to assess the quality of the sequenced libraries and explore the data first hand, without the need of investing a lot of time of computational core analysts in the early stages of analysis.NGseqBasic is an easy-to-use single-command analysis tool for chromatin accessibility (ATAC, DNaseI) and ChIP sequencing data, providing support to also new techniques such as low cell number sequencing and Cut-and-Run. It takes in fastq, fastq.gz or bam files, conducts all quality control, trimming and mapping steps, along with quality control and data processing statistics, and combines all this to a single-click loadable UCSC data hub, with integral statistics html page providing detailed reports from the analysis tools and quality control metrics. The tool is easy to set up, and no installation is needed. A wide variety of parameters are provided to fine-tune the analysis, with optional setting to generate DNase footprint or high resolution ChIP-seq tracks. A tester script is provided to help in the setup, along with a test data set and downloadable example user cases.NGseqBasic has been used in the routine analysis of next generation sequencing (NGS) data in high-impact publications 1,2. The code is actively developed, and accompanied with Git version control and Github code repository. Here we demonstrate NGseqBasic analysis and features using DNaseI-seq data from GSM689849, and CTCF-ChIP-seq data from GSM2579421, as well as a Cut-and-Run CTCF data set GSM2433142, and provide the one-click loadable UCSC data hubs generated by the tool, allowing for the ready exploration of the run results and quality control files generated by the tool.AvailabilityDownload, setup and help instructions are available on the NGseqBasic web site http://userweb.molbiol.ox.ac.uk/public/telenius/NGseqBasicManual/external/Bioconda users can load the tool as library “ngseqbasic”. The source code with Git version control is available in https://github.com/Hughes-Genome-Group/NGseqBasic/releases.Contactjelena.telenius@imm.ox.ac.uk

Transcriptional regulation occurs via changes to the rates of various biochemical processes. Sequencing-based approaches that average together many cells have suggested that polymerase binding and polymerase release from promoter-proximal pausing are two key regulated steps in the transcriptional process. However, single cell studies have revealed that transcription occurs in short, discontinuous bursts, suggesting that transcriptional burst initiation and termination might also be regulated steps. Here, we develop and apply a quantitative framework to connect changes in both Pol II ChIP-seq and single cell transcriptional measurements to changes in the rates of specific steps of transcription. Using a number of global and targeted transcriptional regulatory perturbations, we show that burst initiation rate is indeed a key regulated step, demonstrating that transcriptional activity can be frequency modulated. Polymerase pause release is a second key regulated step, but the rate of polymerase binding is not changed by any of the biological perturbations we examined. Our results establish an important role for transcriptional burst regulation in the control of gene expression.

How chromatin shapes pathways that promote genome–epigenome integrity in response to DNA damage is an issue of crucial importance. We report that human bromodomain (BRD)-containing proteins, the primary “readers” of acetylated chromatin, are vital for the DNA damage response (DDR). We discovered that more than one-third of all human BRD proteins change localization in response to DNA damage. We identified ZMYND8 (zinc finger and MYND [myeloid, Nervy, and DEAF-1] domain containing 8) as a novel DDR factor that recruits the nucleosome remodeling and histone deacetylation (NuRD) complex to damaged chromatin. Our data define a transcription-associated DDR pathway mediated by ZMYND8 and the NuRD complex that targets DNA damage, including when it occurs within transcriptionally active chromatin, to repress transcription and promote repair by homologous recombination. Thus, our data identify human BRD proteins as key chromatin modulators of the DDR and provide novel insights into how DNA damage within actively transcribed regions requires chromatin-binding proteins to orchestrate the appropriate response in concordance with the damage-associated chromatin context.

We report a highly sensitive and cost-effective method for genome-wide identification of chromatin interactions in eukaryotic cells. Combining proximity ligation with chromatin immunoprecipitation and sequencing, the method outperforms the state of art approach in sensitivity, accuracy and ease of operation. Application of the method to mouse embryonic stem cells improves mapping of enhancer-promoter interactions.

In order to efficiently edit eukaryotic genomes, it is critical to test the impact of chromatin dynamics on CRISPR/Cas9 function and develop strategies to adapt the system to eukaryotic contexts. So far, research has extensively characterized the relationship between the CRISPR endonuclease Cas9 and the composition of the RNADNA duplex that mediates the system’s precision. Evidence suggests that chromatin modifications and DNA packaging can block eukaryotic genome editing by custom-built DNA endonucleases like Cas9; however, the underlying mechanism of Cas9 inhibition is unclear. Here, we demonstrate that closed, gene-silencing-associated chromatin is a mechanism for the interference of Cas9-mediated DNA editing. Our assays use a transgenic cell line with a drug-inducible switch to control chromatin states (open and closed) at a single genomic locus. We show that closed chromatin inhibits editing at specific target sites, and that artificial reversal of the silenced state restores editing efficiency. These results provide new insights to improve Cas9-mediated editing in human and other mammalian cells.

Chromatin ImmunoPrecipitation followed by massively parallel sequencing (ChIP-Seq) has been instrumental to our current view of chromatin structure and function, and identifies correlating histone marks, which together demarcate biologically-relevant domains. However, as with most genome-wide assays, ChIP-seq is an ensemble measurement that reports on the average occupancy of individual modifications in a population of cells. Consequently, our understanding of the combinatorial nature of chromatin states relies almost exclusively on spatial correlations. Here, we report the development of a novel protocol, called indexed Combinatorial ChIP (comb-ChIP), which has the power to determine the genome-wide co-occurrence of histone marks at single nucleosome resolution. We show that at regions of overlapping ChIP signals, certain combinations of marks (H3K36me3 and H3K79me3) tend to co-occur on the same nucleosome, while other combinations (H3K4me3 and H3K36me3) do not, reflecting differences in the underlying chromatin pathways. We further use comb-ChIP to detect changes in histone mark co-occurrence upon genetic perturbation, illuminating new aspects of the Set2-RPD3S pathway. Overall, comb-ChIP promises to greatly improve our understanding of the structural and functional complexity of chromatin.

To better understand genome regulation, it is important to uncover the role of transcription factors in the process of chromatin structure establishment and maintenance. Here we present a data-driven approach to systematically characterize transcription factors that are relevant for this process. Our method uses a linear mixed modeling approach to combine data sets of transcription factor binding motif enrichments in open chromatin and gene expression across the same set of cell lines. Applying this approach to the ENCODE data set we confirm already known and imply numerous novel transcription factors in playing a role in the establishment or maintenance of open chromatin.

DNA replication initiates from thousands of start sites throughout the Drosophila genome and must be coordinated with other ongoing nuclear processes such as transcription to ensure genetic and epigenetic inheritance. Considerable progress has been made toward understanding how chromatin modifications regulate the transcription program; in contrast, we know relatively little about the role of the chromatin landscape in defining how start sites of DNA replication are selected and regulated. Here, we describe the Drosophila replication program in the context of the chromatin and transcription landscape for multiple cell lines using data generated by the modENCODE consortium. We find that while the cell lines exhibit similar replication programs, there are numerous cell line-specific differences that correlate with changes in the chromatin architecture. We identify chromatin features that are associated with replication timing, early origin usage, and ORC binding. Primary sequence, activating chromatin marks, and DNA-binding proteins (including chromatin remodelers) contribute in an additive manner to specify ORC-binding sites. We also generate accurate and predictive models from the chromatin data to describe origin usage and strength between cell lines. Multiple activating chromatin modifications contribute to the function and relative strength of replication origins, suggesting that the chromatin environment does not regulate origins of replication as a simple binary switch, but rather acts as a tunable rheostat to regulate replication initiation events.

The Saccharomyces cerevisiae protein Nhp6A is a model for the abundant and multifunctional high-mobility group B (HMGB) family of chromatin-associated proteins. Nhp6A binds DNA in vitro without sequence specificity and bends DNA sharply, but its role in chromosome biology is poorly understood. We show by whole-genome chromatin immunoprecipitation (ChIP) and high-resolution whole-genome tiling arrays (ChIP–chip) that Nhp6A is localized to specific regions of chromosomes that include ∼23% of RNA polymerase II promoters. Nhp6A binding functions to stabilize nucleosomes, particularly at the transcription start site of these genes. Both genomic binding and transcript expression studies point to functionally related groups of genes that are bound specifically by Nhp6A and whose transcription is altered by the absence of Nhp6. Genomic analyses of Nhp6A mutants specifically defective in DNA bending reveal a critical role of DNA bending for stabilizing chromatin and coregulation of transcription but not for targeted binding by Nhp6A. We conclude that the chromatin environment, not DNA sequence recognition, localizes Nhp6A binding, and that Nhp6A stabilizes chromatin structure and coregulates transcription.

The human transcription factor TP53 is a pivotal roadblock against cancer. A key unresolved question is how the p53 protein selects its genomic binding sites in vivo out of a large pool of potential consensus sites. We hypothesized that chromatin may play a significant role in this site-selection process. To test this, we used a custom DNA microarray to measure p53 binding at approximately 2000 sites predicted to possess high-sequence specificity, and identified both strongly bound and weakly bound sites. When placed within a plasmid, weakly bound sites become p53 responsive and regain p53 binding when stably integrated into random genomic locations. Notably, strongly bound sites reside preferentially within genomic regions whose DNA sequence is predicted to encode relatively high intrinsic nucleosome occupancy. Using in vivo nucleosome occupancy measurements under conditions where p53 is inactive, we experimentally confirmed this prediction. Furthermore, upon p53 activation, nucleosomes are partially displaced from a relatively broad region surrounding the bound p53 sites, and this displacement is rapidly reversed upon inactivation of p53. Thus, in contrast to the general assumption that transcription-factor binding is preferred in sites that have low nucleosome occupancy prior to factor activation, we find that p53 binding occurs preferentially within a chromatin context of high intrinsic nucleosome occupancy.

Genome-wide mapping of nucleosomes has revealed a great deal about the relationships between chromatin structure and control of gene expression, and has led to mechanistic hypotheses regarding the rules by which chromatin structure is established. High-throughput sequencing has recently become the technology of choice for chromatin mapping studies, yet analysis of these experiments is still in its infancy. Here, we introduce a pipeline for analyzing deep sequencing maps of chromatin structure and apply it to data from S. cerevisiae. We analyze a digestion series where nucleosomes are isolated from under- and overdigested chromatin. We find that certain classes of nucleosomes are unusually susceptible or resistant to overdigestion, with promoter nucleosomes easily digested and mid-coding region nucleosomes being quite stable. We find evidence for highly sensitive nucleosomes located within “nucleosome-free regions,” suggesting that these regions are not always completely naked but instead are likely associated with easily digested nucleosomes. Finally, since RNA polymerase is the dominant energy-consuming machine that operates on the chromatin template, we analyze changes in chromatin structure when RNA polymerase is inactivated via a temperature-sensitive mutation. We find evidence that RNA polymerase plays a role in nucleosome eviction at promoters and is also responsible for retrograde shifts in nucleosomes during transcription. Loss of RNA polymerase results in a relaxation of chromatin structure to more closely match in vitro nucleosome positioning preferences. Together, these results provide analytical tools and experimental guidance for nucleosome mapping experiments, and help disentangle the interlinked processes of transcription and chromatin packaging.

Transcription factors are most commonly thought of as proteins that regulate expression of specific genes, independently of the order of those genes along the chromosome. By screening genome-wide chromatin immunoprecipitation (ChIP) profiles in yeast, we find that more than 10% of DNA-binding transcription factors concentrate at the subtelomeric regions near to chromosome ends. None of the proteins identified were previously implicated in regulation at telomeres, yet genomic and proteomic studies reveal that a subset of factors show many interactions with established telomere binding complexes. For many factors, the subtelomeric binding pattern is dynamic and undergoes flux toward or away from the telomere as physiological conditions shift. We find that subtelomeric binding is dependent on environmental conditions and correlates with the induction of gene expression in response to stress. Taken together, these results underscore the importance of genome structure in understanding the regulatory dynamics of transcriptional networks.

CHD7 is a member of the chromodomain helicase DNA binding domain family of ATP-dependent chromatin remodeling enzymes. De novo mutation of the CHD7 gene is a major cause of CHARGE syndrome, a genetic disease characterized by a complex constellation of birth defects (Coloboma of the eye, Heart defects, Atresia of the choanae, severe Retardation of growth and development, Genital abnormalities, and Ear abnormalities). To gain insight into the function of CHD7, we mapped the distribution of the CHD7 protein on chromatin using the approach of chromatin immunoprecipitation on tiled microarrays (ChIP-chip). These studies were performed in human colorectal carcinoma cells, human neuroblastoma cells, and mouse embryonic stem (ES) cells before and after differentiation into neural precursor cells. The results indicate that CHD7 localizes to discrete locations along chromatin that are specific to each cell type, and that the cell-specific binding of CHD7 correlates with a subset of histone H3 methylated at lysine 4 (H3K4me). The CHD7 sites change concomitantly with H3K4me patterns during ES cell differentiation, suggesting that H3K4me is part of the epigenetic signature that defines lineage-specific association of CHD7 with specific sites on chromatin. Furthermore, the CHD7 sites are predominantly located distal to transcription start sites, most often contained within DNase hypersensitive sites, frequently conserved, and near genes expressed at relatively high levels. These features are similar to those of gene enhancer elements, raising the possibility that CHD7 functions in enhancer mediated transcription, and that the congenital anomalies in CHARGE syndrome are due to alterations in transcription of tissue-specific genes normally regulated by CHD7 during development.

Identifying genomic targets of transcription factors is fundamental for understanding transcriptional regulatory networks. Current technology enables identification of all targets of a single transcription factor, but there is no realistic way to achieve the converse: identification of all proteins that bind to a promoter of interest. We have developed a method that promises to fill this void. It employs the yeast retrotransposon Ty5, whose integrase interacts with the Sir4 protein. A DNA-binding protein fused to Sir4 directs insertion of Ty5 into the genome near where it binds; the Ty5 becomes a “calling card” the DNA-binding protein leaves behind in the genome. We constructed customized calling cards for seven transcription factors of yeast by including in each Ty5 a unique DNA sequence that serves as a “molecular bar code.” Ty5 transposition was induced in a population of yeast cells, each expressing a different transcription factor–Sir4 fusion and its matched, bar-coded Ty5, and the calling cards deposited into selected regions of the genome were identified, revealing the transcription factors that visited that region of the genome. In each region we analyzed, we found calling cards for only the proteins known to bind there: In the GAL1–10 promoter we found only calling cards for Gal4; in the HIS4 promoter we found only Gcn4 calling cards; in the PHO5 promoter we found only Pho4 and Pho2 calling cards. We discuss how Ty5 calling cards might be implemented for mapping all targets of all transcription factors in a single experiment.

Sequence motifs that are potentially recognized by DNA-binding proteins occur far more often in genomic DNA than do observed in vivo protein–DNA interactions. To determine how chromatin influences the utilization of particular DNA-binding sites, we compared the in vivo genome-wide binding location of the yeast transcription factor Leu3 to the binding location observed on the same genomic DNA in the absence of any protein cofactors. We found that the DNA-sequence motif recognized by Leu3 in vitro and in vivo was functionally indistinguishable, but Leu3 bound different genomic locations under the two conditions. Accounting for nucleosome occupancy in addition to DNA-sequence motifs significantly improved the prediction of protein–DNA interactions in vivo, but not the prediction of sites bound by purified Leu3 in vitro. Use of histone modification data does not further improve binding predictions, presumably because their effect is already manifest in the global histone distribution. Measurements of nucleosome occupancy in strains that differ in Leu3 genotype show that low nucleosome occupancy at loci bound by Leu3 is not a consequence of Leu3 binding. These results permit quantitation of the epigenetic influence that chromatin exerts on DNA binding-site selection, and provide evidence for an instructive, functionally important role for nucleosome occupancy in determining patterns of regulatory factor targeting genome-wide.

We have developed a new method for determining the DNA-binding specificity of proteins. In DIP-chip (DNA immunoprecipitation with microarray detection), protein·DNA complexes are isolated from an in vitro mixture of purified protein and naked genomic DNA. Whole-genome DNA microarrays are used to identify the protein-bound DNA fragments, and the sequence of the identified fragments is used to derive binding-site descriptions. Using objective criteria for assessing the accuracy of DNA-binding motifs, and using yeast Leu3p as a model, we demonstrate that motifs determined by DIP-chip are as effective at predicting the location of bound proteins in vivo as are motifs determined by conventional low-throughput in vitro methods.

An increasing number of studies indicate a central role for chromatin remodeling in the regulation of gene expression. Current methods for high-resolution studies of the relationship between chromatin accessibility and transcription are low throughput, making a genome-wide study impractical. To enable the simultaneous measurement of the global chromatin accessibility state at the resolution of single genes, we developed the Chromatin Array technique, in which chromatin is separated by its condensation state using either the solubility differences of mono- and oligonucleosomes in specific buffers or controlled DNase I digestion and selection of the large refractory (condensed) DNA fragments. By probing with a comparative genomic hybridization style microarray, we can determine the condensation state of thousands of individual loci and correlate this with transcriptional activity. Applying this technique to the breast tumor model cell line, MCF7, we found that when the condensation is homogeneous in the population of cells, expression is inversely proportional to the level of accessibility and the two methods of accessibility-based target selection correlate well. Using functional annotation and comparative genomic hybridization data, we have begun to decipher the possible biological implications of the relationship between chromatin accessibility and expression.

Previously, identification of promoters regulated by mammalian transcription factors has relied upon overexpression studies. Here we present the identification of a large set of promoters that are bound by E2F in physiological conditions. Probing a human CpG microarray with chromatin immunoprecipitated using an antibody to E2F4, we have identified 68 unique target loci; 15% are bidirectional promoters and 25% recruit E2F via a mechanism distinct from the defined consensus site. Interestingly, although E2F has been shown previously to regulate genes involved in cell cycle progression, many of the new E2F target genes encode proteins involved in DNA repair or recombination. We suggest that human CpG microarrays, in combination with chromatin immunoprecipitation, will allow rapid identification of target promoters for many mammalian transcription factors.

Cells respond to DNA damage by activating a network of signaling pathways that control cell cycle progression and DNA repair. Genetic studies in yeast suggested that several checkpoint proteins, including the RFC-related Rad17 protein, and the PCNA-related Rad1–Rad9–Hus1 protein complex might function as sensors of DNA damage. In this study, we show that the human Rad17 protein recruits the Rad9 protein complex onto chromatin after damage. Rad17 binds to chromatin prior to damage and is phosphorylated by ATR on chromatin after damage but Rad17's phosphorylation is not required for Rad9 loading onto chromatin. The chromatin associations of Rad17 and ATR are largely independent, which suggests that they localize to DNA damage independently. Furthermore, the phosphorylation of Rad17 requires Hus1, suggesting that the Rad1–Rad9–Hus1 complex recruited by Rad17 enables ATR to recognize its substrates. Our data are consistent with a model in which multiple checkpoint protein complexes localize to sites of DNA damage independently and interact to trigger the checkpoint-signaling cascade.

Background: Recent chromatin immunoprecipitation (ChIP) experiments in fly, mouse, and human have revealed the existence of high-occupancy target (HOT) regions or “hotspots” that show enrichment across many assayed DNA-binding proteins. Similar co-enrichment observed in yeast so far has been treated as artifactual, and has not been fully characterized. Results: Here we reanalyze ChIP data from both array-based and sequencing-based experiments to show that in the yeast S. cerevisiae, the collective enrichment phenomenon is strongly associated with proximity to noncoding RNA genes and with nucleosome depletion. DNA sequence motifs that confer binding affinity for the proteins are largely absent from these hotspots, suggesting that protein-protein interactions play a prominent role. The hotspots are condition-specific, suggesting that they reflect a chromatin state or protein state, and are not a static feature of underlying sequence. Additionally, only a subset of all assayed factors is associated with these loci, suggesting that the co-enrichment cannot be simply explained by a chromatin state that is universally more prone to immunoprecipitation. Conclusions: Together our results suggest that the co-enrichment patterns observed in yeast represent transcription factor co-occupancy. More generally, they make clear that great caution must be used when interpreting ChIP enrichment profiles for individual factors in isolation, as they will include factor-specific as well as collective contributions.

Centromere protein A (CENP-A) is an essential histone H3-related protein that constitutes the specialized chromatin of an active centromere. It has been suggested that this protein plays a key role in the epigenetic marking and transformation of noncentromeric genomic DNA into functional neocentromeres. Neocentromeres have been identified on more than two-thirds of the human chromosomes, presumably involving different noncentromeric DNA sequences, but it is unclear whether some generalized sequence properties account for these neocentromeric sites. Using a novel method combining chromatin immunoprecipitation and genomic array hybridization, we have identified a 460-kb CENP-A-binding DNA domain of a neocentromere derived from the 20p12 region of an invdup (20p) human marker chromosome. Detailed sequence analysis indicates that this domain contains no centromeric α-satellite, classical satellites, or other known pericentric repetitive sequence motifs. Putative gene loci are detected, suggesting that their presence does not preclude neocentromere formation. The sequence is not significantly different from surrounding non-CENP-A-binding DNA in terms of the prevalence of various interspersed repeats and binding sites for DNA-interacting proteins (Topoisomerase II and High-Mobility-Group protein I). Notable variations include a higher AT content similar to that seen in human α-satellite DNA and a reduced prevalence of long terminal repeats (LTRs), short interspersed repeats (SINEs), and Alus. The significance of these features in neocentromerization is discussed.

The binding of p53 protein to DNA is stimulated by its interaction with covalent as well as noncovalent modifiers. We report the identification of a factor from HeLa nuclear extracts that activates p53 DNA binding. This factor was purified to homogeneity and identified as the high mobility group protein, HMG-1. HMG-1 belongs to a family of highly conserved chromatin-associated nucleoproteins that bend DNA and facilitate the binding of various transcription factors to their cognate DNA sequences. We demonstrate that recombinant His-tagged HMG-1 enhances p53 DNA binding in vitro and also that HMG-1 and p53 can interact directly in vitro. Unexpectedly, HMG-1 also stimulates DNA binding by p53Delta30, a carboxy-terminally deleted form of the protein that is considered to be constitutively active, suggesting that HMG-1 stimulates p53 by a mechanism that is distinct from other known activators of p53. Finally, using transient transfection assays we show that HMG-1 can increase p53 and p53Delta30-mediated transactivation in vivo. HMG-1 promotes the assembly of higher order p53 nucleoprotein structures, and these data, along with the fact that HMG-1 is capable of bending DNA, suggest that HMG-1 may activate p53 DNA binding by a novel mechanism involving a structural change in the target DNA.

Purified, reconstituted chromatin templates containing regular, physiological nucleosome spacing were transcribed in vitro by RNA polymerase II along with the Gal4-VP16 activator. When Gal4-VP16 was prebound to DNA before reconstitution of either H1-deficient or H1-containing chromatin, the resulting templates were transcribed with a similar efficiency. Under such conditions, we observed long-range (1000 bp) activation of transcription in vitro with H1-containing chromatin, but not naked DNA templates. When Gal4-VP16 was added to preassembled chromatin, the H1-deficient chromatin was transcriptionally active, whereas the H1-containing chromatin, which possessed properties similar to native chromatin, was transcriptionally inert. We then mimicked DNA replication and chromatin assembly at a replication fork and found that Gal4-VP16 could potentiate transcription during, but not after, replication and assembly of histone H1-containing chromatin. These experiments provide biochemical data that support a DNA replication-dependent mechanism for reconfiguration of chromatin structure and activation of transcription by Gal4-VP16.

T h e expression of eukaryotic genes is mediated at the transcriptional level by a complex interplay between trans-acting factors and cis-acting elements. Interactions of proteins with specific DNA sequences are most commonly studied in vitro by applying the popular gel mobility-shift assays and footprinting techniques. These experiments are usually done using nuclear extracts from cell lines or tissues of interest. Beyond this type of analysis, it is desirable to have evidence that a protein which complexes with a particular DNA sequence in the test tube is actually bound to the same sequences under in vivo conditions. Within the cell, the accessibility of a sequence to a trans-acting factor may be restricted by a particular chromat in environment. Dynamic interactions of protein factors with their DNA target may be accompanied by changes in chromat in structure. These complex situations and processes are difficult to reproduce in vitro. To investigate protein-DNA contacts in living cells, in vivo footprinting strategies have been devised. (1'2) These methods are based on direct genomic sequencing (3) where total genomic DNA is subjected to base-specific chemical cleavage and the sequence of a particular single-copy gene is then revealed without conventional cloning procedures. For complex genornes, such as those of mammals , special techniques are required to accomplish this goal. The most sensitive method uses the ligation-mediated polymerase chain reaction (LMPCR). ~4-6) This method significantly increases the sensitivity and specificity of the original genomic sequencing methods. All of the component fragments of a DNA sequence ladder are s imultaneously amplified by PCR and displayed as sequencing gel autoradiograms. LMPCR preserves the quantitative representation of each fragment of the sequence ladder in the original population of cleaved molecules and each sequence ladder can be amplified very reproducibly. This has made possible its use for in vivo footprinting, ~4) genomic sequencing, ~s) methylat ion analysis, ~s'7'8) and studies on chromat in structure at single-nucleotide resolution. {9)

Chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-Seq) is a powerful method used to identify genome-wide binding patterns of transcription factors and distribution of various histone modifications associated with different chromatin states. In most published studies, ChIP-Seq has been performed on cultured cells grown under controlled conditions, allowing generation of large amounts of material in a homogeneous biological state. Although such studies have provided great insight into the dynamic landscapes of animal genomes, they do not allow the examination of transcription factor binding and chromatin states in adult tissues, developing embryonic structures, or tumors. Such knowledge is critical to understanding the information required to create and maintain a complex biological tissue and to identify noncoding regions of the genome directly involved in tissues affected by complex diseases such as autism. Studying these tissue types with ChIP-Seq can be challenging due to the limited availability of tissues and the lack of complex biological states able to be achieved in culture. These inherent differences require alterations of standard cross-linking and chromatin extraction typically used in cell culture. Here we describe a general approach for using small amounts of animal tissue to perform ChIP-Seq directed at histone modifications and transcription factors. Tissue is homogenized before treatment with formaldehyde to ensure proper cross-linking, and a two-step nuclear isolation is performed to increase extraction of soluble chromatin. Small amounts of soluble chromatin are then used for immunoprecipitation (IP) and prepared for multiplexed high-throughput sequencing.

Chromatin immunoprecipitation (ChIP) combined with genomic analysis techniques provide a global snapshot of protein–DNA interactions in the context of chromatin, yielding insights into which genomic loci might be regulated by the DNA-associated protein under investigation. This protocol describes how to perform ChIP on intact or dissected Xenopus embryos. The ChIP-isolated DNA fragments are suitable for high-throughput sequencing (ChIP-Seq) or for quantitative PCR (ChIP-qPCR). In this protocol, embryonic tissue is harvested from Xenopus tropicalis or Xenopus laevis at the developmental stage of interest, and DNA-associated proteins are immobilized to their endogenous genomic binding sites with formaldehyde. Nuclei are extracted from embryos and subjected to sonication so as to shear the chromatin to a size that allows sufficient positional resolution of protein binding to genomic DNA. Chromatin fragments bound by the protein of interest are immunoprecipitated using antibody-coupled beads, washed under high-stringency conditions, and stripped from the beads with anionic detergents. The chemical cross-links are reversed, and the coimmunoprecipitated DNA is purified. The resulting DNA fragments can be analyzed by qPCR or used to create a ChIP-Seq library. General advice for qPCR and for making ChIP-Seq libraries is offered, and approaches for analyzing ChIP-Seq data are outlined.

The DNA context of nucleotides that a protein recognizes can influence the strength of the protein-DNA interaction. Moreover, in prokaryotes, understanding the quantitative differences in binding affinities that result in part from the DNA context is often important in describing regulatory mechanisms. Nevertheless, these issues have not been a major focus yet for the investigation of protein-DNA interactions in eukaryotes. In this study, we explored the binding specificity and the range of affinities that the BPV-1 E2 transcriptional activator has for DNA. Because E2 binding sites are positioned near several different BPV-1 promoters, such quantitative information may be important to understand transcriptional regulatory mechanisms in BPV-1. Gel retardation assays and DNA footprinting were used to quantitate the affinities of the E2 binding sites in the viral genome. In the process, five sites were discovered, which, on the basis of sequence, had not been predicted previously to interact with the E2 protein. Equilibrium and kinetic studies show that the range of E2 affinities of the 17 sites varied over 300-fold. The sequence elements responsible for E2 recognition of DNA were determined by missing contact analysis of several sites and a point mutation analysis of one site. The results presented show that the affinity of an E2 binding site is to a large extent determined by the availability of specific contacts, but the data also strongly suggest that DNA structure plays an important role.