Altered Treg Infiltration after Discoidin Domain Receptor 1 (DDR1) Inhibition and Knockout Promotes Tumor Growth in Lung Adenocarcinoma.

Lung cancer is the leading cause of cancer-related death worldwide. Discoidin domain receptor 1 (DDR1), a tyrosine kinase receptor, has been associated with poor prognosis in patients with non-small cell lung cancer (NSCLC). However, its role in tumorigenesis remains poorly understood. This work aimed to explore the impact of DDR1 expression on immune cell infiltration in lung adenocarcinoma. Pharmacological inhibition and knockout of DDR1 were used in an immunocompetent mouse model of KRAS/p53-driven lung adenocarcinoma (LUAD). Tumor cells were engrafted subcutaneously, after which tumors were harvested for investigation of immune cell composition via flow cytometry. The Cancer Genome Atlas (TCGA) cohort was used to perform gene expression analysis of 509 patients with LUAD. Pharmacological inhibition and knockout of DDR1 increased the tumor burden, with DDR1 knockout tumors showing a decrease in CD8+ cytotoxic T cells and an increase in CD4+ helper T cells and regulatory T cells. TCGA analysis revealed that low-DDR1-expressing tumors showed higher FoxP3 (regulatory T-cell marker) expression than high-DDR1-expressing tumors. Our study showed that under certain conditions, the inhibition of DDR1, a potential therapeutic target in cancer treatment, might have negative effects, such as inducing a pro-tumorigenic tumor microenvironment. As such, further investigations are necessary.

Inflammatory bone marrow signaling in pediatric acute myeloid leukemia distinguishes patients with poor outcomes.

High levels of the inflammatory cytokine IL-6 in the bone marrow are associated with poor outcomes in pediatric acute myeloid leukemia (pAML), but its etiology remains unknown. Using RNA-seq data from pre-treatment bone marrows of 1489 children with pAML, we show that > 20% of patients have concurrent IL-6, IL-1, IFNα/ß, and TNFα signaling activity and poorer outcomes. Targeted sequencing of pre-treatment bone marrow samples from affected patients (n = 181) revealed 5 highly recurrent patterns of somatic mutation. Using differential expression analyses of the most common genomic subtypes (~60% of total), we identify high expression of multiple potential drivers of inflammation-related treatment resistance. Regardless of genomic subtype, we show that JAK1/2 inhibition reduces receptor-mediated inflammatory signaling by leukemic cells in-vitro. The large number of high-risk pAML genomic subtypes presents an obstacle to the development of mutation-specific therapies. Our findings suggest that therapies targeting inflammatory signaling may be effective across multiple genomic subtypes of pAML.

Oncogenic role of a developmentally regulated NTRK2 splice variant.

Temporally regulated alternative splicing choices are vital for proper development, yet the wrong splice choice may be detrimental. Here, we highlight a previously unidentified role for the neurotrophin receptor splice variant TrkB.T1 in neurodevelopment, embryogenesis, transformation, and oncogenesis across multiple tumor types in humans and mice. TrkB.T1 is the predominant NTRK2 isoform across embryonic organogenesis, and forced overexpression of this embryonic pattern causes multiple solid and nonsolid tumors in mice in the context of tumor suppressor loss. TrkB.T1 also emerges as the predominant NTRK isoform expressed in a wide range of adult and pediatric tumors, including those harboring tropomyosin receptor kinase fusions. Affinity purification-mass spectrometry proteomic analysis reveals distinct interactors with known developmental and oncogenic signaling pathways such as Wnt, transforming growth factor-ß, Sonic Hedgehog, and Ras. From alterations in splicing factors to changes in gene expression, the discovery of isoform specific oncogenes with embryonic ancestry has the potential to shape the way we think about developmental systems and oncology.

A B-cell developmental gene regulatory network is activated in infant AML.

