Rank-statistics based enrichment-site prediction algorithm developed for chromatin immunoprecipitation on chip experiments

Background: High density oligonucleotide tiling arrays are an effective and powerful platform for conducting unbiased genome-wide studies. The ab initio probe selection method employed in tiling arrays is unbiased, and thus ensures consistent sampling across coding and non-coding regions of the genome. Tiling arrays are increasingly used in chromatin immunoprecipitation (IP) experiments (ChIP on chip). ChIP on chip facilitates the generation of genome-wide maps of in-vivo interactions between DNA-associated proteins including transcription factors and DNA. Analysis of the hybridization of an immunoprecipitated sample to a tiling array facilitates the identification of ChIP-enriched segments of the genome. These enriched segments are putative targets of antibody assayable regulatory elements. The enrichment response is not ubiquitous across the genome. Typically 5 to 10% of tiled probes manifest some significant enrichment. Depending upon the factor being studied, this response can drop to less than 1%. The detection and assessment of significance for interactions that emanate from non-canonical and/or un-annotated regions of the genome is especially challenging. This is the motivation behind the proposed algorithm. Results: We have proposed a novel rank and replicate statistics-based methodology for identifying and ascribing statistical confidence to regions of ChIP-enrichment. The algorithm is optimized for identification of sites that manifest low levels of enrichment but are true positives, as validated by alternative biochemical experiments. Although the method is described here in the context of ChIP on chip experiments, it can be generalized to any treatment-control experimental design. The results of the algorithm show a high degree of concordance with independent biochemical validation methods. The sensitivity and specificity of the algorithm have been characterized via quantitative PCR and independent computational approaches. Conclusion: The algorithm ranks all enrichment sites based on their intra-replicate ranks and inter-replicate rank consistency. Following the ranking, the method allows segmentation of sites based on a meta p-value, a composite array signal enrichment criterion, or a composite of these two measures. The sensitivities obtained subsequent to the segmentation of data using a meta p-value of 10(-5), an array signal enrichment of 0.2 and a composite of these two values are 88%, 87% and 95%, respectively

NF-Y co-associates with FOS at promoters, enhancers, repetitive elements, and inactive chromatin regions, and is stereo-positioned with growth-controlling transcription factors.

NF-Y, a trimeric transcription factor (TF) composed of two histone-like subunits (NF-YB and NF-YC) and a sequence-specific subunit (NF-YA), binds to the CCAAT motif, a common promoter element. Genome-wide mapping reveals 5,000-15,000 NF-Y binding sites depending on the cell type, with the NF-YA and NF-YB subunits binding asymmetrically with respect to the CCAAT motif. Despite being characterized as a proximal promoter TF, only 25% of NF-Y sites map to promoters. A comparable number of NF-Y sites are located at enhancers, many of which are tissue specific, and nearly half of the NF-Y sites are in select subclasses of HERV LTR repeats. Unlike most TFs, NF-Y can access its target DNA motif in inactive (non-modified) or polycomb-repressed chromatin domains. Unexpectedly, NF-Y extensively co-localizes with FOS in all genomic contexts, and this often occurs in the absence of JUN and the AP-1 motif. NF-Y also co-associates with a select cluster of growth-controlling and oncogenic TFs, consistent with the abundance of CCAAT motifs in the promoters of genes overexpressed in cancer. Interestingly, NF-Y and several growth-controlling TFs bind in a stereo-specific manner, suggesting a mechanism for cooperative action at promoters and enhancers. Our results indicate that NF-Y is not merely a commonly-used, proximal promoter TF, but rather performs a more diverse set of biological functions, many of which are likely to involve co-association with FOS

Sequential organizing activities of engrailed, hedgehog and decapentaplegic in the Drosophila wing.

The Drosophila wing is formed by two cell populations, the anterior and posterior compartments, which are distinguished by the activity of the selector gene engrailed (en) in posterior cells. Here, we show that en governs growth and patterning in both compartments by controlling the expression of the secreted proteins hedgehog (hh) and decapentaplegic (dpp) as well as the response of cells to these signaling molecules. First, we demonstrate that en activity programs wing cells to express hh whereas the absence of en activity programs them to respond to hh by expressing dpp. As a consequence, posterior cells secrete hh and induce a stripe of neighboring anterior cells across the compartment boundary to secrete dpp. Second, we demonstrate that dpp can exert a long-range organizing influence on surrounding wing tissue, specifying anterior or posterior pattern depending on the compartmental provenance, and hence the state of en activity, of the responding cells. Thus, dpp secreted by anterior cells along the compartment boundary has the capacity to organize the development of both compartments. Finally, we report evidence suggesting that dpp may exert its organizing influence by acting as a gradient morphogen in contrast to hh which appears to act principally as a short range inducer of dpp

Genome-wide mapping indicates that p73 and p63 Co-occupy target sites and have similar DNA-binding profiles in vivo

Background: The p53 homologs, p63 and p73, share, ~85% amino acid identity in their DNA-binding domains, but they have distinct biological functions. Principal Findings: Using chromatin immunoprecipitation and high-resolution tiling arrays covering the human genome, we identify p73 DNA binding sites on a genome-wide level in ME180 human cervical carcinoma cells. Strikingly, the p73 binding profile is indistinguishable from the previously described binding profile for p63 in the same cells. Moreover, the p73:p63 binding ratio is similar at all genomic loci tested, suggesting that there are few, if any, targets that are specific for one of these factors. As assayed by sequential chromatin immunoprecipitation, p63 and p73 co-occupy DNA target sites in vivo, suggesting that p63 and p73 bind primarily as heterotetrameric complexes in ME180 cells. Conclusions: The observation that p63 and p73 associate with the same genomic targets suggest that their distinct biological functions are due to cell-type specific expression and/or protein domains that involve functions other than DNA binding. © 2010 Yang et al

