Trimming ancient genomic DNA without existing alignment - leeHom or Adapterremoval?
0
0
Entering edit mode
5.8 years ago
new15b ▴ 30

Do any of you have experience with either leeHom or Adapterremoval? If you were working with a denovo assembly, would it make more sense to use one or the other for very ancient DNA? As someone very new to bioinformatics, I appreciate the thorough documentation for Adapterremoval, plus I don't have any priori assemblies, so I question whether leeHom would provide any advantages.

If you guys have any tips, I'd greatly appeciate it.

genome sequencing • 1.5k views
ADD COMMENT
1
Entering edit mode

leeHom seems to be specifically designed for ancient DNA where as the other program seems to be a general adapter removal one (there are plenty of others in that category).

You could also look at BBMerge from BBMap suite which can do read merging taking into account errors. You can trim the adapters off as a part of that process or using bbduk afterwards..

ADD REPLY
0
Entering edit mode

I found this paper about Adapterremoval where they claim version 2 is very suitable for ancient DNA: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4751634/

What is your opinion about this?

I will look into BBMerge, thank you!

ADD REPLY
0
Entering edit mode

In that case you may have to try all three tools on your data to see which works best. Sometimes the application depends on data one has.

ADD REPLY
0
Entering edit mode

leeHom has the advantage of computing quality scores accurately for ancient DNA fragments and being the most accurate in terms of sequence length. I recently added to bioconda, let me know if you have problems using it.

ADD REPLY

Login before adding your answer.

Traffic: 1822 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6