Infant Acute Myeloid Leukemia (AML) is a poorly-addressed, heterogeneous malignancy distinguished by surprisingly few mutations per patient but accompanied by myriad age-specific translocations. These characteristics make treatment of infant AML challenging. While infant AML is a relatively rare disease, it has enormous impact on families, and in terms of life-years-lost and life limiting morbidities. To better understand the mechanisms that drive infant AML, we performed integrative analyses of genome-wide mRNA, miRNA, and DNA-methylation data in diagnosis-stage patient samples. Here, we report the activation of an onco-fetal B-cell developmental gene regulatory network in infant AML. AML in infants is genomically distinct from AML in older children/adults in that it has more structural genomic aberrations and fewer mutations. Differential expression analysis of ~1500 pediatric AML samples revealed a large number of infant-specific genes, many of which are associated with B cell development and function. 18 of these genes form a well-studied B-cell gene regulatory network that includes the epigenetic regulators BRD4 and POU2AF1, and their onco-fetal targets LIN28B and IGF2BP3. All four genes are hypo-methylated in infant AML. Moreover, micro-RNA Let7a-2 is expressed in a mutually exclusive manner with its target and regulator LIN28B. These findings suggest infant AML may respond to bromodomain inhibitors and immune therapies targeting CD19, CD20, CD22, and CD79A.

Neural G0: a quiescent-like state found in neuroepithelial-derived cells and glioma.

Single-cell RNA sequencing has emerged as a powerful tool for resolving cellular states associated with normal and maligned developmental processes. Here, we used scRNA-seq to examine the cell cycle states of expanding human neural stem cells (hNSCs). From these data, we constructed a cell cycle classifier that identifies traditional cell cycle phases and a putative quiescent-like state in neuroepithelial-derived cell types during mammalian neurogenesis and in gliomas. The Neural G0 markers are enriched with quiescent NSC genes and other neurodevelopmental markers found in non-dividing neural progenitors. Putative glioblastoma stem-like cells were significantly enriched in the Neural G0 cell population. Neural G0 cell populations and gene expression are significantly associated with less aggressive tumors and extended patient survival for gliomas. Genetic screens to identify modulators of Neural G0 revealed that knockout of genes associated with the Hippo/Yap and p53 pathways diminished Neural G0 in vitro, resulting in faster G1 transit, down-regulation of quiescence-associated markers, and loss of Neural G0 gene expression. Thus, Neural G0 represents a dynamic quiescent-like state found in neuroepithelial-derived cells and gliomas.

A kinase-deficient NTRK2 splice variant predominates in glioma and amplifies several oncogenic signaling pathways.

Independent scientific achievements have led to the discovery of aberrant splicing patterns in oncogenesis, while more recent advances have uncovered novel gene fusions involving neurotrophic tyrosine receptor kinases (NTRKs) in gliomas. The exploration of NTRK splice variants in normal and neoplastic brain provides an intersection of these two rapidly evolving fields. Tropomyosin receptor kinase B (TrkB), encoded NTRK2, is known for critical roles in neuronal survival, differentiation, molecular properties associated with memory, and exhibits intricate splicing patterns and post-translational modifications. Here, we show a role for a truncated NTRK2 splice variant, TrkB.T1, in human glioma. TrkB.T1 enhances PDGF-driven gliomas in vivo, augments PDGF-induced Akt and STAT3 signaling in vitro, while next generation sequencing broadly implicates TrkB.T1 in the PI3K signaling cascades in a ligand-independent fashion. These TrkB.T1 findings highlight the importance of expanding upon whole gene and gene fusion analyses to include splice variants in basic and translational neuro-oncology research.

Variability in estimated gene expression among commonly used RNA-seq pipelines.

RNA-sequencing data is widely used to identify disease biomarkers and therapeutic targets using numerical methods such as clustering, classification, regression, and differential expression analysis. Such approaches rely on the assumption that mRNA abundance estimates from RNA-seq are reliable estimates of true expression levels. Here, using data from five RNA-seq processing pipelines applied to 6,690 human tumor and normal tissues, we show that nearly 88% of protein-coding genes have similar gene expression profiles across all pipelines. However, for >12% of protein-coding genes, current best-in-class RNA-seq processing pipelines differ in their abundance estimates by more than four-fold when applied to exactly the same samples and the same set of RNA-seq reads. Expression fold changes are similarly affected. Many of the impacted genes are widely studied disease-associated genes. We show that impacted genes exhibit diverse patterns of discordance among pipelines, suggesting that many inter-pipeline differences contribute to overall uncertainty in mRNA abundance estimates. A concerted, community-wide effort will be needed to develop gold-standards for estimating the mRNA abundance of the discordant genes reported here. In the meantime, our list of discordantly evaluated genes provides an important resource for robust marker discovery and target selection.

Comprehensive Transcriptome Profiling of Cryptic CBFA2T3-GLIS2 Fusion-Positive AML Defines Novel Therapeutic Options: A COG and TARGET Pediatric AML Study.

PURPOSE: A cryptic inv(16)(p13.3q24.3) encoding the CBFA2T3-GLIS2 fusion is associated with poor outcome in infants with acute megakaryocytic leukemia. We aimed to broaden our understanding of the pathogenesis of this fusion through transcriptome profiling. EXPERIMENTAL DESIGN: Available RNA from children and young adults with de novo acute myeloid leukemia (AML; N = 1,049) underwent transcriptome sequencing (mRNA and miRNA). Transcriptome profiles for those with the CBFA2T3-GLIS2 fusion (N = 24) and without (N = 1,025) were contrasted to define fusion-specific miRNAs, genes, and pathways. Clinical annotations defined distinct fusion-associated disease characteristics and outcomes. RESULTS: The CBFA2T3-GLIS2 fusion was restricted to infants <3 years old (P < 0.001), and the presence of this fusion was highly associated with adverse outcome (P < 0.001) across all morphologic classifications. Further, there was a striking paucity of recurrent cooperating mutations, and transduction of cord blood stem cells with this fusion was sufficient for malignant transformation. CBFA2T3-GLIS2 positive cases displayed marked upregulation of genes with cell membrane/extracellular matrix localization potential, including NCAM1 and GABRE. Additionally, miRNA profiling revealed significant overexpression of mature miR-224 and miR-452, which are intronic miRNAs transcribed from the GABRE locus. Gene-set enrichment identified dysregulated Hippo, TGFß, and hedgehog signaling, as well as NCAM1 (CD56) interaction pathways. Therapeutic targeting of fusion-positive leukemic cells with CD56-directed antibody-drug conjugate caused significant cytotoxicity in leukemic blasts. CONCLUSIONS: The CBFA2T3-GLIS2 fusion defines a highly refractory entity limited to infants that appears to be sufficient for malignant transformation. Transcriptome profiling elucidated several highly targetable genes and pathways, including the identification of CD56, providing a highly plausible target for therapeutic intervention.

Neutrophil content predicts lymphocyte depletion and anti-PD1 treatment failure in NSCLC.

Immune checkpoint inhibitor (ICI) treatment has recently become a first-line therapy for many non-small cell lung cancer (NSCLC) patients. Unfortunately, most NSCLC patients are refractory to ICI monotherapy, and initial attempts to address this issue with secondary therapeutics have proven unsuccessful. To identify entities precluding CD8+ T cell accumulation in this process, we performed unbiased analyses on flow cytometry, gene expression, and multiplexed immunohistochemical data from a NSCLC patient cohort. The results revealed the presence of a myeloid-rich subgroup, which was devoid of CD4+ and CD8+ T cells. Of all myeloid cell types assessed, neutrophils were the most highly associated with the myeloid phenotype. Additionally, the ratio of CD8+ T cells to neutrophils (CD8/PMN) within the tumor mass optimally distinguished between active and myeloid cases. This ratio was also capable of showing the separation of patients responsive to ICI therapy from those with stable or progressive disease in 2 independent cohorts. Tumor-bearing mice treated with a combination of anti-PD1 and SX-682 (CXCR1/2 inhibitor) displayed relocation of lymphocytes from the tumor periphery into a malignant tumor, which was associated with induction of IFN-γ-responsive genes. These results suggest that neutrophil antagonism may represent a viable secondary therapeutic strategy to enhance ICI treatment outcomes.

