Question: [ChIP-seq. signal plot] How to align the based line of different samples?
1
gravatar for chiefcat
2.3 years ago by
chiefcat140
chiefcat140 wrote:

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 COMMENTlink modified 6 months ago by ricardo388920 • written 2.3 years ago by chiefcat140

I have a similar issue with my data, but the baseline are remarkably different between my samples. what is the best approach for dealing with this if the goal is to compare occupancy levels across different conditions?

plotProfile

ADD REPLYlink modified 6 months ago • written 6 months ago by ricardo388920

How did you normalize the data in the first place? This reads like your signal/noise levels are remarkably different.

ADD REPLYlink modified 6 months ago • written 6 months ago by ATpoint31k
2
gravatar for Devon Ryan
2.2 years ago by
Devon Ryan94k
Freiburg, Germany
Devon Ryan94k wrote:

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 COMMENTlink written 2.2 years ago by Devon Ryan94k

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 REPLYlink written 2.2 years ago by chiefcat140
1

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

ADD REPLYlink written 2.2 years ago by Devon Ryan94k
0
gravatar for Dave Gerrard
2.2 years ago by
Dave Gerrard190
Manchester
Dave Gerrard190 wrote:

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 COMMENTlink written 2.2 years ago by Dave Gerrard190

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

ADD REPLYlink written 2.2 years ago by Devon Ryan94k
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: 2159 users visited in the last hour