Deconvoluting multi-person biological mixtures and accurate characterization and identification of separated contributors using non-targeted single-cell DNA sequencing

The genetic characterization and identification of individuals who contributed to biological mixtures are complex and mostly unresolved tasks. These tasks are relevant in various fields, particularly in forensic investigations, which frequently encounters crime scene stains generated by more than one person. Currently, forensic mixture deconvolution is mostly performed subsequent to forensic DNA profiling at the level of the mixed DNA profiles, which comes with several limitations. Some previous studies attempted at separating single cells prior to forensic DNA profiling. However, these approaches are biased at selection of the cells and, due to their targeted DNA analysis on low template DNA, provide incomplete and unreliable forensic DNA profiles. We recently demonstrated the feasibility of performing mixture deconvolution prior to forensic DNA profiling through the utilization of a non-targeted single-cell transcriptome sequencing (scRNA-seq). In addition to individual-specific mixture deconvolution, this approach also allowed accurate characterisation of biological sex, biogeographic ancestry and individual identification of the separated mixture contributors. However, RNA has the forensic disadvantage of being prone to degradation, and sequencing RNA - focussing on coding regions - limits the number of single nucleotide polymorphisms (SNPs) utilized for genetic mixture deconvolution, characterization, and identification. These limitations can be overcome by performing single-cell sequencing on the level of DNA instead of RNA. Here, for the first time, we applied non-targeted single-cell DNA sequencing (scDNA-seq) by applying the scATAC-seq (Assay for Transposase-Accessible Chromatin with sequencing) technique to address the challenges of mixture deconvolution in the forensic context. We demonstrated that scATAC-seq, together with our recently developed De-goulash data analysis pipeline, is capable of deconvoluting complex blood mixtures of five individuals from both sexes with varying biogeographic ancestries. We further showed that our approach achieved correct genetic characterization of the biological sex and the biogeographic ancestry of each of the separated mixture contributors and established their identity. Furthermore, by analysing in-silico generated scATAC-seq data mixtures, we demonstrated successful individual-specific mixture deconvolution of i) highly complex mixtures of 11 individuals, ii) balanced mixtures containing as few as 20 cells (10 per each individual), and iii) imbalanced mixtures with a ratio as low as 1:80. Overall, our proof-of-principle study demonstrates the general feasibility of scDNA-seq in general, and scATAC-seq in particular, for mixture deconvolution, genetic characterization and individual identification of the separated mixture contributors. Furthermore, it shows that compared to scRNA-seq, scDNA-seq detects more SNPs from fewer cells, providing higher sensitivity, that is valuable in forensic genetics.

scAbsolute: measuring single-cell ploidy and replication status

Cancer cells often exhibit DNA copy number aberrations and can vary widely in their ploidy. Correct estimation of the ploidy of single-cell genomes is paramount for downstream analysis. Based only on single-cell DNA sequencing information, scAbsolute achieves accurate and unbiased measurement of single-cell ploidy and replication status, including whole-genome duplications. We demonstrate scAbsolute's capabilities using experimental cell multiplets, a FUCCI cell cycle expression system, and a benchmark against state-of-the-art methods. scAbsolute provides a robust foundation for single-cell DNA sequencing analysis across different technologies and has the potential to enable improvements in a number of downstream analyses.

Inferring single-cell copy number profiles through cross-cell segmentation of read counts

BACKGROUND

Copy number alteration (CNA) is one of the major genomic variations that frequently occur in cancers, and accurate inference of CNAs is essential for unmasking intra-tumor heterogeneity (ITH) and tumor evolutionary history. Single-cell DNA sequencing (scDNA-seq) makes it convenient to profile CNAs at single-cell resolution, and thus aids in better characterization of ITH. Despite that several computational methods have been proposed to decipher single-cell CNAs, their performance is limited in either breakpoint detection or copy number estimation due to the high dimensionality and noisy nature of read counts data.

RESULTS

By treating breakpoint detection as a process to segment high dimensional read count sequence, we develop a novel method called DeepCNA for cross-cell segmentation of read count sequence and per-cell inference of CNAs. To cope with the difficulty of segmentation, an autoencoder (AE) network is employed in DeepCNA to project the original data into a low-dimensional space, where the breakpoints can be efficiently detected along each latent dimension and further merged to obtain the final breakpoints. Unlike the existing methods that manually calculate certain statistics of read counts to find breakpoints, the AE model makes it convenient to automatically learn the representations. Based on the inferred breakpoints, we employ a mixture model to predict copy numbers of segments for each cell, and leverage expectation-maximization algorithm to efficiently estimate cell ploidy by exploring the most abundant copy number state. Benchmarking results on simulated and real data demonstrate our method is able to accurately infer breakpoints as well as absolute copy numbers and surpasses the existing methods under different test conditions. DeepCNA can be accessed at: https://github.com/zhyu-lab/deepcna .

