Resolving sub-clonal heterogeneity within cell-line growths by single cell sequencing genomic DNA

ABSTRACT We performed shallow single-cell sequencing of genomic DNA across 1,475 cells from a well-studied cell-line, COLO829, to resolve overall tumor complexity and clonality. This melanoma tumor-line has been previously characterized by multiple technologies and provides a benchmark for evaluating somatic alterations, though has exhibited conflicting and indeterminate copy number states. We identified at least four major sub-clones by discriminant analysis of principal components (DAPC) of single cell copy number data. Break-point and loss of heterozygosity (LOH) analysis of aggregated data from sub-clones revealed a complex rearrangement of chromosomes 1, 10 and 18 that was maintained in all but two sub-clones. Likewise, two of the sub-clones were distinguished by loss of 1 copy of chromosome 8. Re-analysis of previous spectral karyotyping data and bulk sequencing data recapitulated these sub-clone hallmark features and explains why the original bulk sequencing experiments generated conflicting copy number results. Overall, our results demonstrate how shallow copy number profiling together with clustering analysis of single cell sequencing can uncover significant hidden insights even in well studied cell-lines.