Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Metheor only calculate methylation diversity based on CpG site? #16

Open
feihongloveworld opened this issue May 5, 2023 · 5 comments

Comments

@feihongloveworld
Copy link

Metheor only calculate methylation diversity based on CpG site?

@dohlee
Copy link
Owner

dohlee commented May 5, 2023

Yes, it considers only CpG sites.

@feihongloveworld
Copy link
Author

but MHL was designed on MHB which contain more than 3 CpG sits. so I'm wondering weather MHLs based on a CpG sits are different with MHLs based on MHB

@dohlee
Copy link
Owner

dohlee commented May 5, 2023

That's right. MHL was designed for a genomic region in the original paper, but CpG-wise MHL was then proposed in Scherer et al., 2020, which is the benchmark paper for Metheor. Please refer to https://doi.org/10.1093/nar/gkaa120. Basically, you can think that a set of sequencing reads covering each CpG defines the genomic region considered for MHL calculation.

@feihongloveworld
Copy link
Author

I read Scherer's paper just now. whether PM and ME only are calculated based on the 4-CpG window? if the CpG number is more than 4 in a block, should I use overlapped windows with step 1 to get PM and ME?

@dohlee
Copy link
Owner

dohlee commented May 6, 2023

Yes, PM and ME values are assigned for four consecutive CpGs (i.e. CpG quartets). If you want to compute PM and ME in certain genomic block, you can simply average all PM/ME values for CpG quartets overlapping with that genomic block.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants