Principal component analysis for predicting transcription-factor binding motifs from array-derived data

BACKGROUND: The responses to interleukin 1 (IL-1) in human chondrocytes constitute a complex regulatory mechanism, where multiple transcription factors interact combinatorially to transcription-factor binding motifs (TFBMs). In order to select a critical set of TFBMs from genomic DNA information and an array-derived data, an efficient algorithm to solve a combinatorial optimization problem is required. Although computational approaches based on evolutionary algorithms are commonly employed, an analytical algorithm would be useful to predict TFBMs at nearly no computational cost and evaluate varying modelling conditions. Singular value decomposition (SVD) is a powerful method to derive primary components of a given matrix. Applying SVD to a promoter matrix defined from regulatory DNA sequences, we derived a novel method to predict the critical set of TFBMs. RESULTS: The promoter matrix was defined to establish a quantitative relationship between the IL-1-driven mRNA alteration and genomic DNA sequences of the IL-1 responsive genes. The matrix was decomposed with SVD, and the effects of 8 potential TFBMs (5'-CAGGC-3', 5'-CGCCC-3', 5'-CCGCC-3', 5'-ATGGG-3', 5'-GGGAA-3', 5'-CGTCC-3', 5'-AAAGG-3', and 5'-ACCCA-3') were predicted from a pool of 512 random DNA sequences. The prediction included matches to the core binding motifs of biologically known TFBMs such as AP2, SP1, EGR1, KROX, GC-BOX, ABI4, ETF, E2F, SRF, STAT, IK-1, PPARγ, STAF, ROAZ, and NFκB, and their significance was evaluated numerically using Monte Carlo simulation and genetic algorithm. CONCLUSION: The described SVD-based prediction is an analytical method to provide a set of potential TFBMs involved in transcriptional regulation. The results would be useful to evaluate analytically a contribution of individual DNA sequences

High Resolution Models of Transcription Factor-DNA Affinities Improve In Vitro and In Vivo Binding Predictions

Accurately modeling the DNA sequence preferences of transcription factors (TFs), and using these models to predict in vivo genomic binding sites for TFs, are key pieces in deciphering the regulatory code. These efforts have been frustrated by the limited availability and accuracy of TF binding site motifs, usually represented as position-specific scoring matrices (PSSMs), which may match large numbers of sites and produce an unreliable list of target genes. Recently, protein binding microarray (PBM) experiments have emerged as a new source of high resolution data on in vitro TF binding specificities. PBM data has been analyzed either by estimating PSSMs or via rank statistics on probe intensities, so that individual sequence patterns are assigned enrichment scores (E-scores). This representation is informative but unwieldy because every TF is assigned a list of thousands of scored sequence patterns. Meanwhile, high-resolution in vivo TF occupancy data from ChIP-seq experiments is also increasingly available. We have developed a flexible discriminative framework for learning TF binding preferences from high resolution in vitro and in vivo data. We first trained support vector regression (SVR) models on PBM data to learn the mapping from probe sequences to binding intensities. We used a novel -mer based string kernel called the di-mismatch kernel to represent probe sequence similarities. The SVR models are more compact than E-scores, more expressive than PSSMs, and can be readily used to scan genomics regions to predict in vivo occupancy. Using a large data set of yeast and mouse TFs, we found that our SVR models can better predict probe intensity than the E-score method or PBM-derived PSSMs. Moreover, by using SVRs to score yeast, mouse, and human genomic regions, we were better able to predict genomic occupancy as measured by ChIP-chip and ChIP-seq experiments. Finally, we found that by training kernel-based models directly on ChIP-seq data, we greatly improved in vivo occupancy prediction, and by comparing a TF's in vitro and in vivo models, we could identify cofactors and disambiguate direct and indirect binding

Dynamic telomerase gene suppression via network effects of GSK3 inhibition

<b>Background</b>: Telomerase controls telomere homeostasis and cell immortality and is a promising anti-cancer target, but few small molecule telomerase inhibitors have been developed. Reactivated transcription of the catalytic subunit hTERT in cancer cells controls telomerase expression. Better understanding of upstream pathways is critical for effective anti-telomerase therapeutics and may reveal new targets to inhibit hTERT expression. <b>Methodology/Principal Findings</b>: In a focused promoter screen, several GSK3 inhibitors suppressed hTERT reporter activity. GSK3 inhibition using 6-bromoindirubin-3′-oxime suppressed hTERT expression, telomerase activity and telomere length in several cancer cell lines and growth and hTERT expression in ovarian cancer xenografts. Microarray analysis, network modelling and oligonucleotide binding assays suggested that multiple transcription factors were affected. Extensive remodelling involving Sp1, STAT3, c-Myc, NFκB, and p53 occurred at the endogenous hTERT promoter. RNAi screening of the hTERT promoter revealed multiple kinase genes which affect the hTERT promoter, potentially acting through these factors. Prolonged inhibitor treatments caused dynamic expression both of hTERT and of c-Jun, p53, STAT3, AR and c-Myc. <b>Conclusions/Significance</b>: Our results indicate that GSK3 activates hTERT expression in cancer cells and contributes to telomere length homeostasis. GSK3 inhibition is a clinical strategy for several chronic diseases. These results imply that it may also be useful in cancer therapy. However, the complex network effects we show here have implications for either setting