High-throughput sequencing reveals a simple model of nucleosome energetics

We use nucleosome maps obtained by high-throughput sequencing to study sequence specificity of intrinsic histone-DNA interactions. In contrast with previous approaches, we employ an analogy between a classical one-dimensional fluid of finite-size particles in an arbitrary external potential and arrays of DNA-bound histone octamers. We derive an analytical solution to infer free energies of nucleosome formation directly from nucleosome occupancies measured in high-throughput experiments. The sequence-specific part of free energies is then captured by fitting them to a sum of energies assigned to individual nucleotide motifs. We have developed hierarchical models of increasing complexity and spatial resolution, establishing that nucleosome occupancies can be explained by systematic differences in mono- and dinucleotide content between nucleosomal and linker DNA sequences, with periodic dinucleotide distributions and longer sequence motifs playing a secondary role. Furthermore, similar sequence signatures are exhibited by control experiments in which genomic DNA is either sonicated or digested with micrococcal nuclease in the absence of nucleosomes, making it possible that current predictions based on high-throughput nucleosome positioning maps are biased by experimental artifacts.Comment: 36 pages, 13 figure

Direct and long-range action of a DPP morphogen gradient.

During development of the Drosophila wing, the decapentaplegic (dpp) gene is expressed in a stripe of cells along the anteroposterior compartment boundary and gives rise to a secreted protein that exerts a long-range organizing influence on both compartments. Using clones of cells that express DPP, or in which DPP receptor activity has been constitutively activated or abolished, we show that DPP acts directly and at long range on responding cells, rather than by proxy through the short-range induction of other signaling molecules. Further, we show that two genes, optomotor-blind and spalt are transcriptionally activated at different distances from DPP-secreting cells and provide evidence that these genes respond to different threshold concentrations of DPP protein. We propose that DPP acts as a gradient morphogen during wing development

Extensive chromatin fragmentation improves enrichment of protein binding sites in chromatin immunoprecipitation experiments

Extensive sonication of formaldehyde-crosslinked chromatin can generate DNA fragments averaging 200 bp in length (range 75–300 bp). Fragmentation is largely random with respect to genomic region and nucleosome position. ChIP experiments employing such extensively fragmented samples show 2- to 4-fold increased enrichment of protein binding sites over control genomic regions, when compared to samples sonicated to a more conventional size range (300–500 bp). The basis of improved fold enrichments is that immunoprecipitation of protein-bound regions is unaffected by fragment size, whereas immunoprecipitation of control genomic regions decreases progressively along with reduced fragment size due to fewer nonspecific binding sites. The use of extensively sonicated samples improves mapping of protein binding sites, and it extends the dynamic range for quantitative measurements of histone density. We show that many yeast promoter regions are virtually devoid of histones

NF-Y coassociates with FOS at promoters, enhancers, repetitive elements, and inactive chromatin regions, and is stereo-positioned with growth-controlling transcription factors

NF-Y, a trimeric transcription factor (TF) composed of two histone-like subunits (NF-YB and NF-YC) and a sequencespecific subunit (NF-YA), binds to the CCAAT motif, a common promoter element. Genome-wide mapping reveals 5000-15,000 NF-Y binding sites depending on the cell type, with the NF-YA and NF-YB subunits binding asymmetrically with respect to the CCAAT motif. Despite being characterized as a proximal promoter TF, only 25% of NF-Y sites map to promoters. A comparable number of NF-Y sites are located at enhancers, many of which are tissue specific, and nearly half of the NF-Y sites are in select subclasses of HERV LTR repeats. Unlike most TFs, NF-Y can access its target DNA motif in inactive (nonmodified) or polycomb-repressed chromatin domains. Unexpectedly, NF-Y extensively colocalizes with FOS in all genomic contexts, and this often occurs in the absence of JUN and the AP-1 motif. NF-Y also coassociates with a select cluster of growth-controlling and oncogenic TFs, consistent with the abundance of CCAAT motifs in the promoters of genes overexpressed in cancer. Interestingly, NF-Y and several growth-controlling TFs bind in a stereo-specific manner, suggesting a mechanism for cooperative action at promoters and enhancers. Our results indicate that NF-Y is not merely a commonly used proximal promoter TF, but rather performs a more diverse set of biological functions, many of which are likely to involve coassociation with FOS

Tomato protoplast DNA transformation: physical linkage and recombination of exogenous DNA sequences

Tomato protoplasts have been transformed with plasmid DNA's, containing a chimeric kanamycin resistance gene and putative tomato origins of replication. A calcium phosphate-DNA mediated transformation procedure was employed in combination with either polyethylene glycol or polyvinyl alcohol. There were no indications that the tomato DNA inserts conferred autonomous replication on the plasmids. Instead, Southern blot hybridization analysis of seven kanamycin resistant calli revealed the presence of at least one kanamycin resistance locus per transformant integrated in the tomato nuclear DNA. Generally one to three truncated plasmid copies were found integrated into the tomato nuclear DNA, often physically linked to each other. For one transformant we have been able to use the bacterial ampicillin resistance marker of the vector plasmid pUC9 to 'rescue' a recombinant plasmid from the tomato genome. Analysis of the foreign sequences included in the rescued plasmid showed that integration had occurred in a non-repetitive DNA region. Calf-thymus DNA, used as a carrier in transformation procedure, was found to be covalently linked to plasmid DNA sequences in the genomic DNA of one transformant. A model is presented describing the fate of exogenously added DNA during the transformation of a plant cell. The results are discussed in reference to the possibility of isolating DNA sequences responsible for autonomous replication in tomato.