News for hackers. Stuff that matters.
Re:A farewell to bioinformatics (2012) by sampo on Sunday January 27, 2013 @05:28AM (#5123836)
Genome assembly relies heavily on sequence alignment. So: Is genome assembly hard just because sequence alignment is hard? Or would genome assembly present separate algorithmic problems even if there was a super-efficient solution to sequence alignment?
Re:A farewell to bioinformatics (2012) by alephnil on Sunday January 27, 2013 @11:19AM (#5124651)
It is far more difficult than sequence alignment. Sequence alignment has quadratic complexity, while fragment assembly is NP-hard. Se for example

http://scholar.google.com/scholar?cluster=131745416915434219...

Re:A farewell to bioinformatics (2012) by dalke on Sunday January 27, 2013 @06:15PM (#5126111)
Yes, for pairwise sequence alignment. The globally optimized multiple sequence alignment problem is NP-complete.
Re:A farewell to bioinformatics (2012) by dalke on Sunday January 27, 2013 @06:10PM (#5126087)
These are different sorts of alignments, with different sorts of math behind them.

Genome assembly is the shortest common super sequence problem. It involves finding the best rearrangement and overlap of reads which minimize the overall sequence, given the expected errors in the read technology. It would still be hard even if all of the reads were perfect.

Sequence alignment looks at two or more sequences in their entirety, and does a best fit alignment using a given model of how substitutions and gaps can occur. This model may be based on chemical or evolutionary knowledge.

A "super-efficient solution to sequence alignment" doesn't lead to a way to tell how the reads should be assembled into a single large sequence, even ignoring possible read errors.

Re:A farewell to bioinformatics (2012) by ufo on Sunday January 27, 2013 @10:27AM (#5124505)
An extra difficulty with genome assembly is that DNA often has lots and lots of repeated junk sequences that can confuse the algorithms. I don't work with bioinformatics to know how they usually get around this though.
Re:A farewell to bioinformatics (2012) by njbooher on Sunday January 27, 2013 @11:31AM (#5124692)
Repeats aren't necessarily junk (e.g. TAL Effectors http://en.wikipedia.org/wiki/TAL_effector#DNA_recognition). Resolving them requires long reads. PacBio is currently of interest as an alternative to Sanger sequencing for this, although the error rate of PacBio reads is a bit of an issue.
Re:A farewell to bioinformatics (2012) by apathy on Wednesday January 30, 2013 @08:52PM (#5142908)
pacbio is dead, they just don't know it yet. BGI (or somebody, doesn't matter, BGI is just the obvious candidate) would need to buy 50 SMART sequencers a year just for PacBio to stay in business. That seems unlikely given the lower cost and complexity of Illumina and Life sequencers