Abstract 213: Phylogenetic reconstruction across 303 metastatic tumor samples using Ultima whole-genome sequencing dramatically increases subclonal resolution
Authors: Hess Julian , Martin Elizabeth , Soifer Ilya , Benjamin Hila , Miller Mendy , Cibulskis Carrie , Danysh Brian P. , Coole Matthew , Gabriel Stacey , Juric Dejan , Lipson Doron , Getz Gad 
Published: April 2023

Comprehensive characterization of tumor evolution is essential for understanding the drivers of metastasis and treatment resistance, which are still largely unknown. Studying evolution and resistance requires large cohorts of patients; for each patient, a comprehensive phylogenetic tree is used to follow clones over time and space. Deep whole-genome sequencing (WGS) enables the creation of more robust phylogenies, since (i) each clone is characterized by far more clone-specific mutations (compared to whole-exome or panel sequencing); and (ii) WGS yields more accurate clonal and sub-clonal allelic copy-number estimates, enabling more precise estimation of the fraction of cancer cells that harbor each mutation.

Due to recent development of more affordable WGS by Ultima Genomics, we set out to test our ability to construct phylogenetic trees based on Ultima WGS. Ultima sequencing has demonstrated high sensitivity and specificity for detecting germline polymorphisms, but its performance on somatic mutation calling has not yet been rigorously studied.

First, we adjusted our somatic mutation detection pipeline to work with Ultima-generated WGS data and tested its performance using well-studied cell line tumor/normal pairs. We assembled a ground truth set of somatic mutations based on deep sequencing data from multiple platforms for the breast cancer cell line pairs (HCC1954/HCC1954-BL and HCC1143/HCC1143-BL) and the melanoma tumor/normal cell line pair (COLO829/COLO829-BL). We found that mutation detection performance was on-par with other platforms for single nucleotide variations.

Next, we identified 32 patients from various cancer types (breast, lung, cholangiocarcinoma, and melanoma) treated with a variety of therapies. We collected between 4 to 14 samples per patient at autopsy (total of 303), and sequenced them with Ultima to an average coverage of 63x (range 48x-91x). We identified mutations using our Ultima-adjusted WGS pipeline and constructed phylogenetic trees using our PhylogicNDT suite of tools. We then compared the results to trees previously constructed from mostly whole-exomes (and a few whole-genomes) sequenced on Illumina. We demonstrate that the trees generated from Ultima WGS data are consistent with the previously generated trees, and proved to be much more detailed.

At commercialization, we anticipate that WGS for projects like this will be performed at ~1/6th the cost of current NGS offerings (on par with current WES costs). Therefore, it will enable larger studies and more accurate phylogenetic reconstruction, which will advance the study of tumor evolution and resistance.

Citation Format: Julian Hess, Elizabeth Martin, Ilya Soifer, Hila Benjamin, Mendy Miller, Carrie Cibulskis, Brian P. Danysh, Matthew Coole, Stacey Gabriel, Dejan Juric, Doron Lipson, Gad Getz. Phylogenetic reconstruction across 303 metastatic tumor samples using Ultima whole-genome sequencing dramatically increases subclonal resolution [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 213.