An integrated single-cell RNA-seq map of human neuroblastoma tumors and preclinical models uncovers divergent mesenchymal-like gene expression programs.

Neuroblastoma is a common pediatric cancer, where preclinical studies suggest that a mesenchymal-like gene expression program contributes to chemotherapy resistance. However, clinical outcomes remain poor, implying we need a better understanding of the relationship between patient tumor heterogeneity and preclinical models. Here, we generated single-cell RNA-seq maps of neuroblastoma cell lines, patient-derived xenograft models (PDX), and a genetically engineered mouse model (GEMM). We developed an unsupervised machine learning approach ('automatic consensus nonnegative matrix factorization' (acNMF)) to compare the gene expression programs found in preclinical models to a large cohort of patient tumors. We confirmed a weakly expressed, mesenchymal-like program in otherwise adrenergic cancer cells in some pre-treated high-risk patient tumors, but this appears distinct from the presumptive drug-resistance mesenchymal programs evident in cell lines. Surprisingly however, this weak-mesenchymal-like program was maintained in PDX and could be chemotherapy-induced in our GEMM after only 24 hours, suggesting an uncharacterized therapy-escape mechanism. Collectively, our findings improve the understanding of how neuroblastoma patient tumor heterogeneity is reflected in preclinical models, provides a comprehensive integrated resource, and a generalizable set of computational methodologies for the joint analysis of clinical and pre-clinical single-cell RNA-seq datasets.

Conditional c-MYC activation in catecholaminergic cells drives distinct neuroendocrine tumors: neuroblastoma vs somatostatinoma.

The MYC proto-oncogenes (c-MYC, MYCN , MYCL ) are among the most deregulated oncogenic drivers in human malignancies including high-risk neuroblastoma, 50% of which are MYCN -amplified. Genetically engineered mouse models (GEMMs) based on the MYCN transgene have greatly expanded the understanding of neuroblastoma biology and are powerful tools for testing new therapies. However, a lack of c-MYC-driven GEMMs has hampered the ability to better understand mechanisms of neuroblastoma oncogenesis and therapy development given that c-MYC is also an important driver of many high-risk neuroblastomas. In this study, we report two transgenic murine neuroendocrine models driven by conditional c-MYC induction in tyrosine hydroxylase (Th) and dopamine ß-hydroxylase (Dbh)-expressing cells. c-MYC induction in Th-expressing cells leads to a preponderance of Pdx1 + somatostatinomas, a type of pancreatic neuroendocrine tumor (PNET), resembling human somatostatinoma with highly expressed gene signatures of δ cells and potassium channels. In contrast, c-MYC induction in Dbh-expressing cells leads to onset of neuroblastomas, showing a better transforming capacity than MYCN in a comparable C57BL/6 genetic background. The c-MYC murine neuroblastoma tumors recapitulate the pathologic and genetic features of human neuroblastoma, express GD2, and respond to anti-GD2 immunotherapy. This model also responds to DFMO, an FDA-approved inhibitor targeting ODC1, which is a known MYC transcriptional target. Thus, establishing c-MYC-overexpressing GEMMs resulted in different but related tumor types depending on the targeted cell and provide useful tools for testing immunotherapies and targeted therapies for these diseases.

Genome-wide mapping of cancer dependency genes and genetic modifiers of chemotherapy in high-risk hepatoblastoma.

A lack of relevant genetic models and cell lines hampers our understanding of hepatoblastoma pathogenesis and the development of new therapies for this neoplasm. Here, we report an improved MYC-driven hepatoblastoma-like murine model that recapitulates the pathological features of embryonal type of hepatoblastoma, with transcriptomics resembling the high-risk gene signatures of the human disease. Single-cell RNA-sequencing and spatial transcriptomics identify distinct subpopulations of hepatoblastoma cells. After deriving cell lines from the mouse model, we map cancer dependency genes using CRISPR-Cas9 screening and identify druggable targets shared with human hepatoblastoma (e.g., CDK7, CDK9, PRMT1, PRMT5). Our screen also reveals oncogenes and tumor suppressor genes in hepatoblastoma that engage multiple, druggable cancer signaling pathways. Chemotherapy is critical for human hepatoblastoma treatment. A genetic mapping of doxorubicin response by CRISPR-Cas9 screening identifies modifiers whose loss-of-function synergizes with (e.g., PRKDC) or antagonizes (e.g., apoptosis genes) the effect of chemotherapy. The combination of PRKDC inhibition and doxorubicin-based chemotherapy greatly enhances therapeutic efficacy. These studies provide a set of resources including disease models suitable for identifying and validating potential therapeutic targets in human high-risk hepatoblastoma.

Dynamics of Age- versus Therapy-Related Clonal Hematopoiesis in Long-term Survivors of Pediatric Cancer.

We present the first comprehensive investigation of clonal hematopoiesis (CH) in 2,860 long-term survivors of pediatric cancer with a median follow-up time of 23.5 years. Deep sequencing over 39 CH-related genes reveals mutations in 15% of the survivors, significantly higher than the 8.5% in 324 community controls. CH in survivors is associated with exposures to alkylating agents, radiation, and bleomycin. Therapy-related CH shows significant enrichment in STAT3, characterized as a CH gene specific to survivors of Hodgkin lymphoma, and TP53. Single-cell profiling of peripheral blood samples revealed STAT3 mutations predominantly present in T cells and contributed by SBS25, a mutational signature associated with procarbazine exposure. Serial sample tracking reveals that larger clone size is a predictor for future expansion of age-related CH clones, whereas therapy-related CH remains stable decades after treatment. These data depict the distinct dynamics of these CH subtypes and support the need for longitudinal monitoring to determine the potential contribution to late effects. SIGNIFICANCE: This first comprehensive CH analysis in long-term survivors of pediatric cancer presents the elevated prevalence and therapy exposures/diagnostic spectrum associated with CH. Due to the contrasting dynamics of clonal expansion for age-related versus therapy-related CH, longitudinal monitoring is recommended to ascertain the long-term effects of therapy-induced CH in pediatric cancer survivors. See related commentary by Collord and Behjati, p. 811. This article is highlighted in the In This Issue feature, p. 799.

Cell type identification in spatial transcriptomics data can be improved by leveraging cell-type-informative paired tissue images using a Bayesian probabilistic model.

Spatial transcriptomics technologies have recently emerged as a powerful tool for measuring spatially resolved gene expression directly in tissues sections, revealing cell types and their dysfunction in unprecedented detail. However, spatial transcriptomics technologies are limited in their ability to separate transcriptionally similar cell types and can suffer further difficulties identifying cell types in slide regions where transcript capture is low. Here, we describe a conceptually novel methodology that can computationally integrate spatial transcriptomics data with cell-type-informative paired tissue images, obtained from, for example, the reverse side of the same tissue section, to improve inferences of tissue cell type composition in spatial transcriptomics data. The underlying statistical approach is generalizable to any spatial transcriptomics protocol where informative paired tissue images can be obtained. We demonstrate a use case leveraging cell-type-specific immunofluorescence markers obtained on mouse brain tissue sections and a use case for leveraging the output of AI annotated H&E tissue images, which we used to markedly improve the identification of clinically relevant immune cell infiltration in breast cancer tissue. Thus, combining spatial transcriptomics data with paired tissue images has the potential to improve the identification of cell types and hence to improve the applications of spatial transcriptomics that rely on accurate cell type identification.

Simultaneous monitoring of disease and microbe dynamics through plasma DNA sequencing in pediatric patients with acute lymphoblastic leukemia.

Treatment of acute lymphoblastic leukemia (ALL) necessitates continuous risk assessment of leukemic disease burden and infections that arise in the setting of immunosuppression. This study was performed to assess the feasibility of a hybrid capture next-generation sequencing panel to longitudinally measure molecular leukemic disease clearance and microbial species abundance in 20 pediatric patients with ALL throughout induction chemotherapy. This proof of concept helps establish a technical and conceptual framework that we anticipate will be expanded and applied to additional patients with leukemia, as well as extended to additional cancer types. Molecular monitoring can help accelerate the attainment of insights into the temporal biology of host-microbe-leukemia interactions, including how those changes correlate with and alter anticancer therapy efficacy. We also anticipate that fewer invasive bone marrow examinations will be required, as these methods improve with standardization and are validated for clinical use.

The chemotherapeutic CX-5461 primarily targets TOP2B and exhibits selective activity in high-risk neuroblastoma.

