Question: [ChIP-seq. signal plot] How to align the based line of different samples?
0
Entering edit mode

Hi all,

Does anyone have an idea how to make the baseline of different samples to the same level in a ChIP-seq. signal plot?

I've been using deepTools for normalization and plotting (have tried different normalization method provided in deepTools e.g. 1x normalize, log2 ratio, SES normalization but still get the same pattern as the plot below).

Any suggestion is welcome, thanks!!

enter image description here

ADD COMMENTlinkeditmoderate 2.1 years ago chiefcat • 100 • updated 9 months ago Biostar 20
0
Entering edit mode

There's no way to do this directly in deepTools, you'll need to save the data underlying the plot and fudge it in R (I have a feeling you're trying to over-interpret slight differences in ensemble enrichment over neighboring background between time points...be careful about doing that).

ADD COMMENTlinkeditmoderate 2.1 years ago Devon Ryan 90k
Entering edit mode
0

Thanks Devon. Yes, I also think it is dangourous to interpret the slight differences (I will try to convince my PI not doing that as well. I think it is much worth to spend time on searching real differences).

By the way, do you mean using --outFileNameData to save the data underpaying the plot? Which package can fudge the profiles? Thanks!

ADD REPLYlinkeditmoderate 2.1 years ago
chiefcat
• 100
Entering edit mode
1

Yes, --outFileNameData, you can modify the data in base R (you're just adding/removing a small offset).

ADD REPLYlinkeditmoderate 2.1 years ago
Devon Ryan
90k
0
Entering edit mode

I would suggest NOT creating lines on the same plot because for many of these kinds of differences it is hard to normalise out the library differences. In my opinion, and difference is better shown with two (or more) side by side heatplots running through the transcription start site. I've not used deepTools recently but know that ngsplot can produce side-side-side heatplots of two or more samples with the genes/transcripts in the same order on each (and optionally clustered and/or using ranks). These plots show better the co-variation across reference points and are relatively robust to differences in library size or distribution.

ADD COMMENTlinkeditmoderate 2.1 years ago Dave Gerrard • 190
Entering edit mode
0

What you're describing is the default output of plotProfile and plotHeatmap from deepTools.

ADD REPLYlinkeditmoderate 2.1 years ago
Devon Ryan
90k

Login before adding your answer.

Powered by the version 2.0