Genomic distinctions between metastatic lower and upper tract urothelial carcinoma revealed through rapid autopsy.

BACKGROUND: Little is known about the genomic differences between metastatic urothelial carcinoma (LTUC) and upper tract urothelial carcinoma (UTUC). We compare genomic features of primary and metastatic UTUC and LTUC tumors in a cohort of patients with end stage disease. METHODS: We performed whole exome sequencing on matched primary and metastatic tumor samples (N=37) from 7 patients with metastatic UC collected via rapid autopsy. Inter- and intra-patient mutational burden, mutational signatures, predicted deleterious mutations, and somatic copy alterations (sCNV) were analyzed. RESULTS: We investigated 3 patients with UTUC (3 primary samples, 13 metastases) and 4 patients with LTUC (4 primary samples, 17 metastases). We found that sSNV burden was higher in metastatic LTUC compared to UTUC. Moreover, the APOBEC mutational signature was pervasive in metastatic LTUC and less so in UTUC. Despite a lower overall sSNV burden, UTUC displayed greater inter- and intra-individual genomic distances at the copy number level between primary and metastatic tumors than LTUC. Our data also indicate that metastatic UTUC lesions can arise from small clonal populations present in the primary cancer. Importantly, putative druggable mutations were found across patients with the majority shared across all metastases within a patient. CONCLUSIONS: Metastatic UTUC demonstrated a lower overall mutational burden but greater structural variability compared to LTUC. Our findings suggest that metastatic UTUC displays a greater spectrum of copy number divergence from LTUC. Importantly, we identified druggable lesions shared across metastatic samples, which demonstrate a level of targetable homogeneity within individual patients.

Publisher Correction: The molecular landscape of pediatric acute myeloid leukemia reveals recurrent structural alterations and age-specific mutational interactions.

In the version of this article originally published, the color key in Fig. 1a was wrong. In the Cytogenetics key, the box over t(8;21) originally was green. It should have been red, matching the color of the sections of the pie graphs below the key that were labeled with 15% and 19%.

Crebbp Loss Drives Small Cell Lung Cancer and Increases Sensitivity to HDAC Inhibition.

CREBBP, encoding an acetyltransferase, is among the most frequently mutated genes in small cell lung cancer (SCLC), a deadly neuroendocrine tumor type. We report acceleration of SCLC upon Crebbp inactivation in an autochthonous mouse model. Extending these observations beyond the lung, broad Crebbp deletion in mouse neuroendocrine cells cooperated with Rb1/Trp53 loss to promote neuroendocrine thyroid and pituitary carcinomas. Gene expression analyses showed that Crebbp loss results in reduced expression of tight junction and cell adhesion genes, including Cdh1, across neuroendocrine tumor types, whereas suppression of Cdh1 promoted transformation in SCLC. CDH1 and other adhesion genes exhibited reduced histone acetylation with Crebbp inactivation. Treatment with the histone deacetylase (HDAC) inhibitor Pracinostat increased histone acetylation and restored CDH1 expression. In addition, a subset of Rb1/Trp53/Crebbp-deficient SCLC exhibited exceptional responses to Pracinostat in vivo Thus, CREBBP acts as a potent tumor suppressor in SCLC, and inactivation of CREBBP enhances responses to a targeted therapy.Significance: Our findings demonstrate that CREBBP loss in SCLC reduces histone acetylation and transcription of cellular adhesion genes, while driving tumorigenesis. These effects can be partially restored by HDAC inhibition, which exhibited enhanced effectiveness in Crebbp-deleted tumors. These data provide a rationale for selectively treating CREBBP-mutant SCLC with HDAC inhibitors. Cancer Discov; 8(11); 1422-37. ©2018 AACR. This article is highlighted in the In This Issue feature, p. 1333.

Copy number profiling across glioblastoma populations has implications for clinical trial design.