CONCLUSIONS

Profiling single-cell CNAs based on deep learning is becoming a new paradigm of scDNA-seq data analysis, and DeepCNA is an enhancement to the current arsenal of computational methods for investigating cancer genomics.

ConDoR: tumor phylogeny inference with a copy-number constrained mutation loss model

A tumor contains a diverse collection of somatic mutations that reflect its past evolutionary history and that range in scale from single nucleotide variants (SNVs) to large-scale copy-number aberrations (CNAs). However, no current single-cell DNA sequencing (scDNA-seq) technology produces accurate measurements of both SNVs and CNAs, complicating the inference of tumor phylogenies. We introduce a new evolutionary model, the constrained k-Dollo model, that uses SNVs as phylogenetic markers but constrains losses of SNVs according to clusters of cells. We derive an algorithm, ConDoR, that infers phylogenies from targeted scDNA-seq data using this model. We demonstrate the advantages of ConDoR on simulated and real scDNA-seq data.

Novel scheme for defining the clinical implications of TP53 mutations in myeloid neoplasia

BACKGROUND

TP53 mutations (TP53MT) occur in diverse genomic configurations. Particularly, biallelic inactivation is associated with poor overall survival in cancer. Lesions affecting only one allele might not be directly leukemogenic, questioning the presence of cryptic biallelic subclones in cases with dismal prognosis.

METHODS

We have collected clinical and molecular data of 7400 patients with myeloid neoplasms and applied a novel model by identifying an optimal VAF cutoff using a statistically robust strategy of sampling-based regression on survival data to accurately classify the TP53 allelic configuration and assess prognosis more precisely.

RESULTS

Overall, TP53MT were found in 1010 patients. Following the traditional criteria, 36% of the cases were classified as single hits, while 64% exhibited double hits genomic configuration. Using a newly developed molecular algorithm, we found that 579 (57%) patients had unequivocally biallelic, 239 (24%) likely contained biallelic, and 192 (19%) had most likely monoallelic TP53MT. Interestingly, our method was able to upstage 192 out of 352 (54.5%) traditionally single hit lesions into a probable biallelic category. Such classification was further substantiated by a survival-based model built after re-categorization. Among cases traditionally considered monoallelic, the overall survival of those with probable monoallelic mutations was similar to the one of wild-type patients and was better than that of patients with a biallelic configuration. As a result, patients with certain biallelic hits, regardless of the disease subtype (AML or MDS), had a similar prognosis. Similar results were observed when the model was applied to an external cohort. In addition, single-cell DNA studies unveiled the biallelic nature of previously considered monoallelic cases.

CONCLUSION

Our novel approach more accurately resolves TP53 genomic configuration and uncovers genetic mosaicism for the use in the clinical setting to improve prognostic evaluation.

SIEVE: joint inference of single-nucleotide variants and cell phylogeny from single-cell DNA sequencing data

We present SIEVE, a statistical method for the joint inference of somatic variants and cell phylogeny under the finite-sites assumption from single-cell DNA sequencing. SIEVE leverages raw read counts for all nucleotides and corrects the acquisition bias of branch lengths. In our simulations, SIEVE outperforms other methods in phylogenetic reconstruction and variant calling accuracy, especially in the inference of homozygous variants. Applying SIEVE to three datasets, one for triple-negative breast (TNBC), and two for colorectal cancer (CRC), we find that double mutant genotypes are rare in CRC but unexpectedly frequent in the TNBC samples.

Parsimonious Clone Tree Integration in cancer

BACKGROUND

Every tumor is composed of heterogeneous clones, each corresponding to a distinct subpopulation of cells that accumulated different types of somatic mutations, ranging from single-nucleotide variants (SNVs) to copy-number aberrations (CNAs). As the analysis of this intra-tumor heterogeneity has important clinical applications, several computational methods have been introduced to identify clones from DNA sequencing data. However, due to technological and methodological limitations, current analyses are restricted to identifying tumor clones only based on either SNVs or CNAs, preventing a comprehensive characterization of a tumor's clonal composition.

RESULTS

