AmpliconArchitect icon indicating copy to clipboard operation
AmpliconArchitect copied to clipboard

Long read sequencing + AA to call ecDNA?

Open alhafidzhamdan opened this issue 3 years ago • 2 comments

Hi this might be off topic but have you tried using long read data on AA? I'm planning to generate HIFI reads for some of my samples and wonder if it's any good with AA? Probably not as HIFI won't have discordant reads, and if not, what would you recommend, and if you have any insights on this.

Thanks in advance! A

alhafidzhamdan avatar Mar 19 '21 11:03 alhafidzhamdan

Hi,

We have not tested this before. However, because AA uses discordant read pairs to identify the presence of junctions (which are refined by split-reads), the absence of discordant read pairs from a single-end sequencing technology such as PacBio or ONT will likely result in missing breakpoints in the region of interest. Others have alternatively attempted de novo assembly approaches for long-read resolution of ecDNA structures, but keep in mind that if the read does not span multiple breakpoints in a complex amplicon, the assembly may contain errors or incomplete contigs. The incorporation of long-reads into AA is a project which is on our radar.

Best, Jens

jluebeck avatar Mar 21 '21 22:03 jluebeck

Hi Jens,

Thanks for that. I guess what you are saying is that single-end sequencing will miss regions where an amplicon does not span multiple breakpoints hence may miss these. Do keep me updated if you end up developing such a tool to integrate/refine AA with long read data. Thanks Al

alhafidzhamdan avatar Mar 24 '21 10:03 alhafidzhamdan