Discovering Motifs in Ranked Lists of DNA Sequences

Abstract

Computational methods for discovery of sequence elements that are enriched in a target set compared with a background set are fundamental in molecular biology research. One example is the discovery of transcription factor binding motifs that are inferred from ChIP–chip (chromatin immuno-precipitation on a microarray) measurements. Several major challenges in sequence motif discovery still require consideration: (i) the need for a principled approach to partitioning the data into target and background sets; (ii) the lack of rigorous models and of an exact p-value for measuring motif enrichment; (iii) the need for an appropriate framework for accounting for motif multiplicity; (iv) the tendency, in many of the existing methods, to report presumably significant motifs even when applied to randomly generated data. In this paper we present a statistical framework for discovering enriched sequence elements in ranked lists that resolves these four issues. We demonstrate the implementation of this framework in a software application, termed DRIM (discovery of rank imbalanced motifs), which identifies sequence motifs in lists of ranked DNA sequences. We applied DRIM to ChIP–chip and CpG methylation data and obtained the following results. (i) Identification of 50 novel putative transcription factor (TF) binding sites in yeast ChIP–chip data. The biological function of some of them was further investigated to gain new insights on transcription regulation networks in yeast. For example, our discoveries enable the elucidation of the network of the TF ARO80. Another finding concerns a systematic TF binding enhancement to sequences containing CA repeats. (ii) Discovery of novel motifs in human cancer CpG methylation data. Remarkably, most of these motifs are similar to DNA sequence elements bound by the Polycomb complex that promotes histone methylation. Our findings thus support a model in which histone methylation and CpG methylation are mechanistically linked. Overall, we demonstrate that the statistical framework embodied in the DRIM software tool is highly effective for identifying regulatory sequence elements in a variety of applications ranging from expression and ChIP–chip to CpG methylation data. DRIM is publicly available at http://bioinfo.cs.technion.ac.il/drim.

Microtubule disruption synergizes with oncolytic virotherapy by inhibiting interferon translation and potentiating bystander killing

“…3a). On the basis of gene set-enrichment analyses using both GOrilla and DAVID gene set databases [21][22][23][24] , this group of 248 genes, and more specifically the 158 virus-induced genes dampened by colchicine, is significantly enriched in type I IFN-responsive genes (Fig. 3a, lower inset and Fig.…”

Longitudinal muscle gene expression patterns associated with differential intramuscular fat in cattle

“…We submitted this ranked list to GOrilla (Eden et al, 2007) and explored it for functional enrichment. This enriched for fatty acid metabolism with THRSP, PCK1, ADIPOQ, ADIG, SCD, G0S2, ACSM1, PLIN1, CIDEA and TUSC5, all in the top 25 (Table 1; Supplementary Material S1 1 ).…”

Development of Methods for Automatic Extraction of Knowledge From Texts of Scientific Publications for the Creation of a Knowledge Base Solanum Tuberosum

“…Such methods are widely used in known computer systems intended for interpreting experimental data, e.g. DAVID [34], PANTHER [35,36], GORILLA [37,38], etc.…”