CRIT/workflow

From GersteinInfo

Revision as of 21:12, 30 January 2011 by Tara (Talk | contribs)
Jump to: navigation, search

Contents

Transcription Factor Example

Motivation and Problem Set Up

Cis regulatory elements as a means of regulating gene expression have been extensively studied. However, beyond such motifs, are there inherent properties of the targets themselves that make them more or less likely to be regulated by a given class of transcription factors? As an example, do essential transcription factors preferentially regulate essential targets? Are there genome composition features such as GC or codon bias that influence which targets are regulated by which TFs?

Input Data

Here, we use three different datasets as shown.

Data Input Set up

These objects are named as follows in the R dataset:

(1) T: Transcription factors and their associated properties

(2) C: Connector Matrix matching transcription factors to their associated targets

(3) G: Gene targets and their associated properties

T and G are both post processed from:

Y. Xia, E. A. Franzosa, and M. B. Gerstein. Integrated assessment of genomic correlates of protein evolutionary rate. PLoS Comput Biol, 5(6):e1000413–e1000413, 2009.

C is post processed from:

C. T. Harbison, et al. Transcriptional regulatory code of a eukaryotic genome. Nature, 431(7004):99–104, 2004.

As in Harbison et al, p<.001 was used to indicate a TF-gene target. We binarized the matrix such that any TF-gene pair with a pval<.001 had a 1 and anything greater than this had a 0.

Example Code

#Example Code

#Load Data
load(file="TFExample.RData")

#Load CRIT functions
source(file="CRIT.R")

#Generate label for feature of interest - set x for column variable
tLabel<-initializer(T[,x], type="median")

#Determine set of targets sensitive to this feature
DC<-discriminator(C, tLabel, multCorrect=TRUE)

#Generate new label based on sensitivity identified in previous step
gLabel<-labelSlicer(DC, .05)

#Identify features that seem to discriminate between sens/insens targets
DG<-discriminator(t(G), gLabel, multCorrect=TRUE)
Personal tools