Arborist: Prioritizing Bulk DNA Inferred Tumor Phylogenies via Low-pass Single-cell DNA Sequencing Data

Abstract

Cancer arises from an evolutionary process that can be reconstructed from DNA sequencing and modeled by tumor phylogenies. High coverage bulk DNA sequencing (bulk DNA-seq) is widely available, but tumor phylogeny inference requires deconvolution, often resulting in non-uniqueness in the solution space. Single-cell DNA sequencing (scDNA-seq) holds potential to yield higher resolution tumor phylogenies, but the sparsity of emerging low-pass sequencing technologies poses challenges for the study of single-nucleotide variants. Increasing availability of data sequenced with both modalities provides an opportunity to capitalize on the advantages of these technologies. While inference methods exist for bulk DNA-seq and for low-pass scDNA-seq, no joint inference methods currently exist. As a first step, we propose a method named ARBORIST that prioritizes tumor phylogenies inferred via bulk DNA-seq using low-pass scDNA-seq data. ARBORIST takes as input a candidate set of trees with corresponding SNV clustering, along with variant and total read count data from scDNA-seq and uses variational inference to approximate a lower bound on the marginal likelihood of each tree in the candidate set. On simulated data, matching characteristics of current scDNA-seq data, ARBORIST outperforms both bulk and low-pass single-cell reconstruction methods. On a biological dataset, ARBORIST conclusively resolves the evolutionary relationship between different SNV clusters on a malignant peripheral nerve sheath tumor, which is supported by orthogonal validation via a proxy for copy number. ARBORIST provides a principled framework for integrating bulk DNA-seq and low-pass scDNA-seq data, improving confidence in tumor phylogeny reconstruction. Availabilityhttps://github.com/VanLoo-lab/Arborist