de bruijn graph for paired end reads
1
0
Entering edit mode
3.9 years ago
Star70 • 0

I have a set of paired end reads from an unknown species. I want to assemble them to find the main genome. I have a program based on the basic de bruijn graph algorithm which is for single reads. Are there any other classic algorithm which works better in the case of paired end reads.? (Or any other version of de bruijn graph algorithm) Also I have 3 versions of these reads from 3 sample of the same species. Could them help me to improve my assembler?

RNA-Seq • 1.1k views
ADD COMMENT
1
Entering edit mode
3.9 years ago

There are heaps of De Bruijin graph implementations for paired end reads:

Velvet https://www.ebi.ac.uk/~zerbino/velvet/

ABySS http://www.bcgsc.ca/platform/bioinfo/software/abyss

MasuRCA http://www.genome.umd.edu/masurca.html

DISCOVAR (needs 250 bp paired reads, PCR-free) https://software.broadinstitute.org/software/discovar/blog/?page_id=23

SOAPdenovo http://soap.genomics.org.cn/soapdenovo.html

ALLPATHS-LG (needs matepaired data) http://software.broadinstitute.org/allpaths-lg/blog/

and many many many more

As for your three versions from three samples, do you expect the samples to be different? If they're replicates (identical) then you can treat them as three libraries in the assembler and assemble it all together.

ADD COMMENT
1
Entering edit mode

I'd like to also mention SPAdes, which is the only De Bruijin graph assembler I'm aware of that actually makes the graph using paired k-mers (k-bimers) from paired reads.

ADD REPLY
1
Entering edit mode

Thanks, I missed that! SPAdes has given me really good (best?) results with smaller (nonplant) genomes

ADD REPLY
0
Entering edit mode

It consistently gives us the best results with bacteria/archaea. We also routinely use it for metagenomes, but the resource requirements are so much higher than Megahit that it is only usable on low-complexity metagenomes, or metagenomes that have been highly processed (normalized, error-corrected, and low-depth reads removed). I've never tried it on a eukaryote, but it would not surprise me if it did a good job. Particularly on haploid euks like some fungi we assemble.

ADD REPLY

Login before adding your answer.

Traffic: 2435 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