search

iprada / Circle-Map

A method for circular DNA detection based on probabilistic mapping of ultrashort reads
MIT License
62 stars 19 forks source link

ecDNA only from Chr10 and Chr11 #63

Open lanliting opened 3 years ago

lanliting commented 3 years ago

Hi: I am using Circle-Map to detect ecDNA from circle-seq data(paired-end). However, strangely, the output bedfiles show that ecDNAs only exist in Chr10 and Chr11, and another sample only in Chr10,11,12,13,14,15. I don't know if it is the truth or is caused by some errors. My input BAM files contains all the chromosome. And I check the log files, it shows some warnings: 6 Can these warnings contribute to the results above?

iprada commented 3 years ago

Hi,

Based on the time it took to finish, you had very little sequencing depth. Can you elaborate more on your data? Sequencing depth, circle enrichement procedure etc.

best,

Inigo

lanliting commented 3 years ago

Thank you，Inigo Our circle enrichment was based on column chromatography followed by digestion of remaining linear chromosomal DNA with exonuclease. Then we use rolling circle amplification and shear the amplified DNA to an average target peak size of 300 bp for library preparation and sequencing(150bp, paired). We have sequenced 11G bases for each sample.

Sent with GitHawk

lanliting commented 3 years ago

I still did not figure out this issue, could you give some instructions?

iprada commented 3 years ago

Dear Ianliting,

Your experimental set-up looks fine to me. However, the issue you are encountering doesn't seem related to Circle-Map per se. As it does not seem to be a software issue. In case you haven't done it already, I advice you to contact some experts on circular DNA purification that could help troubleshooting your project.

best,

Iñigo