3 ms·
I would not have expected such a misunderstanding of the fruits of the knowledge of the human reference sequence from J Craig Venter. Frankly, it sounds like so
by carbocation 16y ago
I would not have expected such a misunderstanding of the fruits of the knowledge of the human reference sequence from J Craig Venter. Frankly, it sounds like sour grapes. For some reason, I feel inclined to try to explain some of the technical reasons why the public project's approach was necessary for the completion of the human genome.
There is broad agreement (perhaps "universal - 1") that Venter's sequencing approach would never have been able to achieve the level of completeness that we have now thanks to the Human Genome Project. This is because of duplicate DNA in the human genome, and the nature of short read lengths.
Let's try this with colors. Say you have a terrible sequencer that only reads 7bp at a time, max. You come across the following 2 strands of DNA:
GREENgattacaGREEN
REDgattacaRED
They're on 2 entirely different chromosomes, but you don't know this since there exists no reference sequence - you're mapping de novo, remember? Good. So you're get readouts like "GREENgattac" and "gattaca" and "attacaRED" - but you don't have any reads >7bp, so you never see what is on both sides of gattaca at once.
Thus, since you can never span gattaca, you will never know if your DNA in this region should be GREENgattacaRED + REDgattacaGREEN, or REDgattacaRED + GREENgattacaGREEN. For de novo alignment with short reads, this type of shortcoming is deadly.
Because the public consortium / Human Genome Project used much longer read lengths and a method involving bacterial artificial chromosomes, they were able to map the more difficult regions of the genome having repeats.
Edit I should have addressed the public consortium's approach more thoroughly, though to be clear I am less familiar with the approach and have never used it myself. Essentially, you chop up the 3 gigabase genome into a bunch of 150 kilobase segments (in this case, bacterial artificial chromosomes). You can map these back to their source chromosomes, of which we already knew we had 23. You then chop these BACs up into smaller pieces and sequence. In the worst case, you might be unsure of where some of your reads from each BAC map within the 150kb window that they come from, but hey, at least you know which chromosome and more or less which portion of the chromosome they come from. In other words, in the situation described above, the RED and the GREEN never co-occur, since they are on different BACs and are not sequenced together. Thus you know that RED/RED and GREEN/GREEN are the right pairings (and you know which belongs on which chromosome).
- ramanujan 16y agoSure, but Venter and Celera do deserve props for introducing whole genome shotgun, which became the de facto standard for getting most of the genome. This doesn't invalidate anything you're saying, because finishing is important, but the whole trend towards short reads and pushing more of the job to the computer was begun by Celera. Eugene Myers is/was an amazing computer scientist, and the public project would never have finished as quickly without the push from Celera. [beat] Also, all of this really deals with the technical aspects of obtaining a sequence, which is not really what Venter is talking about. He's saying that a good sequence isn't that useful! I'm not sure what he means with this, or whether he's just trying to be controversial.
- carbocation 16y agoI agree with you on all points. Like I said, I've never actually used the BAC approach but I've definitely had my hands all over short read data. I think that underscores your point about the popularity of short read length methods. It's good stuff - but definitely aided by the presence of a relatively complete reference. It is true that my post is somewhat off topic, too; I didn't want to dwell on Venter's naysaying.
- dasht 16y agocarbocation: Can you say a word or two about the current view of short reads for re-sequencing? I worked for a time on software for a particular lab to do alignment of short(-ish) reads like 13-bp+gap-of-100-200+7bp+gap-of-5-8+7bp+gap-of-100-200+13-bp. The idea was that you give me millions or more of those reads, and a reference like the human genome, and I give you back all of the alignments including partial matches (a few bp off, here and there). Interesting programming problem. This aren't the actual bp counts and gap sizes I was working with - they were in that range. I got concerned that the reads we were working with were too ambiguous -- similar to the problem you describe with de novo. I was having trouble doing the combinatorics to prove myself wrong and started asking: hey, where's the mathematical analysis that shows reads of this size can possibly work? Nobody could produce it and this led to some internal strife. That was a few years back. What's current thinking on what kind of reads plausibly work for re-sequencing? (If you happen to know. And where can I read something about the math that proves it.)