how to use Circle-map to get counts for a specific ecDNA? about circle-map

[iprada]

Hi,

Could you elaborate a bit more?

Do you wanna get the read counts for every ecDNA in your sample? If so, do you program in python?

[shaka-emperor]

exactly, I wanna get the read counts for every ecDNA in my sample, but I don't know how to do it...

[iprada]

Do you program in python?

[shaka-emperor]

sorry... I only know it's a programming language. if it helps, I'd love to learn it

[iprada]

If you program in python, you could use pysam (https://pysam.readthedocs.io/en/latest/api.html) to get the counts. I also thing that "bedtools coverage" does exactly what you want. If you do not know what I am talking about, I encourage you to find some bioinformatics collaborator that could help you.

Best,

Iñigo

[shaka-emperor]

Thanks a lot~~~ these clues are really valuable ~~~

[iprada]

You are very welcome.

Sorry I cannot be more helpful. I am finishing my PhD and I am quite busy. Feel free to open a new issue if you have more questions

[shaka-emperor]

Dear iprada
Hi, I used bedtools coverage to get read counts for a specific ecDNA. But how to use these reads to trace back the copy number of a specific ecDNA?
Besides, is it appropriate to directly use Circle-Map realignmy_unknown_circle.bed file (split reads + discordant reads, or just split reads) to indicate ecDNA copy number? Because my understanding is that these reads contain breakpoints. For instance, if this putative ecDNA have 5 split reads, then I can say the raw data contain 5 copies of this ecDNA?

[iprada]

Hi @shaka-emperor,

This outside my scope and it heavily depends on your experimental setting. If you do not have anybody with a rigurous statistics/quantitative bioinformatics background, I encourage you to find a collaborator that could help you there.

I am sorry I cannot help more on this.

Best,

Iñigo

[shaka-emperor]

Thank you anyway, you have helped me a lot~~