English

Minimum Segmentation for Pan-genomic Founder Reconstruction in Linear Time

Data Structures and Algorithms 2019-01-09 v2

Abstract

Given a threshold LL and a set R={R1,,Rm}\mathcal{R} = \{R_1, \ldots, R_m\} of mm haplotype sequences, each having length nn, the minimum segmentation problem for founder reconstruction is to partition the sequences into disjoint segments R[i1+1,i2],R[i2+1,i3],,R[ir1+1,ir]\mathcal{R}[i_1{+}1,i_2], \mathcal{R}[i_2{+}1, i_3], \ldots, \mathcal{R}[i_{r-1}{+}1, i_r], where 0=i1<<ir=n0 = i_1 < \cdots < i_r = n and R[ij1+1,ij]\mathcal{R}[i_{j-1}{+}1, i_j] is the set {R1[ij1+1,ij],,Rm[ij1+1,ij]}\{R_1[i_{j-1}{+}1, i_j], \ldots, R_m[i_{j-1}{+}1, i_j]\}, such that the length of each segment, ijij1i_j - i_{j-1}, is at least LL and K=maxj{R[ij1+1,ij]}K = \max_j\{ |\mathcal{R}[i_{j-1}{+}1, i_j]| \} is minimized. The distinct substrings in the segments R[ij1+1,ij]\mathcal{R}[i_{j-1}{+}1, i_j] represent founder blocks that can be concatenated to form KK founder sequences representing the original R\mathcal{R} such that crossovers happen only at segment boundaries. We give an optimal O(mn)O(mn) time algorithm to solve the problem, improving over earlier O(mn2)O(mn^2). This improvement enables to exploit the algorithm on a pan-genomic setting of haplotypes being complete human chromosomes, with a goal of finding a representative set of references that can be indexed for read alignment and variant calling.

Keywords

Cite

@article{arxiv.1805.03574,
  title  = {Minimum Segmentation for Pan-genomic Founder Reconstruction in Linear Time},
  author = {Tuukka Norri and Bastien Cazaux and Dmitry Kosolobov and Veli Mäkinen},
  journal= {arXiv preprint arXiv:1805.03574},
  year   = {2019}
}
R2 v1 2026-06-23T01:49:47.480Z