Question: Assign BigWig scores to Bed file (gene list): ChIP-seq
0
gravatar for bradbenjamin526
2.3 years ago by
bradbenjamin5260 wrote:

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 • 971 views
ADD COMMENTlink modified 2.3 years ago by Devon Ryan90k • written 2.3 years ago by bradbenjamin5260
0
gravatar for Devon Ryan
2.3 years ago by
Devon Ryan90k
Freiburg, Germany
Devon Ryan90k wrote:

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.

ADD COMMENTlink written 2.3 years ago by Devon Ryan90k

I'm trying this now, will report back with results!

ADD REPLYlink written 2.3 years ago by bradbenjamin5260
0
gravatar for Alex Reynolds
2.3 years ago by
Alex Reynolds28k
Seattle, WA USA
Alex Reynolds28k wrote:

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
ADD COMMENTlink modified 2.3 years ago • written 2.3 years ago by Alex Reynolds28k
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: 1521 users visited in the last hour