Question: Clustering of enhancer sites!!
0
gravatar for asharmatelome
5.2 years ago by
United States
asharmatelome60 wrote:

Hi, I have list of 8 TF's and they bind to various enhancers in the genome.

I have this data in excel and now I want to make a matrix where I want to see

TF A+B= X number of enhancers,

TF A+B+C= X +/- n number of enhancers,

Finally I want to cluster everything

can anyone help?

Thanks

 

R blog assembly genome • 1.4k views
ADD COMMENTlink modified 5.2 years ago • written 5.2 years ago by asharmatelome60

Thanks,

I tried Cluster3, but I dont I believe I have to download UniConvertor to get the data. I tried downloading this but I couldn't because it require some additional software package sk1.

I am trying with intersectBed but looks like it runs on 'R' and I have to accept that this is the first time I am trying to use this 'code' language. In my PhD I spend 6 yrs on bench and now I need to do these analysis.

 

ADD REPLYlink written 5.2 years ago by asharmatelome60

intersectbed will not run on R. it will work on unix command line. (http://bedtools.readthedocs.org/en/latest/)

Cluster 3.0 is easy to download as a GUI application on Windows/Mac. Why do you have to download Uniconvertor ? ( ps. i dont know what Uniconvertor is)

ps. this thread -How To Make Visual Graphs To Represent Common Transcription Factor Binding Sites In Different Enhancers? might help if you can get R working for you.

ADD REPLYlink modified 5.2 years ago • written 5.2 years ago by aditi.qamra260

It would be easier if you would show us how your data looks like. Can you update your question?

ADD REPLYlink written 5.2 years ago by PoGibas4.8k
2
gravatar for Sukhdeep Singh
5.2 years ago by
Sukhdeep Singh9.8k
Netherlands
Sukhdeep Singh9.8k wrote:

Once you have information in the form of the co-ordinates (which enhancer sites, TF A or B binds), use [intersectBed](http://bedtools.readthedocs.org/en/latest/content/tools/intersect.html) or multiIntersectBed, to get this information in the form of table, then you can make a heatmap after that.

Else use [Cluster3](http://bonsai.hgc.jp/~mdehoon/software/cluster/) for clustering with the data, where rows contain information for each enhancer and columns are different TF's. The values could me how strongly they bind or how much overlap or how close they are, you can decide that.

ADD COMMENTlink written 5.2 years ago by Sukhdeep Singh9.8k
1
gravatar for aditi.qamra
5.2 years ago by
aditi.qamra260
Toronto
aditi.qamra260 wrote:

one quick way would be to use intersectBed with your combined list of enhancer file and different combinations of TFs. (eg. A+B, A+B+C..etc)

Im sure there is a more elegant solution with bedtools though.

For clustering - depends on what your question is. maybe you can explain further.

 

 

ADD COMMENTlink written 5.2 years ago by aditi.qamra260
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 930 users visited in the last hour