How to use methyKit with two platform data with different coverage distribution

[liuyang2006]

Dear author,

I have two platform generated methylation data. The coverage of Platform1 is very high >25x, however, Platform2's data coverage is very low ~3-5x.

Each platform I have one sample in ctrl and one sample in test group.

I want to call DMC in each platform ctrl/test group, and compare the called CpGs between the two platform. In this case, can I directly use the Fisher test in methyKit, no matter the coverage difference of the two platform? Thank you for your help!

Best,
Yang

[alexg9010]

Hi @liuyangzzu,

In general, when processing methylation data, we advise a minimum coverage threshold of 10 reads per base, thus many sites from your second platform's data may already be filtered out. IMHO, given this strong difference in coverage, it would not be reasonable to perform the logistic regression-based tests.

To answer your question, when there is only one sample per sample group, calculateDiffMeth() automatically applies Fisher’s exact test.

Best
Alex

[al2na]

Yes, you can !! Please update us about the results. This is very interesting, would be good to know how concordant Rrbs and ONT

On Thu, Jul 7, 2022 at 9:53 PM Yang Liu ***@***.***> wrote: Dear author, I have two platform generated data, RRBS and ONT. The coverage of RRBS is very high >25x, however, ONT data coverage is very low ~3-5x. Each platform I have one sample in ctrl and one sample in test group. I want to call DMC in each platform ctrl/test group, and compare the CpG difference between the two platform. In this case, can I directly use the Fisher test in methyKit, no matter the coverage difference of the two platform? Thank you for your help! Best, Yang — Reply to this email directly, view it on GitHub <#260>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AAE32EKE5PHQUNJKTUJ72M3VS4YVJANCNFSM526TB43A> . You are receiving this because you are subscribed to this thread.Message ID: ***@***.***>

-- Sent from mobile, excuse the brevity