The signal elements in promoter sequences are not well characterized. We developed statistical tests to find nucleotide words (generally of length 8) that appear localized relative to TSSs (transcription start site). These words constituted """"""""seeds"""""""" for expansion to develop PSSMs (position-specific scoring matrices) characterizing systems of co-regulated genes. To this end, Dr. Marino-Ramirez collected a database of about 4700 sequences around the TSS of human genes, later increasing the size of the database by about a factor of 2. The database was exceptionally well characterized, and ideal for our statistical study. We used a Poisson scan statistic to determine whether occurrences of a given 8-letter DNA word are clustered unusually relative to the TSS. The Poisson scan statistic also identified clusters of significant words. About 80 of these words occurred in two or three clusters. By validating our results with microarray data and gene ontology information, we were therefore able to show that the same 8-letter word could have two different biological functions, depending on its position with respect to the TSS. Although this kind of positional dependency is a known phenomenon, our study showed that it is widespread in the human genome. ? ? In addition, with gold standard datasets and rigorous statistical tests, we showed that Markov models and positional information improve TFBS prediction significantly (although not yet to practical accuracies). Moreover, we showed that the Markov models used in extant TFBS programs is inferior, both theoretically and practically, to the theoretically correct Markov model we proposed. Positional information and the theoretically sound Markov models have been incorporated into the publicly available motif-finding program A-GLAM.
| Kim, Nak-Kyeong; Tharakaraman, Kannan; Spouge, John L (2006) Adding sequence context to a Markov background model improves the identification of regulatory elements. Bioinformatics 22:2870-5 |
| Tharakaraman, Kannan; Marino-Ramirez, Leonardo; Sheetlin, Sergey L et al. (2006) Scanning sequences after Gibbs sampling to find multiple occurrences of functional elements. BMC Bioinformatics 7:408 |
| Tharakaraman, Kannan; Marino-Ramirez, Leonardo; Sheetlin, Sergey et al. (2005) Alignments anchored on genomic landmarks can aid in the identification of regulatory elements. Bioinformatics 21 Suppl 1:i440-8 |
| Marino-Ramirez, Leonardo; Spouge, John L; Kanga, Gavin C et al. (2004) Statistical analysis of over-represented words in human promoter sequences. Nucleic Acids Res 32:949-58 |
| Frith, Martin C; Hansen, Ulla; Spouge, John L et al. (2004) Finding functional sequence elements by multiple local alignment. Nucleic Acids Res 32:189-200 |
