Single-cell open chromatin profiling via scATAC-seq has become a mainstream measurement of open chromatin in single-cells. Here we present epiAneufinder, an algorithm that exploits the read count information from scATAC-seq data to extract genome-wide copy number alterations (CNAs) for individual cells, allowing the study of CNA heterogeneity present in a sample at the single-cell level. Using different cancer scATAC-seq datasets, we show that epiAneufinder can identify intratumor clonal heterogeneity in populations of single cells based on their CNA profiles. We demonstrate that these profiles are concordant with the ones inferred from single-cell whole genome sequencing data for the same samples. EpiAneufinder allows the inference of single-cell CNA information from scATAC-seq data, without the need of additional experiments, unlocking a layer of genomic variation which is otherwise unexplored.
GrantsHelmholtz Association Impuls-und Vernetzungsfonds of the Helmholtz-Gemeinschaft We would like to thank Dr. Antonio Scialdone for discussions about calling CNAs in scRNA-seq data, and Anna Danese for help in embedding of scATAC-seq data. We thank Dr. Andres Castellanos for the discussions regarding the BCC tumor composition. We thank