SOX10 mediates glioblastoma cell-state plasticity.

Phenotypic plasticity is a cause of glioblastoma therapy failure. We previously showed that suppressing the oligodendrocyte-lineage regulator SOX10 promotes glioblastoma progression. Here, we analyze SOX10-mediated phenotypic plasticity and exploit it for glioblastoma therapy design. We show that low SOX10 expression is linked to neural stem-cell (NSC)-like glioblastoma cell states and is a consequence of temozolomide treatment in animal and cell line models. Single-cell transcriptome profiling of Sox10-KD tumors indicates that Sox10 suppression is sufficient to induce tumor progression to an aggressive NSC/developmental-like phenotype, including a quiescent NSC-like cell population. The quiescent NSC state is induced by temozolomide and Sox10-KD and reduced by Notch pathway inhibition in cell line models. Combination treatment using Notch and HDAC/PI3K inhibitors extends the survival of mice carrying Sox10-KD tumors, validating our experimental therapy approach. In summary, SOX10 suppression mediates glioblastoma progression through NSC/developmental cell-state transition, including the induction of a targetable quiescent NSC state. This work provides a rationale for the design of tumor therapies based on single-cell phenotypic plasticity analysis.

T-bet suppresses proliferation of malignant B cells in chronic lymphocytic leukemia.

ABSTRACT: The T-box transcription factor T-bet is known as a master regulator of the T-cell response but its role in malignant B cells has not been sufficiently explored. Here, we conducted single-cell resolved multi-omics analyses of malignant B cells from patients with chronic lymphocytic leukemia (CLL) and studied a CLL mouse model with a genetic knockout of Tbx21. We found that T-bet acts as a tumor suppressor in malignant B cells by decreasing their proliferation rate. NF-κB activity, induced by inflammatory signals provided by the microenvironment, triggered T-bet expression, which affected promoter-proximal and distal chromatin coaccessibility and controlled a specific gene signature by mainly suppressing transcription. Gene set enrichment analysis identified a positive regulation of interferon signaling and negative control of proliferation by T-bet. In line, we showed that T-bet represses cell cycling and is associated with longer overall survival of patients with CLL. Our study uncovered a novel tumor suppressive role of T-bet in malignant B cells via its regulation of inflammatory processes and cell cycling, which has implications for the stratification and therapy of patients with CLL. Linking T-bet activity to inflammation explains the good prognostic role of genetic alterations in the inflammatory signaling pathways in CLL.

Branched-chain amino acid transaminase 1 regulates glioblastoma cell plasticity and contributes to immunosuppression.

BACKGROUND: Glioblastoma is the most common malignant brain tumor in adults. Cellular plasticity and the poorly differentiated features result in a fast relapse of the tumors following treatment. Moreover, the immunosuppressive microenvironment proved to be a major obstacle to immunotherapeutic approaches. Branched-chain amino acid transaminase 1 (BCAT1) was shown to drive the growth of glioblastoma and other cancers;however, its oncogenic mechanism remains poorly understood. METHODS: Using human tumor data, cell line models and orthotopic immuno-competent and -deficient mouse models, we investigated the phenotypic and mechanistic effects of BCAT1 on glioblastoma cell state and immunomodulation. RESULTS: Here, we show that BCAT1 is crucial for maintaining the poorly differentiated state of glioblastoma cells and that its low expression correlates with a more differentiated glioblastoma phenotype. Furthermore, orthotopic tumor injection into immunocompetent mice demonstrated that the brain microenvironment is sufficient to induce differentiation of Bcat1-KO tumors in vivo. We link the transition to a differentiated cell state to the increased activity of ten-eleven translocation demethylases and the hypomethylation and activation of neuronal differentiation genes. In addition, the knockout of Bcat1 attenuated immunosuppression, allowing for an extensive infiltration of CD8+ cytotoxic T-cells and complete abrogation of tumor growth. Further analysis in immunodeficient mice revealed that both tumor cell differentiation and immunomodulation following BCAT1-KO contribute to the long-term suppression of tumor growth. CONCLUSIONS: Our study unveils BCAT1's pivotal role in promoting glioblastoma growth by inhibiting tumor cell differentiation and sustaining an immunosuppressive milieu. These findings offer a novel therapeutic avenue for targeting glioblastoma through the inhibition of BCAT1.