Survival in high-risk pediatric neuroblastoma has remained around 50% for the last 20 years, with immunotherapies and targeted therapies having had minimal impact. Here, we identify the small molecule CX-5461 as selectively cytotoxic to high-risk neuroblastoma and synergistic with low picomolar concentrations of topoisomerase I inhibitors in improving survival in vivo in orthotopic patient-derived xenograft neuroblastoma mouse models. CX-5461 recently progressed through phase I clinical trial as a first-in-human inhibitor of RNA-POL I. However, we also use a comprehensive panel of in vitro and in vivo assays to demonstrate that CX-5461 has been mischaracterized and that its primary target at pharmacologically relevant concentrations, is in fact topoisomerase II beta (TOP2B), not RNA-POL I. This is important because existing clinically approved chemotherapeutics have well-documented off-target interactions with TOP2B, which have previously been shown to cause both therapy-induced leukemia and cardiotoxicity-often-fatal adverse events, which can emerge several years after treatment. Thus, while we show that combination therapies involving CX-5461 have promising anti-tumor activity in vivo in neuroblastoma, our identification of TOP2B as the primary target of CX-5461 indicates unexpected safety concerns that should be examined in ongoing phase II clinical trials in adult patients before pursuing clinical studies in children.

Accurate genomic variant detection in single cells with primary template-directed amplification.

Improvements in whole genome amplification (WGA) would enable new types of basic and applied biomedical research, including studies of intratissue genetic diversity that require more accurate single-cell genotyping. Here, we present primary template-directed amplification (PTA), an isothermal WGA method that reproducibly captures >95% of the genomes of single cells in a more uniform and accurate manner than existing approaches, resulting in significantly improved variant calling sensitivity and precision. To illustrate the types of studies that are enabled by PTA, we developed direct measurement of environmental mutagenicity (DMEM), a tool for mapping genome-wide interactions of mutagens with single living human cells at base-pair resolution. In addition, we utilized PTA for genome-wide off-target indel and structural variant detection in cells that had undergone CRISPR-mediated genome editing, establishing the feasibility for performing single-cell evaluations of biopsies from edited tissues. The improved precision and accuracy of variant detection with PTA overcomes the current limitations of accurate WGA, which is the major obstacle to studying genetic diversity and evolution at cellular resolution.

Integrative genomic analyses reveal mechanisms of glucocorticoid resistance in acute lymphoblastic leukemia.

Identification of genomic and epigenomic determinants of drug resistance provides important insights for improving cancer treatment. Using agnostic genome-wide interrogation of mRNA and miRNA expression, DNA methylation, SNPs, CNAs and SNVs/Indels in primary human acute lymphoblastic leukemia cells, we identified 463 genomic features associated with glucocorticoid resistance. Gene-level aggregation identified 118 overlapping genes, 15 of which were confirmed by genome-wide CRISPR screen. Collectively, this identified 30 of 38 (79%) known glucocorticoid-resistance genes/miRNAs and all 38 known resistance pathways, while revealing 14 genes not previously associated with glucocorticoid-resistance. Single cell RNAseq and network-based transcriptomic modelling corroborated the top previously undiscovered gene, CELSR2. Manipulation of CELSR2 recapitulated glucocorticoid resistance in human leukemia cell lines and revealed a synergistic drug combination (prednisolone and venetoclax) that mitigated resistance in mouse xenograft models. These findings illustrate the power of an integrative genomic strategy for elucidating genes and pathways conferring drug resistance in cancer cells.

Evaluation of Plasma Microbial Cell-Free DNA Sequencing to Predict Bloodstream Infection in Pediatric Patients With Relapsed or Refractory Cancer.

