Abstract:Over the past two decades, a series of works have aimed at studying the problem of genome assembly: the process of reconstructing a genome from sequence reads. An early formulation of the genome assembly problem showed that genome reconstruction is NP-hard when framed as finding the shortest sequence that contains all observed reads. Although this original formulation is very simplistic and does not allow for mate-pair information, subsequent formulations have also proven to be NP-hard, and/or may not be guaranteed to return a correct assembly. In this paper, we provide an alternate perspective on the genome assembly problem by showing genome assembly is easy when provided with sufficient mate-pair information. Moreover, we quantify the number of mate-pair libraries necessary and sufficient for accurate genome assembly, in terms of the length of the longest repetitive region within a genome. In our analysis, we consider an idealized sequencing model where each mate-pair library generates pairs of error free reads with a fixed and known insert size at each position in the genome. Even in this idealized model, we show that accurate genome reconstruction cannot be guaranteed in the worst case unless at least roughly R/2L mate-pair libraries are produced, where R is the length of the longest repetitive region in the genome and L is the length of each read. On the other hand, if (R/L)+1 mate-pair libraries are provided, then a simple algorithm can be used to find a correct genome assembly easily in polynomial time. Although (R/L)+1 mate-pair libraries can be too much to produce in practice, the previous bounds only hold in the worst case. In our last result, we show that if additional conditions hold on a genome, a correct assembly can be guaranteed with only O(log (R/L)) mate-pair libraries.

Phase transition in the computational complexity of the shortest common superstring and genome assembly

Dna Directed Model for Shortest Superstring Problem

Computational Genome Analysis

Maximum likelihood genome assembly

A Case Study in Genome-Level Fragment Assembly

Theoretical Bounds on Mate-Pair Information for Accurate Genome Assembly

Complex genome assembly based on long-read sequencing

Bayesian genome assembly and assessment by Markov Chain Monte Carlo sampling

A Parallel DNA Fragment Assembly Algorithm Based on Eulerian Superpath Approach

Modeling biological problems in computer science: a case study in genome assembly

Complexity and Enumeration in Models of Genome Rearrangement

Folding@Home and Genome@Home: Using distributed computing to tackle previously intractable problems in computational biology

A Parallel DNA Fragment Assembly Algorithm Based on Eulerian Superpath Problem

Phased nanopore assembly with Shasta and modular graph phasing with GFAse

Computational complexity and pragmatic solutions for flexible tile based DNA self-assembly

CuMen: Clustering Sequences Based on Maximal Frequent Sequential Pattern and its Application in Genome Sequence Assembly

Parallel String Graph Construction and Transitive Reduction for De Novo Genome Assembly

Quantum Algorithms for the Shortest Common Superstring and Text Assembling Problems

Phasing or purging: tackling the genome assembly of a highly heterozygous animal species in the era of high-accuracy long reads

DBG2OLC: Efficient Assembly of Large Genomes Using Long Erroneous Reads of the Third Generation Sequencing Technologies

Genome assembly using quantum and quantum-inspired annealing