Time-resolved, integrated analysis of clonally evolving genomes.

Clonal genome evolution is a key feature of asexually reproducing species and human cancer development. While many studies have described the landscapes of clonal genome evolution in cancer, few determine the underlying evolutionary parameters from molecular data, and even fewer integrate theory with data. We derived theoretical results linking mutation rate, time, expansion dynamics, and biological/clinical parameters. Subsequently, we inferred time-resolved estimates of evolutionary parameters from mutation accumulation, mutational signatures and selection. We then applied this framework to predict the time of speciation of the marbled crayfish, an enigmatic, globally invasive parthenogenetic freshwater crayfish. The results predict that speciation occurred between 1986 and 1990, which is consistent with biological records. We also used our framework to analyze whole-genome sequencing datasets from primary and relapsed glioblastoma, an aggressive brain tumor. The results identified evolutionary subgroups and showed that tumor cell survival could be inferred from genomic data that was generated during the resection of the primary tumor. In conclusion, our framework allowed a time-resolved, integrated analysis of key parameters in clonally evolving genomes, and provided novel insights into the evolutionary age of marbled crayfish and the progression of glioblastoma.

scSNPdemux: a sensitive demultiplexing pipeline using single nucleotide polymorphisms for improved pooled single-cell RNA sequencing analysis.

BACKGROUND: Here we present scSNPdemux, a sample demultiplexing pipeline for single-cell RNA sequencing data using natural genetic variations in humans. The pipeline requires alignment files from Cell Ranger (10× Genomics), a population SNP database and genotyped single nucleotide polymorphisms (SNPs) per sample. The tool works on sparse genotyping data in VCF format for sample identification. RESULTS: The pipeline was tested on both single-cell and single-nuclei based RNA sequencing datasets and showed superior demultiplexing performance over the lipid-based CellPlex and Multi-seq sample multiplexing technique which incurs additional single cell library preparation steps. Specifically, our pipeline demonstrated superior sensitivity and specificity in cell-identity assignment over CellPlex, especially on immune cell types with low RNA content. CONCLUSIONS: We designed a streamlined pipeline for single-cell sample demultiplexing, aiming to overcome common problems in multiplexing samples using single cell libraries which might affect data quality and can be costly.

Simultaneous Inhibition of PI3Kgamma and PI3Kdelta Deteriorates T-cell Function With Implications for Chronic Lymphocytic Leukemia.

Chronic lymphocytic leukemia (CLL) is a common and incurable B-cell malignancy. Recent therapeutic approaches that target the B-cell receptor signaling pathway include inhibition of phosphatidylinositol-3-kinase (PI3K). The PI3K isoform delta is constitutively active in CLL, making it an attractive therapeutic target. However, the expression of PI3K isoforms is not exclusive to leukemic cells, as other immune cells in the tumor microenvironment also rely on PI3K activity. Subsequently, therapeutic inhibition of PI3K causes immune-related adverse events (irAEs). Here, we analyzed the impact of the clinically approved PI3Kδ inhibitors idelalisib and umbralisib, the PI3Kγ inhibitor eganelisib, and the dual-γ and -δ inhibitor duvelisib on the functional capacity of T cells. All investigated inhibitors reduced T-cell activation and proliferation in vitro, which is in line with PI3K being a crucial signaling component of the T-cell receptor signaling. Further, dual inhibition of PI3Kγ and PI3Kδ showed strong additive effects suggesting a role also for PI3Kγ in T cells. Extrapolation of this data to a clinical setting could provide an explanation for the observed irAEs in CLL patients undergoing treatment with PI3K inhibitors. Consequently, this highlights the need for a close monitoring of patients treated with PI3K inhibitors, and particularly duvelisib, due to their potentially increased risk of T-cell deficiencies and associated infections.

Sphingosine-1-Phosphate Recruits Macrophages and Microglia and Induces a Pro-Tumorigenic Phenotype That Favors Glioma Progression.

