Identifying cancer cells from calling single-nucleotide variants in scRNA-seq data.

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.

Late-onset NPM1 mutation in a MYC-amplified relapsed / refractory acute myeloid leukemia patient treated with gemtuzumab ozogamicin and glasdegib.

Not available.

Somatic epimutations enable single-cell lineage tracing in native hematopoiesis across the murine and human lifespan.

Current approaches to lineage tracing of stem cell clones require genetic engineering or rely on sparse somatic DNA variants, which are difficult to capture at single-cell resolution. Here, we show that targeted single-cell measurements of DNA methylation at single-CpG resolution deliver joint information about cellular differentiation state and clonal identities. We develop EPI-clone, a droplet-based method for transgene-free lineage tracing, and apply it to study hematopoiesis, capturing hundreds of clonal trajectories across almost 100,000 single-cells. Using ground-truth genetic barcodes, we demonstrate that EPI-clone accurately identifies clonal lineages throughout hematopoietic differentiation. Applied to unperturbed hematopoiesis, we describe an overall decline of clonal complexity during murine ageing and the expansion of rare low-output stem cell clones. In aged human donors, we identified expanded hematopoietic clones with and without genetic lesions, and various degrees of clonal complexity. Taken together, EPI-clone enables accurate and transgene-free single-cell lineage tracing at scale.

Integrative single-cell expression and functional studies unravels a sensitization to cytarabine-based chemotherapy through HIF pathway inhibition in AML leukemia stem cells.

Relapse remains a major challenge in the clinical management of acute myeloid leukemia (AML) and is driven by rare therapy-resistant leukemia stem cells (LSCs) that reside in specific bone marrow niches. Hypoxia signaling maintains cells in a quiescent and metabolically relaxed state, desensitizing them to chemotherapy. This suggests the hypothesis that hypoxia contributes to the chemoresistance of AML-LSCs and may represent a therapeutic target to sensitize AML-LSCs to chemotherapy. Here, we identify HIFhigh and HIFlow specific AML subgroups (inv(16)/t(8;21) and MLLr, respectively) and provide a comprehensive single-cell expression atlas of 119,000 AML cells and AML-LSCs in paired diagnostic-relapse samples from these molecular subgroups. The HIF/hypoxia pathway signature is attenuated in AML-LSCs compared with more differentiated AML cells but is more expressed than in healthy hematopoietic cells. Importantly, chemical inhibition of HIF cooperates with standard-of-care chemotherapy to impair AML growth and to substantially eliminate AML-LSCs in vitro and in vivo. These findings support the HIF pathway in the stem cell-driven drug resistance of AML and unravel avenues for combinatorial targeted and chemotherapy-based approaches to specifically eliminate AML-LSCs.

A global genetic interaction network by single-cell imaging and machine learning.

Cellular and organismal phenotypes are controlled by complex gene regulatory networks. However, reference maps of gene function are still scarce across different organisms. Here, we generated synthetic genetic interaction and cell morphology profiles of more than 6,800 genes in cultured Drosophila cells. The resulting map of genetic interactions was used for machine learning-based gene function discovery, assigning functions to genes in 47 modules. Furthermore, we devised Cytoclass as a method to dissect genetic interactions for discrete cell states at the single-cell resolution. This approach identified an interaction of Cdk2 and the Cop9 signalosome complex, triggering senescence-associated secretory phenotypes and immunogenic conversion in hemocytic cells. Together, our data constitute a genome-scale resource of functional gene profiles to uncover the mechanisms underlying genetic interactions and their plasticity at the single-cell level.

Clonally resolved single-cell multi-omics identifies routes of cellular differentiation in acute myeloid leukemia.

Inter-patient variability and the similarity of healthy and leukemic stem cells (LSCs) have impeded the characterization of LSCs in acute myeloid leukemia (AML) and their differentiation landscape. Here, we introduce CloneTracer, a novel method that adds clonal resolution to single-cell RNA-seq datasets. Applied to samples from 19 AML patients, CloneTracer revealed routes of leukemic differentiation. Although residual healthy and preleukemic cells dominated the dormant stem cell compartment, active LSCs resembled their healthy counterpart and retained erythroid capacity. By contrast, downstream myeloid progenitors constituted a highly aberrant, disease-defining compartment: their gene expression and differentiation state affected both the chemotherapy response and leukemia's ability to differentiate into transcriptomically normal monocytes. Finally, we demonstrated the potential of CloneTracer to identify surface markers misregulated specifically in leukemic cells. Taken together, CloneTracer reveals a differentiation landscape that mimics its healthy counterpart and may determine biology and therapy response in AML.