SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
short read length on XL+ vlee2 454 Pyrosequencing 26 04-23-2012 06:43 AM
short read assembly john@wurbio De novo discovery 1 05-31-2011 09:18 AM
The best short read aligner Deutsche Bioinformatics 4 04-14-2011 07:12 PM
New Short Read Aligner sparks Bioinformatics 48 08-26-2009 08:01 AM

Reply
 
Thread Tools
Old 02-06-2014, 10:25 AM   #1
lan_ub
Junior Member
 
Location: New york state

Join Date: Jun 2013
Posts: 2
Question How short an Illumina read needs to be discard?

I'm working on ancient DNA mitochondrial genome mapping. We extracted genomic DNA from "fossil" bones. Then we did paired-end Illumina Hiseq shotgun sequencing using single stranded library. (Based on protocol from this paper: http://www.nature.com/nprot/journal/....2013.038.html)
The problem for ancient DNA is that, most of the DNA fragments we extracted are exogenous DNA (from microbes or human or other organisms). What we want to do next is to try to map reads to a reference mitochondrial genome.

At first, I trimmed the adapters. Then I discarded the reads less than 20bp. After using BWA, only a few hundred reads can be mapped to the reference. If I keep all the reads, no matter how short they are, more than 10,000 reads can be mapped. So, my question is, what is the minimum lenth for an Illumina read that can be used? If there's any paper discussed this question that would be fantastic. Thanks in advance!
lan_ub is offline   Reply With Quote
Old 02-06-2014, 10:54 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,059
Default

Have you checked to see what other groups have done for difficult samples (e.g. neandrathals, denisovans etc) such as yours?

If only something smaller than 20 bp is hybridizing then where did the rest of the DNA in that read originate from (artifact)?
GenoMax is offline   Reply With Quote
Old 02-09-2014, 02:50 PM   #3
jimmybee
Senior Member
 
Location: Adelaide, Australia

Join Date: Sep 2010
Posts: 119
Default

Quote:
Originally Posted by lan_ub View Post
I'm working on ancient DNA mitochondrial genome mapping. We extracted genomic DNA from "fossil" bones. Then we did paired-end Illumina Hiseq shotgun sequencing using single stranded library. (Based on protocol from this paper: http://www.nature.com/nprot/journal/....2013.038.html)
The problem for ancient DNA is that, most of the DNA fragments we extracted are exogenous DNA (from microbes or human or other organisms). What we want to do next is to try to map reads to a reference mitochondrial genome.

At first, I trimmed the adapters. Then I discarded the reads less than 20bp. After using BWA, only a few hundred reads can be mapped to the reference. If I keep all the reads, no matter how short they are, more than 10,000 reads can be mapped. So, my question is, what is the minimum lenth for an Illumina read that can be used? If there's any paper discussed this question that would be fantastic. Thanks in advance!
Generally its 25bp. A 20mer will map to a variety of different genomes, if given the chance, whereas 25mers are more stringent are more likely to be "real" when mapped to your closest reference.
jimmybee is offline   Reply With Quote
Reply

Tags
ancient dna, illumina reads, minimum length, mitochondial genome

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 09:37 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO