DESeq2 analyze MeRIP-seq
1
0
Entering edit mode
BioEpi • 0
@11b02720
Last seen 4 hours ago
United States

Would you like to tell me whether DESeq2 can analyze m6A-seq to detect differential m6A peaks? Other papers firstly used input reads to normalize m6A-seq reads, and then performed fisher’s exact test. It seems that DESeq2 can not accept normalized reads.

DESeq2 MeRIP-seq • 267 views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 1 hour ago
United States

In DESeq2 instead of pre-scaling counts, you would just include the assay type as a design variable. Then you are essentially doing a paired analysis with assay type as the condition effect. E.g. ~sample + assay will compare, within each sample, the enriched vs input LFC. I would look at MA plots though as a diagnostic, and consider if you have features could be used to define controlGenes, but we don't have any automated workflows to do this. See ?estimateSizeFactors for description of controlGenes, and this step can then be run before DESeq if you want to specify that certain features be used for estimating the size factors.

ADD COMMENT
0
Entering edit mode

Hello Michal, Thanks for your reply. Maybe what I describled was not clearly. I have 12 samples:
case group:
case_input: (case1_input.bam,case2_input.bam,case3_input.bam)
case_m6A: (case1_IP.bam,case2_IP.bam,case3_IPbam)
control group:
control_input: (control1_input.bam,control2_input.bam,control3_input.bam)
control_m6A: (control1_IP.bam,control2_IP.bam,control3_IPbam)
I want to analyze m6A-seq to detect differential m6A peaks between case and control groups. For m6A-seq, input samples must be used to normalize/control IP samples. So I want to know whether DESeq2 can regress out/control input samples. Thanks again.

ADD REPLY
0
Entering edit mode

You can accommodate additional designs than just enriched vs input by using interaction terms. E.g. you are asking with your comment whether you can see enrich / input changes across another variable. That's just adding an interaction term. Do you have a local statistician that you could work with on choosing a design and interpreting results?

ADD REPLY
0
Entering edit mode

do you mean design (~assay(IP/input) + condition (case/control))?

ADD REPLY
0
Entering edit mode

Sure, that’s one way to code this up. It would be good to discuss an analysis plan with a statistician but we and many groups have used this approach to analyze enrichment or differential enrichment within DESeq2.

ADD REPLY

Login before adding your answer.

Traffic: 221 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.6