Glioblastoma is the most aggressive brain tumor in adults. Treatment failure is predominantly caused by its high invasiveness and its ability to induce a supportive microenvironment. As part of this, a major role for tumor-associated macrophages/microglia (TAMs) in glioblastoma development was recognized. Phospholipids are important players in various fundamental biological processes, including tumor-stroma crosstalk, and the bioactive lipid sphingosine-1-phosphate (S1P) has been linked to glioblastoma cell proliferation, invasion, and survival. Despite the urgent need for better therapeutic approaches, novel strategies targeting sphingolipids in glioblastoma are still poorly explored. Here, we showed that higher amounts of S1P secreted by glioma cells are responsible for an active recruitment of TAMs, mediated by S1P receptor (S1PR) signaling through the modulation of Rac1/RhoA. This resulted in increased infiltration of TAMs in the tumor, which, in turn, triggered their pro-tumorigenic phenotype through the inhibition of NFkB-mediated inflammation. Gene set enrichment analyses showed that such an anti-inflammatory microenvironment correlated with shorter survival of glioblastoma patients. Inhibition of S1P restored a pro-inflammatory phenotype in TAMs and resulted in increased survival of tumor-bearing mice. Taken together, our results establish a crucial role for S1P in fine-tuning the crosstalk between glioma and infiltrating TAMs, thus pointing to the S1P-S1PR axis as an attractive target for glioma treatment.

Molecular characterization of Richter syndrome identifies de novo diffuse large B-cell lymphomas with poor prognosis.

Richter syndrome (RS) is the transformation of chronic lymphocytic leukemia (CLL) into aggressive lymphoma, most commonly diffuse large B-cell lymphoma (DLBCL). We characterize 58 primary human RS samples by genome-wide DNA methylation and whole-transcriptome profiling. Our comprehensive approach determines RS DNA methylation profile and unravels a CLL epigenetic imprint, allowing CLL-RS clonal relationship assessment without the need of the initial CLL tumor DNA. DNA methylation- and transcriptomic-based classifiers were developed, and testing on landmark DLBCL datasets identifies a poor-prognosis, activated B-cell-like DLBCL subset in 111/1772 samples. The classification robustly identifies phenotypes very similar to RS with a specific genomic profile, accounting for 4.3-8.3% of de novo DLBCLs. In this work, RS multi-omics characterization determines oncogenic mechanisms, establishes a surrogate marker for CLL-RS clonal relationship, and provides a clinically relevant classifier for a subset of primary "RS-type DLBCL" with unfavorable prognosis.

The coming decade in precision oncology: six riddles.

High-throughput methods to investigate tumour omic landscapes have quickly catapulted cancer specialists into the precision oncology era. The singular lesson of precision oncology might be that, for it to be precise, treatment must be personalized, as each cancer's complex molecular and immune landscape differs from patient to patient. Transformative therapies include those that are targeted at the sequelae of molecular abnormalities or at immune mechanisms, and, increasingly, pathways previously thought to be undruggable have become druggable. Critical to applying precision medicine is the concept that the right combination of drugs must be chosen for each patient and used at the right stage of the disease. Multiple puzzles remain that complicate therapy choice, including evidence that deleterious mutations are common in normal tissues and non-malignant conditions. The host's role is also likely to be key in determining treatment response, especially for immunotherapy. Indeed, maximizing the impact of immunotherapy will require omic analyses to match the right immune-targeted drugs to the individualized patient and tumour setting. In this Perspective, we discuss six key riddles that must be solved to optimize the application of precision oncology to otherwise lethal malignancies.

Focal structural variants revealed by whole genome sequencing disrupt the histone demethylase KDM4C in B-cell lymphomas.

