Abstract:Modern pangenome graphs are built using haplotype-resolved genome assemblies. During read mapping to a pangenome graph, prioritizing alignments that are consistent with the known haplotypes has been shown to improve genotyping accuracy. However, the existing rigorous formulations for sequence-to-graph co-linear chaining and alignment problems do not consider the haplotype paths in a pangenome graph. This often leads to spurious read alignments to those paths that are unlikely recombinations of the known haplotypes. In this paper, we develop novel formulations and algorithms for haplotype-aware sequence alignment to an acyclic pangenome graph. We consider both sequence-to-graph chaining and sequence-to-graph alignment problems. Drawing inspiration from the commonly used models for genotype imputation, we assume that a query sequence is an imperfect mosaic of the reference haplotypes. Accordingly, we extend previous chaining and alignment formulations by introducing a recombination penalty for a haplotype switch. First, we solve haplotype-aware sequence-to-graph alignment in (| | | | |ℋ|) time, where is the query sequence, is the set of edges, and ℋ is the set of haplotypes represented in the graph. To complement our solution, we prove that an algorithm significantly faster than (| | | | |ℋ|) is impossible under the Strong Exponential Time Hypothesis (SETH). Second, we propose a haplotype-aware chaining algorithm that runs in (|ℋ| log |ℋ| ) time after graph preprocessing, where is the count of input anchors. We then establish that a chaining algorithm significantly faster than (|ℋ| ) is impossible under SETH. As a proof-of-concept of our algorithmic solutions, we implemented the chaining algorithm in the Minichain aligner ( ). We demonstrate the advantage of the algorithm by aligning sequences sampled from human major histocompatibility complex (MHC) to a pangenome graph of 60 MHC haplotypes. The proposed algorithm offers better consistency with ground-truth recombinations when compared to a haplotype-agnostic algorithm.

Graphasing: phasing diploid genome assembly graphs with single-cell strand sequencing

Phasing Diploid Genome Assembly Graphs with Single-Cell Strand Sequencing

Haplotype-resolved de novo assembly with phased assembly graphs

Haplotype-resolved de novo assembly using phased assembly graphs with hifiasm

Chromosome-scale, haplotype-resolved assembly of human genomes

Fully phased human genome assembly without parental data using single-cell strand sequencing and long reads

Learning to Untangle Genome Assembly with Graph Convolutional Networks

Simultaneous de novo calling and phasing of genetic variants at chromosome-scale using NanoStrand-seq

GraphPlas: Refined Classification of Plasmid Sequences Using Assembly Graphs

Phased nanopore assembly with Shasta and modular graph phasing with GFAse

Scalable telomere-to-telomere assembly for diploid and polyploid genomes with double graph

Pangenome graph construction from genome alignments with Minigraph-Cactus

Robust haplotype-resolved assembly of diploid individuals without parental data

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

New algorithms for accurate and efficient de novo genome assembly from long DNA sequencing reads

Gapless assembly of complete human and plant chromosomes using only nanopore sequencing

Haplotype-aware sequence alignment to pangenome graphs

Determination of complete chromosomal haplotypes by bulk DNA sequencing

Haplotype-resolved assembly of diploid genomes without parental data

ComHapDet: A Spatial Community Detection Algorithm for Haplotype Assembly

Efficient whole genome haplotyping and high-throughput single molecule phasing with barcode-linked reads