To overcome these challenges, we formulate the identification of clones in terms of both SNVs and CNAs as a integration problem while accounting for uncertainty in the input SNV and CNA proportions. We thus characterize the computational complexity of this problem and we introduce PACTION (PArsimonious Clone Tree integratION), an algorithm that solves the problem using a mixed integer linear programming formulation. On simulated data, we show that tumor clones can be identified reliably, especially when further taking into account the ancestral relationships that can be inferred from the input SNVs and CNAs. On 49 tumor samples from 10 prostate cancer patients, our integration approach provides a higher resolution view of tumor evolution than previous studies.

CONCLUSION

PACTION is an accurate and fast method that reconstructs clonal architecture of cancer tumors by integrating SNV and CNA clones inferred using existing methods.

Copy Number Variation Detection by Single-Cell DNA Sequencing with SCOPE

Whole-genome single-cell DNA sequencing (scDNA-seq) enables the characterization of copy number profiles at the cellular level. This circumvents the averaging effects associated with bulk-tissue sequencing and has increased resolution yet decreased ambiguity in deconvolving cancer subclones and elucidating cancer evolutionary history. ScDNA-seq data is, however, sparse, noisy, and highly variable even within a homogeneous cell population, due to the biases and artifacts that are introduced during the library preparation and sequencing procedure. Here, we describe SCOPE, a normalization and copy number estimation method for scDNA-seq data. We give an overview of the methodology and illustrate SCOPE with step-by-step demonstrations.

Somatic Single-Nucleotide Variant Calling from Single-Cell DNA Sequencing Data Using SCAN-SNV

SCAN-SNV is a recent computational tool for somatic single-nucleotide variant (SNV) identification from the single-cell DNA sequencing data. The workflow of the SCAN-SNV package is as follows. First, candidate somatic SNVs and credible heterozygous single-nucleotide polymorphisms (hSNP) are obtained by analyzing single-cell and matched bulk sequencing data, respectively. Subsequently, SCAN-SNV estimates genome-wide allele-specific amplification balance (AB) at any position of DNA sequencing data using a probabilistic spatial statistical model. Finally, candidate somatic SNVs that are likely artifacts according to the AB predictions are further removed to obtain putative mutations. This chapter provides a step-by-step practical guide of the package by explaining how to install and use the variance caller in a real-world example.

Distinguishing linear and branched evolution given single-cell DNA sequencing data of tumors

Background

Cancer arises from an evolutionary process where somatic mutations give rise to clonal expansions. Reconstructing this evolutionary process is useful for treatment decision-making as well as understanding evolutionary patterns across patients and cancer types. In particular, classifying a tumor’s evolutionary process as either linear or branched and understanding what cancer types and which patients have each of these trajectories could provide useful insights for both clinicians and researchers. While comprehensive cancer phylogeny inference from single-cell DNA sequencing data is challenging due to limitations with current sequencing technology and the complexity of the resulting problem, current data might provide sufficient signal to accurately classify a tumor’s evolutionary history as either linear or branched.

Results

We introduce the Linear Perfect Phylogeny Flipping (LPPF) problem as a means of testing two alternative hypotheses for the pattern of evolution, which we prove to be NP-hard. We develop Phyolin, which uses constraint programming to solve the LPPF problem. Through both in silico experiments and real data application, we demonstrate the performance of our method, outperforming a competing machine learning approach.

Conclusion

Phyolin is an accurate, easy to use and fast method for classifying an evolutionary trajectory as linear or branched given a tumor’s single-cell DNA sequencing data.

Dysregulated KRAS gene-signaling axis and abnormal chromatin remodeling drive therapeutic resistance in heterogeneous-sized circulating tumor cells in gastric cancer patients

The mechanism by which heterogeneous-sized circulating tumor cells (CTCs) in gastric cancer (GC) patients are resistant to the targeted therapy and/or chemotherapy remains unclear. This study investigated prognostic value and genomic variations of size-heterogenous CTCs, in an attempt to unravel the molecular mechanisms underlying the therapeutic resistance, which is relevant to poor prognosis in GC. Aneuploid CTCs, detected in 111 advanced GC patients, were categorized into small (≤white blood cell [WBC], 25.54%) and large (>WBC, 74.46%) cells. Pre-treatment patients possessing ≥3 baseline small CTCs with trisomy 8 (_SCTCs^tri) or ≥6 large multiploid CTCs (_LCTCs^multi) showed an inferior median progression-free survival. Moreover, the cut-off value of ≥6 _LCTCs^multi was also an effective prognosticator for poor median overall survival. Single cell-based DNA sequencing of 50 targeted CTCs indicated that _SCTCs^tri and _LCTCs^multi harbored distinct gene variations respectively. Mutations in the KRAS and Rap1 pathway were remarkably abundant in _SCTCs^tri, whereas several unique mutations in the MET/PI3K/AKT pathway and SMARCB1 gene were identified in _LCTCs^multi. Obtained results suggested that _SCTCs^tri and _LCTCs^multi exhibited different mechanisms to therapy resistance and correlated with patients' poor outcome.

