Abstract:Abstract Motivation Exponential growth in sequencing databases has motivated scalable De Bruijn graph-based (DBG) indexing for searching these data, using annotations to label nodes with sample IDs. Low-depth sequencing samples correspond to fragmented subgraphs, complicating finding the long contiguous walks required for alignment queries. Aligners that target single-labelled subgraphs reduce alignment lengths due to fragmentation, leading to low recall for long reads. While some (e.g. label-free) aligners partially overcome fragmentation by combining information from multiple samples, biologically irrelevant combinations in such approaches can inflate the search space or reduce accuracy. Results We introduce a new scoring model, ‘multi-label alignment’ (MLA), for annotated DBGs. MLA leverages two new operations: To promote biologically relevant sample combinations, ‘Label Change’ incorporates more informative global sample similarity into local scores. To improve connectivity, ‘Node Length Change’ dynamically adjusts the DBG node length during traversal. Our fast, approximate, yet accurate MLA implementation has two key steps: a single-label seed-chain-extend aligner (SCA) and a multi-label chainer (MLC). SCA uses a traditional scoring model adapting recent chaining improvements to assembly graphs and provides a curated pool of alignments. MLC extracts seed anchors from SCAs alignments, produces multi-label chains using MLA scoring, then finally forms multi-label alignments. We show via substantial improvements in taxonomic classification accuracy that MLA produces biologically relevant alignments, decreasing average weighted UniFrac errors by 63.1%–66.8% and covering 45.5%–47.4% (median) more long-read query characters than state-of-the-art aligners. MLAs runtimes are competitive with label-combining alignment and substantially faster than single-label alignment. Availability and implementation The data, scripts, and instructions for generating our results are available at https://github.com/ratschlab/mla.

Near-optimal labeling schemes for nearest common ancestors

Nearest Common Ancestors: Universal Trees and Improved Labeling Schemes

A 6-approximation algorithm for computing smallest common aon-supertree with application to the reconstruction of glycan trees

Branch Code: A Labeling Scheme for Efficient Query Answering on Trees

An Improved Lower Bound on the Largest Common Subtree of Random Leaf-Labeled Binary Trees

Optimal L(3,2,1)-labeling of trees

Simplifying and Characterizing DAGs and Phylogenetic Networks via Least Common Ancestor Constraints

Tree height and the asymptotic mean of the Colijn-Plazzotta rank of unlabeled binary rooted trees

Listing, Verifying and Counting Lowest Common Ancestors in DAGs: Algorithms and Fine-Grained Lower Bounds

Adjacency labeling schemes and induced-universal graphs

On Two Measures of Distance between Fully-Labelled Trees

Two Metrics on Rooted Unordered Trees with Labels

Label-guided seed-chain-extend alignment on annotated De Bruijn graphs

On the correctness of Maximum Parsimony for data with few substitutions in the NNI neighborhood of phylogenetic trees

A 3-factor approximation algorithm for a Maximum Acyclic Agreement Forest on k rooted, binary phylogenetic trees

On the enumeration of leaf-labelled increasing trees with arbitrary node-degree

Achieving DNA Labeling Capacity with Minimum Labels through Extremal de Bruijn Subgraphs

Size-constrained Weighted Ancestors with Applications

Tree Learning: Optimal Algorithms and Sample Complexity

On the optimality of the neighbor-joining algorithm

Near-Optimal Resilient Labeling Schemes