pHapCompass: Probabilistic Assembly and Uncertainty Quantification of Polyploid Haplotype Phase

TL;DR

pHapCompass introduces a probabilistic approach for assembling haplotypes in polyploid genomes, explicitly modeling read ambiguity and providing uncertainty quantification, which improves accuracy and realism over synthetic benchmarks.

Contribution

It develops a novel probabilistic algorithm for polyploid haplotype assembly that explicitly models read assignment ambiguity and quantifies uncertainty, addressing limitations of prior methods.

Findings

pHapCompass achieves competitive accuracy across various polyploid complexities.

The method provides reliable uncertainty quantification for haplotype phases.

Benchmarking shows improved performance over existing assemblers.

Abstract

Computing haplotypes from sequencing data, i.e. haplotype assembly, is an important component of molecular and population genetics problems, including interpreting the effects of genetic variation on complex traits and reconstructing genealogical relationships. Assembling the haplotypes of polyploid genomes remains a significant challenge due to the exponential search space of haplotype phasings and read assignment ambiguity; the latter challenge is particularly difficult for haplotype assemblers since the information contained within the observed sequence reads is often insufficient for unambiguous haplotype assignment in polyploid genomes. We present pHapCompass, probabilistic haplotype assembly algorithms for diploid and polyploid genomes that explicitly model and propagate read assignment ambiguity to compute a distribution over polyploid haplotype phasings. We develop graph theoretic algorithms to enable statistical inference and uncertainty quantification despite an exponential space of possible phasings. Since prior work evaluates polyploid haplotype assembly on synthetic genomes that do not reflect the realistic genomic complexity of polyploidy organisms, we develop a computational workflow for simulating genomes and DNA-seq for auto- and allopolyploids. Additionally, we generalize the vector error rate and minimum error correction evaluation criteria for partially phased haplotypes. Benchmarking of pHapCompass and several existing polyploid haplotype assemblers shows that pHapCompass yields competitive performance across varying genomic complexities and polyploid structures while retaining an accurate quantification of phase uncertainty. The source code for pHapCompass, simulation scripts, and datasets are freely available at https://github.com/bayesomicslab/pHapCompass.