Cancer stem cells drive disease progression and relapse in many types of cancer. Despite this, a thorough characterization of these cells remains elusive and with it the ability to eradicate cancer at its source. In acute myeloid leukemia (AML), leukemic stem cells (LSCs) underlie mortality but are difficult to isolate due to their low abundance and high similarity to healthy hematopoietic stem cells (HSCs). Here, we demonstrate that LSCs, HSCs, and pre-leukemic stem cells can be identified and molecularly profiled by combining single-cell transcriptomics with lineage tracing using both nuclear and mitochondrial somatic variants. While mutational status discriminates between healthy and cancerous cells, gene expression distinguishes stem cells and progenitor cell populations. Our approach enables the identification of LSC-specific gene expression programs and the characterization of differentiation blocks induced by leukemic mutations. Taken together, we demonstrate the power of single-cell multi-omic approaches in characterizing cancer stem cells.

Outcomes are poor in patients with higher-risk myelodysplastic syndromes (HR MDS) and frontline treatment options are limited. This phase 1b study investigated safety and efficacy of venetoclax, a selective B-cell lymphoma 2 inhibitor, at the recommended phase 2 dose (RP2D; 400 mg for 14 days per 28-day cycle), in combination with azacitidine (75 mg/m2 for 7 days per 28-day cycle) for treatment-naive HR MDS. Safety was the primary outcome, and complete remission (CR) rate was the primary efficacy outcome. Secondary outcomes included rates of modified overall response (mOR), hematologic improvement (HI), overall survival (OS), and time to next treatment (TTNT). As of May 2023, 107 patients received venetoclax and azacitidine combination at the RP2D. Best response of CR or marrow CR was observed in 29.9% and 50.5% (mOR, 80.4%), respectively. Median OS was 26.0 months, with 1- and 2-year survival estimates of 71.2% and 51.3%, respectively. Among 59 patients with baseline red blood cell and/or platelet transfusion–dependence, 24 (40.7%) achieved transfusion independence on study, including 11 (18.6%) in CR. Fifty-one (49.0%) of 104 evaluable patients achieved HI. Median TTNT excluding transplantation was 13.4 months. Adverse events reflected known safety profiles for venetoclax and azacitidine, including constipation (53.3%), nausea (49.5%), neutropenia (48.6%), thrombocytopenia (44.9%), febrile neutropenia (42.1%), and diarrhea (41.1%). Overall, venetoclax plus azacitidine at the RP2D was well tolerated and had favorable outcomes. A phase 3 study (NCT04401748) is ongoing to confirm survival benefit of this combination. This trial was registered at www.clinicaltrials.gov as #NCT02942290.

Somatic structural variants (SVs) are widespread in cancer, but their impact on disease evolution is understudied due to a lack of methods to directly characterize their functional consequences. We present a computational method, scNOVA, which uses Strand-seq to perform haplotype-aware integration of SV discovery and molecular phenotyping in single cells by using nucleosome occupancy to infer gene expression as a readout. Application to leukemias and cell lines identifies local effects of copy-balanced rearrangements on gene deregulation, and consequences of SVs on aberrant signaling pathways in subclones. We discovered distinct SV subclones with dysregulated Wnt signaling in a chronic lymphocytic leukemia patient. We further uncovered the consequences of subclonal chromothripsis in T cell acute lymphoblastic leukemia, which revealed c-Myb activation, enrichment of a primitive cell state and informed successful targeting of the subclone in cell culture, using a Notch inhibitor. By directly linking SVs to their functional effects, scNOVA enables systematic single-cell multiomic studies of structural variation in heterogeneous cell populations.

Acute myeloid leukemia (AML) is a hematological malignancy with an undefined heritable risk. Here we perform a meta-analysis of three genome-wide association studies, with replication in a fourth study, incorporating a total of 4018 AML cases and 10488 controls. We identify a genome-wide significant risk locus for AML at 11q13.2 (rs4930561; P = 2.15 × 10⁻⁸; KMT5B). We also identify a genome-wide significant risk locus for the cytogenetically normal AML sub-group (N = 1287) at 6p21.32 (rs3916765; P = 1.51 × 10⁻¹⁰; HLA). Our results inform on AML etiology and identify putative functional genes operating in histone methylation (KMT5B) and immune function (HLA).

Clinicians in molecular tumor boards (MTB) are confronted with a growing amount of genetic high-throughput sequencing data. Today, at German university hospitals, these data are usually handled in complex spreadsheets from which clinicians have to obtain the necessary information. The aim of this work was to gather a comprehensive list of requirements to be met by cBioPortal to support processes in MTBs according to clinical needs. Therefore, oncology experts at nine German university hospitals were surveyed in two rounds of interviews. To generate an interview guideline a scoping review was conducted. For visual support in the second round, screenshot mockups illustrating the requirements from the first round were created. Requirements that cBioPortal already meets were skipped during the second round. In the end, 24 requirements with sometimes several conceivable options were identified and 54 screenshot mockups were created. Some of the identified requirements have already been suggested to the community by other users or are currently being implemented in cBioPortal. This shows, that the results are in line with the needs expressed by various disciplines. According to our findings, cBioPortal has the potential to significantly improve the processes and analyses of an MTB after the implementation of the identified requirements.

Somatic mutations in genes coding for splicing factors, e.g., SF3B1, U2AF1, SRSF2, and others are found in approximately 50% of patients with myelodysplastic syndromes (MDS). These mutations have been predicted to frequently occur early in the mutational hierarchy of the disease, therefore, making them particularly attractive potential therapeutic targets. Recent studies in cell lines engineered to carry splicing factor mutations have revealed a strong association with elevated levels of DNA:RNA intermediates (R-loops) and a dependency on proper ATR function. However, data confirming this hypothesis in a representative cohort of primary MDS patient samples have so far been missing. Using CD34+ cells isolated from MDS patients with and without splicing factor mutations as well as healthy controls we show that splicing factor mutation- associated R-loops lead to elevated levels of replication stress and ATR pathway activation. Moreover, splicing factor mutated CD34+ cells are more susceptible to pharmacological inhibition of ATR resulting in elevated levels of DNA damage, cell cycle blockade, and cell death. This can be enhanced by combination treatment with the low-dose splicing modulatory compound Pladienolide B. We further confirm the direct association between R-loops and ATR sensitivity and the presence of a splicing factor mutation using lentiviral overexpression of wild-type and mutant SRSF2 P95H in cord blood CD34+ cells. Collectively, our results from n=53 MDS patients identify replication stress and associated ATR signaling to be critical pathophysiological mechanisms in primary MDS CD34+ cells carrying splicing factor mutations, and provide a preclinical rationale for targeting ATR signaling in these patients.