Background: Copy number alterations form prognostic molecular subtypes of glioblastoma with clear differences in median overall survival. In this study, we leverage molecular data from several glioblastoma cohorts to define the distribution of copy number subtypes across random cohorts as well as cohorts with selection biases for patients with inherently better outcome. Methods: Copy number subtype frequency was established for 4 glioblastoma patient cohorts. Two randomly selected cohorts include The Cancer Genome Atlas (TCGA) and the German Glioma Network (GGN). Two more selective cohorts include the phase II trial ARTE in elderly patients with newly diagnosed glioblastoma and a multi-institutional cohort focused on paired resected initial/recurrent glioblastoma. The paired initial/recurrent cohort also had exome data available, which allowed for evaluation of multidimensional scaling analysis. Results: Smaller selective glioblastoma cohorts are enriched for copy number subtypes that are associated with better survival, reflecting the selection of patients who do well enough to enter a clinical trial or who are deemed well enough to undergo resection at recurrence. Adding exome data to copy number data provides additional data reflective of outcome. Conclusions: The overall outcome for diffuse glioma patients is predicted by DNA structure at initial tumor resection. Molecular signature shifts across glioblastoma populations reflect the inherent bias of patient selection toward longer survival in clinical trials. Therefore it may be important to include molecular profiling, including copy number, when enrolling patients for clinical trials in order to balance arms and extrapolate relevance to the general glioblastoma population.

A De Novo Mouse Model of C11orf95-RELA Fusion-Driven Ependymoma Identifies Driver Functions in Addition to NF-κB.

The majority of supratentorial ependymomas (ST-ependymomas) have few mutations but frequently display chromothripsis of chromosome 11q that generates a fusion between C11orf95 and RELA (RELAFUS). Neural stem cells transduced with RELAFUSex vivo form ependymomas when implanted in the brain. These tumors display enhanced NF-κB signaling, suggesting that this aberrant signal is the principal mechanism of oncogenesis. However, it is not known whether RELAFUS is sufficient to drive de novo ependymoma tumorigenesis in the brain and, if so, whether these tumors also arise from neural stem cells. We show that RELAFUS drives ST-ependymoma formation from periventricular neural stem cells in mice and that RELAFUS-induced tumorigenesis is likely dependent on a series of cell signaling pathways in addition to NF-κB.

Erratum: The molecular landscape of pediatric acute myeloid leukemia reveals recurrent structural alterations and age-specific mutational interactions.

This corrects the article DOI: 10.1038/nm.4439.

The molecular landscape of pediatric acute myeloid leukemia reveals recurrent structural alterations and age-specific mutational interactions.

We present the molecular landscape of pediatric acute myeloid leukemia (AML) and characterize nearly 1,000 participants in Children's Oncology Group (COG) AML trials. The COG-National Cancer Institute (NCI) TARGET AML initiative assessed cases by whole-genome, targeted DNA, mRNA and microRNA sequencing and CpG methylation profiling. Validated DNA variants corresponded to diverse, infrequent mutations, with fewer than 40 genes mutated in >2% of cases. In contrast, somatic structural variants, including new gene fusions and focal deletions of MBNL1, ZEB2 and ELF1, were disproportionately prevalent in young individuals as compared to adults. Conversely, mutations in DNMT3A and TP53, which were common in adults, were conspicuously absent from virtually all pediatric cases. New mutations in GATA2, FLT3 and CBL and recurrent mutations in MYC-ITD, NRAS, KRAS and WT1 were frequent in pediatric AML. Deletions, mutations and promoter DNA hypermethylation convergently impacted Wnt signaling, Polycomb repression, innate immune cell interactions and a cluster of zinc finger-encoding genes associated with KMT2A rearrangements. These results highlight the need for and facilitate the development of age-tailored targeted therapies for the treatment of pediatric AML.

Multidimensional scaling of diffuse gliomas: application to the 2016 World Health Organization classification system with prognostically relevant molecular subtype discovery.