Histone methylation-modifiers, such as EZH2 and KMT2D, are recurrently altered in B-cell lymphomas. To comprehensively describe the landscape of alterations affecting genes encoding histone methylation-modifiers in lymphomagenesis we investigated whole genome and transcriptome data of 186 mature B-cell lymphomas sequenced in the ICGC MMML-Seq project. Besides confirming common alterations of KMT2D (47% of cases), EZH2 (17%), SETD1B (5%), PRDM9 (4%), KMT2C (4%), and SETD2 (4%), also identified by prior exome or RNA-sequencing studies, we here found recurrent alterations to KDM4C in chromosome 9p24, encoding a histone demethylase. Focal structural variation was the main mechanism of KDM4C alterations, and was independent from 9p24 amplification. We also identified KDM4C alterations in lymphoma cell lines including a focal homozygous deletion in a classical Hodgkin lymphoma cell line. By integrating RNA-sequencing and genome sequencing data we predict that KDM4C structural variants result in loss-offunction. By functional reconstitution studies in cell lines, we provide evidence that KDM4C can act as a tumor suppressor. Thus, we show that identification of structural variants in whole genome sequencing data adds to the comprehensive description of the mutational landscape of lymphomas and, moreover, establish KDM4C as a putative tumor suppressive gene recurrently altered in subsets of B-cell derived lymphomas.

BCAT1 redox function maintains mitotic fidelity.

The metabolic enzyme branched-chain amino acid transaminase 1 (BCAT1) drives cell proliferation in aggressive cancers such as glioblastoma. Here, we show that BCAT1 localizes to mitotic structures and has a non-metabolic function as a mitotic regulator. Furthermore, BCAT1 is required for chromosome segregation in cancer and induced pluripotent stem cells and tumor growth in human cerebral organoid and mouse syngraft models. Applying gene knockout and rescue strategies, we show that the BCAT1 CXXC redox motif is crucial for controlling cysteine sulfenylation specifically in mitotic cells, promoting Aurora kinase B localization to centromeres, and securing accurate chromosome segregation. These findings offer an explanation for the well-established role of BCAT1 in promoting cancer cell proliferation. In summary, our data establish BCAT1 as a component of the mitotic apparatus that safeguards mitotic fidelity through a moonlighting redox functionality.

The HHIP-AS1 lncRNA promotes tumorigenicity through stabilization of dynein complex 1 in human SHH-driven tumors.

Most lncRNAs display species-specific expression patterns suggesting that animal models of cancer may only incompletely recapitulate the regulatory crosstalk between lncRNAs and oncogenic pathways in humans. Among these pathways, Sonic Hedgehog (SHH) signaling is aberrantly activated in several human cancer entities. We unravel that aberrant expression of the primate-specific lncRNA HedgeHog Interacting Protein-AntiSense 1 (HHIP-AS1) is a hallmark of SHH-driven tumors including medulloblastoma and atypical teratoid/rhabdoid tumors. HHIP-AS1 is actively transcribed from a bidirectional promoter shared with SHH regulator HHIP. Knockdown of HHIP-AS1 induces mitotic spindle deregulation impairing tumorigenicity in vitro and in vivo. Mechanistically, HHIP-AS1 binds directly to the mRNA of cytoplasmic dynein 1 intermediate chain 2 (DYNC1I2) and attenuates its degradation by hsa-miR-425-5p. We uncover that neither HHIP-AS1 nor the corresponding regulatory element in DYNC1I2 are evolutionary conserved in mice. Taken together, we discover an lncRNA-mediated mechanism that enables the pro-mitotic effects of SHH pathway activation in human tumors.

Modeling the B-cell receptor signaling on single cell level reveals a stable network circuit topology between nonmalignant B cells and chronic lymphocytic leukemia cells and between untreated cells and cells treated with kinase inhibitors.