The Genomic Context and Corecruitment of SP1 Affect ERRα Coactivation by PGC-1α in Muscle Cells

The peroxisome proliferator-activated receptor-γ coactivator 1α (PGC-1α) coordinates the transcriptional network response to promote an improved endurance capacity in skeletal muscle, eg, by coactivating the estrogen-related receptor-α (ERRα) in the regulation of oxidative substrate metabolism. Despite a close functional relationship, the interaction between these 2 proteins has not been studied on a genomic level. We now mapped the genome-wide binding of ERRα to DNA in a skeletal muscle cell line with elevated PGC-1α and linked the DNA recruitment to global PGC-1α target gene regulation. We found that, surprisingly, ERRα coactivation by PGC-1α is only observed in the minority of all PGC-1α recruitment sites. Nevertheless, a majority of PGC-1α target gene expression is dependent on ERRα. Intriguingly, the interaction between these 2 proteins is controlled by the genomic context of response elements, in particular the relative GC and CpG content, monomeric and dimeric repeat-binding site configuration for ERRα, and adjacent recruitment of the transcription factor specificity protein 1. These findings thus not only reveal a novel insight into the regulatory network underlying muscle cell plasticity but also strongly link the genomic context of DNA-response elements to control transcription factor-coregulator interactions

Artificial ants deposit pheromone to search for regulatory DNA elements

BACKGROUND: Identification of transcription-factor binding motifs (DNA sequences) can be formulated as a combinatorial problem, where an efficient algorithm is indispensable to predict the role of multiple binding motifs. An ant algorithm is a biology-inspired computational technique, through which a combinatorial problem is solved by mimicking the behavior of social insects such as ants. We developed a unique version of ant algorithms to select a set of binding motifs by considering a potential contribution of each of all random DNA sequences of 4- to 7-bp in length. RESULTS: Human chondrogenesis was used as a model system. The results revealed that the ant algorithm was able to identify biologically known binding motifs in chondrogenesis such as AP-1, NFκB, and sox9. Some of the predicted motifs were identical to those previously derived with the genetic algorithm. Unlike the genetic algorithm, however, the ant algorithm was able to evaluate a contribution of individual binding motifs as a spectrum of distributed information and predict core consensus motifs from a wider DNA pool. CONCLUSION: The ant algorithm offers an efficient, reproducible procedure to predict a role of individual transcription-factor binding motifs using a unique definition of artificial ants

Identification of co-regulated candidate genes by promoter analysis.

EThOS - Electronic Theses Online ServiceGBUnited Kingdo

A framework to identify epigenome and transcription factor crosstalk

While changes in chromatin are integral to transcriptional reprogramming during cellular differentiation, it is currently unclear how chromatin modifications are targeted to specific loci. To systematically identify transcription factors (TFs) that can direct chromatin changes during cell fate decisions, we model the genome-wide dynamics of chromatin marks in terms of computationally predicted TF binding sites. By applying this computational approach to a time course of Polycomb-mediated H3K27me3 marks during neuronal differentiation of murine stem cells, we identify several motifs that likely regulate dynamics of this chromatin mark. Among these, the motifs bound by REST and by the SNAIL family of TFs are predicted to transiently recruit H3K27me3 in neuronal progenitors. We validate these predictions experimentally and show that absence of REST indeed causes loss of H3K27me3 at target promoters in trans, specifically at the neuronal progenitor state. Moreover, using targeted transgenic insertion, we show that promoter fragments containing REST or SNAIL binding sites are sufficient to recruit H3K27me3 in cis, while deletion of these sites results in loss of H3K27me3. These findings illustrate that the occurrence of TF binding sites can determine chromatin dynamics. Local determination of Polycomb activity by Rest and Snail motifs exemplifies such TF based regulation of chromatin. Furthermore, our results show that key TFs can be identified ab initio through computational modeling of epigenome datasets using a modeling approach that we make readily accessible