Abstract:Broad copy number alterations (CNAs) at chromosomal band resolution have been used as cancer diagnostic and prognostic biomarkers for decades. These events have been characterized by cytogenetic imaging, an approach which is powerful for assessing heterogeneity but limited for locus fine mapping. While CNA detection by next-generation sequencing has become a standard analysis, existing methods rarely model CNA heterogeneity in bulk tumor samples and often require a paired normal sample for control. To overcome these limitations, we developed Seq2Karyotype (S2K, https://github.com/chenlab-sj/Seq2Karyotype), a new algorithm for in-silico karyotyping using single-sample whole-genome sequencing (WGS) data. S2K performs joint modeling of read-depth and allelic imbalance (AI) of high-quality heterozygous SNPs to identify reference diploid regions, followed by modeling and segmenting the deviation from the reference. Empirical coverage and AI of segmented regions are fitted to models of single and admixed CNAs to estimate clonality. The final karyotyping considers both the model fitness and minimization of evolution steps.To evaluate S2K's performance, we analyzed two cell lines commonly used for benchmark test: COLO829 (melanoma) and HCT1395 (breast cancer); both had single-cell (sc) WGS for validation. For COLO829, S2K replicated the four populations detected by scWGS but derived a different clonality estimate. The predominant clone, defined by loss of 1p, 10p, and chr18, was estimated to have 67% cellular fraction (CF) by bulk sample analysis of S2K in contrast to the 10% CF by scWGS analysis. In HCT1395, S2K identified three new CNAs present in >50% of the cells of the matching germline sample which were subsequently validated by FISH, karyotyping and SKY mapping.To demonstrate S2K's utility, we analyzed three diverse data sets: 17 neuroblastoma cell lines, 24 pediatric AML samples with karyotyping data, and two blood samples from children with myelodysplastic syndromes (MDS) known to harbor mosaic uniparental disomy (UPD). CNA-based intra-tumor heterogeneity was detected in 88% (15/17) of the neuroblastoma cell lines, comprised of 2-4 distinct populations with diverse ranges of CF (~10%-90%) and varying CNA patterns (e.g. admixture of tetraploid and diploid cells), which were validated by cytogenetics or scWGS. In the patient AML samples, S2K detected >95% of the previously reported cytogenetic events and 30% of additional copy-neutral loss-of-heterozygosity events. The mosaic UPD events in MDS patients were detected with projected clonality of 60% and 25%, respectively. These results not only demonstrate the accuracy of in-silico karyotyping performed by S2K but also reveal the dynamic intra-tumor heterogeneity in cancer cell lines, which may impact the design and interpretation of future experiments using these cell lines. Citation Format: Limeng Pu, Karol Szlachta, Virginia Valentine, Xiaolong Chen, Jian Wang, Dennis Kennetz, Daniel Putnam, Sivaraman Natarajan, Li Dong, Thomas Look, Marcin Wlodarski, Lu Wang, Steven Burden, John Easton, Xiang Chen, Jinghui Zhang. Seq2Karyotype (S2K): A method for deconvoluting heterogeneity of copy number alterations using single-sample whole-genome sequencing data [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular s); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl) nr 7419.

Single Nucleotide Polymorphism (SNP) and Antibody-based Cell Sorting (SNACS): A tool for demultiplexing single-cell DNA sequencing data

demuxSNP: supervised demultiplexing single-cell RNA sequencing using cell hashing and SNPs

Demuxalot: scaled up genetic demultiplexing for single-cell sequencing

coherent genetic demultiplexing in single-cell and single-nuclei experiments

Polled Digital Cell Sorter (p-DCS): Automatic identification of hematological cell types from single cell RNA-sequencing clusters

Demuxafy: improvement in droplet assignment by integrating multiple single-cell demultiplexing and doublet detection methods

A risk-reward examination of sample multiplexing reagents for single cell RNA-Seq

SCITO-seq: single-cell combinatorial indexed cytometry sequencing

DemoTape: Computational demultiplexing of targeted single-cell sequencing data

A statistical learning method for simultaneous copy number estimation and subclone clustering with single-cell sequencing data

MULTI-seq: sample multiplexing for single-cell RNA sequencing using lipid-tagged indices

Accurate identification of locally aneuploid cells by incorporating cytogenetic information in single cell data analysis

Accurate single-cell genotyping utilizing information from the local genome territory

SCCNAInfer: a robust and accurate tool to infer the absolute copy number on scDNA-seq data

msCNVS: medium throughput single cell copy number variation sequencing with barcoded library construction free of preamplification toward clinical implementation

Improved allele-specific single-cell copy number estimation in low-coverage DNA-sequencing

CASB: a concanavalin A‐based sample barcoding strategy for single‐cell sequencing

deMULTIplex2: robust sample demultiplexing for scRNA-seq

Ultrasensitive Multiplex Detection of Single Nucleotide Polymorphisms Based on Short-Chain Hybridization Combined with Online Preconcentration of Capillary Electrophoresis

Abstract 7419: Seq2Karyotype (S2K): A method for deconvoluting heterogeneity of copy number alterations using single-sample whole-genome sequencing data

scSNV-seq: high-throughput phenotyping of single nucleotide variants by coupled single-cell genotyping and transcriptomics