Subtracting Background Signal Before DESeq Analysis
0
Entering edit mode
vanbelj • 10
@vanbelj-21216
Last seen 1 hour ago

I am looking for advice on how to appropriately subtract background signal from true signal before performing analysis with DESeq2.

I've performed a cut-and-run experiment, which is akin to ChIP-Seq and reveals protein-genome interactions. The technique results in background signal at locations where DNA is less protected, and this background can obscure signal from bonafide interactions - such as when performing DE analysis (I say DE, but we are talking about differential binding here).

My protein is recruited to DNA downstream of other proteins. To quantify the amount of non-specific background signal, I degraded an upstream protein and then performed cut-and-run. Signal remained throughout the genome - this represents non-specific background signal. Background is often large relative to a gain in signal in response to treatment, and so these changes are often obscured when performed DE analysis . Luckily, the amount of signal appears to be reproducible (low variation between replicates).

I quantified gene-counts with Salmon. For each gene, I intend to subtract the average background signal away from average total signal. However, I am unsure if this would be appropriate because Salmon does more than a simple read count.

After subtracting out background signal, I intend to use the counts data in DESeq2 analysis.

Thanks for any advice.

Salmon DESeq2 • 105 views
ADD COMMENTlink
0
Entering edit mode
@mikelove
Last seen 1 day ago
United States

So there are some threads here about contrasting against background:

DESeq2 testing ratio of ratios (RIP-Seq, CLIP-Seq, ribosomal profiling)

The idea is to test a ratio of ratios, but don't modify the counts themselves.

Does this help to get you going? Feel free to post followup questions.

ADD COMMENTlink
0
Entering edit mode

Hi Michael,

In the linked thread, four samples are mentioned: Input for A, Input for B, IP for A, IP for B

And a ratio of ratios if taken: (IP for B / Input for B) / (IP for A / Input for A)

However, what I have is more akin to: IP for A, IP for B, IP for C

, where C is the background. Would the following ratio of ratios be appropriate in this case? (A/C) / (B/C)

Intuitively I think this ratio is an appropriate way of taking the background into account, but I'm not sure if its valid for the math that is going on behind the scenes of DESeq2.

Thanks!

ADD REPLYlink
0
Entering edit mode

This has been asked before on the support site a few times, but I don't know how you'd find those threads.

If you have three groups, (A/C) / (B/C) is just a pairwise comparison of A and B, C drops away. There is no such thing as, 'comparison of B and A with respect to a group C'.

ADD REPLYlink
0
Entering edit mode

Hi Michael,

In comparing A and B directly, would DESeq2 identify the following as DE? For a given gene, assume counts of:

Sample A: 500 Sample B: 525 Sample C: 500

I'm trying to get an idea of the sensitivity of DESeq2 when background signal is high. Luckily I have low variation between replicates, which my aid in detection. I'm unsure if comparing 500 vs. 525 will flag as DE? My original thought was to subtract out average counts from background, to get:

Sample A: 0 Sample B: 25

, which I would expect DESeq2 to identify as DE (as long as counts pre-filter is set low enough).

ADD REPLYlink
0
Entering edit mode

With DESeq2 you cannot modify the counts. If you want to restrict the DE results between A and B with respect to C, you can do that by requiring that A > C and B > C with pairwise tests, but the above doesn't strike me as fitting into the DESeq2 paradigm.

ADD REPLYlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Traffic: 242 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.4