Abstract:Motivation: Some economically important plants including wheat and cotton have more than two copies of each chromosome. With the decreasing cost and increasing read length of next-generation sequencing technologies, reconstructing the multiple haplotypes of a polyploid genome from its sequence reads becomes practical. However, the computational challenge in polyploid haplotyping is much greater than that in diploid haplotyping, and there are few related methods.Results: This article models the polyploid haplotyping problem as an optimal poly-partition problem of the reads, called the Polyploid Balanced Optimal Partition model. For the reads sequenced from a k-ploid genome, the model tries to divide the reads into k groups such that the difference between the reads of the same group is minimized while the difference between the reads of different groups is maximized. When the genotype information is available, the model is extended to the Polyploid Balanced Optimal Partition with Genotype constraint problem. These models are all NP-hard. We propose two heuristic algorithms, H-PoP and H-PoPG, based on dynamic programming and a strategy of limiting the number of intermediate solutions at each iteration, to solve the two models, respectively. Extensive experimental results on simulated and real data show that our algorithms can solve the models effectively, and are much faster and more accurate than the recent state-of-the-art polyploid haplotyping algorithms. The experiments also show that our algorithms can deal with long reads and deep read coverage effectively and accurately. Furthermore, H-PoP might be applied to help determine the ploidy of an organism.Availability and Implementation: https://github.com/MinzhuXie/H-PoPGContact: xieminzhu@hotmail.comSupplementary information: Supplementary data are available at Bioinformatics online.

A generalized framework of AMOVA with any number of hierarchies and any level of ploidies

A Generalized Framework for AMOVA with Multiple Hierarchies and Ploidies

Polygene: Population Genetics Analyses for Autopolyploids Based on Allelic Phenotypes

Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data.

A mechanistic model of linkage analysis in allohexaploids

Inference of Individual Ploidy Level Using Codominant Markers

Performing Parentage Analysis in the Presence of Inbreeding and Null Alleles.

Vcfpop: Performing Population Genetics Analyses for Autopolyploids and Aneuploids Based on Next-Generation Sequencing Datasets.

Performing Parentage Analysis for Polysomic Inheritances Based on Allelic Phenotypes

VCFPOP: performing population genetics analyses for polyploids and anisoploids based on next-generation sequencing variant calling dataset

Estimating Pairwise Relatedness Between Individuals with Different Levels of Ploidy

The Analysis of Polyploid Genetic Data

A Multivalent Three-Point Linkage Analysis Model of Autotetraploids

A Maximum-Likelihood Estimation of Pairwise Relatedness for Autopolyploids

Multilocus Tetrasomic Linkage Analysis Using Hidden Markov Chain Model

A Unifying Framework for Bivalent Multilocus Linkage Analysis of Allotetraploids.

H-Pop and H-Popg: Heuristic Partitioning Algorithms for Single Individual Haplotyping of Polyploids

Genotyping single nucleotide polymorphisms and inferring ploidy by amplicon sequencing for polyploid, ploidy‐variable organisms

Model and Algorithm for Linkage Disequilibrium Analysis in a Non-Equilibrium Population.

Analysis of distributional variation through multi-scale Beta-Binomial modeling

A Pairwise Relatedness Estimator for Polyploids