The restriction scaffold problem

Amir Ben-Dor, Richard M. Karp, Benno Schwikowski, Ron Shamir

Research output: Contribution to conferencePaperpeer-review

Abstract

Most shotgun sequencing projects undergo a long and costly phase of finishing, in which a partial assembly forms several contigs whose order, orientation and relative distance is unknown. We propose here a new technique that supplements the shotgun assembly data by cheap and simple complete restriction digests of the target. By computationally combining information from the contig sequences and the fragment sizes measured for several different enzymes, we seek to form a "scaffold" on which the contigs will be placed in their correct orientation, order and distance. We give a heuristic search algorithm for solving the problem and report on promising preliminary simulation results. The key to the success of the search scheme is the very rapid solution of two time-critical subproblems that are solved precisely in linear time. Our simulations indicate that with noise levels of some 3% relative error in measuring fragment sizes, using five enzymes, most datasets of 20 contigs can be correctly ordered, and the remaining ones have most of their pairs of neighboring contigs correct. Hence, the technique has a potential to provide real help to finishing. Even when the target clone remains unfinished, the ability to order and orient the contigs correctly makes the partial assembly both more accessible and more useful for biologists.

Original languageEnglish
Pages58-66
Number of pages9
DOIs
StatePublished - 2002
EventRECOMB 2002: Proceedings of the Sixth Annual International Conference on Computational Biology - Washington, DC, United States
Duration: 18 Apr 200221 Apr 2002

Conference

ConferenceRECOMB 2002: Proceedings of the Sixth Annual International Conference on Computational Biology
Country/TerritoryUnited States
CityWashington, DC
Period18/04/0221/04/02

Fingerprint

Dive into the research topics of 'The restriction scaffold problem'. Together they form a unique fingerprint.

Cite this