Importance: Bloodstream infection (BSI) is a common, life-threatening complication of treatment for cancer. Predicting BSI before onset of clinical symptoms would enable preemptive therapy, but there is no reliable screening test. Objective: To estimate sensitivity and specificity of plasma microbial cell-free DNA sequencing (mcfDNA-seq) for predicting BSI in patients at high risk of life-threatening infection. Design, Setting, and Participants: A prospective pilot cohort study of mcfDNA-seq for predicting BSI in pediatric patients (<25 years of age) with relapsed or refractory cancers at St Jude Children's Research Hospital, a specialist quaternary pediatric hematology-oncology referral center. Remnant clinical blood samples were collected during chemotherapy and hematopoietic cell transplantation. Samples collected during the 7 days before and at onset of BSI episodes, along with negative control samples from study participants, underwent blinded testing using a mcfDNA-seq test in a Clinical Laboratory Improvement Amendments/College of American Pathologists-approved laboratory. Main Outcomes and Measures: The primary outcomes were sensitivity of mcfDNA-seq for detecting a BSI pathogen during the 3 days before BSI onset and specificity of mcfDNA-seq in the absence of fever or infection in the preceding or subsequent 7 days. Results: Between August 9, 2017, and June 4, 2018, 47 participants (27 [57%] male; median age [IQR], 10 [5-14] years) were enrolled; 19 BSI episodes occurred in 12 participants, and predictive samples were available for 16 episodes, including 15 bacterial BSI episodes. In the 3 days before the onset of infection, predictive sensitivity of mcfDNA-seq was 75% for all BSIs (12 of 16; 95% CI, 51%-90%) and 80% (12 of 15; 95% CI, 55%-93%) for bacterial BSIs. The specificity of mcfDNA-seq, evaluated on 33 negative control samples from enrolled participants, was 82% (27 of 33; 95% CI, 66%-91%) for any bacterial or fungal organism and 91% (30 of 33; 95% CI, 76%-97%) for any common BSI pathogen, and the concentration of pathogen DNA was lower in control than predictive samples. Conclusions and Relevance: A clinically relevant pathogen can be identified by mcfDNA-seq days before the onset of BSI in a majority of episodes, potentially enabling preemptive treatment. Clinical application appears feasible pending further study. Trial Registration: ClinicalTrials.gov identifier: NCT03226158.

Single-cell RNA sequencing reveals the impact of chromosomal instability on glioblastoma cancer stem cells.

BACKGROUND: Intra-tumor heterogeneity stems from genetic, epigenetic, functional, and environmental differences among tumor cells. A major source of genetic heterogeneity comes from DNA sequence differences and/or whole chromosome and focal copy number variations (CNVs). Whole chromosome CNVs are caused by chromosomal instability (CIN) that is defined by a persistently high rate of chromosome mis-segregation. Accordingly, CIN causes constantly changing karyotypes that result in extensive cell-to-cell genetic heterogeneity. How the genetic heterogeneity caused by CIN influences gene expression in individual cells remains unknown. METHODS: We performed single-cell RNA sequencing on a chromosomally unstable glioblastoma cancer stem cell (CSC) line and a control normal, diploid neural stem cell (NSC) line to investigate the impact of CNV due to CIN on gene expression. From the gene expression data, we computationally inferred large-scale CNVs in single cells. Also, we performed copy number adjusted differential gene expression analysis between NSCs and glioblastoma CSCs to identify copy number dependent and independent differentially expressed genes. RESULTS: Here, we demonstrate that gene expression across large genomic regions scales proportionally to whole chromosome copy number in chromosomally unstable CSCs. Also, we show that the differential expression of most genes between normal NSCs and glioblastoma CSCs is largely accounted for by copy number alterations. However, we identify 269 genes whose differential expression in glioblastoma CSCs relative to normal NSCs is independent of copy number. Moreover, a gene signature derived from the subset of genes that are differential expressed independent of copy number in glioblastoma CSCs correlates with tumor grade and is prognostic for patient survival. CONCLUSIONS: These results demonstrate that CIN is directly responsible for gene expression changes and contributes to both genetic and transcriptional heterogeneity among glioblastoma CSCs. These results also demonstrate that the expression of some genes is buffered against changes in copy number, thus preserving some consistency in gene expression levels from cell-to-cell despite the continuous change in karyotype driven by CIN. Importantly, a gene signature derived from the subset of genes whose expression is buffered against copy number alterations correlates with tumor grade and is prognostic for patient survival that could facilitate patient diagnosis and treatment.

Murine hematopoietic stem cell activity is derived from pre-circulation embryos but not yolk sacs.

