How can I quantify degradation of DNA present in the sample after sequencing?

1

Entering edit mode

23 months ago

ilykos ▴ 10

What tools are there to approximately quantify the DNA degradation in sequenced data post-alignment (.bam) or on read level (.fq)?

The goal is not to have an absolute metric of degradation, but to compare and stratify different samples in relation to each other. The less interpretation burden the better.

An example of a tool like this would be mapDamage, and I am looking for other options.

qc sequencing DNA • 738 views

ADD COMMENT • link updated 23 months ago by Ram 44k • written 23 months ago by ilykos ▴ 10

1

Entering edit mode

Are you using Illumina short read sequencing? If you are, I'm not sure you will be able to infer the degree of damage, as short DNA fragments can persist in nature for a long time. You could look into softwares used for ancient DNA analysis.

ADD REPLY • link 23 months ago by mark.ziemann ★ 1.9k

0

Entering edit mode

Thanks, Mark. Yes, all my data is PE150. At the moment I am looking at the output from pyDamage and mapDamage. The latter manages to distinguish between normal-embedded-old, tumor-embedded-old, normal-embedded-recent, and fresh samples, which is promising.

Could you explain further the connection between persistence of DNA fragments and the ability to evaluate the level of degradation?

ADD REPLY • link 23 months ago by ilykos ▴ 10

Login before adding your answer.