Comprehensive description of the current breast cancer microenvironment advancements via single-cell analysis

Breast cancer is a heterogeneous disease with a complex microenvironment consisting of tumor cells, immune cells, fibroblasts and vascular cells. These cancer-associated cells shape the tumor microenvironment (TME) and influence the progression of breast cancer and the therapeutic responses in patients. The exact composition of the intra-tumoral cells is mixed as the highly heterogeneous and dynamic nature of the TME. Recent advances in single-cell technologies such as single-cell DNA sequencing (scDNA-seq), single-cell RNA sequencing (scRNA-seq) and mass cytometry have provided new insights into the phenotypic and functional diversity of tumor-infiltrating cells in breast cancer. In this review, we have outlined the recent progress in single-cell characterization of breast tumor ecosystems, and summarized the phenotypic diversity of intra-tumoral cells and their potential prognostic relevance.

Pancreatic circulating tumor cell detection by targeted single-cell next-generation sequencing

BACKGROUND AND AIMS

Single-cell next-generation sequencing (scNGS) technology has been widely used in genomic profiling, which relies on whole-genome amplification (WGA). However, WGA introduces errors and is especially less accurate when applied to single nucleotide variant (SNV) analysis. Targeted scNGS for SNV without WGA has not been described. We aimed to develop a method to detect circulating tumor cells (CTCs) with DNA SNVs.

METHODS

We tested this targeted scNGS method with three driver mutant genes (KRAS/TP53/SMAD4) on one pancreatic cancer cell line AsPC-1 and then applied it to patients with metastatic PDAC for the validation.

RESULTS

All single-cell of AsPC-1 and spiked-in AsPC-1 cells in healthy donor blood, which were isolated by the filtration with size or by flow cytometry, were detected by targeted scNGS method. All blood samples from six patients with metastatic PDAC, for the validation of target scNGS method, showed CTCs with SNVs of KRAS/TP53/SMAD4 and the positive confirmation of immunofluorescent stainings with Pan-CK/Vimentin/CD45. Four patients with early stage disease, one patient with benign pancreatic cyst and a healthy control sample all showed concordant results between targeted scNGS and CTC enumeration.

CONCLUSIONS

The novel technique of targeted scNGS for SNV analysis, without pre-amplification, is a promising method for identifying and characterizing circulating tumor cells.

Methods for copy number aberration detection from single-cell DNA-sequencing data

Copy number aberrations (CNAs), which are pathogenic copy number variations (CNVs), play an important role in the initiation and progression of cancer. Single-cell DNA-sequencing (scDNAseq) technologies produce data that is ideal for inferring CNAs. In this review, we review eight methods that have been developed for detecting CNAs in scDNAseq data, and categorize them according to the steps of a seven-step pipeline that they employ. Furthermore, we review models and methods for evolutionary analyses of CNAs from scDNAseq data and highlight advances and future research directions for computational methods for CNA detection from scDNAseq data.

SCSIM: Jointly simulating correlated single-cell and bulk next-generation DNA sequencing data

Background

Recently, it has become possible to collect next-generation DNA sequencing data sets that are composed of multiple samples from multiple biological units where each of these samples may be from a single cell or bulk tissue. Yet, there does not yet exist a tool for simulating DNA sequencing data from such a nested sampling arrangement with single-cell and bulk samples so that developers of analysis methods can assess accuracy and precision.

Results

We have developed a tool that simulates DNA sequencing data from hierarchically grouped (correlated) samples where each sample is designated bulk or single-cell. Our tool uses a simple configuration file to define the experimental arrangement and can be integrated into software pipelines for testing of variant callers or other genomic tools.

Conclusions

The DNA sequencing data generated by our simulator is representative of real data and integrates seamlessly with standard downstream analysis tools.

Conbase: a software for unsupervised discovery of clonal somatic mutations in single cells through read phasing

Accurate variant calling and genotyping represent major limiting factors for downstream applications of single-cell genomics. Here, we report Conbase for the identification of somatic mutations in single-cell DNA sequencing data. Conbase leverages phased read data from multiple samples in a dataset to achieve increased confidence in somatic variant calls and genotype predictions. Comparing the performance of Conbase to three other methods, we find that Conbase performs best in terms of false discovery rate and specificity and provides superior robustness on simulated data, in vitro expanded fibroblasts and clonal lymphocyte populations isolated directly from a healthy human donor.