The embryonic site of definitive hematopoietic stem cell (dHSC) origination has been debated for decades. Although an intra-embryonic origin is well supported, the yolk sac (YS) contribution to adult hematopoiesis remains controversial. The same developmental origin makes it difficult to identify specific markers that discern between an intraembryonic versus YS-origin using a lineage trace approach. Additionally, the highly migratory nature of blood cells and the inability of pre-circulatory embryonic cells (i.e., 5-7 somite pairs (sp)) to robustly engraft in transplantation, even after culture, has precluded scientists from properly answering these questions. Here we report robust, multi-lineage and serially transplantable dHSC activity from cultured 2-7sp murine embryonic explants (Em-Ex). dHSC are undetectable in 2-7sp YS explants. Additionally, the engraftment from Em-Ex is confined to an emerging CD31+CD45+c-Kit+CD41- population. In sum, our work supports a model in which the embryo, not the YS, is the major source of lifelong definitive hematopoiesis.

LC3-Associated Phagocytosis in Myeloid Cells Promotes Tumor Immune Tolerance.

Targeting autophagy in cancer cells and in the tumor microenvironment are current goals of cancer therapy. However, components of canonical autophagy play roles in other biological processes, adding complexity to this goal. One such alternative function of autophagy proteins is LC3-associated phagocytosis (LAP), which functions in phagosome maturation and subsequent signaling events. Here, we show that impairment of LAP in the myeloid compartment, rather than canonical autophagy, induces control of tumor growth by tumor-associated macrophages (TAM) upon phagocytosis of dying tumor cells. Single-cell RNA sequencing (RNA-seq) analysis revealed that defects in LAP induce pro-inflammatory gene expression and trigger STING-mediated type I interferon responses in TAM. We found that the anti-tumor effects of LAP impairment require tumor-infiltrating T cells, dependent upon STING and the type I interferon response. Therefore, autophagy proteins in the myeloid cells of the tumor microenvironment contribute to immune suppression of T lymphocytes by effecting LAP.

Notch1-induced brain tumor models the sonic hedgehog subgroup of human medulloblastoma.

While activation of the Notch pathway is observed in many human cancers, it is unknown whether elevated Notch1 expression is sufficient to initiate tumorigenesis in most tissues. To test the oncogenic potential of Notch1 in solid tumors, we expressed an activated form of Notch1 (N1ICD) in the developing mouse brain. N1ICD;hGFAP-cre mice were viable but developed severe ataxia and seizures, and died by weaning age. Analysis of transgenic embryo brains revealed that N1ICD expression induced p53-dependent apoptosis. When apoptosis was blocked by genetic deletion of p53, 30% to 40% of N1ICD;GFAP-cre;p53(+/-) and N1ICD;GFAP-cre;p53(-/-) mice developed spontaneous medulloblastomas. Interestingly, N1ICD-induced medulloblastomas most closely resembled the sonic hedgehog subgroup of human medulloblastoma at the molecular level. Surprisingly, N1ICD-induced tumors do not maintain high levels of the Notch pathway gene expression, except for Notch2, showing that initiating oncogenic events may not be decipherable by analyzing growing tumors in some cases. In summary, this study shows that Notch1 has an oncogenic potential in the brain when combined with other oncogenic hits, such as p53 loss, and provides a novel mouse model of medulloblastoma. Cancer Res; 73(17); 5381-90. ©2013 AACR.

The synaptic ribbon is a site of phosphatidic acid generation in ribbon synapses.

Ribbon synapses continuously transmit graded membrane potential changes into changes of synaptic vesicle exocytosis and rely on intense synaptic membrane trafficking. The synaptic ribbon is considered central to this process. In the present study we asked whether tonically active ribbon synapses are associated with the generation of certain lipids, specifically the highly active signaling phospholipid phosphatidic acid (PA). Using PA-sensor proteins, we demonstrate that PA is enriched at mouse retinal ribbon synapses in close vicinity to the synaptic ribbon in situ. As shown by heterologous expression, RIBEYE, a main component of synaptic ribbons, is responsible for PA binding at synaptic ribbons. Furthermore, RIBEYE is directly involved in the synthesis of PA. Using various independent substrate binding and enzyme assays, we demonstrate that the B domain of RIBEYE possesses lysophosphatidic acid (LPA) acyltransferase (LPAAT) activity, which leads to the generation of PA from LPA. Since an LPAAT-deficient RIBEYE mutant does not recruit PA-binding proteins to artificial synaptic ribbons, whereas wild-type RIBEYE supports PA binding, we conclude that the LPAAT activity of the RIBEYE(B) domain is a physiologically relevant source of PA generation at the synaptic ribbon. We propose that PA generated at synaptic ribbons likely facilitates synaptic vesicle trafficking.