This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Assign BigWig scores to Bed file (gene list): ChIP-seq

Hey all,

I imagine this is an easy question, but I haven't been able to find the answer! I have a bed file list of all HG19 genes. I also have BigWig files of MACS2 peak calls from ChIP-seq data. I would like to iterate through my list of genes, and assign the "score" value of the bed file as the signal from the BigWig at the bed file's indicated "start" value.

I have used DeepTools to achieve a similar goal, but the computeMatrix output, "sorted/filtered regions," does not have a true "score".

Thank you!

chip-seq
• 4,124 views • link
updated by Devon Ryan Devon Ryan • written by bradbenjamin526 • bradbenjamin526

2 answers

It would have made a bit more sense to use bedtools intersect with the BED file from MACS2 and simply intersect that with the genes BED file.

• 0 views • link
written by Devon Ryan Devon Ryan

Convert BigWig to Wig via Kent Tools bigWigToWig and from Wig to BED via BEDOPS wig2bed. Then run BEDOPS bedmap to map scores to genes.

It's fast and you can use Unix pipes to make a clean workflow, once you get out of BigWig world and into BED:

$ wig2bed < in.wig | bedmap --echo --echo-map-score genes.bed - > answer.bed
• 0 views • link
written by Alex Reynolds Alex Reynolds

Log in to answer this question.

More posts like this