This is a test version of Biostars. For the public version, visit https://www.biostars.org.
macs2 "Effective genome size" for repetitive genomes

Hi everyone

I'm working on a chip-seq experiment in Wheat, which has a very large and repeptitive genome.

I'm a bit baffled by the "effective genome size" parameter in macs2. I understand it is related to the repetitiveness of the genome but I'm not sure how to calculate it. I've tried GEM but it gave me an error, so in parallel to trying to solve the GEM problem, maybe someone has an alternative?

Secondly, if I'm looking for peaks in repeptitive as well as non-repetitve regions of the genome, I thought maybe I should use the full length rather than the mappable length. Am I correct?

Finally - if I have a control sample (no antibody), can that be used to estimate the mappability of the genome?

Thanks!

chip-seq
• 5,344 views • link
updated by GouthamAtla GouthamAtla • written by Menachem Sklarz • Menachem Sklarz

1 answer

Effective genome size is after removing the repetetive elements in the genome. So you need to get the uniqely mappable regions.

Though directly not related, few options are given here https://github.com/fidelram/deepTools/wiki/General-deepTools-FAQs#effGenomeSize

Of which GEM-Mappability Calculator would be useful for you.

You may not be able to use the control sample for calculating the mappability as it would not cover the entire genome.

• 1 views • link
written by GouthamAtla GouthamAtla

Thanks for the input and for the links.

Do you think the "effective genome size" should be calculated the same way I'm doing the mapping? For example, if I'm retaining only uniquely mapped reads, then I should calculate the mappability as uniquely mappable regions but if I'm retaining also reads that mapped twice or three times then maybe I should determine the mappability as the regions that are mappable two- or three times?

Is there a site you know of that explains the statistics behind the effective size?

Thanks

• 1 views • link
updated by Ram Ram • written by Menachem Sklarz • Menachem Sklarz

Log in to answer this question.

More posts like this