Cancer subclone detection based on DNA copy number in single cell and spatial omic sequencing data

Abstract

In cancer, somatic mutations such as copy number alterations (CNAs) accumulate during disease progression and lead to functional intra-tumor heterogeneity that can influence the efficacy of cancer therapy. Therefore, studying the functional characteristics and spatial distribution of genetically distinct subclones is crucial to the understanding of tumor evolution and the design of cancer treatment. Here, we present Clonalscope, a method for subclone detection using copy number profiles that can be applied to spatial transcriptomics (ST) data and data from single-cell sequencing platforms such as scRNA-seq and scATAC-seq. Clonalscope implements a nested Chinese restaurant process to identify de novo subclones within one or multiple samples from the same patient. Clonalscope incorporates prior information from paired whole-genome or whole-exome sequencing (WGS/WES) data to achieve more reliable subclone detection and malignant cell labeling. On scRNA-seq and scATAC-seq data from four gastrointestinal tumor samples, Clonalscope successfully labeled malignant cells and identified genetically different subclones, which were validated in detail using matched scDNA-seq data. On ST data from a squamous cell carcinoma and two invasive ductal carcinoma samples, Clonalscope successfully labelled malignant spots, traced subclones between associated datasets, and identified spatially segregated subclones expressing genes associated with drug resistance and survival.