Search

Motivation Single-cell RNA sequencing (scRNA-seq) data are widely used to study cancer cell states and their heterogeneity. However, the tumour microenvironment is usually a mixture of healthy and cancerous cells and it can be difficult to fully separate these two populations based on transcriptomics alone. If available, somatic single-nucleotide variants (SNVs) observed in the scRNA-seq data could be used to identify the cancer population and match that information with the single cells' expression profile. However, calling somatic SNVs in scRNA-seq data is a challenging task, as most variants seen in the short-read data are not somatic, but can instead be germline variants, RNA edits or transcription, sequencing, or processing errors. In addition, only variants present in actively transcribed regions for each individual cell will be seen in the data. Results To address these challenges, we develop CCLONE (Cancer Cell Labelling On Noisy Expression), an interpretable tool adapted to handle the uncertainty and sparsity of SNVs called from scRNA-seq data. CCLONE jointly identifies cancer clonal populations, and their associated variants. We apply CCLONE on two acute myeloid leukaemia datasets and one lung adenocarcinoma dataset and show that CCLONE captures both genetic clones and somatic events for multiple patients. These results show how CCLONE can be used to gather insight into the course of the disease and the origin of cancer cells in scRNA-seq data. Availability and implementation Source code is available at github.com/HaghverdiLab/CCLONE. [ABSTRACT FROM AUTHOR]

A few years ago, it was proposed to use the simultaneous quantification of unspliced and spliced messenger RNA (mRNA) to add a temporal dimension to high-throughput snapshots of single cell RNA sequencing data. This concept can yield additional insight into the transcriptional dynamics of the biological systems under study. However, current methods for inferring cell state velocities from such data (known as RNA velocities) are afflicted by several theoretical and computational problems, hindering realistic and reliable velocity estimation. We discuss these issues and propose new solutions for addressing some of the current challenges in consistency of data processing, velocity inference and visualisation. We translate our computational conclusion in two velocity analysis tools: one detailed method κ-velo and one heuristic method eco-velo, each of which uses a different set of assumptions about the data. Author summary: Single cell transcriptomics has been used to study dynamical biological processes such as cell differentiation or disease progression. An ideal study of these systems would track individual cells in time but this is not directly feasible since cells are destroyed as part of the sequencing protocol. Because of asynchronous progression of cells, single cell snapshot datasets often capture cells at different stages of progression. The challenge is to infer both the overall direction of progression (pseudotime) as well as single cell specific variations in the progression. Computational methods development for inference of the overall direction are well advanced but attempts to address the single cell level variations of the dynamics are newer. Simultaneous measurement of abundances of new (unspliced) and older (spliced) mRNA in the same single cell adds a temporal dimension to the data which can be used to infer the time derivative of single cells progression through the dynamical process. State-of-the-art methods for inference of cell state velocities from RNA-seq data (also known as RNA velocity) have multiple unaddressed issues. In this manuscript, we discuss these issues and propose new solutions. In previous works, agreement of RNA velocity estimations with pseudotime has been used as validation. We show that this in itself is not proof that the method works reliably and the overall direction of progression has to be distinguished from individual cells' behaviour. We propose two new methods (one detailed and one cost efficient heuristic) for estimation and visualisation of RNA velocities and show that our methods faithfully capture the single-cell variances and overall trend on simulation. We further apply the methods to different datasets and show how the method can help us gain biological insight from real data. [ABSTRACT FROM AUTHOR]