Recent updating of the World Health Organization (WHO) classification of central nervous system (CNS) tumors in 2016 demonstrates the first organized effort to restructure brain tumor classification by incorporating histomorphologic features with recurrent molecular alterations. Revised CNS tumor diagnostic criteria also attempt to reduce interobserver variability of histological interpretation and provide more accurate stratification related to clinical outcome. As an example, diffuse gliomas (WHO grades II-IV) are now molecularly stratified based upon isocitrate dehydrogenase 1 or 2 (IDH) mutational status, with gliomas of WHO grades II and III being substratified according to 1p/19q codeletion status. For now, grading of diffuse gliomas is still dependent upon histological parameters. Independent of WHO classification criteria, multidimensional scaling analysis of molecular signatures for diffuse gliomas from The Cancer Genome Atlas (TCGA) has identified distinct molecular subgroups, and allows for their visualization in 2-dimensional (2D) space. Using the web-based platform Oncoscape as a tool, we applied multidimensional scaling-derived molecular groups to the 2D visualization of the 2016 WHO classification of diffuse gliomas. Here we show that molecular multidimensional scaling of TCGA data provides 2D clustering that represents the 2016 WHO classification of diffuse gliomas. Additionally, we used this platform to successfully identify and define novel copy-number alteration-based molecular subtypes, which are independent of WHO grading, as well as predictive of clinical outcome. The prognostic utility of these molecular subtypes was further validated using an independent data set of the German Glioma Network prospective glioblastoma patient cohort.

Mutant IDH1 regulates the tumor-associated immune system in gliomas.

Gliomas harboring mutations in isocitrate dehydrogenase 1/2 (IDH1/2) have the CpG island methylator phenotype (CIMP) and significantly longer patient survival time than wild-type IDH1/2 (wtIDH1/2) tumors. Although there are many factors underlying the differences in survival between these two tumor types, immune-related differences in cell content are potentially important contributors. In order to investigate the role of IDH mutations in immune response, we created a syngeneic pair mouse model for mutant IDH1 (muIDH1) and wtIDH1 gliomas and demonstrated that muIDH1 mice showed many molecular and clinical similarities to muIDH1 human gliomas, including a 100-fold higher concentration of 2-hydroxygluratate (2-HG), longer survival time, and higher CpG methylation compared with wtIDH1. Also, we showed that IDH1 mutations caused down-regulation of leukocyte chemotaxis, resulting in repression of the tumor-associated immune system. Given that significant infiltration of immune cells such as macrophages, microglia, monocytes, and neutrophils is linked to poor prognosis in many cancer types, these reduced immune infiltrates in muIDH1 glioma tumors may contribute in part to the differences in aggressiveness of the two glioma types.

miR-155 expression and correlation with clinical outcome in pediatric AML: A report from Children's Oncology Group.

BACKGROUND: Aberrant expression of microRNA-155 (miR-155) has been implicated in acute myeloid leukemia (AML) and associated with clinical outcome. PROCEDURE: We evaluated miR-155 expression in 198 children with normal karyotype AML (NK-AML) enrolled in Children's Oncology Group (COG) AML trial AAML0531 and correlated miR-155 expression levels with disease characteristics and clinical outcome. Patients were divided into quartiles (Q1-Q4) based on miR-155 expression level, and disease characteristics were then evaluated and correlated with miR-155 expression. RESULTS: MiR-155 expression varied over 4-log10-fold range relative to its expression in normal marrow with a median expression level of 0.825 (range 0.043-25.630) for the entire study cohort. Increasing miR-155 expression was highly associated with the presence of FLT3/ITD mutations (P < 0.001) and high-risk disease (P < 0.001) and inversely associated with standard-risk (P = 0.008) and low-risk disease (P = 0.041). Patients with highest miR-155 expression had a complete remission (CR) rate of 46% compared with 82% in low expressers (P < 0.001) with a correspondingly lower event-free (EFS) and overall survival (OS) (P < 0.001 and P = 0.002, respectively). In a multivariate model that included molecular risk factors, high miR-155 expression remained a significant independent predictor of OS (P = 0.022) and EFS (0.019). CONCLUSIONS: High miR-155 expression is an adverse prognostic factor in pediatric NK-AML patients. Specifically, high miR-155 expression not only correlates with FLT3/ITD mutation status and high-risk disease but it is also an independent predictor of worse EFS and OS.