Predicting clone genotypes from tumor bulk sequencing of multiple samples

Motivation: Analyses of data generated from bulk sequencing of tumors have revealed extensive ge-nomic heterogeneity within patients. Many computational methods have been developed to enable the inference of genotypes of tumor cell populations (clones) from bulk sequencing data. However, the relative and absolute accuracy of available computational methods in estimating clone counts and clone genotypes is not yet known. Results: We have assessed the performance of nine methods, including eight previously-published and one new method (CloneFinder), by analyzing computer simulated datasets. CloneFinder, LICHeE, CITUP, and cloneHD inferred clone genotypes with low error (<5% per clone) for a majority of datasets in which the tumor samples contained evolutionarily-related clones. Computational methods did not perform well for datasets in which tumor samples contained mixtures of clones from different clonal lineages. Generally, the number of clones was underestimated by cloneHD and overestimated by Phy-loWGS, and BayClone2, Canopy, and Clomial required prior information regarding the number of clones. AncesTree and Canopy did not produce results for a large number of datasets. Conclusions: Deconvolution of clone genotypes from single nucleotide variant (SNV) frequency differ-ences among tumor samples remains challenging, so there is a need to develop more accurate compu-tational methods and robust software for clone genotype inference.