B-cell receptor (BCR) signaling is central for the pathomechanism of chronic lymphocytic leukemia (CLL), and inhibitors of BCR signaling have substantially improved treatment options. To model malignant and nonmalignant BCR signaling, we quantified five components of BCR signaling (ZAP70/SYK, BTK, PLCγ2, AKT, ERK1/2) in single cells from primary human leukemic cells and from nonmalignant tissue. We measured signaling activity in a time-resolved manner after stimulation with BCR crosslinking by anti-IgM and/or anti-CD19 and with or without inhibition of phosphatases with H2 O2 . The phosphorylation of BCR signaling components was increased in malignant cells compared to nonmalignant cells and in IGHV unmutated CLL cells compared to IGHV mutated CLL cells. Intriguingly, inhibition of phosphatases with H2 O2 led to higher phosphorylation levels of BCR components in CLL cells with mutated IGHV compared to unmutated IGHV. We modeled the connectivity of the cascade components by correlating signal intensities across single cells. The network topology remained stable between malignant and nonmalignant cells. To additionally test for the impact of therapeutic compounds on the network topology, we challenged the BCR signaling cascade with inhibitors for BTK (ibrutinib), PI3K (idelalisib), LYN (dasatinib) and SYK (entospletinib). Idelalisib treatment resulted in similar effects in malignant and nonmalignant cells, whereas ibrutinib was mostly active on CLL cells. Idelalisib and ibrutinib had complementary effects on the BCR signaling cascade whose activity was further reduced upon dasatinib and entospletinib treatment. The characterization of the molecular circuitry of leukemic BCR signaling will allow a more refined targeting of this Achilles heel.

A synergistic interaction between HDAC- and PARP inhibitors in childhood tumors with chromothripsis.

Chromothripsis is a form of genomic instability characterized by the occurrence of tens to hundreds of clustered DNA double-strand breaks in a one-off catastrophic event. Rearrangements associated with chromothripsis are detectable in numerous tumor entities and linked with poor prognosis in some of these, such as Sonic Hedgehog medulloblastoma, neuroblastoma and osteosarcoma. Hence, there is a need for therapeutic strategies eliminating tumor cells with chromothripsis. Defects in DNA double-strand break repair, and in particular homologous recombination repair, have been linked with chromothripsis. Targeting DNA repair deficiencies by synthetic lethality approaches, we performed a synergy screen using drug libraries (n = 375 compounds, 15 models) combined with either a PARP inhibitor or cisplatin. This revealed a synergistic interaction between the HDAC inhibitor romidepsin and PARP inhibition. Functional assays, transcriptome analyses and in vivo validation in patient-derived xenograft mouse models confirmed the efficacy of the combinatorial treatment.

Association of mutation signature effectuating processes with mutation hotspots in driver genes and non-coding regions.

Cancer driving mutations are difficult to identify especially in the non-coding part of the genome. Here, we present sigDriver, an algorithm dedicated to call driver mutations. Using 3813 whole-genome sequenced tumors from International Cancer Genome Consortium, The Cancer Genome Atlas Program, and a childhood pan-cancer cohort, we employ mutational signatures based on single-base substitution in the context of tri- and penta-nucleotide motifs for hotspot discovery. Knowledge-based annotations on mutational hotspots reveal enrichment in coding regions and regulatory elements for 6 mutational signatures, including APOBEC and somatic hypermutation signatures. APOBEC activity is associated with 32 hotspots of which 11 are known and 11 are putative regulatory drivers. Somatic single nucleotide variants clusters detected at hypermutation-associated hotspots are distinct from translocation or gene amplifications. Patients carrying APOBEC induced PIK3CA driver mutations show lower occurrence of signature SBS39. In summary, sigDriver uncovers mutational processes associated with known and putative tumor drivers and hotspots particularly in the non-coding regions of the genome.

Clonal evolution in chronic lymphocytic leukemia is scant in relapsed but accelerated in refractory cases after chemo(immune) therapy.

Clonal evolution is involved in the progression of chronic lymphocytic leukemia (CLL). In order to link evolutionary patterns to different disease courses, we performed a long-term longitudinal mutation profiling study of CLL patients. Tracking somatic mutations and their changes in allele frequency over time and assessing the underlying cancer cell fraction revealed highly distinct evolutionary patterns. Surprisingly, in long-term stable disease and in relapse after long-lasting clinical response to treatment, clonal shifts are minor. In contrast, in refractory disease major clonal shifts occur although there is little impact on leukemia cell counts. As this striking pattern in refractory cases is not linked to a strong contribution of known CLL driver genes, the evolution is mostly driven by treatment-induced selection of sub-clones, underlining the need for novel, non-genotoxic treatment regimens.