Comprehensive single-cell genome analysis at nucleotide resolution using the PTA Analysis Toolbox

Abstract

Detection of somatic mutations in single cells has been severely hampered by technical limitations of whole genome amplification. Novel technologies including primary template-directed amplification (PTA) significantly improved the accuracy of single-cell whole genome sequencing (WGS), but still generate hundreds of artefacts per amplification reaction. We developed a comprehensive bioinformatic workflow, called the PTA Analysis Toolkit (PTATO), to accurately detect single base substitutions, small insertions and deletions (indels) and structural variants in PTA-based WGS data. PTATO includes a machine learning approach to distinguish PTA-artefacts from true mutations with high sensitivity (up to 90% for base substitution and 95% for indels), outperforming existing bioinformatic approaches. Using PTATO, we demonstrate that many hematopoietic stem and progenitor cells of patients with Fanconi anemia, which cannot be analyzed using regular WGS technologies, have normal somatic single base substitution burdens, but increased numbers of deletions. Our results show that PTATO enables studying somatic mutagenesis in the genomes of single cells with unprecedented sensitivity and accuracy.