Delineating copy number and clonal substructure in human tumors from single-cell transcriptomes

Single-cell transcriptomic analysis is widely used to study human tumors. However it remains challenging to distinguish normal cell types in the tumor microenvironment from malignant cells and to resolve clonal substructure within the tumor. To address these challenges, we developed an integrative Bayesian segmentation approach called CopyKAT (Copynumber Karyotyping of Aneuploid Tumors) to estimate genomic copy number profiles at an average genomic resolution of 5Mb from read depth in high-throughput scRNA-seq data. We applied CopyKAT to analyze 46,501 single cells from 21 tumors, including triple-negative breast cancer, pancreatic ductal adenocarcinomas, anaplastic thyroid cancer, invasive ductal carcinoma and glioblastoma to accurately (98%) distinguish cancer cells from normal cell types. In three breast tumors, CopyKAT resolved clonal subpopulations that differed in the expression of cancer genes such as KRAS and signatures including EMT, DNA repair, apoptosis and hypoxia. These data show that CopyKAT can aid the analysis of scRNA-seq data in a variety of solid human tumors. We calculated single cell copy number profiles and inferred clonal structures of tumor cells from scRNAseq data using COPYKAT.